Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • what is wrong with my MAQ?

    [lemon@localhost Desktop]$ maq map flc_N.map flc.bfa /var/SEQ/Ningyou7_solexa.bfq
    -- maq-0.6.8
    [ma_load_reads] loading reads...
    [ma_load_reads] set length of the first read as 35.
    [ma_load_reads] 21952885*2 reads loaded.
    [ma_longread2read] encoding reads... 43905770 sequences processed.
    [ma_match] set the minimum insert size as 36.
    [match_core] round 1/3...
    [match_core] making index...
    [match_core] processing sequence sjh (1100 bp)...
    [match_core] round 2/3...
    [match_core] making index...
    [match_core] processing sequence sjh (1100 bp)...
    [match_core] round 3/3...
    [match_core] making index...
    [match_core] processing sequence sjh (1100 bp)...
    [match_core] sorting the hits and dumping the results...
    [ma_load_reads] loading reads...
    [ma_load_reads] 21952885*2 reads loaded.
    [mapping_count_single] 16, 16, 16, 16
    已杀死
    ??
    Last edited by lemon027; 09-11-2008, 05:17 PM. Reason: sloved

  • #2
    Maq man page will suggest you put 2 million (instead of 20 million) reads in a batch to run maq map. Probably the process gets killed due to insufficient memory.

    You can also use maq.pl easyrun if you works with small genomes. It will split the input for you.

    Comment


    • #3
      lemon:

      To follow up on lh3's comment, use the fastq2bfq command with the -n option (say -n 2000000, or 10000000) to automatically split the reads.

      Then run a map command on each file, then use mapmerge.

      Comment


      • #4
        my data is too large? but for anther maq, it had done well

        maq map /home/lemon/Desktop/flc_T.map /home/lemon/Desktop/flc.bfa Tapidor_solexa.bfq
        -- maq-0.6.8
        [ma_load_reads] loading reads...
        [ma_load_reads] set length of the first read as 35.
        [ma_load_reads] 19002103*2 reads loaded.
        [ma_longread2read] encoding reads... 38004206 sequences processed.
        [ma_match] set the minimum insert size as 36.
        [match_core] round 1/3...
        [match_core] making index...
        [match_core] processing sequence sjh (1100 bp)...
        [match_core] round 2/3...
        [match_core] making index...
        [match_core] processing sequence sjh (1100 bp)...
        [match_core] round 3/3...
        [match_core] making index...
        [match_core] processing sequence sjh (1100 bp)...
        [match_core] sorting the hits and dumping the results...
        [ma_load_reads] loading reads...
        [ma_load_reads] 19002103*2 reads loaded.
        [mapping_count_single] 16, 16, 16, 16
        [maq_indel_pe] the indel detector only works with short-insert mate-pair reads.
        [match_data2mapping] 62 out of 38004206 raw reads are mapped with 0 in pairs.
        -- (total, isPE, mapped, paired) = (19002103, 0, 62, 0)

        Comment


        • #5
          thank you I am trying your suggestion

          Comment


          • #6
            Yes ,
            your right~
            the error never come out again!
            but I had to merge the results together?
            can i use one line command to get the all the result?

            Comment


            • #7
              Yes, once you issue the maq mapmerge command (with all the proper filenames of course...see the man page)...you will have one .map file that is all of your data combined!

              Comment

              Latest Articles

              Collapse

              • seqadmin
                Strategies for Sequencing Challenging Samples
                by seqadmin


                Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
                03-22-2024, 06:39 AM
              • seqadmin
                Techniques and Challenges in Conservation Genomics
                by seqadmin



                The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

                Avian Conservation
                Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
                03-08-2024, 10:41 AM

              ad_right_rmr

              Collapse

              News

              Collapse

              Topics Statistics Last Post
              Started by seqadmin, Yesterday, 06:37 PM
              0 responses
              10 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, Yesterday, 06:07 PM
              0 responses
              9 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 03-22-2024, 10:03 AM
              0 responses
              50 views
              0 likes
              Last Post seqadmin  
              Started by seqadmin, 03-21-2024, 07:32 AM
              0 responses
              67 views
              0 likes
              Last Post seqadmin  
              Working...
              X