search

kamimrcht / REINDEER

REINDEER REad Index for abuNDancE quERy
GNU Affero General Public License v3.0
56 stars 6 forks source link

Output contains only * with no counts #26

Open tbenavi1 opened 2 weeks ago

tbenavi1 commented 2 weeks ago

Hello,

I created a Reindeer index on the T2T genome, and then queried the genome against the index in order to get the counts of all the kmers. The commands I ran were: 

bcalm -in chm13v2.0.fasta -kmer-size 31 -abundance-min 1 -nb-cores 32
Reindeer --index -f fof_unitigs.txt -o quick_out -t 32
Reindeer --query -q chm13v2.0.fasta -l quick_out -o quick_out/result.txt

However, for most of the chromosomes the result file only had "*". I have uploaded the reindeer index folder and result file here: https://drive.google.com/drive/folders/1P3ABKTFboLuiQ_YjWfsmz-Q3kLXWCfpL?usp=drive_link

Chromosomes 17, 18, and M have the expected output. I did run a query on a subset of chromosome 1, and it did seem to produce the correct output. I'm not sure if there is some error with really long sequence lengths. Thanks for any suggestions.