Splitting paired end bam files
Hello all,
Any efficient way to split the bam file to two separate files, say one file with only singletons(mate unmapped) and other with paired reads(no matter if it is discordant).
I have already used.
(i) samtools view -f 8 -F 4 -b foo.bam > foo.singletons.bam
(ii) samtools view -f 1 -F 12 -b foo.bam > foo.paired.bam
Is this correct solution to go ahead.
PS: This file has already been filtered for uniquely mapped reads using NH:i:1 tags.
• 4,829 views
•
link
updated
by
Dan D
•
written
by
tharveshliyakat
•
written
by
tharveshliyakat
1 answer
That looks good to me. To save space, you might also try outputting to a named pipe:
Log in to answer this question.
More posts like this
-
Samtools mapped and unmapped flags, what is the difference between the flags?
written by DNAngel
-
Obtaining mapped reads after extracting unmapped reads .
written by shrivadeepak •
-
Samtools SORT or bedtools bamToFastq giving error when converting reads to fastq.
written by DNAngel
-
Error when trying to convert BAM file to FASTQ.
written by DNAngel
-
Apply conditions to forward/reverse strands on BAM file (RNA-Seq)
written by ek699 •
-
Filter bam files using a bed file: Why is the mate missing?
written by komal.rathi
-
Request to validate my RNAseq workflow
written by umn_bist •
-
Thorough Documentation Of Samtools Flagstat
written by Biomonika (Noolean)
-
Structural Variation; Discordant Pairs?
written by Rm
-
How To Remove The Qc-Failure Reads?
written by Junfeng
if you -f 1 to the second command, you should really -f 9 to the first to be fair :)
Even though its unlikely youll ever come across a bam file with both single reads and paired reads, the first command would match both. The single reads it would match would also be pretty random - not all of them.