Input file format in methylkit package

[pj]

Hi,

This my first time using the methylkit package. I have a bed file output from bsmap consisting of chr, start, end, methylation ratio, coverage but the strand is split by coverage (eg. 36 sites split into 31 on + strand and 5 on – strand). So how can I prepare the strand column for methylkit read function? Has anyone encounter this issue or I can prepare the input file without strand information, is it going to work?