Question: Extraction of aligned read using samtools
gravatar for Neu
4 months ago by
Neu10 wrote:

Hi, I have aligned paired end RNA-Seq data and have the BAM file. I want to get the aligned reads using samtools. My command is : samtools view -b -F 4 filename.bam > mapped_out.bam Is it right for the paired end data? Also, when I am using the command, the generated file is found to be bigger than the input file. I would be grateful if anyone can help me to understand this.

sequencing rna-seq next-gen • 153 views
ADD COMMENTlink written 4 months ago by Neu10

I think this answer can help you figure out how to do the job. how to extract pair end reads from bam file using samtool

ADD REPLYlink written 4 months ago by anicet.ebou140

Please use the search function, this question has been asked very often before. File size is not a good measure for anything. Count reads with samtools flagstat instead as file size depends on data composition and compression level.

ADD REPLYlink modified 4 months ago • written 4 months ago by ATpoint17k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1604 users visited in the last hour