Help: Phred quality for ab1 files

Hello,

I'm using the biojava legacy release of 1.8.1 (API here: http://www.biojava.org/docs/api1.8/)

I want to get the Phred quality values for single bases in an ab1 chromatogram. Using methods in the ABITools class, I was able to create a SymbolList of the sequence.

I'm now using the PhredTools class, but I don't know where to start.

Thanks

phred chromatogram abif biojava sequence

• 3,599 views • link

updated 4.1 years ago by Ram

• written 10.5 years ago by rabidbirdman4 •

1 answer

Indigo analyses ab1 files and aligns them to a reference. The command-line version available here also dumps the chromatogram including the base-calls and the Phred quality values so this might be what you are looking for.

• 0 views • link

written 9.4 years ago by trausch

More posts like this

Using RSEM after Cuffcompare
written 2.0 years ago by ZZOYEA •

Hello, I am currently using the Tuxedo pipeline to identify lncRNAs. I obtained transcript assemblies for each sample using cufflinks, combined them into a single …
single cell ATAC-seq cicero gene activities in Seurat
written 6.0 years ago by el24 •

Hi everyone! I am loading my single-cell ATACseq data using cicero, and then I create gene activity scores according to the [tutorial][1]. My goal is …
How to write a BGZF file with Biopython?
written 6.8 years ago by mrals89

Hi everyone, I've read the documentation for [`Bio.bgzf`](https://biopython.org/DIST/docs/api/Bio.bgzf-module.html) and I'm sufficiently confused by it. To create an example bgzf file, I've tried the following. The …
Capillary sequence trimming algorithm
written 7.7 years ago by tospo •

I need to implement a quality trimming algorithm for capillary sequence reads (yes, they still exist!) in a Perl package of my own, so I …
CWL: default values for custom types in a workflow
written 9.1 years ago by thomas.e

Hello, How do I set the default value for a custom type in a workflow. Something like: steps: mystep: run: mything.cwl in: important_arg: default: magicValue …
Determining inhouse cut offs
written 10.0 years ago by skbrimer

Hello group, My boss wants a depth of coverage vs quality of data (i.e. Q10=100x, Q20 =50, Q30=25, etc...) so I'm not sure how to …
Phred Score/Base quality score in Pysam
written 10.9 years ago by Lylthera •

Hi, I'm examining the phred quality scores for bases of certain snps, so I wonder how to get this information - I'm using pysam and …
Using Pbase bioconductor package to map protein substitutions to chromosome positions
written 11.0 years ago by Eli Korvigo

I'm using trying to use [this tutorial][1] to map certain protein positions to genomic positions. The whole tutorial is focused on using the `Proteins` class …
What is difference between Ab1 PCON values and phd.1 file quality scores
written 11.0 years ago by crivenster

I am new to sangar sequencing files. I am using the ABIF parser from biojava legacy package to load the raw data and process the …
LOH and CNV data from VCF files
written 11.7 years ago by GarF •

I was asked to find both the CNV and the LOH datas from 3 VCF files (all from the same patient but 2 of them …