SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Bioinformatics (http://seqanswers.com/forums/forumdisplay.php?f=18)
-   -   alignment and reference genome assembling (http://seqanswers.com/forums/showthread.php?t=10085)

slny 03-15-2011 11:57 AM

alignment and reference genome assembling
 
When I do the visualization of my sequence alignment data, I noticed that I could get the consensus. Then, what's the difference between sequence alignment and reference genome assembling?

Zigster 03-15-2011 01:29 PM

yes you have made the correct observation, they are basically one in the same except in the latter you walk away a new reference from your consensus

swbarnes2 03-15-2011 04:59 PM

Building your sequence de novo generally requires higher coverage, and you'll end up with lots of contigs that can't be bridged due to repetitive or non-complex sequence. Aligning to a very closely related genome will allow you to easily find discrepancies, even if coverage is not so great. If your coverage is, say, an even 5x, you could call a lot of SNPs. I've found plenty of real, sanger-verified SNPs with 3x coverage. But no de novo assembler will do anything with that.

Also, for a lot of applications, what you really want are the discrepancies between your sample and a reference, so aligning to that reference is a much easier way to find them, rather than building a de novo sample reference, and comparing that to something. Also, de novo won't be totally accurate if your sample is heterozygous, or a mix of samples; at least not Velvet, which is a popular one, and the one I'm most familiar with. Whereas with an aligner, it isn't hard to see that so many reads show one allele at a locus, and so many reads show another.

slny 03-16-2011 07:49 AM

Very detailed explanation. Thanks for all the responses.

I'm also using Velvet for sequence assembling. How to visualize the Velvet results? There are a lot of result files from Velvet assembling. How to deal with these files?


All times are GMT -8. The time now is 07:14 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2022, vBulletin Solutions, Inc.