fake reads

Ask questions about sequence improvement / finishing D. mojavensis projects here.
Post Reply
chauser
Posts: 47
Joined: Wed Jan 31, 2007 3:34 pm

fake reads

Post by chauser » Wed Mar 10, 2010 2:21 pm

Chris,

We have a project that still has a gap DGA38L22 (we'll see if reads today close it).
We decided to try to make a fake read to either design primers from, or to close the gap itself.
-blasted 2 contigs vs grimshawii (flybase)
-identified grimshawii scaffold coordinates that span the gap
-pulled up the scaffold sequence (image below) (flybase)
-make fake reads

We actually managed to make fake reads which was pretty cool using the sequence identified and it matches reads in our assembly (good).

However, I believe n's in a scaffold assembly indicate a gap - is that right?
If so, are we working off of different assemblies?
Is this a valid approach?
The digests give us an estimate for gap size, and we are working based on these estimates.

Chuck


Image

Post Reply