TideHunter
TideHunter copied to clipboard
[abpoa_gen_cons] "Not enough sequences to perform msa"
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