Skip to content
Snippets Groups Projects
Commit d1430a51 authored by Graham McVicker's avatar Graham McVicker
Browse files

modified call to genomic2cdna in get_all_cdna_SNPs b/c slice arg should no longer be necessary

parent 65fdbccf
No related branches found
No related tags found
No related merge requests found
......@@ -13,8 +13,6 @@ Transcript - gene transcript object
=head1 SYNOPSIS
Give standard usage here
=head1 DESCRIPTION
Contains details of coordinates of all exons that make
......@@ -37,6 +35,8 @@ Manipulation:
Email questions to the ensembl developer mailing list <ensembl-dev@ebi.ac.uk>
=head1 METHODS
=cut
package Bio::EnsEMBL::Transcript;
......@@ -979,8 +979,8 @@ sub get_all_cdna_SNPs {
foreach my $type (@cdna_types) {
$snp_hash{$type} = [];
foreach my $snp (@{$all_snps->{$type}}) {
my @coords = $transcript->genomic2cdna($snp->start, $snp->end,
$snp->strand, $slice);
my @coords = $transcript->genomic2cdna($snp->start, $snp->end,
$snp->strand);
#skip snps that don't map cleanly (possibly an indel...)
if(scalar(@coords) != 1) {
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment