Mercurial > repos > gianmarco_piccinno > project_rm
comparison codon_usage.py @ 31:3cb2af2435d3 draft
Uploaded
author | gianmarco_piccinno |
---|---|
date | Wed, 12 Dec 2018 08:48:22 -0500 |
parents | e85a1feaaf38 |
children | 6e6b1fd6eda1 |
comparison
equal
deleted
inserted
replaced
30:e12061ce9601 | 31:3cb2af2435d3 |
---|---|
7 from pprint import pprint | 7 from pprint import pprint |
8 import argparse as ap | 8 import argparse as ap |
9 import sys | 9 import sys |
10 import os | 10 import os |
11 import pandas as pd | 11 import pandas as pd |
12 #from BCBio import GFF | |
13 | 12 |
14 | 13 |
15 def read_input(data = "example.fna", type_ = "fasta"): | 14 def read_input(data = "example.fna", type_ = "fasta"): |
16 | 15 |
17 """ | 16 """ |
141 args = vars(parser.parse_args()) | 140 args = vars(parser.parse_args()) |
142 | 141 |
143 seqs = read_input(data=args['input'], type_=args['type']) | 142 seqs = read_input(data=args['input'], type_=args['type']) |
144 out = codon_usage(seqs, args['codonTable']) | 143 out = codon_usage(seqs, args['codonTable']) |
145 | 144 |
146 with open(args['output'], "w") as outf: | 145 with open(args['output']+".csv", "w") as outf: |
147 out["Table"].to_csv(outf, sep="\t", index_label=["AA", "Codon"]) | 146 out["Table"].to_csv(outf, sep="\t", index_label=["AA", "Codon"]) |
148 | 147 |
149 | 148 |