search

dpryan79 / MethylDackel

A (mostly) universal methylation extractor for BS-seq experiments.
MIT License
164 stars 44 forks source link

How to index genome file for MethylDackel? #150

Open harshraj1919 opened 1 year ago

harshraj1919 commented 1 year ago

I am running following script on my indexed and sorted bam file, but it is giving error on genome file. I indexed my genome file using BWA index.

MethylDackel extract -@ 20 --chunkSize 30000000 \ --methylKit Genome/Cow.fasta C10_BWA_sorted.bam.bai

Error message,

Couldn't load the index for Genome/Cow.fasta, will attempt to build it. Couldn't build the index for Genome/Cow.fasta! File corrupted?

Note : The file is not corrupted, I already used this file for mapping successfully.

Could you please help me to solve the issue?. How to index genome file for MethylDackel

Thanks, Harshraj

Loopy-fox commented 11 months ago

Hi ,did your error solved ? I meet the same issue.