Mercurial > repos > estrain > seqsero2_v1_0_1
changeset 24:4e1596ee5c94
Uploaded
author | estrain |
---|---|
date | Sat, 07 Sep 2019 13:56:48 -0400 |
parents | f8c06e3b0fb9 |
children | 55d2ec45737f |
files | SeqSero2_package.py |
diffstat | 1 files changed, 6 insertions(+), 0 deletions(-) [+] |
line wrap: on
line diff
--- a/SeqSero2_package.py Sat Sep 07 13:51:19 2019 -0400 +++ b/SeqSero2_package.py Sat Sep 07 13:56:48 2019 -0400 @@ -1215,15 +1215,21 @@ def judge_subspecies(fnameA,dirpath): #seqsero2 -a; judge subspecies on just forward raw reads fastq salmID_output=subprocess.Popen("./SalmID.py -i "+fnameA,shell=True,stdout=subprocess.PIPE,stderr=subprocess.PIPE) + salmID_output2=subprocess.Popen("ls",shell=True,stdout=subprocess.PIPE,stderr=subprocess.PIPE) out, err = salmID_output.communicate() + out2, err2 = salmID_output2.communicate() out=out.decode("utf-8") + out2=out2.decode("utf-8") err=err.decode("utf-8") + err2=err2.decode("utf-8") file=open("data_log.txt","a") file.write(out) file.close() print("dirpath " + str(dirpath)) print("SalmID File" + str(out)) + print("SalmID File" + str(out2)) print("SalmID Err File" + str(err)) + print("SalmID Err File" + str(err2)) print("SalmID 1 "+str(out.split("\n")[1])) print("SalmID 0 "+str(out.split("\n")[0])) salm_species_scores=out.split("\n")[1].split("\t")[6:]