Dear all,
I have a bamfile which contains some reads that have supplementary alignments within them, as shown in the screenshot.
I tried to filter the bamfile with samtools using flag 2048. Luckily the supplementary alignments belonging to a specific read were gone. However, I want the whole read to be filtered out whenever it contains supplementary alignments from the bamfile and not only its supplementary alignments. is there a possible way to do this?
Best
•
written
by
mokbel73 •
3 answers
add -e '![SA]' to your samtools view command ( http://www.htslib.org/doc/samtools.html#FILTER_EXPRESSIONS )
In one step, you could remove every alignments with a mapping quality of 0. With most aligners, this corresponds to the reads for which there is ambiguous mapping. Note that this will also remove secondary alignments and the primary alignments from reads that also have secondary/supplementary alignments.
samtools view -q 1 -b input.bam -o output.bam
Log in to answer this question.
-
Does sam 2048 flag discriminate secondary from supplementary alignments?
written by resug •
-
Extracting complete supplementary alignments for long reads from a bam file
written by aacomeau •
-
IGV grey vs white lines
written by Mo •
-
[pysam] How to find both both primary and supplementary read？
written by octpus616 •
-
Remove supplementary alignments from a bam file
written by jeni •
-
Multiple mapping with longreads - Get primary mapping
written by ricardoguerreiro2121
-
Samtools view -f "OR" operator possible?
written by QVINTVS_FABIVS_MAXIMVS
-
How to retain SEQ and QUAL fields in BAM file after filtering?
written by toddknutson
-
How do I tell bwa to map only unmapped portion of the read for supplementary alignments?
written by Matteo Schiavinato
-
What criteria are used to determine supplementary alignment with bwa-mem (flag 2048)
written by annawilliford •
for what it's worth, "supplementary alignments" can be quite interesting. they often indicate things like SVs