What criteria are used to determine supplementary alignment with bwa-mem (flag 2048)
Greetings,
Supplementary alignment (chimeric reads, flag 2048) is set when a read aligns to 2 or more different genome locations. What is the minimum number of bases that must align to a different location to be considered a chimeric read? I have reads marked as SA when >25 bases align elsewhere, but not 8, for example. Does it depend on read length? I have bwa-mem alignments with 45-bp paired-end reads and none of the reads are set as SA reads (although I know they are), while alignments using 75-bp reads (same genome) have plenty of SA reads. BWA versions :0.7.7/0.7.12
Thank you!
Anna
• 6,074 views
•
link
updated
by
Ram
•
written
by
annawilliford •
•
written
by
annawilliford •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Extracting complete supplementary alignments for long reads from a bam file
written by aacomeau •
-
The effect of heterozygosity on genome alignment
written by Iseez •
-
Total number of trimmed reads less than number of mapped reads
written by kspata
-
Get read start site using bedtools from bam file after merging paired end data
written by snishtala03 •
-
How do I tell bwa to map only unmapped portion of the read for supplementary alignments?
written by Matteo Schiavinato
-
Can a part of a read be in the alignment and in the supplementary alignment (-F 0x800) as well?
written by Matteo Schiavinato
-
Difference between chimeric alignments and multiple mapping
written by Vanilla
-
What are the sematics of containing bwa split read alignments?
written by d-cameron
-
Determine relative position (start and end) of chimeric secondary alignment within original long re…
written by mattia.bosio85 •
-
insert size discrepencies using bwa backtrack and bwa mem
written by urmi208 •
The odds are good that you'll need to go through the code to find the answer to this. The paper describing bwa mem is only a page long, so details like this will be omitted.
you could also try tagging your post bwa and hope that lh3 replies