@@ -59,24 +59,24 @@ def inDelphiScore(guideDict,genome_fa,twoBitToFa_path,genome_2bit,tempfiles_path
59
59
#TODO: NEED TO REPLACE NEXT LINE
60
60
# the last term should be replaced with ~20 bps after the cutsite
61
61
# get this via the get_sequence code
62
- print (guide ['guide_seq' ])
62
+ # print(guide['guide_seq'])
63
63
64
64
cleavePos = int (str (cleaves )[2 :- 2 ])
65
- print (cleavePos )
65
+ # print(cleavePos)
66
66
#get strand
67
- print (guide ['strand' ])
67
+ # print(guide['strand'])
68
68
#apply strand logic to coordinates before submitting to get_sequence.py
69
69
if guide ['strand' ] == "+" :
70
70
coords = str ("%s:%s-%s" % (guide ['pam_chrom' ],guide ['guide_genomic_start' ],(guide ['pam_genomic_start' ]+ cleavePos + 20 )))
71
71
elif guide ['strand' ] == "-" :
72
72
coords = str ("%s:%s-%s" % (guide ['pam_chrom' ],(guide ['pam_genomic_start' ]- len (guide ['guide_seq' ])+ 3 ),(guide ['pam_genomic_start' ]+ 1 + cleavePos + 22 )))
73
73
# coords=str("%s:%s-%s"%(guide['pam_chrom'],guide['guide_genomic_start'],(guide['pam_genomic_start']+cleavePos+20)))
74
74
75
- print (coords )
75
+ # print(coords)
76
76
77
77
timestr = time .strftime ("%Y%m%d-%H%M%S" )
78
78
getseq = str (get_sequence .fetch_sequence (twoBitToFa_path , coords , genome_2bit , os .path .join (tempfiles_path ,timestr + '_out.fa' )))
79
- print (os .path .join (tempfiles_path ,timestr + '_out.fa' ))
79
+ # print(os.path.join(tempfiles_path,timestr+'_out.fa'))
80
80
81
81
#get flanking sequence before and after cutsite
82
82
fullSeq = ''
@@ -99,7 +99,7 @@ def inDelphiScore(guideDict,genome_fa,twoBitToFa_path,genome_2bit,tempfiles_path
99
99
fullSeq = "" .join (list (reversed (fullSeq ))).translate (translate_code )
100
100
101
101
cutsite = len (guide ['guide_seq' ]) + cleavePos
102
- print (fullSeq )
102
+ # print(fullSeq)
103
103
104
104
pred_df , stats = inDelphi .predict (fullSeq , cutsite )
105
105
guide ['Precision' ] = round (stats ['Precision' ],2 )
0 commit comments