Is there any flag to remove broken (ONT reads split and map to different locations) reads from bam file.
Note: The ONT reads mapped with gap should not remove.
•
written
by
BioGeek
1 answer
As Pierre Lindenbaum said in the comment, samtools view -F 2048 will remove all those reads that have supplementary alignment. Be sure that the alignment program you use specifies those (like bwa mem does, for example).
If you also want to remove reads that have secondary alignments, there is no flag to do so but you can work around it with awk. Combining it with the one above:
samtools view -h -F 2048 file.bam | awk 'substr($1, 0, 1)=="@" || $0 !~ /ZS:/' | samtools view -h -b > filtered_file.bam
Be careful: ZS: is a tag of HISAT2. If you're using another program, check which tag is used for secondary alignments. I think it is AS: in Bowtie2, and I don't know what it is in bwa. Usually, it is specified in the manual.
Log in to answer this question.
-
Total Alignment: 0 /results of featureCounts RNA-Seq
written by santos48 •
-
Salmon Quant /SubRead featureCounts Results Total Aligments:0
written by santos48 •
-
How to choose the best alignment from mapping with bowtie2?
written by anikcropscience •
-
Reads that map to two different genomes, but with different map coordinates
written by Susmita Mandal
-
PacBio and ONT reads difference !
written by BioGeek
-
How to calculate fraction of RNA-Seq reads mapped to a customised genomic locations?
written by biorepine
-
Where to find PCR duplicate reads in bam file?
written by bioinforesearchquestions
-
flagstats results show 0.0 % properly paired reads
written by SOHAIL
-
bam processing before variant calling
written by user230613
-
How Should I Compute Coverage When Some Reads Map To Multiple Locations?
written by Ryan Thompson
samtools view with flag 2048 ?