Question: Assembly of amplicons
1
gravatar for Picasa
21 months ago by
Picasa350
Picasa350 wrote:

Hi,

I have a pcr amplicon data sequencing of 5 species and 2 genes markers.

My reads were sequenced with HiSeq; so I got 2 x PE of 125 bp.

After sorting my reads per individual, I would like now to assemble my reads in order to retrieve the 2 genes marker (length of 1000 bp aproximately)

1) Anybody have already done that kind of stuff ? Do you know a specific assembler for this job ?

2) The coverage should be different (due to pcr amplication), iam looking for an assembly who can give me the coverage for the contig.

Thanks for your help

pcr amplicon assembly • 987 views
ADD COMMENTlink modified 9 days ago by doctor.dee00550 • written 21 months ago by Picasa350

Try doing following steps: 1. use PEAR to merge forward and reverse reads. 2. Then use any denovo assembler to assemble reads. (If u know reference it would be nice) 3. Try to map to reference to validate (if neccessary).

ADD REPLYlink written 9 days ago by doctor.dee00550
1
gravatar for Brice Sarver
21 months ago by
Brice Sarver2.4k
United States
Brice Sarver2.4k wrote:

Try one of these:

  1. De novo assemble all of your reads. If your amplicons are in high coverage, you ought to be able to recover them. Any standard assembler will do (ABySS, Velvet, Newblet, SOAPdenovo, etc.)

  2. Since you know what you are looking for, you can map your reads to those targets. Extract the reads from the BAM file and assemble them de novo. Have a look at ARC.

ADD COMMENTlink written 21 months ago by Brice Sarver2.4k

Thanks Brice.

I've never head about ARC.

Did you have any experience on that ? in which way could it be "better" than the standard ones.

ADD REPLYlink written 21 months ago by Picasa350

I've used it for a few things, and it's great. I'm one of the authors on the manuscript :)

What do you mean by "standard ones?" Part of the motivation for developing it was that there weren't any approaches that did exactly what we needed.

ADD REPLYlink written 21 months ago by Brice Sarver2.4k
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: 1024 users visited in the last hour