diff scripts/S02b_study_seq_composition_nuc.py @ 9:04a9ada73cc4 draft

planemo upload for repository htpps://github.com/abims-sbr/adaptearch commit f1ba8d136e0129f3e8435b25a95f70f697d51464-dirty
author abims-sbr
date Tue, 03 Jul 2018 10:55:46 -0400
parents f1e24200e5ae
children
line wrap: on
line diff
--- a/scripts/S02b_study_seq_composition_nuc.py	Fri Mar 16 05:11:04 2018 -0400
+++ b/scripts/S02b_study_seq_composition_nuc.py	Tue Jul 03 10:55:46 2018 -0400
@@ -95,13 +95,18 @@
 os.makedirs(Path_IN_loci_NUC)
 os.makedirs(outpath)
 
-infiles = str.split(sys.argv[1], ",")
+infiles = []
+with open(sys.argv[2], 'r') as f:
+    for line in f.readlines():
+        infiles.append(line.strip('\n'))
+
+#infiles = str.split(sys.argv[1], ",")
 for file in infiles:
     os.system("cp %s %s" %(file, Path_IN_loci_NUC))
 
 ## 1 ## List taxa
 LT=[]
-cmd="grep '>' %s" % sys.argv[2]
+cmd="grep '>' {}".format(sys.argv[1])
 result = subprocess.check_output(cmd, shell=True)
 result=result.split('\n')
 for i in result:   
@@ -110,11 +115,7 @@
         LT.append(sp)
 print LT
 
-#LT = ["Ap", "Ac", "Pg"]
-#os.system("grep '>' %s" %(sys.argv[1]))
-
 ## 2 ## PathIN
-# fileIN_properties = open("01_AminoAcid_Properties2.csv", "r")
 Lloci_NUC = os.listdir(Path_IN_loci_NUC)