Question: Extraction of aligned read using samtools
0
gravatar for Neu
8 days ago by
Neu10
India
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 • 60 views
ADD COMMENTlink written 8 days 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 8 days 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 2 days ago • written 2 days ago by ATpoint13k
Please log in to add an answer.

Help
Access

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