This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Mapsplice Reads Format Error

I'm trying to use mapsplice,I used an SRA database file ERR015688.fastq , against hg19 reference genome and bowtie index into the same directory

my command:

$ python bin/mapsplice_segments.py -Q fq -u /storage/gwajnberg/ERR015688/ERR015688.fastq -c /storage/data/gabriel/ucsc/hg19/ -B /storage/data/gabriel/ucsc/hg19/hg19 -o /storage/gwajnberg/mapsplice_out/

what happened:

[Wed Aug 31 16:07:31 2011] Preparing output location /storage/gwajnberg/mapsplice_out//


[Wed Aug 31 16:07:32 2011] Beginning Mapsplice run (v1.15.2)

-----------------------------------------------

bin directory: [/storage/app/MapSplice_1.15.2/bin/] 

[Wed Aug 31 16:07:32 2011] Checking for chromosomes files or directory

base name of chromosome file not consistent with chromosome name

base chromosome file name:hg19

chromosome name:chr10

[Wed Aug 31 16:07:32 2011] format chromosome files

[Wed Aug 31 16:09:43 2011] Checking for Bowtie index files

[Wed Aug 31 16:09:43 2011] Indexing chromosome sequences

[Wed Aug 31 19:09:30 2011] check reads format

read name contain blank space or tab

the 1th read in /storage/gwajnberg/ERR015688/ERR015688.fastq

@ERR015688.1 F2JR09301DGOLY length=270

    [FAILED]

Error: check reads format failed

Does anyone know what happened??? thanks

read error fastq
• 5,234 views • link
updated by Istvan Albert Istvan Albert • written by Gabrielw • Gabrielw

can you add the the result of "head -n 20 /storage/gwajnberg/ERR015688/ERR015688.fastq" ?

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

I 've discovered what happened example: @ERR015688.1 F2JR09301DGOLY length=270 It doesn t accept blank spaces...so I had to remove them...then it worked... now the program stopped at another point!

• 0 views • link
written by Gabrielw • Gabrielw

I had the same error... you can use ´sed 's/ /_/g' ERR015688.fastq´

• 0 views • link
written by Geparada Geparada

I had the same error... and I use: sed 's/ /_/g' reads.fastq and then works nice.

• 0 views • link
written by Geparada Geparada

1 answer

from the error message it sounds like the sequence for each chromosome needs to be in a separate file that is named after the chromosome.

See the manual: http://www.netlab.uky.edu/p/bioinfo/MapSpliceManual

• 41 views • link
written by Istvan Albert Istvan Albert

Log in to answer this question.

More posts like this