I am trying to do this
samtools view -f 0x2 -b XX.bam | samtools view -b -F 0x200 - | samtools view -b -F 0x400 - |genomeCoverageBed
as I am not sure of how to combine bitwise flag.
Unfortunately this 'lazy' way consumes 2 extra threads.
There's a python script here http://biostar.stackexchange.com/questions/13079/flag-in-sam-format
But I am not sure how to change it so that I might be able to filter in one samtools step? (or two?)
•
written
by
Kevin
2 answers
Seems a bit silly to answer my own question but can anyone verify if what I can do is
samtools view -f 0x2 -b XX.bam | samtools view -b -F 1536|genomeCoverageBed
combined the last two commands bit flag here http://picard.sourceforge.net/explain-flags.html
Hi Kevin,
That looks right to me. 0x200 = 2^9 = 512, 0x400 = 2^10 = 1024, 512 + 1024 = 1536.
Log in to answer this question.
-
How to plot alignment statistics; how many reads are mapped to the genome?
written by anamaria
-
Samtools SORT or bedtools bamToFastq giving error when converting reads to fastq.
written by DNAngel
-
Variants calling by bcftools
written by Mbillah
-
Error in bcftools view -bvcg command
written by hafiz.talhamalik
-
Using Diffbind R-studio
written by dimitrischat
-
Tutorial: How to select reads with specific flag using samtools
written by imlituan •
-
Tutorial: Convert bam files to fastq in numbers as indicated in samtools flagstat stats.
written by kirannbishwa01
-
Extracting mapped paired end FastQ reads from BAM file using Samtools - trying to understand -f 0x0…
written by acxcrew007 •
-
simulating samtools -f option in pysam
written by sa9
-
Mapping shorts reads using Bowtie
written by bioinfo
