def main(argv): identifiers = files.read_line_of_words(argv[0]) for identifier in identifiers: write_to_file(identifier) print pairwise_local_alignment(identifiers[0], identifiers[1])
def main(argv): ids = files.read_line_of_words(argv[0]) Entrez.email = "*****@*****.**" handle = Entrez.efetch(db='nucleotide', id=ids, rettype='fasta') records = SeqIO.parse(handle, 'fasta') shortest = min(records, key=len) print shortest.format('fasta')
def main(argv): leaves = files.read_line_of_words(argv[0]) print '\n'.join('(%s)%s;' % (enum, leaves[0]) for enum in tree.enumerate_unrooted(leaves[1:]))
def main(argv): words = files.read_line_of_words(argv[0]) freq = arrays.frequency_table(words) print '\n'.join('%s %s' % (w, f) for w, f in freq.iteritems())
def main(argv): id1, id2 = files.read_line_of_words(argv[0]) print pairwise_global_alignment(id1, id2)