Question: Fastq File Without Index?
0
gravatar for lsvijfhuizen
6.5 years ago by
lsvijfhuizen90
lsvijfhuizen90 wrote:

Dear All,

I have a fastq file which failed to align with bowtie. The error message is: reads file does not look like a FASTQ file. I don't know what is going wrong. Is it possible that there is no index in the fastq file? This is how my fastq file looks like:

@D5N3XBQ1:149:C0TCLACXX:2:1101:1342:2245 1:N:0:

TTAATAAAAGCTATTGAT

+

BB**?BHGG9?0?8888B

@D5N3XBQ1:149:C0TCLACXX:2:1101:1583:2228 1:N:0:

ATCGGGCGCAGGGACAGA

+

E99)?0?@6;'5>(;BE@

@D5N3XBQ1:149:C0TCLACXX:2:1101:1688:2238 1:N:0:

GCTGCCCTCCACCATATC

Thanks for all your ideas,

Lisanne

fastq bowtie index • 3.7k views
ADD COMMENTlink modified 6.5 years ago by Dave Richardson370 • written 6.5 years ago by lsvijfhuizen90
5
gravatar for Dave Richardson
6.5 years ago by
Cambridge, UK
Dave Richardson370 wrote:

Lisanne,

I don't normally see blank lines in a FASTQ file, I'd expect it to look like this:

@D5N3XBQ1:149:C0TCLACXX:2:1101:1342:2245 1:N:0:
TTAATAAAAGCTATTGAT
+
BB**?BHGG9?0?8888B
@D5N3XBQ1:149:C0TCLACXX:2:1101:1583:2228 1:N:0:
ATCGGGCGCAGGGACAGA
+
E99)?0?@6;'5>(;BE@
@D5N3XBQ1:149:C0TCLACXX:2:1101:1688:2238 1:N:0:
GCTGCCCTCCACCATATC

It also looks like the 3rd read is truncated, as the quality scores are missing.

However, a quick test with bowtie suggests that the blank lines aren't a problem. Is your fastq file compressed? (e.g. with gzip) That would give the error you're describing:

bowtie index.ebwt test.fastq.gz 
Error: reads file does not look like a FASTQ file
Command: bowtie index.ebwt test.fastq.gz

Piping uncompressed fastq into bowtie works

zcat test.fastq.gz | bowtie index.ebwt -
D5N3XBQ1:149:C0TCLACXX:2:1101:1342:2245 1:N:0:    +    XIV    548234    TTAATAAAAGCTATTGAT    BB**?BHGG9?0 8888B    0    5:G>A,14:G>T
D5N3XBQ1:149:C0TCLACXX:2:1101:1583:2228 1:N:0:    -    XI    611698    TCTGTCCCTGCGCCCGAT    @EB;(>5';6@?0?)99E    0    8:T>G,13:C>T
# reads processed: 2
# reads with at least one reported alignment: 2 (100.00%)
# reads that failed to align: 0 (0.00%)
Reported 2 alignments to 1 output stream(s)
ADD COMMENTlink modified 6.5 years ago • written 6.5 years ago by Dave Richardson370
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1138 users visited in the last hour