<div dir="ltr"><div><div><div><div><div><div><div>I'm writing a plugin and am running into an issue with strand. The strand I get from Bio::EnsEMBL::Feature::strand doesn't match the strand I'm seeing in the output or the strand that is in the line hash provided to my run method. Specifically I'm having an issue with the negative strand.<br><br></div>This is my test plugin:<br><br>package StrandTest;<br>use strict;<br>use warnings;<br>use Bio::EnsEMBL::Variation::Utils::BaseVepPlugin;<br>use base qw(Bio::EnsEMBL::Variation::Utils::BaseVepPlugin);<br><br>sub run {<br>    my ($self, $tva, $line_hash) = @_;<br><br>    my $strand = $tva->variation_feature->strand;<br>    my $extra_strand = $line_hash->{Extra}->{STRAND};<br><br>    print("Strand: $strand/$extra_strand\n");<br><br>    return {};<br>}<br><br>1;<br><br></div>My output shows this: "Strand: 1/-1"<br><br></div>The strand written to my output file is -1, and the strand in the line hash is -1, but the strand from the VariationFeature is 1.<br><br></div>My VCF:<br><br>##fileformat=VCFv4.0<br>#CHROM    POS    ID    REF    ALT    QUAL    FILTER    INFO    FORMAT    PAT<br>6    133078809    .    C    A    60    PASS    .    GT:VR:RR:DP:GQ    0/1:12:15:27:.<br><br></div>I'm running VEP with these options: --plugin StrandTest --offline --vcf --force_overwrite --pick -i test.vcf -o test.vep<br><br></div>Any idea what's going on? For the time being I'm going to use "$line_hash->{Extra}->{STRAND}" for my strand info. Is that safe to do? Are there going to be cases where that hash is undef?<br><br></div>Thanks,<br>Matt Wood<br><div><div><div><div><div><br><br></div></div></div></div></div></div>