This is a test version of Biostars. For the public version, visit https://www.biostars.org.
WGS variant filtering for exome

Hello everyone, I have a set of WGS data that I have processed until variant calling with freebayes, but I am only interested in conserving the mutations that are found in exonic regions, which would be the most reliable way of accomplishing this?

Thanks in advanced

wes snp strelka freebayes wgs
• 696 views • link
updated by Pierre Lindenbaum Pierre Lindenbaum • written by isMM • isMM

1 answer 

gunzip -c /path/to/your.gtf.gz | awk -F '\t' '($3=="exon") {printf("%s\t%d\t%s\n",$1,int($4)-1,$5);}' | LC_ALL=C sort -T . -t $'\t' -k1,1 -k2,2n | bedtools merge > exons.bed

bcftools view --regions-file exons.bed indexed.vcf.gz
• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Log in to answer this question.

More posts like this