Question: Using Phrap Assembly For Recombinining Sequences
gravatar for hadasa
9.3 years ago by
hadasa1.0k wrote:

What are the best phrap parameters for assembling sequence reads(forward and reverse) in cloned parasite isolates that are known to undergo recombination?

assembly recombination • 2.0k views
ADD COMMENTlink written 9.3 years ago by hadasa1.0k

What parasite and what level of recombination?

ADD REPLYlink written 9.3 years ago by Jarretinha3.3k

I agree with Jarretinha, knowing the scale of recombination would be important to know

ADD REPLYlink written 9.3 years ago by Phis1.0k

the parasite is Plasmodium falciparum. i do not know exactly what the full scale or recombination is. but let us say high recombination rates among the surface antigens.

ADD REPLYlink written 9.3 years ago by hadasa1.0k
gravatar for Lee Katz
8.8 years ago by
Lee Katz2.9k
Atlanta, GA
Lee Katz2.9k wrote:

I'm not sure if recombination rates have anything to do with it. You've sequenced a locus with two reads (F/R) and want to assemble those reads. Whatever level of recombination, those reads are true data. Just assemble with normal parameters.

phred *.ab1 -sd . -qd . -pd .
phd2fasta -os locus.fna -oq locus.fna.qual
phrap locus.fna -ace 2>phrap.err 1>phrap.out
# output is in locus.fna.contigs and locus.fna.contigs.qual

If these forward and reverse sequences do not assemble, consider the possibility that there is some physical obstruction to their overlap. Maybe a recombination event created an insertion between the forward and reverse sequences, such that you need to sequence the insert before you can assemble all sequences. Or, your forward or reverse primer locus was recombined to a different place in the genome. To overcome this obstruction, you may need to primer walk.

ADD COMMENTlink modified 10 months ago by RamRS22k • written 8.8 years ago by Lee Katz2.9k
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: 1571 users visited in the last hour