This is a test version of Biostars. For the public version, visit https://www.biostars.org.
mummer results merging

Hi, I splited my assembly into 200 parts (part1.fasta - part199.fasta) and ran nucmer -p test pseudochromosome-1-11.fasta part1.fasta. As result I got the following extensions:

  • *.delta
  • *.ntref
  • *.mgaps

How should I merge the 200 results into one alignment?

Thank you in advance.

mummer nucmer
• 2,326 views • link
written by Ric Ric

I think that is the solution: sed 1,2d BAFB.contigs.clean.*.delta > BAFB.contigs.clean.delta

• 0 views • link
written by Ric Ric

I got ERROR: Could not parse delta file, BAFB.contigs.clean.delta error no: 400 with show-coords BAFB.contigs.clean.delta > Scaffold.coords

• 0 views • link
written by Ric Ric

0 answers

No answers yet.

Log in to answer this question.

More posts like this
  • RNA-Seq expression analysis
    written by Ric Ric
    Hi, I ran successfully [nf-core/rnaseq][1] pipeline in the following way: ``` nextflow run nf-core/rnaseq -profile singularity -r 3.14.0 --input samplesheet.csv --fasta /work/waterhouse_team/NB/LAB360/NbLab360.genome.fasta --save_reference --gff /work/waterhouse_team/NB/LAB360/NbLab360.v103.gff.gz …
  • bcftools for plant SNPs
    written by Ric Ric
    Hi, I ran bcftools in the following way: bcftools mpileup -O b -o ${3}-${basename}.bcf -f $asm $bam bcftools call --ploidy 2 -m -o ${3}-${basename}.vcf ${3}-${basename}.bcf …
  • Merging contigs: Minimus2 error
    written by fiona.newberry fiona.newberry
    I am using minimus2 to merge overlapping contigs to improve my assembly. I am getting the following error from minimus2 in .runAmos.log CONSTRUCTIONTIME /opt/software/conda_env/opt/mummer-3.23/mummer PHE_1_CDHIT.ntref …
  • Retraining Augustus with BUSCO
    written by Ric Ric
    Hi, I tried to follow this [instruction][1] on how to retrain Augustus with BUSCO. I ran BUSCO: run_BUSCO.py -i NbV1Ch07.fasta -o augustus_busco -l /assembly/busco-embryophyta_odb9/ -m …
  • splitting a SAM file
    written by Ric Ric
    I wrote the below script which takes a SAM file as input and creates for each contig a SAM and FASTA file: #!/usr/bin/env python3 import …
  • Nucmer and delta-filter not symmetric for reference and query
    written by Fabio Marroni Fabio Marroni
    Dear all, I was surprised by the fact that running the following two sets of commands I obtained different results **Command 1:** nucmer --prefix=Ref_1_vs_Ref_2 Ref_1.fasta …
  • reducing memory consumption of circoletto
    written by Ric Ric
    Hi, 1. I splitted my query genome in 200 parts with pyfasta and I ran each of them separetly with blastall: **blastall -p blastn -i …
  • Aligning de-novo assembly to reference strain genome
    written by aindap aindap
    Dear Biostars Community - I generated a de-novo assembly of my viral strain. I wanted to analyze differences of my denovo result to a reference …
  • Mitochondrial denovo assembly
    written by mandar.bobade • mandar.bobade
    Dear all, I am trying to generate assembly of expected genome size (~800Kb) and have done lot of things so far. Following are the details …
  • Need A Public Canonical Dna Alignment Problem With Known Good Solution
    written by daniel.wilkerson • daniel.wilkerson
    <p>I am implementing a DNA alignment tool and I want to test it for correctness and performance against the best known tools. The basic alignment …