This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Problems With Mirdeep2

I want to find miRNAs on a fastq file with mirDeep2. After adapter clipping and alignment (with mapper.pl from mirDeep2) against human_g1k_v37 (indexed with bowtie) I get following error message:

Error: Genome file /PathToFile/human_g1k_v37.fasta has not allowed whitespaces in its first identifier

The first line of the reference file looks like this:

>1 dna:chromosome chromosome:GRCh37:1:1:249250621:1

If I delete all white spaces in the identifiers, and rerun the mapping and mirDeep2, I get an other error message:

The mapped reference id 1 from file exmaple_collapsed.arf is not an id of the genome file /PathToFile/human_g1k_v37.fasta

Any idea how i have to convert the reference file? Thank you!

hg19 reference
• 9,806 views • link
updated by gzl • gzl • written by pmuench • pmuench

1 answer

You need to rename the sequences in your reference fasta, just remove everything after the space:

 perl -plane 's/\s+.+$//' < genome.fa > new_genome.fa

You will need to create a new index.

• 0 views • link
written by JC JC

Hi,

I know this post has been out for a while but I'd appreciate some help. I tried a few data sets on mirdeep2 and I always got this error:

First line of FASTA reads file is not in accordance with the fasta format specifications
Please make sure your file is in accordance with the fasta format specifications and does not contain whitespace in IDs or sequences

***** Please check if the option you used (options c) designates the correct format of the supplied reads file earthwormShort1.fa *****

I have tried to remove whitespaces several times to ensure my fasta file is OK but the error keeps returning. Please, kindly help.

• 0 views • link
updated by Ram Ram • written by Gabe Anderson • Gabe Anderson

Log in to answer this question.

More posts like this
  • mirdeep2 mapped reference error
    written by jomagrax • jomagrax
    Hi everyone! after succesfully indexed the reference genome with bowti-build and mapped with mapper.pl mapper.pl /home/superuser/jose/TFG/rawdata/control_1_trimmed.fastq -e -h -i -j -m -p /home/superuser/CAP-miRSEQ/ref/pig_genome.fa -s pig_control_1_collapsed.fa …
  • allow mismatch in miRDeep2
    written by anthihongthuy • anthihongthuy
    I used mirDeep2 for my analysis with the following command lines: **Method 1:** 1. mapper.pl (without option q) <code bash> mapper.pl trimmed_data.fq -e -h -i …
  • samtools flagstat returning 0% align reads
    written by Bioinformatist Newbie Bioinformatist Newbie
    Hello, I am using PASH 3.0 for aligning the .fastq file against GRCh37 genome (Using PASH is mandatory here). I am using the following command: …
  • mirdeep2 step-by-step tutorial
    written by robm9119 robm9119
    Hi folks, Does anyone know of any tutorial for running mirdeep2? I tried mirdeep2's own manual: https://www.mdc-berlin.de/36105849/en/research/research_teams/systems_biology_of_gene_regulatory_elements/projects/miRDeep/documentation but couldn't make sense of the directions. If …
  • how to setup miRDeep2 on osx 10.10 (Yosemite)
    written by khblack • khblack
    Hello I want use miRDeep2 on my mac (Yoseimite). So I downloaded miRDeep2 [here][1] [1]: https://www.mdc-berlin.de/8551903/en/ I open iterm and change dir to Download/mirdeep2 type …
  • Reading colon delimited description from Fasta file possible?
    written by bfeeny • bfeeny
    I am working with human_g1k_v37.fasta which is found on the 1000genomes site, specifically: ftp://ftp.1000genomes.ebi.ac.uk/vol1/ftp/technical/reference/ I have parsed this into files of individual chromosomes. The header …
  • Mapping statistics error from mirdeep2
    written by Jordan Jordan
    Hi, I'm using mirdeep2 for finding miRNAs on mouse data. I ran mapper.pl as a first step. It runs well and give me output. But …
  • Iupac Bases In The Hg19 1000G Reference ?
    written by Gabriel R. Gabriel R.
    <p>If you download the reference from 1000g :</p> <p><a href="ftp://ftp.1000genomes.ebi.ac.uk/vol1/ftp/technical/reference/human_g1k_v37.fasta.gz">ftp://ftp.1000genomes.ebi.ac.uk/vol1/ftp/technical/reference/human_g1k_v37.fasta.gz</a></p> <p>Why is there a double "RR" in chr 3 ?</p> <pre><code>zcat human_g1k_v37.fasta.gz.1 | grep -n …
  • How Do I Interpret The Output Generated By Variationhunter?
    written by Bioscientist Bioscientist
    <p>I'm using Variationhunter to analyze CNV; and when I used hg18 reference genome sometime ago, I usually got the output as:</p> <pre><code>IL5_286:1:311:10:40 chr1 142188880 142188915 …
  • 1000 Genomes Human Reference Assembly File With Non Actgn Chars In Chr 3, Why?
    written by Pablo Marin-Garcia Pablo Marin-Garcia
    <p>I am using the 1000 genomes human assembly reference file </p> <p>wget <a href="ftp://ftp-trace.ncbi.nih.gov/1000genomes/ftp/technical/reference/human_g1k_v37.fasta.gz">ftp://ftp-trace.ncbi.nih.gov/1000genomes/ftp/technical/reference/human_g1k_v37.fasta.gz</a></p> <p>and It was working OK with bwa but not with other …