Get Phred Score/Base Quality for BAM File Per base
Hi all,
I have a BAM file I'm try to get some base quality scores from in the form of a bed file. I have used QualityScoreDistribution from GATK, and that works at giving overall scores, but not per base metric. To explain this a little better:
I have a BAM file and an associated BED file for the sequences I want to get phred scores from. If I could somehow use these two files, and then generate a file that has Chromosome Pos and Phred Score.
If anybody could point me in the right direction? Thanks
• 1,757 views
•
link
updated
by
Pierre Lindenbaum
•
written
by
gtasource
•
written
by
gtasource
1 answer
Use samtools mpileup to get all the qualities at a given position.
• 0 views
•
link
written
by
Pierre Lindenbaum
Log in to answer this question.
More posts like this
-
How to make bam from start pos of read 1 and end pos of read 2 from paired end bam files
written by genomics_student •
-
Is it okay to convert bam files to fastq and then getting seqkit results?
written by Baran •
-
Issue with bed to bam conversion - sequence in bam headers
written by kperlejewski •
-
Using Samtools depth to get number of forward and reverse reads per base from a .bam file
written by jspainhour3 •
-
Annotate pairwise alignment with query sequence quality scores
written by Thomas
-
Get Introns from a gft annotation file
written by xd_d
-
Phred Score/Base quality score in Pysam
written by Lylthera •
-
Find reference base for a position in a given read from a BAM file
written by Andrew •
-
Samtools mpileup quality scores all zero
written by fwuffy
-
Ngs Question ~ Consensus
written by Dirkjan Schokker •
