Question: flagstat interpretation from a bam file
gravatar for A
5.4 years ago by
A3.9k wrote:

sorry friends,

i had a bam file but my adviser believes that because of some reason the results seem not significant, i used flagstat (below)

5541714 + 0 in total (QC-passed reads + QC-failed reads)
0 + 0 secondary
0 + 0 supplementary
0 + 0 duplicates
5541714 + 0 mapped (100.00%:-nan%)
0 + 0 paired in sequencing
0 + 0 read1
0 + 0 read2
0 + 0 properly paired (-nan%:-nan%)
0 + 0 with itself and mate mapped
0 + 0 singletons (-nan%:-nan%)
0 + 0 with mate mapped to a different chr
0 + 0 with mate mapped to a different chr (mapQ>=5)

may you please tell me how is the alignment??

thank you

samtools bam flagstat • 2.6k views
ADD COMMENTlink modified 5.3 years ago by Biostar ♦♦ 20 • written 5.4 years ago by A3.9k

Well, based only on this flagstat, you've got perfect alignment (all the reads are aligned perfectly). I will suspect that this is a post-processed bam file (not the raw bam alignment file). So we cannot tell yo anything about it. Also, what do you mean by not significant? Are you trying to do exome sequencing? RNA Sequencing? Significant in terms of what? Is it paired end data? 


ADD REPLYlink modified 5.4 years ago • written 5.4 years ago by Sam3.3k

thanks Sam,

I am doing ribo-seq and this is bam produced of footprints. but when trimming the adapters (illumina TruSeq Universal adapter), the rate of trimmed reads were too low, then supervisor told maybe the results are influenced of bad trimming or multi alignment..and asked me to check the bowtie2 output but from biostars I knew that by flagstat I can check the quality

now I don know if the reads aligned properly or not

ADD REPLYlink modified 14 months ago by _r_am32k • written 5.4 years ago by A3.9k

So let me be clear, what you have done is:

1. Trimming

2. Alignment (bowtie2)

3. Flagstat

Have you run fastQC to see how the reads quality looks? Also, when you mean the rate of trimmed reads were too low, you mean only a few reads need to be trimmed?

ADD REPLYlink written 5.4 years ago by Sam3.3k


i mean in the results, reads with adapter were just 4 percent

in fastqc result overrepresentation sequennces were vary based on the file

yeah all right i did which you mentioned

ADD REPLYlink written 5.4 years ago by A3.9k

Actually, you should try performing the fastQC before and after the trimming. That should let you know whether if there is any over-representative sequences in the start of your reads. If there is only a small amount of over-representative sequence (or adaptors, which fastQC can sometimes detect), then it is normal for you to only trim 4% of all your reads.

Usually, when we perform Exome Sequencing and RNA Seq, if the data is nice (no read-through), we can have nice fastq files that doesn't contain the adaptor sequences, therefore there is no need of trimming. However, as I have never played with Ribo-seq data, I am not sure if that is the case for you. 

ADD REPLYlink written 5.4 years ago by Sam3.3k

thank you Sam

ADD REPLYlink written 5.4 years ago by A3.9k
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: 1552 users visited in the last hour