Automated liftover of NM_001261456.1:c.1762A>G (rs509749) to NM_001261457.1 via GRCh37

Automatically project variant from one transcript to another via common reference. http://www.ncbi.nlm.nih.gov/projects/SNP/snp_ref.cgi?rs=509749



In [1]:

    
import hgvs.parser
hgvsparser = hgvs.parser.Parser()
var_c1 = hgvsparser.parse_hgvs_variant('NM_001261456.1:c.1762A>G')



In [2]:

    
import hgvs.dataproviders.uta
hdp = hgvs.dataproviders.uta.connect()



In [3]:

    
import hgvs.projector
pj = hgvs.projector.Projector(hdp=hdp,
                              alt_ac='NC_000001.10',
                              src_ac=var_c1.ac,
                              dst_ac='NM_001261457.1')



In [4]:

    
pj.project_variant_forward(var_c1)









    Out[4]:





SequenceVariant(ac=NM_001261457.1, type=c, posedit=1534A>G)



In [4]: