CNVkit bwa-mem parameters
2
0
Entering edit mode
4 months ago
atariw • 0

Which are the correct BWA-mem parameters suggested/required/tolerated by CNVKit ? Are -M and -a parameters ok ? Thx a lot

CNVkit BWA-mem • 230 views
ADD COMMENT
0
Entering edit mode
4 months ago
Eric T. ★ 2.6k

Referring to the BWA MEM docs for posterity: http://bio-bwa.sourceforge.net/bwa.shtml

Yes, secondary alignments created with -a and marked as such with -M will be handled appropriately by CNVkit, i.e. excluded from the read depth count. CNVkit calculates depth using pysam, which uses htslib internally, and will therefore handle alignments and pileups the same way as other standard-compliant tools like samtools and bwa.

ADD COMMENT
0
Entering edit mode
3 months ago
atariw • 0

Thx. May I also run BWA without any option ? (i.e. would it work to detect CNV by off-target reads ? or -a is required ? )

ADD COMMENT
0
Entering edit mode

Yes, it works either way. CNVkit will ignore secondary alignments if they are present.

ADD REPLY

Login before adding your answer.

Traffic: 1032 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