TempReadAlignment::TempReadAlignment(AlignmentPtr alignment,
char* path)
: _path(path)
{
alignment->close();
alignment->open(_path, true);
_alignment = alignment;
}
void GenomeUpdateTest::createCallBack(AlignmentPtr alignment)
{
hal_size_t alignmentSize = alignment->getNumGenomes();
CuAssertTrue(_testCase, alignmentSize == 0);
Genome* ancGenome = alignment->addRootGenome("AncGenome", 0);
vector<Sequence::Info> seqVec(1);
seqVec[0] = Sequence::Info("Sequence", 1000000, 5000, 700000);
ancGenome->setDimensions(seqVec);
alignment->close();
alignment->open(_createPath, false);
ancGenome = alignment->openGenome("AncGenome");
seqVec[0] = Sequence::Info("Sequence", 10000005, 14000, 2000001);
ancGenome->setDimensions(seqVec);
}