Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
PubMed: Rapid hybrid de novo assembly of a microbial genome using only short reads: C Newsbot! Literature Watch 0 10-19-2011 11:40 PM
4-dimensional hybrid cckandmd Sample Prep / Library Generation 0 04-01-2011 08:04 PM
Discussion about MIRA hybrid assembly of 454 reads with Illumina unpaired data edge De novo discovery 5 11-16-2009 01:17 AM
de novo hybrid assembly glacerda Bioinformatics 6 08-06-2008 09:52 PM
Hybrid selection dismith Events / Conferences 0 08-03-2008 01:07 PM

Thread Tools
Old 06-07-2010, 12:31 PM   #1
Junior Member
Location: baltimore

Join Date: Feb 2009
Posts: 9
Default Hybrid reads


I want to align nextgen seq reads that should only partially align to one reference seq. I'm most interested in the part of the read that doesnt align to the reference, but need the aligned portion to filter for it. What aligner would work best? BFAST? I'm working on a Mac OSX to start.

Thanks much!
veena is offline   Reply With Quote
Old 06-07-2010, 01:34 PM   #2
Senior Member
Location: San Diego

Join Date: May 2008
Posts: 912

It's very old fashioned, and probably too slow for mammalian genomes, but SOAP has an option to iteratively trim from the 3' end until the remnent aligns well enough.

From the SOAP output, you could get the names of reads which had to be severely trimmed to ailgn, and how far they had to be trimmed. Then pull out those reads, and look at the sequence after the aligned portion.
swbarnes2 is offline   Reply With Quote
Old 06-07-2010, 09:36 PM   #3
Senior Member
Location: Australia

Join Date: Sep 2008
Posts: 136

Another option is to use Bowtie, and iteratively map data that doesnt align clipping the 3' or 5' end of the read (whichever end you want to remove) a few times... hopefully bowtie is fast enough that mapping over and over isnt too much of a problem
frozenlyse is offline   Reply With Quote

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 07:27 PM.

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