#!/bin/perl -w
    use strict;
    
    my $site_file=$ARGV[0];

    #input the methylation site info
    my %site_info=();
    open(IN_site,"<$site_file");
    while(<IN_site>)
    {
	    chomp;
	    my @arr=split(/\t/,$_,4);
	    if($arr[0] eq "MULTI"){next;}
	    my @chr_temp=split //,$arr[0];
	    for(1..3)
	    {
		    shift(@chr_temp);
	    }
	    $arr[0]=join("",@chr_temp);
	    $site_info{$arr[0]}{$arr[1]}{"strand"}=$arr[2];
	    $site_info{$arr[0]{$arr[1]}}{"value"}=$arr[3];
    }
    close IN_site;