# HG changeset patch
# User cpt
# Date 1704433866 0
# Node ID 673d1776d3b98e7d1454e0235b31d54494da9d70
# Parent fd70980a516b2db8b444cc15577cbb06234cc421
planemo upload commit f33bdf952d796c5d7a240b132af3c4cbd102decc
diff -r fd70980a516b -r 673d1776d3b9 findSpanin.py
--- a/findSpanin.py Mon Jun 05 02:42:01 2023 +0000
+++ b/findSpanin.py Fri Jan 05 05:51:06 2024 +0000
@@ -19,6 +19,7 @@
#### INPUT : putative_isp.fa & putative_osp.fa (in that order)
#### PARAMETERS :
+
###############################################################################
def write_output(candidates):
"""output file function...maybe not needed"""
@@ -132,7 +133,6 @@
spanin_data["positive"]["amt_positive"] = 0
spanin_data["positive"]["pos_amt_unique"] = 0
if spanin_data["negative"]["coords"]:
-
# do something else...
# print('in other function')
# print(spanin_data['negative']['coords'])
@@ -196,7 +196,6 @@
if __name__ == "__main__":
-
# Common parameters for both ISP / OSP portion of script
parser = argparse.ArgumentParser(
@@ -330,8 +329,7 @@
for k, v in SPANIN_TYPES["EMBEDDED"].items():
# print(k)
f.writelines(
- ""
- + str(k)
+ +str(k)
+ " ==> Amount of corresponding candidate o-spanins(s): "
+ str(len(v))
+ "\n"
@@ -339,8 +337,7 @@
if each_spanin_type == "SEPARATED":
for k, v in SPANIN_TYPES["SEPARATED"].items():
f.writelines(
- ""
- + str(k)
+ +str(k)
+ " ==> Amount of corresponding candidate o-spanins(s): "
+ str(len(v))
+ "\n"
@@ -348,8 +345,7 @@
if each_spanin_type == "OVERLAPPED":
for k, v in SPANIN_TYPES["OVERLAPPED"].items():
f.writelines(
- ""
- + str(k)
+ +str(k)
+ " ==> Amount of corresponding candidate o-spanins(s): "
+ str(len(v))
+ "\n"
@@ -414,15 +410,13 @@
if each_posp[6] == "+":
if each_posp[2] in pair_dict["pairs"]["pair_number"].keys():
f.write(
- ""
- + str(pair_dict["pairs"]["pair_number"][each_posp[2]])
+ +str(pair_dict["pairs"]["pair_number"][each_posp[2]])
+ "\n"
)
elif each_posp[6] == "-":
if each_posp[1] in pair_dict["pairs"]["pair_number"].keys():
f.write(
- ""
- + str(pair_dict["pairs"]["pair_number"][each_posp[1]])
+ +str(pair_dict["pairs"]["pair_number"][each_posp[1]])
+ "\n"
)
else:
@@ -491,15 +485,13 @@
if each_posp[2] in pair_dict["pairs"]["pair_number"].keys():
# print('ovl ; +')
f.write(
- ""
- + str(pair_dict["pairs"]["pair_number"][each_posp[2]])
+ +str(pair_dict["pairs"]["pair_number"][each_posp[2]])
+ "\n"
)
elif each_posp[6] == "-":
if each_posp[1] in pair_dict["pairs"]["pair_number"].keys():
f.write(
- ""
- + str(pair_dict["pairs"]["pair_number"][each_posp[1]])
+ +str(pair_dict["pairs"]["pair_number"][each_posp[1]])
+ "\n"
)
else:
@@ -565,15 +557,13 @@
if each_posp[6] == "+":
if each_posp[2] in pair_dict["pairs"]["pair_number"].keys():
f.write(
- ""
- + str(pair_dict["pairs"]["pair_number"][each_posp[2]])
+ +str(pair_dict["pairs"]["pair_number"][each_posp[2]])
+ "\n"
)
elif each_posp[6] == "-":
if each_posp[1] in pair_dict["pairs"]["pair_number"].keys():
f.write(
- ""
- + str(pair_dict["pairs"]["pair_number"][each_posp[1]])
+ +str(pair_dict["pairs"]["pair_number"][each_posp[1]])
+ "\n"
)
else:
diff -r fd70980a516b -r 673d1776d3b9 findSpanin.xml
--- a/findSpanin.xml Mon Jun 05 02:42:01 2023 +0000
+++ b/findSpanin.xml Fri Jan 05 05:51:06 2024 +0000
@@ -2,7 +2,7 @@
With the outputs from the ISP and OSP candidate tools, cull the list down to candidate pairs
macros.xml
- cpt-macros.xml
+