-
Notifications
You must be signed in to change notification settings - Fork 70
Description
I'm having trouble getting results from fastANI when comparing two fungal genomes. Here's the command I ran and the output I received:
fastANI -r /path/to/reference/fungi/GCA_004299765.1.fna -q /path/to/query/fungi/GCA_002105135.1.fna -o $PWD/test.2.csvThe output shows:
>>>>>>>>>>>>>>>>>>
Reference = [/path/to/reference/fungi/GCA_004299765.1.fna]
Query = [/path/to/query/fungi/GCA_002105135.1.fna]
Kmer size = 16
Fragment length = 3000
Threads = 1
ANI output file = /path/to/output/test.2.csv
Sanity Check = 0
>>>>>>>>>>>>>>>>>>
INFO [thread 0], skch::main, Count of threads executing parallel_for : 1
INFO [thread 0], skch::Sketch::build, window size for minimizer sampling = 24
INFO [thread 0], skch::Sketch::build, minimizers picked from reference = 2279545
INFO [thread 0], skch::Sketch::index, unique minimizers = 2234242
INFO [thread 0], skch::Sketch::computeFreqHist, Frequency histogram of minimizers = (1, 2196715) ... (70, 1)
INFO [thread 0], skch::Sketch::computeFreqHist, consider all minimizers during lookup.
INFO [thread 0], skch::main, Time spent sketching the reference : 2.4358 sec
INFO [thread 0], skch::main, Start Map 1
INFO [thread 0], skch::main, Time spent mapping fragments in query #1 : 2.31659 sec
INFO [thread 0], skch::main, Time spent post mapping : 6.8426e-05 sec
INFO [thread 0], skch::main, ready to exit the loop
INFO, skch::main, parallel_for execution finished
After running the command, the output file (test.2.csv) is empty:
$ ll
total 16
drwxrwxr-x 2 user group 4096 Jul 17 07:10 ./
drwxrwxr-x 3 user group 4096 Jul 17 06:56 ../
-rw-rw-r-- 1 user group 391 Jul 17 07:07 test.1.csv
-rw-rw-r-- 1 user group 0 Jul 17 07:14 test.2.csv
-rw-rw-r-- 1 user group 391 Jul 17 07:05 test.csvOther similar commands have produced non-empty results (e.g., test.1.csv and test.csv), but this specific comparison between these two fungal genomes is resulting in an empty output file. Both input files are valid FASTA files containing genomic sequences from different fungal species.
I've checked the input files, and they seem to be in the correct format. The reference genome is from Talaromyces funiculosus, and the query genome is from Syncephalastrum racemosum.
Has anyone encountered this issue before or have any suggestions on what might be causing this? Any help would be greatly appreciated!