Question: .Ab1 Trace File Quality Threshold
2
gravatar for Nicolas Rosewick
8.0 years ago by
Belgium, Brussels
Nicolas Rosewick8.4k wrote:

Hi,

I've several .ab1 files from a sanger sequencing experiment. What is a good threshold for the quality of the sequencing ? And is it a way to export them into a fastq file. It would simplier to assemble them into a longer contig.

Thanks,

N.

fastq quality • 4.3k views
ADD COMMENTlink modified 8.0 years ago by lh331k • written 8.0 years ago by Nicolas Rosewick8.4k
4
gravatar for ALchEmiXt
8.0 years ago by
ALchEmiXt1.9k
The Netherlands
ALchEmiXt1.9k wrote:

In addition to the answer by Casey you could use "old class" assemblers that handle qualityscores directly like (academic free) PHRED/PHRAP/CONSED and their derivatives. Alternatively use some DNA sequence analysis solution as provided (commercially) by for instance DNASTAR or CLC bio.... But most likely there will be free substitutes for those as well. We do all our Sanger reads this way.

ADD COMMENTlink written 8.0 years ago by ALchEmiXt1.9k

+1. Sanger ABI reads predate fastq by some years and there are lots of free tools available which may be "old", but still work perfectly well.

ADD REPLYlink written 8.0 years ago by Neilfws48k

Another old-school assembler that's easy to use and works well on fasta + quals is PCAP: http://seq.cs.iastate.edu/

ADD REPLYlink written 8.0 years ago by Casey Bergman18k
1
gravatar for Casey Bergman
8.0 years ago by
Casey Bergman18k
Athens, GA, USA
Casey Bergman18k wrote:

For a single Sanger read, you should be using at least Q>20 and preferable Q>30, although there is no hard and fast rule. Read more about what Q-scores mean here: http://en.wikipedia.org/wiki/Phred_quality_score

Galaxy has a fasta+qual -> fastq converter under NGS: QC and manipulation -> Roche-454 data -> Combine FASTA and QUAL into FASTQ. You could use this on the web or pull this tool out of the galaxy distribution and run it locally.

ADD COMMENTlink written 8.0 years ago by Casey Bergman18k
1
gravatar for lh3
8.0 years ago by
lh331k
United States
lh331k wrote:

See io_lib from Staden.

As to trimming, you should consider the algorithm used by phrap/phred, which is later used in bwa. Setting a simple per-base threshold is not the best way to trim reads. The bwa manual gives the algorithm in the -q option of the aln command. It is quite simple once you understand the math. There is a seqanswers thread explains what it does, but I cannot find it now.

ADD COMMENTlink written 8.0 years ago by lh331k
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: 1847 users visited in the last hour