Changeset 6141 for trunk/STAT/ST_quality.cxx
- Timestamp:
- 14/08/09 16:29:27 (3 years ago)
- Files:
-
- 1 modified
-
trunk/STAT/ST_quality.cxx (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
trunk/STAT/ST_quality.cxx
r5855 r6141 131 131 } 132 132 AWT_dna_base base = awt_dna_table.char_to_enum(source_sequence[pos]); 133 if (base != ST_UNKNOWN && base != ST_GAP) { // don t count gaps133 if (base != ST_UNKNOWN && base != ST_GAP) { // don't count gaps 134 134 double val = max / (0.0001 + vec->b[base]); 135 135 double log_val = log(val);
