ref_coord.py 文件源码

def READ_FASTA_ENTRY(file_name):
    fasta_sequences=[]
    sequence_name=[]
    full_names_dict={}

    sequences_dict={}


    if os.stat(file_name)[6]!=0: #not empty

        fh = open(file_name, "r")
        for record in SeqIO.parse(fh, "fasta"):
            short_name=str(record.id).split(' ')[0]
            sequences_dict[short_name]=str(record.seq)


    return sequences_dict



#------------------------------------------------------------------------------------------