search

GuoliangLi-HZAU / BatMeth2

BS-seq analysis pipeline
38 stars 13 forks source link

genome index #17

Open fanghe0720 opened 4 years ago

fanghe0720 commented 4 years ago

Hi I'm trying to index human genome hg38 using BatMeth2 but it gave me error messages. I tried to index your genome.fa file from your google drive and it gave me the same error. Could you help to take a look?

Parsing FASTA file.. ParseFASTToPacked() : FASTA file does not begin with '>'! Deleting auxilliary files... mv: cannot stat /genome.fa.non.fmv': No such file or directory mv: cannot stat/genome.fa.non.bwt': No such file or directory mv: cannot stat /genome.fa.non.sa': No such file or directory rm: cannot remove/genome.fa.non.sai': No such file or directory BWTLoad() : cannot open bwtCodeFile!

fanghe0720 commented 4 years ago

I just noticed I didn't paste the full error message.

BatMeth2: build index rrbs
restriction enzyme digestion sites: CCGG
build_index_rrbs: malloc.c:2395: sysmalloc: Assertion `(old_top == initial_top (av) && old_size == 0) || ((unsigned long) (old_size) >= MINSIZE && prev_inuse (old_top) && ((unsigned long) old_end & (pagesize - 1 )) == 0)' failed.
~/BatMeth2/BatMeth2/bin/build_all: line 18: 88960 Aborted (core dumped) $CPATH/build_index_rrbs $BASEX