Splitting BAM file from paired end into forward and reverse strands?
Hi,
I have splitted BAM file from single end reads into forward and reverse strands using the following command
samtools view -b -F 16 *for forward strand
samtools view -b -f 16 *for reverse strand
As I understand the sum of these two should be equal to the main BAM file. So my questing is for paired end reads, will it be different or the same code will work for both single and paired end reads? Help me understand.
Thanks,
Susmita
• 2,561 views
•
link
written
by
Susmita Mandal
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
RNA-Seq - Strand-specific coverage for a paired library
written by npb27 •
-
Strand bias in supplementary alignments from Nanopore
written by dariober
-
Viral positive and negative strand with paired sequencing and bowtie
written by Lluís R.
-
Split paired-end reads BAM file into Forward and Reverse Strands
written by GA •
-
Reads mapped on wrong strand
written by qstefano •
-
Split BAM file from paired end read into forward and reverse strand w.r.t reference genome
written by Susmita Mandal
-
SAM flag and select reads that map uniquely
written by epigene
-
BAM Flag From Bismark Alignment: Forward and Reverse Chain
written by Shicheng Guo
-
Single-end strand specific Rna-Seq: how to separate alignments by strand
written by termanini •
-
Tutorial: How To Separate Illumina Based Strand Specific Rna-Seq Alignments By Strand
written by Istvan Albert
That code will pay no attention to any of the flags having to do with paired reads. It also won't pay attention to whether or not the read mapped.
Can you be little more elaborative? I did not get it.