Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
mosaik-aligner installation problem on redhat 6 joachim.jacob Bioinformatics 5 05-10-2013 08:36 AM
Ion Torrent PGM data with Mosaik-aligner Magnus Bioinformatics 0 05-19-2011 03:59 AM
Aligner/SNP caller for Illumina data sklages Bioinformatics 2 10-06-2010 04:48 AM
Powerful PC for MOSAIK aligner? neni Bioinformatics 4 04-30-2010 02:45 AM
Mosaik Aligner problem with SOLiD sequences g_solid Bioinformatics 0 02-07-2010 05:30 PM

Thread Tools
Old 12-22-2010, 01:29 AM   #1
Location: India

Join Date: Oct 2010
Posts: 59
Default Mosaik-aligner for SNP detection, Need Help.

Need help!

I am using Mosaik-aligner on ubuntu-9.10 for SNP detection, it operates without showing any errors. I have 454 sequencing data and i am using contigs formed by gsAssembler software that is available with 454 titanium system, formed from reads,

Steps for using Mosaik-aligner are,

./MosaikBuild (.fastq and ref.fasta)

./MosaikAlign (using parameters -hs 14 -mmp 0.05 -act 55 -bw 51 -mhp 100 for 454)


at the end i get *.ace file which i tried to open in Tablet, Gambit etc but it is not showing any alignments as it shows with the test data (c_elegans) given with mosaik. It shows how much sequences has aligned to me ref seq uniquly and non-uniquly but in .ace file no data except ref seq is seen.

*.bed and *.eland files shows contigs/reads aligned to my reference derived using MosaikText.

If any one knows the solution pls help me i am stucked here?

Whether i should use reads (from which contigs are formed) instead of contigs? (which i have did but the result is same, nothing seen)
ketan_bnf is offline   Reply With Quote
Old 12-22-2010, 05:39 AM   #2
Location: Montreal, Quebec, Canada

Join Date: Dec 2009
Posts: 30


when your are building the dat file specifying your input FASTQ?

you shouldn't see the c.elegans file in the ace file you're getting!

Do you see all your reads stats in the post build/assembly report?
ntremblay is offline   Reply With Quote
Old 12-22-2010, 09:17 PM   #3
Location: India

Join Date: Oct 2010
Posts: 59

yes, i am using FASTQ file for dat file formation.

yes, you are right i should not see c.elegans file in the ace file and i am not seeing it as well.

I am not merging c_elegans data with my data, i just meant to say processing of my data with mosaik-aligner dosen't give output as it gives with test data of c_elegans .

I am pasting my stats of post align/assembly report.

1. orf@ubuntu:~/mosaik-aligner/data$ ../bin/MosaikAligner -in sequence_archives/454AllContigs_normal.dat -out sequence_archives/454AllContigs_normal_chr5_aligned.dat -ia reference/chr5.dat -hs 13 -act 55 -mmp 0.05 -bw 51 -mhp 100

MosaikAlign gives output showing,

Alignment statistics (mates):
# filtered out: 3804 ( 92.9 %)
# unique: 123 ( 3.0 %)
# non-unique: 168 ( 4.1 %)
total: 4095
total aligned: 291 ( 7.1 %)

Miscellaneous statistics:
aligned mate bp: 95883
alignment candidates/s: 854.4

2. orf@ubuntu:~/mosaik-aligner/data$ ../bin/MosaikSort -in sequence_archives/454AllContigs_normal_chr5_aligned.dat -out sequence_archives/454AllContigs_normal_chr5_sorted.dat

MosaikSort gives output showing,

Single-end read statistics:
reads alignments
# non-unique: 168 (57.7 %) 20932 (99.4 %)
# unique: 123 (42.3 %) 123 ( 0.6 %)
total: 291 21055

[test@BIOTECH-PC-1 data]$ ../bin/MosaikText -in sequence_archives/454AllContigs_normal_chr5_sorted.dat -ref chr5 -bam bam/454AllContigs_normal_chr5.bam
MosaikText 1.0.1388 2010-02-01
Michael Stromberg Marth Lab, Boston College Biology Department

- converting the alignment archive to the following formats: BAM

Converting alignment archive:
100%[================================] 123.0 alignments/s in 1 s

[test@BIOTECH-PC-1 data]$ ../bin/MosaikAssembler
-in sequence_archives/454AllContigs_normal_chr5_sorted.dat -out assembly/chr5/454AllContigs_normal_chr5_sorted -ia reference/chr5.dat -f ace
./ line 1: #!/bin/sh: No such file or directory
MosaikAssembler 1.0.1388 2010-02-01
Michael Stromberg Marth Lab, Boston College Biology Department

alignment count reference sequence
123 chr5

Processing reference sequence chr5:
- inserting gaps into reference sequence... finished.
- creating ungapped to gapped conversion table... finished.
- writing assembly header... finished.
- locating first read for this reference sequence... finished.

- saving alignments from chr5:
100%[================================] 123.0 alignments/s in 1 s

- appending read data to header data... finished.

Whether i should use contigs or reads for aligning?
ketan_bnf 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 05:21 PM.

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