one exon missing

[EricSHo]

Hi,

I've just installed cruzdb in Python 2.7. I found that the first exon was missed from this gene A4GNT, NM_016161. Here is how I fetched the sequence:
g = Genome(db="hg19")
gene='A4GNT'
gene_obj = g.refGene.filter_by(name2=gene).first()
print(len(gene_obj.mrna_sequence))

got 2, but RefSeq shows 3 exons

print(len(gene_obj.mrna_sequence[0]))

1161

print(len(gene_obj.mrna_sequence[1]))

434

print(len(gene_obj.mrna_sequence[2]))

index out of range error

I cross checked NM_016161 in NCBI, it seems the 1st exon (176 bps long) was not included.

Any help appreciated.

Eric.

[brentp]

what does len(gene_obj.coding_exons) show?

[EricSHo]

len(gene_obj.coding_exons) is 2 instead of 3