Grab hg19 positions in dbSNP that have unknown status

dbsnp_unknown.py

def get_dbsnp_unknown():
    snvs = {}
    with open('/vol/Users/yk336/database/human/dbsnp/hg19/snp.txt') as f:
        for line in f:
            sp = line.strip().split('\t')
            chrom = sp[1]
            pos = sp[3]
            status = sp[12]
            snv_type = sp[11]
            if status == 'unknown' and snv_type == 'single':
                snvs[chrom + ':' + pos] = True
    return snvs