diff --git a/gene_bed.sh b/gene_bed.sh new file mode 100644 index 0000000..798291d --- /dev/null +++ b/gene_bed.sh @@ -0,0 +1,7 @@ +zcat Homo_sapiens.GRCh38.115.gtf.gz \ + | awk '{if ($3 == "gene") {print $0}}' \ + | grep -E 'gene_name \"BRCA1\"|gene_name \"PALB2\"|gene_name \"CHEK2\"|gene_name \"ATM\"' \ + | gff2bed - \ + | sort -k1,1 -k2,2n \ + > test.bed +