Question: How To Use Unmapped Reads To Reconstruct Strains Specific Genome Segments?
gravatar for William
5.8 years ago by
William4.4k wrote:

Are there any common approaches to make use of unmapped reads to reconstruct strain specific genome segments?

Are these just your normal denovo aproaches do you somehow also take the mapped reads and the reference in to account? Maybe for example read pairs in which one reads maps and not the other? Or reference assisted assembly?

And is this at all possible with Solid (fragment and PE ) sequencing data? A denovo assembly or reference assisted assembly?

solid denovo • 1.8k views
ADD COMMENTlink modified 3.1 years ago by Philipp Bayer6.0k • written 5.8 years ago by William4.4k

Both your approaches are reasonable. If you have a relatively small genome, enough CPUs and RAM, and good libraries, it's easier to perform a complete denovo assembly and then identify additional segments in your strain.

ADD REPLYlink written 5.8 years ago by Fabio Marroni2.1k
gravatar for Philipp Bayer
3.1 years ago by
Philipp Bayer6.0k
Philipp Bayer6.0k wrote:

There was a recent paper which did this with a ton of rice accessions:

Exploring the rice dispensable genome using a metagenome-like assembly strategy

We proposed an assembly strategy incorporating SOAPdenovo [21], Meta-IDBA [22], and GICL [23] to assemble the dispensable sequences (“Materials and methods”; Fig. 1a). Both reads of a read pair of which at least one read could not be mapped to the reference genome were extracted to do de novo assembly (“Materials and methods”). In total, 395,891,829 read pairs (7.0 % of all the read pairs) were collected, of which 628,163,124 reads (79.3 %) could not be mapped to the Nipponbare genome

ADD COMMENTlink written 3.1 years ago by Philipp Bayer6.0k
gravatar for fridhackery
5.8 years ago by
fridhackery140 wrote:

I don't know if its common, yet, but certianly doable: A: Collect read pairs where at least one one read is mapped

ADD COMMENTlink written 5.8 years ago by fridhackery140
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: 1070 users visited in the last hour