Yes, it works either way. CNVkit will ignore secondary alignments if they are present.
CNVkit bwa-mem parameters
Which are the correct BWA-mem parameters suggested/required/tolerated by CNVKit ? Are -M and -a parameters ok ? Thx a lot
• 1,766 views
•
link
2 answers
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.
• 0 views
•
link
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 ? )
• 0 views
•
link
• 0 views
•
link
Log in to answer this question.