QAseq-Fusion_v2.1/py_script/test.py

#!/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)