The problems when i map to genome by mapper.pl in miRDeep2
1
0
Entering edit mode
3.1 years ago
51203903042 ▴ 20

When I did the first step with mirdeep2, which is to map the sequencing data to the genome, I encountered a problem:

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.

My sequencing data is like this,

>@A00808:451:HK3 TGDRXX:2 :2101:1380:1000 1:N:0:CCGAAGTA
CCGGGTGCTGTAGGCTT
>@A00808:451:HK3 TGDRXX:2 :2101:1398:1000 1:N:0:CCGAAGTA
TATGGCACTGGTAGAATTCACT
>@A00808:451:HK3 TGDRXX:2 :2101:2013:1000 1:N:0:CCGAAGTA
AAAAGCTGGGTTGAGAGGGCGA

Isn't this data in FASTA format? How can I solve it?

RNA-Seq • 702 views
ADD COMMENT
0
Entering edit mode

I believe mirdeep2 is expecting a FASTA file. What you've provided is not a FASTA file.

ADD REPLY
0
Entering edit mode

Since you said that your data is in fasta format I used the 101010 button in editor window to display your example properly. Otherwise > at beginning of lines is shown as quote by Biostars . If it does not exactly look like above then your data is NOT in fasta format.

ADD REPLY
0
Entering edit mode

I have examed the file ,there are indeed spaces in the middle, which caused this error, thank you everyone!

ADD REPLY
0
Entering edit mode

Please accept @swbarnes2's answer to provide closure to this thread.

ADD REPLY
0
Entering edit mode
3.1 years ago

Did you read the error? Are you absolutely sure that your file does not contain whitespaces?

ADD COMMENT

Login before adding your answer.

Traffic: 2444 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6