rnn_char.py 文件源码

def dna_from_reference(chrom='21'):
    reference_hg19 = '/dsde/data/deep/vqsr/Homo_sapiens_assembly19.fasta'
    record_dict = SeqIO.to_dict(SeqIO.parse(reference_hg19, "fasta"))
    dna = str(record_dict[chrom].seq[35000000:37000000])
    return dna