Skip to content

Commit

Permalink
Add files via upload
Browse files Browse the repository at this point in the history
  • Loading branch information
rwnull committed Jun 3, 2020
1 parent e151ded commit 259bfe5
Show file tree
Hide file tree
Showing 2 changed files with 20 additions and 21 deletions.
24 changes: 12 additions & 12 deletions HCRProbeMaker2_7_v2.py
Original file line number Diff line number Diff line change
Expand Up @@ -6,7 +6,7 @@
__copyright__ = "Copyright 2019-2020, The Ozpolat Lab (http://bduyguozpolat.org/)"
__credits__ = ["B. Duygu Ozpolat"]
__license__ = "GPL 3.0"
__version__ = "2020_2.0"
__version__ = "2020_0.2.0"



Expand Down Expand Up @@ -93,7 +93,7 @@ def amp(ampl):

print ""
print ""
print "HCR3.0_Probe_Maker_Output"
print "HCR3.0 Probe Maker Output"
print ""
print ""

Expand All @@ -105,7 +105,7 @@ def amp(ampl):
print "Figure Layout:"
print ""
print str(amplifier+"_"+str(name)+"_"+str(count)) # Name based on inputs
print "Pair#\t1st_Half_of_Initiator_I1\tSpacer\tProbe\t\tProbe\tSpacer\t2nd_Half_of_Initiator_I1" # HEADER of Output
print "Pair#,1st Half of Initiator I1,Spacer,Probe,,Probe,Spacer,2nd Half of Initiator I1" # HEADER of Output


position=cdna-pause # This controls how far from the 5'end of the mRNA probes begin
Expand All @@ -116,11 +116,11 @@ def amp(ampl):
while position>52: # 52 is the cutoff for fitting an entire pair at the end of the gene. the program will cycle back over the RNA if not limited like this
downstream=str(fullseq[position-25:position])
upstream=str(fullseq[position-52:position-27])
pairlib[pair]=str(str(pair)+"\t"+str(cdna-position+25)+"\t"+str(fullseq[position-25:position])+"\t"+str(cdna-position)+"\t\t"+str(cdna-position+52)+"\t"+str(fullseq[position-52:position-27])+"\t"+str(cdna-position+27))
idtlibu[pair]=str(amplifier+"_"+str(name)+"_"+str(count)+"\t"+upinit+uspc+upstream) # This is a library used for IDT output
idtlibd[pair]=str(amplifier+"_"+str(name)+"_"+str(count)+"\t"+downstream+dspc+dninit)
pairlib[pair]=str(str(pair)+","+str(cdna-position+25)+","+str(fullseq[position-25:position])+","+str(cdna-position)+",,"+str(cdna-position+52)+","+str(fullseq[position-52:position-27])+","+str(cdna-position+27))
idtlibu[pair]=str(amplifier+"_"+str(name)+"_"+str(count)+","+upinit+uspc+upstream) # This is a library used for IDT output
idtlibd[pair]=str(amplifier+"_"+str(name)+"_"+str(count)+","+downstream+dspc+dninit)
position-=54 # 54 is the number of bases covered by one probe set (25bp probe 1, 2bp spacer, 25bp probe 2, 2bp spacer).
print str(pair)+"\t"+upinit+"\t"+uspc+"\t"+upstream+"\t\t"+downstream+"\t"+dspc+"\t"+dninit
print str(pair)+","+upinit+","+uspc+","+upstream+",,"+downstream+","+dspc+","+dninit
if pair<count:
pair+=1
else:
Expand All @@ -134,10 +134,10 @@ def amp(ampl):
print ""
print "Below are the hybridizing sequences and where they align to the cDNA:"
print ""
print "Pair#\tcDNAcoord\tProbe\tcDNAcoord\t\tcDNAcoord\tProbe\tcDNAcoord" # Header for Output
print "Pair#,cDNAcoord,Probe,cDNAcoord,,cDNAcoord,Probe,cDNAcoord" # Header for Output
i=1
while i <= pair:
print(pairlib[i])
print pairlib[i]
i+=1

