Mercurial > repos > rdvelazquez > mmvc
annotate extract_mmvc_consensus @ 3:4ca91bf87f1f draft default tip
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/mmvc/ commit b'35710f754cc8be97c423c95eb91197df01d0520d\n'-dirty
| author | rdvelazquez | 
|---|---|
| date | Wed, 20 Mar 2019 15:58:26 -0400 | 
| parents | c5c1faa433e4 | 
| children | 
| rev | line source | 
|---|---|
| 
0
 
c5c1faa433e4
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/mmvc/ commit b'f226ae444ad694bb8d4f8faa80ed541b41cbc993\n'
 
rdvelazquez 
parents:  
diff
changeset
 | 
1 #!/usr/bin/env python3 | 
| 
 
c5c1faa433e4
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/mmvc/ commit b'f226ae444ad694bb8d4f8faa80ed541b41cbc993\n'
 
rdvelazquez 
parents:  
diff
changeset
 | 
2 | 
| 
 
c5c1faa433e4
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/mmvc/ commit b'f226ae444ad694bb8d4f8faa80ed541b41cbc993\n'
 
rdvelazquez 
parents:  
diff
changeset
 | 
3 import json | 
| 
 
c5c1faa433e4
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/mmvc/ commit b'f226ae444ad694bb8d4f8faa80ed541b41cbc993\n'
 
rdvelazquez 
parents:  
diff
changeset
 | 
4 import sys | 
| 
 
c5c1faa433e4
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/mmvc/ commit b'f226ae444ad694bb8d4f8faa80ed541b41cbc993\n'
 
rdvelazquez 
parents:  
diff
changeset
 | 
5 | 
| 
 
c5c1faa433e4
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/mmvc/ commit b'f226ae444ad694bb8d4f8faa80ed541b41cbc993\n'
 
rdvelazquez 
parents:  
diff
changeset
 | 
6 # first argument is script name; second should be mmvcOutput path; third should be name for fasta header; the fourth argument should be the path to save the fasta at | 
| 
 
c5c1faa433e4
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/mmvc/ commit b'f226ae444ad694bb8d4f8faa80ed541b41cbc993\n'
 
rdvelazquez 
parents:  
diff
changeset
 | 
7 mmvcOutput = sys.argv[1] | 
| 
 
c5c1faa433e4
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/mmvc/ commit b'f226ae444ad694bb8d4f8faa80ed541b41cbc993\n'
 
rdvelazquez 
parents:  
diff
changeset
 | 
8 fastaHeaderName = sys.argv[2] | 
| 
 
c5c1faa433e4
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/mmvc/ commit b'f226ae444ad694bb8d4f8faa80ed541b41cbc993\n'
 
rdvelazquez 
parents:  
diff
changeset
 | 
9 fastaFilePath = sys.argv[3] | 
| 
 
c5c1faa433e4
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/mmvc/ commit b'f226ae444ad694bb8d4f8faa80ed541b41cbc993\n'
 
rdvelazquez 
parents:  
diff
changeset
 | 
10 | 
| 
 
c5c1faa433e4
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/mmvc/ commit b'f226ae444ad694bb8d4f8faa80ed541b41cbc993\n'
 
rdvelazquez 
parents:  
diff
changeset
 | 
11 with open(mmvcOutput) as f: | 
| 
 
c5c1faa433e4
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/mmvc/ commit b'f226ae444ad694bb8d4f8faa80ed541b41cbc993\n'
 
rdvelazquez 
parents:  
diff
changeset
 | 
12 data = json.load(f) | 
| 
 
c5c1faa433e4
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/mmvc/ commit b'f226ae444ad694bb8d4f8faa80ed541b41cbc993\n'
 
rdvelazquez 
parents:  
diff
changeset
 | 
13 consensus = data["consensus"] | 
| 
3
 
4ca91bf87f1f
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/mmvc/ commit b'35710f754cc8be97c423c95eb91197df01d0520d\n'-dirty
 
rdvelazquez 
parents: 
0 
diff
changeset
 | 
14 consensusFasta = ">" + fastaHeaderName + "\n" + consensus + "\n" | 
| 
0
 
c5c1faa433e4
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/mmvc/ commit b'f226ae444ad694bb8d4f8faa80ed541b41cbc993\n'
 
rdvelazquez 
parents:  
diff
changeset
 | 
15 f = open(fastaFilePath, "w") | 
| 
 
c5c1faa433e4
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/mmvc/ commit b'f226ae444ad694bb8d4f8faa80ed541b41cbc993\n'
 
rdvelazquez 
parents:  
diff
changeset
 | 
16 f.write(consensusFasta) | 
