caluculate q30 from bam file with samtools

Hello,

Is it possible to calculate Q30 from a bam with Samtools? I've looked at Samtools stats and didn't see that option.

Thanks

next-gen

• 2,830 views • link

written 6.5 years ago by genya35 •

Q30 for entire dataset or average per read?

• 0 views • link

written 6.5 years ago by GenoMax

One value for the entire bam file please.

• 0 views • link

written 6.5 years ago by genya35 •

Please use ADD REPLY and leave the answer field for answers. That keeps the thread logically organized.

• 0 views • link

written 6.5 years ago by ATpoint

1 answer

See answers in : FASTQ Phred33 average base quality score

• 0 views • link

written 6.5 years ago by GenoMax

More posts like this

Get discordant read pairs from bam
written 2.8 years ago by paulaotero.sanchez •

Hi, I would like to know if there is a way to isolate discordant read pairs from a bam file. For example, when I use …
Calling with samtools mpileup calls fewer SNPs than expected.
written 3.4 years ago by kimgeng •

Hello, I'm master course student, and I'm embarrased that I'm very poor at controlling bam files and samtools I tried to variant_calling with samtools mpileup, …
get % of >= Q30 Bases (PF) per index with InterOp
written 5.5 years ago by genya35 •

Hello, Does anyone know if it's possible to calculate % of >= Q30 Bases (PF) per 'sample_id/index' with InterOp. I'm aware that these are in …
calculating total yield for Illumina sequencing run
written 5.9 years ago by genya35 •

Hello, I'm trying to calculate ChipResultsSummary 'yield' for Hiseq sequencing run, that is listed in BustardSummary.xml file manually, so I can avoid parsing the file. …
blastn output format with full subject's sequence
written 7.4 years ago by genya35 •

Hello, Is there an outfmt option for blastn that will return the full subject's sequence in tubular format? I have tried outfmt 6 but it …
samtools and docker
written 8.2 years ago by genya35 •

Hello, I've pull and image secessfully and now trying to run samtools sort, however, docker is not finding the bam file. Any suggestions? Here is …
How to achieve the same alignment with HiSAT2 command line vs. HiSAT2 Galaxy
written 8.8 years ago by genya35 •

Hello, I've created a small reference fasta file that consists of combined sequences of two exons from two different genes where I know the fusion …
How can I use Samtools to map coverage distribution?
written 9.8 years ago by SaltedPork

I want to use Samtools to map copverage distribution and retrieve the regions with the highest coverage and view it in IGV as a track. …
samtools merge extremely slow
written 10.6 years ago by ehzed •

Hello, For parallel processing I split the original fastq files into smaller ones (1000000 reads per file) and performed trimming, alignment, post processing, etc. Now …
Use samtools index on multiple bam files at the same time from the command line?
written 10.6 years ago by nash.claire

Hi all, Simple question, is it possible to run samtools index (or any of the samtools in fact) on multiple bam files at the same …