SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Assembling against a reference Tools pari_89 Bioinformatics 8 05-12-2014 04:15 AM
D. melanogaster known SNPs and indels id0 Bioinformatics 6 08-26-2013 12:25 PM
Assembling against a reference pari_89 Bioinformatics 1 06-27-2013 07:24 AM
find all snps/indels prbndr Bioinformatics 2 09-20-2011 10:43 AM
454 assembling against reference genome donniemarco 454 Pyrosequencing 2 08-17-2009 06:22 AM

Reply
 
Thread Tools
Old 06-16-2014, 06:26 AM   #1
sdmoore
Member
 
Location: Florida

Join Date: Jun 2014
Posts: 11
Default Won't assembling reads to a reference remove SNPs and InDels?

When aligning reads to a reference to generate a config for SNP/InDel detection, won't the alignment necessarily remove the reads that have the desired information?

I see this as a commonly-listed step before searching for mutations.

Other threads I have read here (and papers referred to in them) refer to software that does the alignments (I have used Bowtie and Ugene), but I clearly see the absence in the resulting contigs of known antibiotic resistance markers that are in the genomes (large insertions). My filtered Illumina sets are paired (which don't seem to work in the alignments), but the individual libraries provide about 10 fold coverage with most reads >200.
sdmoore is offline   Reply With Quote
Old 06-16-2014, 07:04 AM   #2
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

Most of the pipelines aren't intended for finding very large insertions of novel sequence. For those, just use a different method, such as this one.
dpryan is offline   Reply With Quote
Old 06-16-2014, 07:32 AM   #3
sdmoore
Member
 
Location: Florida

Join Date: Jun 2014
Posts: 11
Default

Thanks for that paper, I had it in my folder, so I must have at least skimmed it.
I have not gotten decent assemblies using the paired end libraries (processed), I think they may all have to be the same length, which I didn't realize when I was filtering/trimming each.

The assemblers I have used don't spit out the orphans, which is exactly what I need. I will find a way to get those.

Thanks for the help.
sdmoore is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 01:13 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO