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