How To Use Unmapped Reads To Reconstruct Strains Specific Genome Segments?
2
2
Entering edit mode
8.3 years ago
William ★ 4.9k

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 • 2.3k views
ADD COMMENT
0
Entering edit mode

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 REPLY
3
Entering edit mode
5.6 years ago

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 COMMENT
1
Entering edit mode
8.3 years ago
fridhackery ▴ 150

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

ADD COMMENT

Login before adding your answer.

Traffic: 1686 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6