This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Does Jellyfish support zipped fastq sequences?

I used jellyfish to counting kmers and my input files were zipped fastq. I used the commond :

jellyfish count -m 17 -o 17mer_counts -s 500000000 -t 5 -c 3 read_1.fq.gz read_2.fq.gz

And I got the errors like this:

terminate called after throwing an instance of 'std::runtime_error'
  what():  Unsupported format

I don't know how to solve it.

genome assembly kmer
• 8,765 views • link
updated by GenoMax GenoMax • written by comma1993 • comma1993

Answer to your question is no. Workaround is described by @Sej below.

• 0 views • link
written by GenoMax GenoMax

1 answer

Jellyfish manual: http://www.genome.umd.edu/docs/JellyfishUserGuide.pdf

2.1 How to read compressed files (or other format) Jellyfish only reads FASTA or FASTQ formatted input files. By reading from pipes, jellyfish can read compressed files, like this: zcat *.fastq.gz | jellyfish count /dev/fd/0 ... or by using the '<()' redirection provided by the shell (e.g. bash, zsh): jellyfish count <(zcat file1.fastq.gz) <(zcat file2.fasta.gz) ...

• 0 views • link
updated by WouterDeCoster WouterDeCoster • written by Sej Modha Sej Modha

Thanks for your reply~~It helps a lot.

• 0 views • link
written by comma1993 • comma1993

Log in to answer this question.

More posts like this