fastq_screen error: Failed to calculate read length
1
0
Entering edit mode
18 months ago
firestar ★ 1.6k

Running it like this:

fastq_screen --conf $FASTQ_SCREEN_DATA/fastq_screen.conf --aligner bowtie2 --threads 6 03dpf.filtered.tagged.unmapped-1.fasta.gz

Runs into this error. I cannot seem to find any information about this.

Using fastq_screen v0.15.1
Reading configuration from '/FastQ_Screen_Genomes/fastq_screen.conf'
Adding database Human
Adding database Mouse
Adding database Rat
Adding database Drosophila
Adding database Worm
Adding database Yeast
Adding database Arabidopsis
Adding database Ecoli
Adding database rRNA
Adding database MT
Adding database PhiX
Adding database Lambda
Adding database Vectors
Adding database Adapters
Using 6 threads for searches
Option --subset set to 100000 reads
Processing 03dpf.filtered.tagged.unmapped-1.fasta.gz
Failed to calculate read length from 03dpf.filtered.tagged.unmapped-1.fasta.gz
Processing complete

The output file looks quite empty:

#Fastq_screen version: 0.15.1   #Aligner: bowtie2   #Reads in subset: 100000
Genome  #Reads_processed    #Unmapped   %Unmapped   #One_hit_one_genome %One_hit_one_genome #Multiple_hits_one_genome   %Multiple_hits_one_genome   #One_hit_multiple_genomes   %One_hit_multiple_genomes   Multiple_hits_multiple_genomes  %Multiple_hits_multiple_genomes

I am running this tool for the first time, so I could be missing something obvious.

fastq_screen • 961 views
ADD COMMENT
1
Entering edit mode

Failed to calculate read length from 03dpf.filtered.tagged.unmapped-1.fasta.gz

Based on the name it looks like this is a fasta input? As the name suggests fastq_screen is likely expecting fastq format input.

ADD REPLY
1
Entering edit mode
18 months ago
Trivas ★ 1.8k

From the documentation:

FastQ Screen allows you to screen a library of sequences in FastQ format

You're passing a .fasta file.

ADD COMMENT
0
Entering edit mode

Right. I completely overlooked that. Didn't even think that quality score would be that relevant for what this tool does.

ADD REPLY
1
Entering edit mode

You could convert the fasta file to include fake Q scores if you still want to use it as an input.

ADD REPLY
0
Entering edit mode

I created the .fasta from a BAM, so I am just going to go back to the BAM and convert that to a .fastq and see if that works. But good to know about fake Qs if I ever need that. Thanks :)

ADD REPLY

Login before adding your answer.

Traffic: 1539 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