#!/usr/bin/env python

import IPython

import bdi.sources.uta1

import hgvs.parser
import hgvs.hgvsmapper
import hgvs.validator

alt_aln_method = 'splign'
hgvs_g = 'NC_000007.13:g.36561662C>T'
hgvs_c = 'NM_001637.3:c.1582G>A'
hgvs_p = 'NP_001628.1:p.(Gly528Arg)'

u1 = bdi.sources.uta1.connect()

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

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

var_g_to_c = hgvsmapper.hgvsg_to_hgvsc( var_g, var_c.ac )
var_g_to_r = hgvsmapper.hgvsg_to_hgvsr(var_g, var_c.ac, 'splign')
#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>
