Filtering short reads from .fastq

[samt]

I am getting an error from Bowtie as follows:

Error: Read (SRR015250.2
S0014_20071116_1_ES_EStranscriptome_1_39_146_F3 length=35) is less
than 2 characters long

Is there a program or option in Bowtie to filter these out from the .fastq file?

[maubp]

Has your file already been processed?

It looks like the record header says it is length 35, but Bowtie says it is less than 2. Why not search the file to see what this read looks like? It is possible the file has been corrupted or truncated...

[If you have a valid FASTQ file I could share a short Biopython script to filter out short reads]

[samt]

Hi,

I took this fastq file straight from NCBIs Short Read Archive
the contents look like:
@SRR015253.1 LIZ_20071025_2_GrimmondsMES_SS7747_13_23_38_F3 length=35
T1011122220100230032132.2111111002.1
+SRR015253.1 LIZ_20071025_2_GrimmondsMES_SS7747_13_23_38_F3 length=35
!)+%.*%*+2'0%%%-%+%*5'%!%9+'%+<+0%!%
@SRR015253.2 LIZ_20071025_2_GrimmondsMES_SS7747_13_23_119_F3 length=35
T0101233211103200232333.2111211002.1
+SRR015253.2 LIZ_20071025_2_GrimmondsMES_SS7747_13_23_119_F3 length=35
!,.+'+')'390%%%%%%%'%%%!-<++++<99%!%
@SRR015253.3 LIZ_20071025_2_GrimmondsMES_SS7747_13_23_146_F3 length=35
T0312202213101213131111.1110131102.1
+SRR015253.3 LIZ_20071025_2_GrimmondsMES_SS7747_13_23_146_F3 length=35
!93<*/18+%:9%+075*%:;+6!3<26%/<%-%!%

I'm not sure how to make sense of this but bowtie seems to think the reads are less than 2 bp

[nilshomer]

Originally Posted by samt

Hi,

I took this fastq file straight from NCBIs Short Read Archive
the contents look like:
@SRR015253.1 LIZ_20071025_2_GrimmondsMES_SS7747_13_23_38_F3 length=35
T1011122220100230032132.2111111002.1
+SRR015253.1 LIZ_20071025_2_GrimmondsMES_SS7747_13_23_38_F3 length=35
!)+%.*%*+2'0%%%-%+%*5'%!%9+'%+<+0%!%
@SRR015253.2 LIZ_20071025_2_GrimmondsMES_SS7747_13_23_119_F3 length=35
T0101233211103200232333.2111211002.1
+SRR015253.2 LIZ_20071025_2_GrimmondsMES_SS7747_13_23_119_F3 length=35
!,.+'+')'390%%%%%%%'%%%!-<++++<99%!%
@SRR015253.3 LIZ_20071025_2_GrimmondsMES_SS7747_13_23_146_F3 length=35
T0312202213101213131111.1110131102.1
+SRR015253.3 LIZ_20071025_2_GrimmondsMES_SS7747_13_23_146_F3 length=35
!93<*/18+%:9%+075*%:;+6!3<26%/<%-%!%

I'm not sure how to make sense of this but bowtie seems to think the reads are less than 2 bp

These are color space data. Bowtie does not support ABI SOLiD data. You could try aligners that do, including BFASThttp://genome.ucla.edu/bfast, BWA, MAQ, or SHRiMP.