Files
coriell_samples/gene_bed.sh

8 lines
224 B
Bash

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