INNER CODE UNIT · Python

seqs

wwood/singlem · singlem/main.py:38

def seqs(args):
    from singlem.sequence_classes import SeqReader as SingleMSeqReader
    from singlem.metagenome_otu_finder import MetagenomeOtuFinder

    if args.alignment_type == 'aa':
        is_protein_alignment = True
    elif args.alignment_type == 'dna':
        is_protein_alignment = False
    else:
        raise Exception("Unexpected alignment type '%s'" % args.alignment_type)

    # Read in the fasta Alignment
    protein_alignment = SingleMSeqReader().alignment_from_alignment_file(args.alignment)
    logging.info("Read in %i aligned protein sequences e.g. %s %s" % (
        len(protein_alignment),
        protein_alignment[0].name,
        protein_alignment[0].seq))

View source record →