def test_load_motif(self):
#motif_db_path = '/home/jiawei/MAmotif/motifscan_pkg/MotifScan/test/test_motif_db'
motif_path = '/mnt/MAmotif/motifscan_pkg/MotifScan/motif/hg19/hg19_jaspar_2014_all_1e-4.txt'
motif_list_path = ''
motif_table = motif.load_motif(motif_list_path, motif_path)
def test_motif_scan(self):
peak_path = 'Ese14_H3K27ac_LICR_Rep2_peaks.xls'
p = peak.load_peak(peak_path,'/home/jiawei/.MotifScan/genome/hg19',1000,'macs')
m = motif.load_motif('/mnt/MAmotif/motifscan_pkg/MotifScan/motif_list/jaspar_2014_vertebrates_motif_list.txt','/mnt/MAmotif/motifscan_pkg/MotifScan/motif/mm9/mm9_jaspar_2014_all_1e-4.txt')
core.motif_scan(p,m,np.array([0.29,0.21,0.21,0.29]),'./tmp')
