#!/usr/bin/env python

import IPython

import bdi.sources.uta0

import hgvs.parser
import hgvs.hgvsmapper


ref = 'GRCh37.p10'
hgvs_g = 'NC_000007.13:g.36561662C>T'
hgvs_c = 'NM_001637.3:c.1582G>A'

bdi_uta0 = bdi.sources.uta0.connect()

hgvsparser = hgvs.parser.Parser()
hgvsmapper = hgvs.hgvsmapper.HGVSMapper( bdi_uta0, cache_transcripts = True )

var_g = hgvsparser.parse_hgvs_variant(hgvs_g)
var_c = hgvsparser.parse_hgvs_variant(hgvs_c)

var_g_to_c = hgvsmapper.hgvsg_to_hgvsc( var_g, var_c.ac, ref )
var_c_to_g = hgvsmapper.hgvsc_to_hgvsg( var_c, ref )

IPython.embed()

## <LICENSE>
## Copyright 2014 HGVS Contributors (https://bitbucket.org/invitae/hgvs)
## 
## Licensed under the Apache License, Version 2.0 (the "License");
## you may not use this file except in compliance with the License.
## You may obtain a copy of the License at
## 
##     http://www.apache.org/licenses/LICENSE-2.0
## 
## Unless required by applicable law or agreed to in writing, software
## distributed under the License is distributed on an "AS IS" BASIS,
## WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
## See the License for the specific language governing permissions and
## limitations under the License.
## </LICENSE>
