Don't count sequence gaps towards length

elsiklab · Jul 26, 2017 · 531ce03 · 531ce03
1 parent 6ec95f8
commit 531ce03
Showing 1 changed file with 5 additions and 1 deletion.
diff --git a/js/Store/SeqFeature/GCContent.js b/js/Store/SeqFeature/GCContent.js
@@ -41,17 +41,21 @@ function (
                     var r = f ? residues[i] : residues.slice(i - hw, i + hw);
                     var nc = 0;
                     var ng = 0;
+                    var len = 0;
                     for (var j = 0; j < r.length; j++) {
                         if (r[j] === 'c' || r[j] === 'C') {
                             nc++;
                         } else if (r[j] === 'g' || r[j] === 'G') {
                             ng++;
                         }
+                        if (r[j] != 'N') {
+                            len++;
+                        }
                     }
                     var pos = query.start;
                     var score;
                     if (thisB.gcMode === 'content') {
-                        score = (ng + nc) / r.length;
+                        score = (ng + nc) / len;
                     } else if (thisB.gcMode === 'skew') {
                         score = (ng - nc) / (ng + nc);
                     }