Skip to content
Snippets Groups Projects
Commit 5c5a11da authored by tjc's avatar tjc
Browse files

fix for finding synonymous and non-synonymous variants

git-svn-id: svn+ssh://svn.internal.sanger.ac.uk/repos/svn/pathsoft/artemis/trunk@15629 ee4ac58c-ac51-4696-9907-e4b3aa274f04
parent 1bb3c28e
No related branches found
No related tags found
No related merge requests found
......@@ -356,11 +356,9 @@ class VCFRecord
{
RangeVector ranges = feature.getLocation().getRanges();
intronlength = 0;
for(int j=0; j< ranges.size(); j++)
{
Range range = (Range) ranges.get(j);
if(j > 0)
{
if(feature.isForwardFeature())
......@@ -379,16 +377,14 @@ class VCFRecord
if(feature.isForwardFeature())
{
mod = (basePosition-feature.getRawFirstBase())%3;
codonStart = basePosition-feature.getRawFirstBase()-mod;
mod = (basePosition-feature.getRawFirstBase()-intronlength)%3;
codonStart = basePosition-feature.getRawFirstBase()-intronlength-mod;
}
else
{
mod = (feature.getRawLastBase()-basePosition)%3;
codonStart = feature.getRawLastBase()-basePosition-mod;
mod = (feature.getRawLastBase()-basePosition-intronlength)%3;
codonStart = feature.getRawLastBase()-basePosition-intronlength-mod;
}
codonStart-=intronlength;
try
{
......@@ -397,7 +393,6 @@ class VCFRecord
char codon[] = feature.getBases().substring(codonStart,
codonStart + 3).toLowerCase().toCharArray();
// String oldBase = new String(codon);
char aaRef = AminoAcidSequence.getCodonTranslation(codon[0],
codon[1], codon[2]);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment