Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • mummerplot error

    Good day,

    Help me please with plotting mummer alignment.
    I've created out.mums file with command
    ../sparseMEM/mummer -mum -b reference.fa query.fa >out.mums.
    Both input fasta file are multi-fasta sequences, ~9.000 of records in each. Actually it is a draft assemblies.
    When I am trying to plot it using mummerplot,following the documentation, I always receive an error:

    " set ticscale 0 0
    ^
    "mummer.gp", line 13: warning: Deprecated syntax - please use 'set tics scale' keyword
    "
    Then I've tried to create small test *.mums file with usual mummer command
    ../mummer-3.22/mummer -mum -b -c testReference.fa testQuery.fa >testOut.mums

    But mummerplot could not plot *.mums file, it returns the same error.

    When I used nucmer for aligning, it has created *.delta file, and this delta file is plotting perfectly. But for my real files usual mummer and nucmer always crashes with segmentation fault error.

    I need to create dot-plot, so what you'll suggest me to do?

    Regards,
    Vlasova Anna
    ---------------------------------------------
    Bioinformatica i Genomica
    CRG-Centre for Genomic Regulation
    C/ Dr. Aiguader, 88
    Edif. PRBB
    08003 Barcelona, Spain

  • #2
    Hi,

    I was wondering if you ever got this problem solved. I am now having the same problem and when I searched online found your post. I will be grateful if you can share with me how you resolved this issue or what other alternatives did you find.

    Thanks for reading.

    Rama

    Comment


    • #3
      Originally posted by aleo View Post
      Good day,

      When I used nucmer for aligning, it has created *.delta file, and this delta file is plotting perfectly. But for my real files usual mummer and nucmer always crashes with segmentation fault error.

      I need to create dot-plot, so what you'll suggest me to do?

      Regards,
      Vlasova Anna
      ---------------------------------------------
      Bioinformatica i Genomica
      CRG-Centre for Genomic Regulation
      C/ Dr. Aiguader, 88
      Edif. PRBB
      08003 Barcelona, Spain
      Hi, can you please tell me how you plotted the nucmer.delta file using mummer please?? I mean which commands did you use. I am trying to map contigs to a reference sequence.

      THank you.

      Comment


      • #4
        Hello, everyone.
        Today, I followed the case to learn how to use mummer,

        In 4.1. Aligning two finished sequences
        The most basic use case is the alignment of two contiguous sequences. For all of the one vs. one use cases the mummer program alone, when coupled with mummerplot, may be all that is necessary to visualize a global alignment of the two sequences. This process alone can be very helpful in determining the large scale differences between the two sequences. For a single reference sequence ref.fasta and a single query sequence qry.fasta in FastA format, type:
        mummer -mum -b -c ref.fasta qry.fasta > ref_qry.mums

        till now, it seems every goes well, a file named ref_qry.mums was generated, but when I try to type the next commend,

        mummerplot --postscript --prefix=ref_qry ref_qry.mums

        an error appeared,
        defined(%hash) is deprecated at /usr/bin/mummerplot line 884.
        (Maybe you should just omit the defined()?)
        defined(%hash) is deprecated at /usr/bin/mummerplot line 894.
        (Maybe you should just omit the defined()?)
        defined(%hash) is deprecated at /usr/bin/mummerplot line 981.
        (Maybe you should just omit the defined()?)
        defined(%hash) is deprecated at /usr/bin/mummerplot line 991.
        (Maybe you should just omit the defined()?)
        defined(%hash) is deprecated at /usr/bin/mummerplot line 1034.
        (Maybe you should just omit the defined()?)
        defined(%hash) is deprecated at /usr/bin/mummerplot line 1044.
        (Maybe you should just omit the defined()?)
        gnuplot 4.4 patchlevel 3
        Reading mummer file ref_qry.mums (use mummer -c)
        ERROR: Could not parse ref_qry.mums
        10003036 5474049 20

        what's wrong? Is there anyone can tell me?

        Comment


        • #5
          Originally posted by louis_ly View Post
          Hello, everyone.
          Today, I followed the case to learn how to use mummer,

          In 4.1. Aligning two finished sequences
          The most basic use case is the alignment of two contiguous sequences. For all of the one vs. one use cases the mummer program alone, when coupled with mummerplot, may be all that is necessary to visualize a global alignment of the two sequences. This process alone can be very helpful in determining the large scale differences between the two sequences. For a single reference sequence ref.fasta and a single query sequence qry.fasta in FastA format, type:
          mummer -mum -b -c ref.fasta qry.fasta > ref_qry.mums

          till now, it seems every goes well, a file named ref_qry.mums was generated, but when I try to type the next commend,

          mummerplot --postscript --prefix=ref_qry ref_qry.mums

          an error appeared,
          defined(%hash) is deprecated at /usr/bin/mummerplot line 884.
          (Maybe you should just omit the defined()?)
          defined(%hash) is deprecated at /usr/bin/mummerplot line 894.
          (Maybe you should just omit the defined()?)
          defined(%hash) is deprecated at /usr/bin/mummerplot line 981.
          (Maybe you should just omit the defined()?)
          defined(%hash) is deprecated at /usr/bin/mummerplot line 991.
          (Maybe you should just omit the defined()?)
          defined(%hash) is deprecated at /usr/bin/mummerplot line 1034.
          (Maybe you should just omit the defined()?)
          defined(%hash) is deprecated at /usr/bin/mummerplot line 1044.
          (Maybe you should just omit the defined()?)
          gnuplot 4.4 patchlevel 3
          Reading mummer file ref_qry.mums (use mummer -c)
          ERROR: Could not parse ref_qry.mums
          10003036 5474049 20

          what's wrong? Is there anyone can tell me?
          Hi there, I think you are running the command wrongly. The command in the documentation says this :

          mummerplot -x "[0,275287]" -y "[0,265111]" -postscript -p mummer mummer.mums

          -x and - y are the coordinates for the dot plots, then use mummer for the .mums file and not prefix.

          I hope this helps

          Comment


          • #6
            Originally posted by louis_ly View Post
            mummerplot --postscript --prefix=ref_qry ref_qry.mums

            an error appeared,
            defined(%hash) is deprecated at /usr/bin/mummerplot line 884.
            (Maybe you should just omit the defined()?)
            defined(%hash) is deprecated at /usr/bin/mummerplot line 894.
            (Maybe you should just omit the defined()?)
            defined(%hash) is deprecated at /usr/bin/mummerplot line 981.
            (Maybe you should just omit the defined()?)
            defined(%hash) is deprecated at /usr/bin/mummerplot line 991.
            (Maybe you should just omit the defined()?)
            defined(%hash) is deprecated at /usr/bin/mummerplot line 1034.
            (Maybe you should just omit the defined()?)
            defined(%hash) is deprecated at /usr/bin/mummerplot line 1044.
            (Maybe you should just omit the defined()?)
            gnuplot 4.4 patchlevel 3
            Reading mummer file ref_qry.mums (use mummer -c)
            ERROR: Could not parse ref_qry.mums
            10003036 5474049 20

            what's wrong? Is there anyone can tell me?
            There are a couple of things wrong. The first set of messages are just Perl warnings you can safely ignore (I actually fixed this issue with my local version and hopefully someone will fix the source). The last line is the informative part. It says there was an error reading your mums file. This may or may not mean there is something wrong with your mummer output, but your command looks fine. Though, I'm not sure those coordinates look correct. You may want to re-run mummer and try mummerplot again.

            You can take at look at the mummerplot examples online and also read the command line documentation for how to call each program. Your mummerplot command looks okay also, so it's likely that something is not right with your mummer file, as the error message would indicate.

            Comment


            • #7
              mummerplot --postscript --prefix=ref_qry ref_qry.mums

              an error appeared,
              defined(%hash) is deprecated at /usr/bin/mummerplot line 884.
              (Maybe you should just omit the defined()?)
              defined(%hash) is deprecated at /usr/bin/mummerplot line 894.
              (Maybe you should just omit the defined()?)
              defined(%hash) is deprecated at /usr/bin/mummerplot line 981.
              (Maybe you should just omit the defined()?)
              defined(%hash) is deprecated at /usr/bin/mummerplot line 991.
              (Maybe you should just omit the defined()?)
              defined(%hash) is deprecated at /usr/bin/mummerplot line 1034.
              (Maybe you should just omit the defined()?)
              defined(%hash) is deprecated at /usr/bin/mummerplot line 1044.
              (Maybe you should just omit the defined()?)
              gnuplot 4.4 patchlevel 3
              Reading mummer file ref_qry.mums (use mummer -c)
              ERROR: Could not parse ref_qry.mums
              10003036 5474049 20
              There an error in mummerplot perl script.
              You have to correct two regular expression.
              in a terminal :
              $ vi path/to/mummer/directory/mummerplot
              go to line 565
              replace :
              if ( /^\s+(\d+)\s+(\d+)\s+(\d+)$/ ) {
              by :
              if ( /^\s*(\d+)\s+(\d+)\s+(\d+)$/ ) {

              do it again on line 575

              Quick explanation : it search for a line that begin with one or more space, but when you have a start superior to 9999999, this space do not exist anymore.
              So if you change '+' by '*' it search for 0 or more spaces.

              Comment


              • #8
                Originally posted by Beuss View Post
                There an error in mummerplot perl script.
                You have to correct two regular expression.
                in a terminal :
                $ vi path/to/mummer/directory/mummerplot
                go to line 565
                replace :
                if ( /^\s+(\d+)\s+(\d+)\s+(\d+)$/ ) {
                by :
                if ( /^\s*(\d+)\s+(\d+)\s+(\d+)$/ ) {

                do it again on line 575

                Quick explanation : it search for a line that begin with one or more space, but when you have a start superior to 9999999, this space do not exist anymore.
                So if you change '+' by '*' it search for 0 or more spaces.
                Thanks for the information, that is helpful. By the way, those other messages above the error are just warnings that can be ignored. What the message is saying is that the syntax
                Code:
                if (defined(%hash)) { ... }
                is deprecated. It also tells you exactly how to fix the issue and avoid the warnings.

                It should be replaced with
                Code:
                if (%hash) { ... }
                That is probably clear though.

                Comment


                • #9
                  Thank Beuss.

                  it works!

                  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
                  8 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, Yesterday, 06:07 PM
                  0 responses
                  8 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 03-22-2024, 10:03 AM
                  0 responses
                  49 views
                  0 likes
                  Last Post seqadmin  
                  Started by seqadmin, 03-21-2024, 07:32 AM
                  0 responses
                  66 views
                  0 likes
                  Last Post seqadmin  
                  Working...
                  X