[ensembl-dev] Genomic Sequence
Nick Fankhauser
lists at nyk.ch
Tue Feb 21 12:44:48 GMT 2012
Hi!
I want to fetch the gnomic sequence of genes as reference for aligning
sequencing reads. I'm using the following method:
my $gene_adaptor = $registry->get_adaptor( 'Human', 'Core', 'Gene' );
my $gene=$gene_adaptor->fetch_by_display_label($gene_symbol);
my $seq = $gene->slice->seq();
Now this retrieved a sequence of length 81195211 bp for TP53 as
gene_symbol, but the genomic size of that gene should be 19143 bp.
Is there anything wrong with this method and how could it be improved?
Nick
---
http://www.ccspmd.ethz.ch/people/fniklaus
More information about the Dev
mailing list