Filter out supplementary aligned reads for variant calling or not
0
0
Entering edit mode
4.1 years ago
gb ★ 2.2k

Dear all,

I am using BWA-mem and freebayes to do some variant calling and I was wondering if I should remove the supplementary aligned reads or not. I can not find much info about this.

freebayes BWA snp • 1.1k views
ADD COMMENT
1
Entering edit mode

In general it is good practice to not do any custom filtering unless the manual of the tool you use specifically recommends it. You should have duplicates marked in the BAM. Any proper variant caller will then internally choose the reads that based on its underlying model or strategy are suitable. The README (https://github.com/ekg/freebayes/blob/master/README.md) does not mention any prefiltering so you should not do it.

ADD REPLY
0
Entering edit mode

I think that is a good point. For now I will not worry about this.

ADD REPLY
0
Entering edit mode

If you later during analysis of the results see that something is wrong you can still dig deeper if supp. alns might be a factor but I doubt so. I learned over the years that one should not overthink things and only start to debug if there is evidence for something going wrong. It might cost a lot of time without making a positive difference.

ADD REPLY

Login before adding your answer.

Traffic: 1145 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6