print ""
Expand All @@ -150,11 +150,11 @@ def amp(ampl):
print "Below is in IDT oPool submission format."
print "Copy and Paste this into an XLSX file for submission to IDT."
print ""
print "Pool name\tsequence" #Header for output
print "Pool name,sequence" #Header for output
i=1
while i <= pair:
print(str(idtlibu[i])
print(str(idtlibd[i])
print str(idtlibu[i])
print str(idtlibd[i])
i+=1


Expand Down
17 changes: 8 additions & 9 deletions HCRProbeMaker3_7_v2.py
Original file line number Diff line number Diff line change
Expand Up @@ -5,7 +5,7 @@
__copyright__ = "Copyright 2019-2020, The Ozpolat Lab (http://bduyguozpolat.org/)"
__credits__ = ["B. Duygu Ozpolat"]
__license__ = "GPL 3.0"
__version__ = "2020_2.0"
__version__ = "2020_0.2.0"



Expand Down Expand Up @@ -105,7 +105,7 @@ def amp(ampl):
print("Figure Layout:")
print("")
print(str(amplifier+"_"+str(name)+"_"+str(count))) # Name based on inputs
print("Pair#\t1st_Half_of_Initiator_I1\tSpacer\tProbe\t\tProbe\tSpacer\t2nd_Half_of_Initiator_I1") # HEADER of Output
print("Pair#,1st Half of Initiator I1,Spacer,Probe,,Probe,Spacer,2nd Half of Initiator I1") # HEADER of Output


position=cdna-pause # This controls how far from the 5'end of the mRNA probes begin
Expand All @@ -116,11 +116,11 @@ def amp(ampl):
while position>52: #52 is the cutoff for fitting an entire pair at the end of the gene. the program will cycle back over the RNA if not limited like this
downstream=str(fullseq[position-25:position])
upstream=str(fullseq[position-52:position-27])
pairlib[pair]=str(str(pair)+"\t"+str(cdna-position+25)+"\t"+str(fullseq[position-25:position])+"\t"+str(cdna-position)+"\t\t"+str(cdna-position+52)+"\t"+str(fullseq[position-52:position-27])+"\t"+str(cdna-position+27))
idtlibu[pair]=str(amplifier+"_"+str(name)+"_"+str(count)+"\t"+upinit+uspc+upstream) # This is a library used for IDT output
idtlibd[pair]=str(amplifier+"_"+str(name)+"_"+str(count)+"\t"+downstream+dspc+dninit)
pairlib[pair]=str(str(pair)+","+str(cdna-position+25)+","+str(fullseq[position-25:position])+","+str(cdna-position)+",,"+str(cdna-position+52)+","+str(fullseq[position-52:position-27])+","+str(cdna-position+27))
idtlibu[pair]=str(amplifier+"_"+str(name)+"_"+str(count)+","+upinit+uspc+upstream) # This is a library used for IDT output
idtlibd[pair]=str(amplifier+"_"+str(name)+"_"+str(count)+","+downstream+dspc+dninit)
position-=54 # 54 is the number of bases covered by one probe set (25bp probe 1, 2bp spacer, 25bp probe 2, 2bp spacer).
print(str(pair)+"\t"+upinit+"\t"+uspc+"\t"+upstream+"\t\t"+downstream+"\t"+dspc+"\t"+dninit)
print(str(pair)+","+upinit+","+uspc+","+upstream+",,"+downstream+","+dspc+","+dninit)
if pair<count:
pair+=1
else:
Expand All @@ -134,7 +134,7 @@ def amp(ampl):
print("")
print("Below are the hybridizing sequences and where they align to the cDNA:")
print("")
print("Pair#\tcDNAcoord\tProbe\tcDNAcoord\t\tcDNAcoord\tProbe\tcDNAcoord") # Header for Output
print("Pair#,cDNAcoord,Probe,cDNAcoord,,cDNAcoord,Probe,cDNAcoord") # Header for Output
i=1
while i <= pair:
print(pairlib[i])
Expand All @@ -150,7 +150,7 @@ def amp(ampl):
print("Below is in IDT oPool submission_format.")
print("Copy and Paste this into an XLSX file for submission to IDT.")
print("")
print("Pool_name\tsequence") #Header for output
print("Pool name,sequence") #Header for output
i=1
while i <= pair:
print(str(idtlibu[i]))
Expand All @@ -171,7 +171,6 @@ def amp(ampl):
print(fullseq)


# In[ ]:



Expand Down

0 comments on commit 259bfe5

Please sign in to comment.