Google Groups no longer supports new Usenet posts or subscriptions. Historical content remains viewable.
Dismiss

[Automated-sequencing] phrap contig ,apply contig quality values?

21 views
Skip to first unread message

Umanga

unread,
Feb 25, 2009, 1:45:51 AM2/25/09
to aut...@magpie.bio.indiana.edu

Greetings all,
Firstof all, come from Computer science background ,so I have very little
knowledge in Bioinformatics.
Now my question ;
I assemble AB1 files using phred/phrap.I noticed that after assembling using
'phrap', it generates contig file (say 'seq.contigs') and a contig-quality
file ('seq.contigs.qual').

Now,I want to get the same contig which generated by ChromasPro using phrap.

I noticed that ,if I ignore the config bases with quality value '0' (which
stored in seq.contigs.qual) , I can get almost same sequence as ChromasPro.
Is this correct?
This there a parameter to do this in 'phrap' itself?

Thanks in advance,
Umanga
--
View this message in context: http://www.nabble.com/phrap-contig-%2Capply-contig-quality-values--tp22197180p22197180.html
Sent from the Bio.net - Autoseq mailing list archive at Nabble.com.

Phillip San Miguel

unread,
Feb 25, 2009, 7:37:05 AM2/25/09
to bionet-genome-...@moderators.isc.org
Umanga wrote:
> Greetings all,
> Firstof all, come from Computer science background ,so I have very little
> knowledge in Bioinformatics.
> Now my question ;
> I assemble AB1 files using phred/phrap.I noticed that after assembling using
> 'phrap', it generates contig file (say 'seq.contigs') and a contig-quality
> file ('seq.contigs.qual').
>
> Now,I want to get the same contig which generated by ChromasPro using phrap.
>
> I noticed that ,if I ignore the config bases with quality value '0' (which
> stored in seq.contigs.qual) , I can get almost same sequence as ChromasPro.
> Is this correct?
> This there a parameter to do this in 'phrap' itself?
>
> Thanks in advance,
> Umanga

Hi Umanga,
Yes, you should remove the bases with zero quality values from the ends
of phrap contigs. A quality value of zero indicates the base call is
random noise. I usually trim the ends of contigs back until I hit the
first base with a 20 or higher quality value.
I don't know of a method to force phrap to do this. It would be nice.
I've never used Chromas or ChromasPro. But different assembly engines
may produce different contigs from the same data set.
Phillip

0 new messages