Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Convert Maq's out put to fasta/fastq?

    Hi!

    I would like to know whether it's possible to convert Maq's out put (out.aln.txt file, after running maq mapview or the unmapped reads files, by incorporating the u option in maq map) back to a fasta/fastq format so that it can be re-used in Maq/ other assemblers.

    Thanks

  • #2
    Not totally sure, but I know that the Vancouver Short Read Analysis Package has a tool for converting Maq to BED (http://vancouvershortr.wiki.sourcefo...t/ConvertToBed) and that Aaron Quinlan's BEDTools have a script for generating BED --> fasta.

    Hope that helps! Definitely not the most direct way though for sure. :P

    Comment


    • #3
      R/BioConductor can do this easily, by importing the MAQ output into an AlignedRead object, creating a ShortReadQ object from that, and then exporting a fastq file.
      @1
      NNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNN
      +
      """"""""""""""""""""""""""""""""""""

      Comment


      • #4
        mapview to fastq:

        maq mapview maq_map_file |awk '{print "@"$1"\n"$15"\n+\n"$16}' > fastq_file

        mapview to fasta:

        maq mapview maq_map_file |awk '{print ">"$1"\n"$15}'> fasta_file

        unmapped reads file to fastq

        awk '{print "@"$1"\n"$3"\n+\n"$4}' unmapped_reads_file > fastq_file

        Comment


        • #5
          I've tried your suggetions. I'm able to create the fastq files with awk and convert that to bfq format required for MAQ. But running MAQ with the bfq files gives me a error:
          Assertion failed: (fp_bfq_l), function ma_match.cc, line 557.
          Any suggestions why this might happen?
          Last edited by bea; 06-04-2009, 12:50 AM.

          Comment


          • #6
            If your fastq name is too long, >35 characters, Maq will delete the front exceed characters, so at this situation, you can't get the actual primary fastq back.
            And the primary fastq paired-end reads, every reads in pe1.fastq file should be corresponding to the same line of pe2.fastq file, so you should sort the name of reads in mapview file at first.

            Comment

            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 on Modified Bases...
              Yesterday, 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
            39 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 04-10-2024, 10:19 PM
            0 responses
            41 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 04-10-2024, 09:21 AM
            0 responses
            35 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