SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Mosaik Alignment Problems rahul.m.dhodapkar Bioinformatics 2 09-14-2010 09:26 AM
Mosaik and gigabuild help alessandra85 Bioinformatics 0 09-03-2010 09:17 AM
Problems visualizing Mosaik alignment with Gambit josecolquitt Bioinformatics 1 05-25-2010 09:14 AM
New version of MOSAIK available snownebula Bioinformatics 9 02-02-2010 02:08 AM
Mosaik alignment chienchi Bioinformatics 0 01-12-2010 01:21 PM

Reply
 
Thread Tools
Old 11-09-2009, 08:22 AM   #1
ernesto
Junior Member
 
Location: uk

Join Date: Oct 2008
Posts: 3
Default Mosaik alignment

Dear All,

I am using MosaikAligner with -m all option. This option map a read into all possible locations (each location with a certain alignment quality). Is there a way to get the alignment with the highest score and to discard the rest of the alignments?

Thanks!!!

Ernesto
ernesto is offline   Reply With Quote
Old 11-10-2009, 10:24 AM   #2
amcrisan
Junior Member
 
Location: Canada

Join Date: Oct 2009
Posts: 9
Default

That's difficult to say. The -m unique option will avoid giving you all the possible reads like the -m all option. However, all it does (according to the documentation) is really set aside any reads which can map to multiple locations, it doesn't take the highest scoring alignment out of all the possible alignment locations for that read. Additionally I don't see any indication that it actually would.

The best you can do in my opinion is to use the -m unique option, and then later you can use MosaikSort and tell it to try resolve multiple vs multiple (-rmm) read pairs, and so on. That being said, you may still get things discarded, as it may not be possible to resolve the multiple vs. multiple.

Lastly, using the -m unique option is faster and takes less space. You could also use the -m all option followed by MosaikSort for I think similar outcomes to using the unique option (not 100% sure though..).

Hope that's helpful.
amcrisan is offline   Reply With Quote
Old 11-12-2009, 10:51 PM   #3
zhanxw
Junior Member
 
Location: umich

Join Date: Jun 2009
Posts: 7
Default

From the documentation, -m unique, only place uniquely mapped read.
Does it mean Mosaik will discard the reads that can aligned at multiple locations?
zhanxw is offline   Reply With Quote
Old 11-13-2009, 11:20 AM   #4
amcrisan
Junior Member
 
Location: Canada

Join Date: Oct 2009
Posts: 9
Default

As I understand it, the non-unique reads will not necessarily be discarded. They will be placed asides and may be resolved at later steps using Mosaik sort (the -rmm option). If you don't choose to do this, then I think the reads are essentially discarded.
amcrisan is offline   Reply With Quote
Old 11-13-2009, 01:12 PM   #5
ernesto
Junior Member
 
Location: uk

Join Date: Oct 2008
Posts: 3
Default

Thanks for you replies, but as far as I understand the -rmm option in MosaikSort is for paired-end reads, and my data is single ended.
Do you think this option will work with single-end reads as well?

Thanks,

Ernesto
ernesto is offline   Reply With Quote
Old 02-02-2010, 03:08 AM   #6
shahid.manzoor
Junior Member
 
Location: uppsala, sweden

Join Date: Jun 2009
Posts: 8
Default

hi all I m using Mosaik and create jump database with Mosaikjumb command which is successfully create the jump database but after that when i run MosaikAligner command it give error like this ERROR: A position (1099511627520) was specified that is larger than the jump positions database (113492972).
So any body can help me to about this error.
shahid.manzoor is offline   Reply With Quote
Old 07-25-2010, 06:41 PM   #7
donniemarco
Member
 
Location: USA

Join Date: Aug 2009
Posts: 17
Default same error

hi everyone, i got the same error:
***************
hi all I m using Mosaik and create jump database with Mosaikjumb command which is successfully create the jump database but after that when i run MosaikAligner command it give error like this ERROR: A position (1099511627520) was specified that is larger than the jump positions database (113492972).
***********

Has anyone got an idea about it. I have submitted this error with mosaik error tracking too.
My parameters are:
-hs 15 -mmp .05 -m all -mhp 100 -act 26 -p 8 -j jumpdb.

Thanks!!!
donniemarco is offline   Reply With Quote
Old 07-26-2010, 06:13 PM   #8
sdriscoll
I like code
 
Location: San Diego, CA, USA

Join Date: Sep 2009
Posts: 423
Default

Quote:
Originally Posted by ernesto View Post
Dear All,

I am using MosaikAligner with -m all option. This option map a read into all possible locations (each location with a certain alignment quality). Is there a way to get the alignment with the highest score and to discard the rest of the alignments?

Thanks!!!

Ernesto
I think most of the time the "highest scoring alignment" is going to be a little ambiguous. If you're aligning short reads ( < 50 bp ) then you're gonna get reads that align in lots of places with the same number of mismatches and virtually no way to say that one is better than another.

I'm not sure what's out there that could do this for you but what might work is to let Mosaik just go for it (-m all) and then to do some post sorting on the output alignment file by binning read IDs that show up in multiple rows and sort their alignment scores. I still think what you'll end up with is a list of alignment locations that have the same score instead of single locations with the highest score as you're looking for. You'd want to look at # of mismatches as well as the alignment scores.
sdriscoll is offline   Reply With Quote
Reply

Tags
mosaik

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


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