private long[] getTagsInBits(TagsByTaxaByte aTBT, int tagIndex, int[] reDirect, int anchorTaxa) {
int lgPerSite = (anchorTaxa / 64) + 1;
long[] seq = new long[lgPerSite];
for (int j = 0; j < aTBT.getTaxaCount(); j++) {
if (reDirect[j] < 0) continue;
int index = reDirect[j] / 64;
int offset = reDirect[j] % 64;
if (aTBT.getReadCountForTagTaxon(tagIndex, j) > 0) { // reference alleles
seq[index] = seq[index] | (1L << offset);
}
}
return seq;
}
