Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • mummer plot unable to parse alignment file

    Hello,

    I am currently trying to make a dotplot in order to compare the sequence similarity between a reference genome and an early stage reference genome that my own lab is constructing. The software I am using to do this is called "mummer".

    To get an initial feel for how the mummer software works, I aligned one chromosome of the reference genome to itself since this should in theory produce a perfect diagonal line.

    The code I used for this was:
    mummer Chromosome1.fasta Chromosome1.fasta > Chromosome1.mums

    Then to create a dotplot of this alignment file I used the following code:
    mummerplot -postscript -p OUTPUT_plot Chromosome1.mums

    This code created the perfectly diagonal line dot plot as predicted.

    My issue is that when I try to use the same code to align two different sequences, mummerplot is not able to parse the alignment file.

    Here is the code I am using:
    mummer Chromosome1.fasta Scaffold1.fasta >Chrom1_v_Scaff1.mums
    mummerplot -postscript -p OUTPUT_plot Chrom1_v_Scaff1.mums

    The following error appears every time I align two different sequences and try to plot the alignment:
    gnuplot 5.2 patchlevel 7
    Reading mummer file Chrom1_v_Scaff1.mums (use mummer -c)
    ERROR: Could not parse Chrom1_v_Scaff1.mums
    26869245 25 82

    I have tried all the other alignment programs in mummer including "-mum", "-maxmatch" and "-nucmer" and they all give the same error.

    I would really appreciate a solution to this issue if anyone has one, thank you.

  • #2
    mummrplot

    Hi,
    did you ever figure out a solution here? I'm running into the same issue.

    Comment

    Latest Articles

    Collapse

    • 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
    • seqadmin
      The Impact of AI in Genomic Medicine
      by seqadmin



      Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
      02-26-2024, 02:07 PM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, 03-14-2024, 06:13 AM
    0 responses
    32 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-08-2024, 08:03 AM
    0 responses
    72 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-07-2024, 08:13 AM
    0 responses
    80 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-06-2024, 09:51 AM
    0 responses
    68 views
    0 likes
    Last Post seqadmin  
    Working...
    X