SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
setting parameters for Novoalign 2-6 bp insertions deletions 12jrowley2 RNA Sequencing 2 11-01-2013 09:48 AM
Illumina insertions deletions andreitudor Bioinformatics 2 02-14-2012 07:20 AM
SOAP2 and deletions/insertions joa_ds Bioinformatics 7 09-21-2010 04:59 AM
How to find insertions in the genome? dina Bioinformatics 0 10-24-2009 11:55 AM
Sequencing random transgene insertions nathan hunkapiller Sample Prep / Library Generation 0 09-08-2009 10:22 AM

Reply
 
Thread Tools
Old 06-23-2010, 01:59 PM   #1
WindGirl2008
Junior Member
 
Location: Baltimore

Join Date: Dec 2008
Posts: 2
Default Looking for Aligner for mapping insertions

Hi, all,

Thanks for the help in advance.
I'm looking to alignment program that can handle deep sequencing reads with insertions. Eg. for some of the reads (100bp)
first 50bp is sequence is not from the ref genome, 50bp is sequence from the ref genome.
The break point is variable in each read.

I'm looking for something that can use the mappable sequences in each read to anchor the read to correct region of the genome. This way i'll be able to detect insertions (structural variants) in the sample comparing to the reference genome.

ELANDv2 does the job but it's not free...><

Any suggestions?

Thanks!
WindGirl2008 is offline   Reply With Quote
Old 06-23-2010, 03:47 PM   #2
mkeehan
Member
 
Location: Hamilton NZ

Join Date: Feb 2010
Posts: 13
Default

MosaikAligner is my favourite hammer.
http://bioinformatics.bc.edu/marthlab/Mosaik

But it seems that you want to do overlap based alignment scoring.
It might require a custom script and program...
mkeehan is offline   Reply With Quote
Old 06-23-2010, 03:53 PM   #3
WindGirl2008
Junior Member
 
Location: Baltimore

Join Date: Dec 2008
Posts: 2
Default

What do you mean by overlap based alignment scoring?
I can write some simple code.

What are the parameters you would recommend for running mosaik for these type of reads? Using the default, the chimera reads simply don't get aligned.

Thanks!!!
WindGirl2008 is offline   Reply With Quote
Old 06-23-2010, 06:07 PM   #4
lh3
Senior Member
 
Location: Boston

Join Date: Feb 2008
Posts: 693
Default

also try bwasw.

i think you need to tune mosaik a bit in order to get chimera, according to a post one or two years ago by the mosaik developer. things may have been changed, though.
lh3 is offline   Reply With Quote
Old 06-23-2010, 06:15 PM   #5
mkeehan
Member
 
Location: Hamilton NZ

Join Date: Feb 2010
Posts: 13
Default

I found this on "split-read" mapping

http://seqanswers.com/forums/archive...hp/t-2267.html
mkeehan is offline   Reply With Quote
Old 06-23-2010, 06:25 PM   #6
lh3
Senior Member
 
Location: Boston

Join Date: Feb 2008
Posts: 693
Default

Yes, this is what I was referring to. Thanks, mkeehan. I believe snownebula is the key developer of mosaik.

EDIT: bwasw was published early this year.

Last edited by lh3; 06-23-2010 at 06:32 PM.
lh3 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 10:59 AM.


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