Header Leaderboard Ad

Collapse

Help with using maqview please

Collapse

Announcement

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

  • Help with using maqview please

    I am trying to view the results of a maq analysis from the alignment file output XXX.map.

    I am having trouble understanding the README file that comes with the maqview package.

    I have successfully indexed the file using the maqindex command but I do not understand the instructions for the ./zrio command, or what it is doing?

    Can someone help me understand what to do next please? Thanks.

  • #2
    have you converted the XXX.map file to a text file from this binary format?

    Comment


    • #3
      Yes, I have used the command maq mapview out.map to extract to a text file XXXX.map. I have no idea what to do next. I am quite new to using Linux and the command line so if anyone has used maqview before, please could you let me know what to do next? Thanks.

      Comment


      • #4
        Assuming you already did the following (or similar)
        Code:
        ./maq-0.7.1/maq fastq2bfq /tmp/s_5_1_fastQ.txt /tmp/s_5_1.bfq
        ./maq-0.7.1/maq fastq2bfq /tmp/s_5_2_fastQ.txt /tmp/s_5_2.bfq
        
        ./maq-0.7.1/maq fasta2bfa ./genomes/Enterobacteria_phage_phiX174/NC001422.fna /tmp/phiX.bfa
        
        ./maq-0.7.1/maq map -1 76 -2 76 -u /tmp/unusedReads.txt -N /tmp/out.aln.map /tmp/phiX.bfa /tmp/s_5_1.bfq /tmp/s_5_2.bfq
        
        ./maq-0.7.1/maq assemble /tmp/out.cns /tmp/phiX.bfa /tmp/out.aln.map
        Then you could run the following:
        Code:
        ./maqview-0.2.5/maqindex -i -c /tmp/out.cns /tmp/out.aln.map
        
        ./maqview-0.2.5/maqview -c /tmp/out.cns /tmp/out.aln.map
        This will open the viewer-window, (though in this example VERY slow, as the PhiX coverage here is ~125000x, the maqview code is not build for such a coverage)

        Cheers
        -Jonathan

        Comment


        • #5
          Thank you very much for your help, it is working fine now . The problem was in fact with my maq analysis, not with my understanding of maqview.

          Using the instructions for using maq on the project website the command "maq match" is used to align reads to the reference. I think this command must be obsolete because it produces an output alignment file that does not work with maqview.

          The correct command you gave me to map the reads to the reference is "maq map" which I used instead and it worked perfectly. I should have noticed the "maq match" command is missing from the list of commands given when you type "maq" on the command line.

          Comment


          • #6
            Yes, that is correct.
            Glad it's working now!

            Best
            -Jonathan

            Comment


            • #7
              Maqview crashes with mac

              Hi,

              I have installed Maqview 0.2.5 on my iMac and found it crashing consistently. I correcly indexed my .map file (maqindex -i XXX.map) and then run the maqview command (maqview XXX.map), I can see the main window, but after a couple of clicks it inevitably crashes complaining about a "Bus error".
              Can anyone help me with this?
              Thanks

              Comment

              Latest Articles

              Collapse

              • seqadmin
                A Brief Overview and Common Challenges in Single-cell Sequencing Analysis
                by seqadmin


                ​​​​​​The introduction of single-cell sequencing has advanced the ability to study cell-to-cell heterogeneity. Its use has improved our understanding of somatic mutations1, cell lineages2, cellular diversity and regulation3, and development in multicellular organisms4. Single-cell sequencing encompasses hundreds of techniques with different approaches to studying the genomes, transcriptomes, epigenomes, and other omics of individual cells. The analysis of single-cell sequencing data i...

                01-24-2023, 01:19 PM
              • seqadmin
                Introduction to Single-Cell Sequencing
                by seqadmin
                Single-cell sequencing is a technique used to investigate the genome, transcriptome, epigenome, and other omics of individual cells using high-throughput sequencing. This technology has provided many scientific breakthroughs and continues to be applied across many fields, including microbiology, oncology, immunology, neurobiology, precision medicine, and stem cell research.

                The advancement of single-cell sequencing began in 2009 when Tang et al. investigated the single-cell transcriptomes
                ...
                01-09-2023, 03:10 PM
              • seqadmin
                AVITI from Element Biosciences: Latest Sequencing Technologies—Part 6
                by seqadmin
                Element Biosciences made its sequencing market debut this year when it released AVITI, its first sequencer. The AVITI System uses avidity sequencing, a novel sequencing chemistry that delivers higher quality data, decreases cycle times, and requires lower reagent concentrations. This new instrument reportedly features lower operating and start-up costs while maintaining quality sequencing.

                Read type and length
                AVITI is a short-read benchtop sequencer that also offers an innovative...
                12-29-2022, 10:44 AM

              ad_right_rmr

              Collapse
              Working...
              X