TideHunter icon indicating copy to clipboard operation
TideHunter copied to clipboard

[abpoa_gen_cons] "Not enough sequences to perform msa"

Open santiago-es opened this issue 2 years ago • 10 comments

So after successfully compiling TideHunter and AbPoa I'm running into this error when attempting to generate a consensus sequence:

[abpoa_gen_cons] "Not enough sequences to perform msa"

However, my sequence file is of length: 195588 (as measured by zcat < *.fq.gz | wc -l) so I feel like that should be more than enough sequence?

I am attempting to use a large gzipped fastq file as TideHunter input. As a minimally reproducible example, I will include on of the constitutent fastq files which is much smaller (FAR....12.fastq.gz):

FAR63237_pass_barcode01_9dc2df5e_12.fastq.gz

These fastqs were generated by the default MinKnow basecaller from FAST5 files produced by nanopore sequencing on a MinION 9.3.4 flow cell.

My execute command for tidehunter was:

./TideHunter-v1.5.3.2/bin/TideHunter -f 3 *.fastq.gz

and I also attempted

./TideHunter-v1.5.3.2/bin/TideHunter -f 3 $(zcat < *.fastq.gz)

which has not yet failed with this error but also seems to be taking much longer than anticipated to run.

Thanks for your help

santiago-es avatar Feb 25 '22 16:02 santiago-es