byte base = s.getBase();
freqs[refBase][base]++;
}
}
h.substitutionMatrix = new SubstitutionMatrix(freqs);
} else h.substitutionMatrix = substitutionMatrix ;
HuffmanParamsCalculator calculator = new HuffmanParamsCalculator();
for (CramRecord r : records)
if (r.getReadFeatures() == null)