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