Commit 571d5609 authored by Erika Sallet's avatar Erika Sallet
Browse files

In the Est plugin, if EST are incompatible, for the same number of gaps...

In the Est plugin, if EST are incompatible,  for the same number of gaps select first the smallest covered region.
parent 5a8b25ab
Pipeline #847 failed with stages
in 43 seconds
-------------------------------------------------------------------------------
v4.2a released Dec. 2016
-------------------------------------------------------------------------------
- new plugin 'SignalWAM' to predict spliced|start signals with weight-array matrix
- In the Est plugin, change the sorting function: for the same number of gaps, select first the smallest covered region.
-------------------------------------------------------------------------------
v4.1c released Sept. 2014
-------------------------------------------------------------------------------
Improve GFF3 output format: remove ending tab.
-------------------------------------------------------------------------------
v4.1c released Janv. 2014
-------------------------------------------------------------------------------
......
......@@ -4,8 +4,8 @@ Reading splice site file (NetGene2)...........forward, reverse done
Loading Sensor.Est............done
Reading cDNA hits............18 sequences read
[AV533768]: inconsistent hit [3571-3657]
[AI992875]: no matching splice site
[AI100480]: inconsistent hit [4201-4301]
[AI992875]: no matching splice site
A total of 3/18 sequences rejected
Reading cDNA hits............2 sequences read
Optimal path length = 614.2795
This diff is collapsed.
......@@ -40,12 +40,12 @@ int HitsCompare(const void *A, const void *B)
return -1;
if ((*UA)->NGaps < (*UB)->NGaps)
return 1;
// select the largest region coverage
// select the smallest region coverage
int regionLgA = (*UA)->End - (*UA)->Start + 1;
int regionLgB = (*UB)->End - (*UB)->Start + 1;
if (regionLgA > regionLgB)
return -1;
if (regionLgA < regionLgB)
return -1;
if (regionLgA > regionLgB)
return 1;
if ((*UA)->Length > (*UB)->Length)
......
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