This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How To Use Unmapped Reads To Reconstruct Strains Specific Genome Segments?

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?

denovo solid
• 3,575 views • link
updated by Philipp Bayer Philipp Bayer • written by William William

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.

• 0 views • link
written by Fabio Marroni Fabio Marroni

2 answers

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

• 1 views • link
written by Philipp Bayer Philipp Bayer

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

• 18 views • link
written by fridhackery • fridhackery

Log in to answer this question.

More posts like this