autoencoder.py 文件源码

def dna_from_reference(chrom='9'):
    #reference_hg19 = '/dsde/data/deep/vqsr/Homo_sapiens_assembly19.fasta'
    reference_hg19 = '/Users/sam/vqsr_data/Homo_sapiens_assembly19.fasta'

    record_dict = SeqIO.to_dict(SeqIO.parse(reference_hg19, "fasta"))
    dna = str(record_dict[chrom].seq[10000000:75000000])
    return dna