#!/usr/bin/python # -*- coding:utf-8 -*- from SearchReadsFusion import DoubleReadsFusion as DRF from SearchReadsFusion import SingleReadsFusion as SRF list1=['A00301:652:HHJ75DSX5:2:1677:20329:35884', 99, 'chr2', 34845763, 60, '64S86M', 'TGCTTGTCTCAGCCTCCCAAAGCACTGGAATTACTGGCGTGGGCTACTGCAACTCGCCTTATTATTTTTTTTGAGGTGGAGTCTCGCTCTGTCCCCAGGCTGGAGTGCAGTGGCGCGATCTCCGCTCACCGCAAGCTCCGCCTCCTGGGT'] list2=['A00301:652:HHJ75DSX5:2:1677:20329:35884', 2163, 'chr2', 42305096, 60, '58H85M7H', 'ACAGAGCGAGACTCCACCTCAAAAAAAATAATAAGGCGAGTTGCAGTAGCCCACGCCAGTAATTCCAGTGCTTTGGGAGGCTGAG'] reads_info=[list1,list2] drf =DRF() results=drf.HardClipFusion(reads_info) print(results)