SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
gigaBayes output format ... inconsistent with manual? jnfass Bioinformatics 3 12-28-2013 10:10 PM
The best output format for paired-end alignment using bowtie chedonat Bioinformatics 4 02-01-2012 02:03 PM
Strange/erroneous TopHat output with paired-end data hoho Bioinformatics 7 03-09-2011 01:21 PM
Bowtie output from paired end reads godzilla07 Bioinformatics 0 01-06-2011 12:36 PM
gigabayes no output shaohua.fan Bioinformatics 7 10-16-2010 10:59 AM

Reply
 
Thread Tools
Old 08-24-2010, 09:06 AM   #1
Gianza
Junior Member
 
Location: Italy

Join Date: Jan 2010
Posts: 7
Unhappy Mosaik, gigaBayes and paired-end output

Hello everybody,
I'm facing a weird problem, really hope that somebody can help me.

I have aligned both 454 and Solexa paired-ends reads on a transcriptome reference for SNP mining purposes.

The two alignments have been merged with MosaikMerge, and MosaikAssembler run fine to build contigs. The MosaikAssembler module itself provide a flag to output reads and qualities (FASTA) which can be used together with ACE directly for gigabayes.
(Just note that all single fasta,qual and ace from each contig have been concatenated, and ace modified accordingly to its format).
Everything was PERFECT with only 454 alignment..confirming that I'm not messing up with concatenation.

BUT, using also paired-ends from Illumina, both ACE and FASTA file in output loose their /1 and /2 tags on read names which identify mates.
Thus, when running gigaBayes, it get confused by this, looking for one mates in ACE, finding another in FASTA. Indeed, it stop as soon unpadded read lengths don't match.

It seems absurd to me that the developer of both tools (Marth Lab) didn't considerate this aspect. Therefore I'm wishing to find my stupid-dumb-error.

the two fastq files as input for Mosaik are as following,

mates_1.fastq

@Sample1-0110:5:1:1002:3520#TGACCA/1
GCAAGACAGCAAGTGACAAATCCCCCATCGATAAATTGTGGGTCGTCTAGTGATC
+Sample1-0110:5:1:1002:3520#TGACCA/1
)-58806)*))/1111+757788707733555AAAAA8776523535AAAAA###
@sample1-0110:5:1:1003:10121#TGACCA/1
AACAATGTCATCGATTTCGCTATCGTATCTACAGAGCGCATGCTCCTATTCAAATCTTGGACCACAAAACAAAAGAG
+Sample1-0110:5:1:1003:10121#TGACCA/1
+*.1AAAAAAAAAAAAAAAAAAAAAAAAA7::7:AAAAAAAAAAAAAAAAAAA7AAAAA787:8AAAAAAA2772AA


mates_2.fastq

@Sample1-0110:5:1:1002:3520#TGACCA/2
AGGCTGCGACGAACAGGGCTCGACAAAGCATTCGAGAATATATAGTTCTTCGACGCTGCTTGATCACTAGACGACC
+Sample1-0110:5:1:1002:3520#TGACCA/2
A?A::B=BB??;?<?BB?+8?BB;BAAA.A?8+??B?BA?BABBABBBBB*?;<?AB?B95??*5AB>BBBB>98A
@Sample1-0110:5:1:1003:10121#TGACCA/2
GTGTTTAAACCAAAAGGGCTCATTGGTTTTGATTTGGAAACATTTGTGTTGTTCTCTTTTGTTTTGTGGTCCAAGA
+Sample1-0110:5:1:1003:10121#TGACCA/2

etc...

I will really appreciate any help! I'm not a real bioinformatician, thus almost impossible for me to work around this.
Gianza 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 02:42 AM.


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