Skip to content

Commit

Permalink
resolves "Exception: Unable to read reference sequence base past end …
Browse files Browse the repository at this point in the history
…of current cached sequence"

It's much easier to just load the entire sequence we're working on. This commit
does that.
  • Loading branch information
ekg committed Apr 23, 2015
1 parent 92eb53a commit c5efe19
Showing 1 changed file with 1 addition and 6 deletions.
7 changes: 1 addition & 6 deletions src/AlleleParser.cpp
Expand Up @@ -601,14 +601,9 @@ bool AlleleParser::loadNextPositionWithInputVariant(void) {

// alignment-based method for loading the first bit of our reference sequence
void AlleleParser::loadReferenceSequence(BamAlignment& alignment) {
DEBUG2("loading reference sequence overlapping first alignment");
loadReferenceSequence(referenceIDToName[alignment.RefID]);
currentPosition = alignment.Position;
currentSequenceStart = alignment.Position;
currentSequenceName = referenceIDToName[alignment.RefID];
currentRefID = alignment.RefID;
rightmostHaplotypeBasisAllelePosition = currentPosition;
DEBUG2("reference.getSubSequence("<< currentSequenceName << ", " << currentSequenceStart << ", " << alignment.AlignedBases.length() << ")");
currentSequence = uppercase(reference.getSubSequence(currentSequenceName, currentSequenceStart, alignment.Length));
}

void AlleleParser::loadReferenceSequence(string& seqname) {
Expand Down

0 comments on commit c5efe19

Please sign in to comment.