Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • 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.

Latest Articles

Collapse

  • seqadmin
    Essential Discoveries and Tools in Epitranscriptomics
    by seqadmin




    The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist...
    04-22-2024, 07:01 AM
  • seqadmin
    Current Approaches to Protein Sequencing
    by seqadmin


    Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...
    04-04-2024, 04:25 PM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 04-11-2024, 12:08 PM
0 responses
59 views
0 likes
Last Post seqadmin  
Started by seqadmin, 04-10-2024, 10:19 PM
0 responses
57 views
0 likes
Last Post seqadmin  
Started by seqadmin, 04-10-2024, 09:21 AM
0 responses
51 views
0 likes
Last Post seqadmin  
Started by seqadmin, 04-04-2024, 09:00 AM
0 responses
55 views
0 likes
Last Post seqadmin  
Working...
X