Commit 4e85cae6 authored by sallet's avatar sallet
Browse files

search an uncompatibility between an altest and a gene allowing a margin

parent 8ca00e03
......@@ -550,6 +550,7 @@ int main (int argc, char * argv [])
std::vector <Prediction*> vPred;
int ExonBorderMatchThreshold = PAR.getI("AltEst.ExonBorderMatchThreshold");
int RepredictMargin = PAR.getI("AltEst.RepredictMargin");
int ExonBorderMatchThreshold2 = PAR.getI("AltEst.IncompatibilityExonBorderMatchThreshold");
Prediction* AltPred;
Gene* uncompatibleGene;
......@@ -560,7 +561,7 @@ int main (int argc, char * argv [])
if (AltEstDB->voae_AltEst[altidx].IsToRemove()) continue;
// Look for an overlapping with an inconsistent gene
uncompatibleGene = AltEstDB->voae_AltEst[altidx].GetUncompatibleGene(pred);
uncompatibleGene = AltEstDB->voae_AltEst[altidx].GetUncompatibleGene(pred, ExonBorderMatchThreshold2);
// skip if not reference gene found
if (uncompatibleGene == NULL) continue;
Supports Markdown
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