sort_allele_orientation.py 文件源码

def get_Short (gene):


    newFasta=''
    for allele in SeqIO.parse(gene, "fasta", generic_dna):      
        try: 
            translatedSequence,sequence=translateSeq(allele.seq)
            newFasta+=">"+allele.id+"\n"+str(sequence)+"\n"
        except Exception as e:
            print e
            pass

    with open(gene, "wb") as f:
        f.write(newFasta)


    return  True