search

marbl / meryl

A genomic k-mer counter (and sequence utility) with nice features.
119 stars 14 forks source link

eval/qv.sh meryl command not working #22

Open mas160 opened 3 years ago

mas160 commented 3 years ago

Hi,

I am using meryl pipelined with merqury via the script eval/qv.sh which repeatedly put .meryl files as input files but I am constantly getting this error: Can't interpret '*.0.meryl': not a meryl command, option, or recognized input file.

Is there an error in the script provided or is it an issue with meryl? even when I try using meryl difference with to .meryl files I get the same error: meryl-1.3/bin/meryl difference 1.meryl 2.meryl output all.meryl I get the error: Can't interpret '2.meryl': not a meryl command, option, or recognized input file.

Best,

arangrhie commented 3 years ago

Hello, try meryl statistics 2.meryl | head and see if that gives a reasonable kmer count summary. It seems like 2.meryl is an invalid kmer db.