INNER CODE UNIT · Python

getPlacements

ncbi/dbsnp · tutorials/rsjson_allele_info_demo.py:53

def getPlacements(info):
    '''
    rs genomic positions
    '''
    rs['alleles'] = []  # holder for one or more variant alleles
    for alleleinfo in info:
        # has top level placement (ptlp) and assembly info
        if alleleinfo['is_ptlp'] and len(
                alleleinfo['placement_annot']['seq_id_traits_by_assembly']
                ) > 0:  # get genomic placement and alleles

            for a in alleleinfo['alleles']:
                spdi = a['allele']['spdi']
                if spdi['inserted_sequence'] == spdi['deleted_sequence']:
                    rs['alleles'].append({'allele': spdi['deleted_sequence']})
                    rs['seq_id'] = spdi['seq_id']
                    rs['position'] = spdi['position']
                else:  # spdi['inserted_sequence'] != spdi['deleted_sequence']:

View source record →