search

zyndagj / BSMAPz

Updated and optimized fork of BSMAP
Other
21 stars 6 forks source link

methratio.py -t end repair option does not work #10

Open tAndreani opened 5 years ago

tAndreani commented 5 years ago

Hello,

am trying to extract the methylated cytosines. The program works perfectly but if I want to use the end repair option (-t or "--trim-fillin") it fails. This is the output:

Traceback (most recent call last):
  File "/data/SHARED_SOFTWARE/anaconda3/envs/scMeth/bin/methratio.py", line 571, in <module>
    main()
  File "/data/SHARED_SOFTWARE/anaconda3/envs/scMeth/bin/methratio.py", line 138, in main
    ret = map(chromWorker, argList)
  File "/data/SHARED_SOFTWARE/anaconda3/envs/scMeth/bin/methratio.py", line 283, in chromWorker
    refseq = refCache.fetch(pos, pos2)
  File "/data/SHARED_SOFTWARE/anaconda3/envs/scMeth/bin/methratio.py", line 214, in fetch
    assert(pos >= self.start)
AssertionError

I am using python-2.7.3

Any suggestions? Thanks in advance, Tommaso

zyndagj commented 5 years ago

Hello,

That option shouldn't affect the sequence cache, but I'll try to reproduce and get back to you.

-Greg

tAndreani commented 5 years ago

Thanks, Greg, let me know if you need the bam file I have used as input. Is not that big because is from a single cell.

Best Tommaso

zyndagj commented 4 years ago

I updated the reference cache in version 1.1.2, which should resolve your issue. Please pull it from github or conda and let me know how it goes.