search

kage-genotyper / kage

Alignment-free genotyper for SNPs and short indels, implemented in Python.
GNU General Public License v3.0
48 stars 2 forks source link

Error when running kage index #15

Open zhiyaoma2t opened 3 months ago

zhiyaoma2t commented 3 months ago

Hi,

I'm trying to run the kage index step of kage but I get the following error: TypeError: init() missing 1 required positional argument: 'dtype'

My input VCF is only including SV. The command I used was:

kage index -r primary.fa -v primary_all_norm_sort_SV.vcf -o index -k 31

Best regards!

ivargr commented 3 months ago

Hi!

Thanks for reporting! That seems a bit strange. Do you get a longer stacktrace (output) that you could share?

zhiyaoma2t commented 3 months ago

Thank you , here is my log file. log.txt

Best regards!

ivargr commented 2 months ago

Thanks for sharing the log.

It seems that the problem actually is that you run out of memory because some large memory allocation in the indexing process. The error message you get is just another error when formatting the first error message. I'm not sure why KAGE is trying to allocate that much memory, but I'm looking into it.

Would you by any chance be able to share the VCF and reference you are using with me? That would make it a bit easier to find exactly what the problem is, but I understand if this is data that it is not possible to share :)

zhiyaoma2t commented 2 months ago

Hi,

I change my VCF file and reference to run the kage index step of kage. But it report another error. here is my log file. log.txt

Best

ivargr commented 2 months ago

Hi!

Sorry for the late reply. This seems like a memory issue again, which happens when KAGE tries to fetch the reference genome sequence around your variants.

May I ask how many variants you have in your VCF? Are these SVs?

Again, if you would be able to share your vcf and reference that would be very useful, but understand if not possible :)