Hi, I have one BAM file which contains all alignments (include those not used in variant calling, such as non-PF, non-mapping and duplicate reads) generated for an assembly. How to filter these useless mapping? I know that Picard MarkDuplicates can be used to remove duplicates.Thank you.
2 answers
You can use samtools to do this. e.g. to remove reads that did not align, you can do:
samtools view -F 0x04 -b in.bam > out.aligned.bam
to only include paired reads, use:
-f 0x02
Check the other bitwise flags on this page.
But you might not want to exclude those as they could be used for finding structural variations.
The bamtools package offers a wide range of filters, including user-definable filters defined in JSON notation. It includes filters for reads failing vendor QC, unmapped reads and pre-marked duplicates.
Log in to answer this question.
-
remove duplicates in Nanopore cDNA PCR data
written by 静 •
-
Unable to remove duplicate sequences
written by juanluis0516 •
-
DNA sequencing workflow involving UMIs
written by ww22runner •
-
Can MarkDuplicates of Picard be used for RNA reads?
written by Nemo •
-
Variants in untargeted genes identified after targeted exome sequencing analysis
written by vinayjrao
-
Effect of removing secondary/supplementary alignments from BAM file?
written by cg1440 •
-
Extracting reads from Pacbio Sequel BAM by read name
written by roryC •
-
Picard remove duplicate leads to some unpaired reads
written by billzt •
-
bam processing before variant calling
written by user230613
-
Remove optical duplicates from bam file
written by Anna S
