Shicheng-Guo / bsmap

Automatically exported from code.google.com/p/bsmap
0 stars 0 forks source link

methratio.py error #8

Open GoogleCodeExporter opened 9 years ago

GoogleCodeExporter commented 9 years ago
What steps will reproduce the problem?
1.Run python methratio.py *.BSP > test.BSP 
2.
3.

What is the expected output? What do you see instead?
This script should generate output of summarized methylation calls for C's in 
the input.
----------------------------------------------------------------------
Instead, I get an following error message:
Mon Dec 19 10:44:32 2011: reading SRR019610.BSP ...
Traceback (most recent call last):
  File "/home2/ramamoor/softwares/bsmap-2.2/methratio.py", line 40, in ?
    map_info = get_alignment(line)
  File "/home2/ramamoor/softwares/bsmap-2.2/methratio.py", line 21, in get_alignment
    seq, strand, cr, pos, refseq = col[1], col[6][0], col[4], int(col[5]), col[8]
IndexError: list index out of range

----------------------------------------------------------------------

I closely looked at the .BSP file from bsmap and it looks like there are lines 
with only 4 columns. Here is an example from the BSMAP output file.
----------------------------------------------------------------------
SRR019489.60    AGAGGATTTTTAGAGATGGTGGAGAAATATAT    IIIIIIIIIII?IIIIIIIIIII?G3II<IDF    N
M 
----------------------------------------------------------------------

What version of the product are you using? On what operating system?

BSMAP v2.2 & RHEL 6.1
Please provide any additional information below.

Please let me know if I need to provide any additional information to fix this 
problem.

Thanks, Govind. BCM

Original issue reported on code.google.com by krgovind...@gmail.com on 19 Dec 2011 at 4:51

GoogleCodeExporter commented 9 years ago
NM means the reads is not mapped, this bug has been fixed in BSMAP-2.4

Original comment by xiyuanxin@gmail.com on 12 Jan 2012 at 10:12