1dedbd46c366091fd9d93d961b4b7b91d83ad743 max Mon Oct 2 07:50:55 2023 -0700 adding polygenic risk scores track, refs #31546 diff --git src/hg/makeDb/scripts/prsEmerge/prs2bigBed.sh src/hg/makeDb/scripts/prsEmerge/prs2bigBed.sh new file mode 100644 index 0000000..7a1420e --- /dev/null +++ src/hg/makeDb/scripts/prsEmerge/prs2bigBed.sh @@ -0,0 +1,17 @@ +# written by Zia Truong and slightly adapted by Max and committed +echo "Converting PRS files to BED files..." +mkdir output +wait + +python3 ~/kent/src/hg/makeDb/scripts/prsEmerge/prs2bed.py # iterates over files in data/, writes to output/ + +echo "Converting BED files to bigBED files..." +#looping through files for loop code based on https://stackoverflow.com/a/10523501 +for bed in ./output/*.bed +do + sort -k1,1 -k2,2n "${bed}" -o "${bed}" + bedToBigBed -as=prs.as -type=bed9+2 "${bed}" /hive/data/genomes/hg19/chrom.sizes "${bed}.bb" + rm ${bed} +done + +echo "Done!"