Commit db4e58f0 authored by Koehorst, Jasper's avatar Koehorst, Jasper
Browse files

version eggnog db added

parent fd8cba38
......@@ -47,18 +47,18 @@ def conversion(identifier):
# "-tool", "conversion", "-f", fasta.getAbsolutePath(), "-gff2rdf", "-i", gff.getAbsolutePath(), "-o", output.getAbsolutePath(), "-id", "MGYG-HGUT-00001", "-codon", "11", "-topology", "linear"
command = "java -jar ./MGnifyParser.jar -tool conversion -f ./genome/" + identifier + "/" + identifier + ".fna -gff2rdf -i ./genome/" + identifier + "/" + identifier + ".gff -id " + identifier + " -codon 11 -topology linear -o ./genome/" + identifier + "/" + identifier + ".conversion.ttl"
print(command)
os.system(command)
# os.system(command)
def interproscan(identifier):
command = "java -jar ./MGnifyParser.jar -tool interpro -i ./genome/" + identifier + "/" + identifier + ".conversion.ttl -o ./genome/" + identifier + "/" + identifier + ".interproscan.ttl -tsv ./genome/" +identifier +"/" + identifier + "_InterProScan.tsv -version InterProScan-v5_35-74_0";
print(command)
os.system(command)
# os.system(command)
def eggnog(identifier):
command = "java -jar ./MGnifyParser.jar -tool eggnog -i ./genome/" + identifier + "/" + identifier + ".interproscan.ttl -o ./genome/" + identifier + "/" + identifier + ".interproscan.eggnog.ttl -tsv ./genome/" +identifier +"/" + identifier + "_eggNOG.tsv";
command = "java -jar ./MGnifyParser.jar -tool eggnog -i ./genome/" + identifier + "/" + identifier + ".interproscan.ttl -o ./genome/" + identifier + "/" + identifier + ".interproscan.eggnog.ttl -tsv ./genome/" +identifier +"/" + identifier + "_eggNOG.tsv -v v2-db5.0";
print(command)
os.system(command)
# os.system(command)
if __name__ == "__main__":
for line in open("identifiers.tsv"):
identifier = line.strip().split("\t")[1]
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment