Biostar Beta. Not for public use.
bwa mem without -M option how to continue using picard
0
Entering edit mode
14 months ago
Peter Chung • 80
Hong Kong

I do the alignment without -M option and now I have all the bam file and would like to do the downstream steps using picard and GATK, what can I do to restore the -M option except realign all the fastq because I have more than hundreds fastq files. Thank you

ADD COMMENTlink
0
Entering edit mode
11 months ago
d-cameron ♦ 2.0k
Australia

This can be achieved with software that can edit SAM files (such as samtools and sambamba)

You want to set the secondary alignment flag (0x100) for all reads that have the supplementary alignment flag set. (0x800).

ADD COMMENTlink
0
Entering edit mode
12 months ago
WCIP | Glasgow | UK

Is it still the case that picard (MarkDuplicates?) chokes on split reads? I would check that first. Using -M option is a hack to make picard work but it would be better avoiding it since the shorter split of a read is not a "secondary alignment" and there is nothing necessarily wrong with split reads. They may be the result of structural variants.

ADD COMMENTlink

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.1