Mercurial > repos > abims-sbr > mutcount
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)
