Question: Assembly With Phrap With Some Missing Quality Files.
gravatar for Joseph Hughes
8.7 years ago by
Joseph Hughes2.7k
Scotland, UK
Joseph Hughes2.7k wrote:

I am trying to assemble sequences with phrap, for some of the sequences I have the associated quality file but not all of them. I have created a fasta file with all my sequences and another file with all the quality information but phrap is not happy as I have some missing quality information for a couple of the sequences:

FATAL ERROR: Inconsistency between quality file and sequence file

Is there a way I can phrap sequences together even if I have some missing quality information for some of the sequences? How?

quality assembly fasta • 1.4k views
ADD COMMENTlink modified 5 months ago by RamRS20k • written 8.7 years ago by Joseph Hughes2.7k
gravatar for Darked89
8.7 years ago by
Barcelona, Spain
Darked894.2k wrote:

You may try to experiment with fake quality values for these sequences which miss them. One may start with just putting 4s, 10s or anything in between ( and check if there are big differences in assembly.

More accurate approach (assuming the same sequencing technology/chemistry) would be to calculate average base quality at a given position from reads you have real values, and put these for reads without quality values.

You may also assembly reads with qualities first, then map the remaining ones to existing contigs. Use the mapping quality to estimate quality of the reads (lower the quality for mismatches), but this is probably an overkill.

ADD COMMENTlink written 8.7 years ago by Darked894.2k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1332 users visited in the last hour