import java.util.Collections;
import java.util.Enumeration;
import java.util.List;
+import java.util.ListIterator;
import java.util.Vector;
import com.stevesoft.pat.Regex;
}
/**
- * {@inheritDoc}
- */
- @Override
- public Range findPositions(int fromCol, int toCol)
- {
- if (cursor != null && cursor.sequence == this
- && cursor.token == changeCount)
- {
- return findPositions(fromCol, toCol, cursor);
- }
-
- /*
- * count residues before fromCol
- */
- int j = 0;
- int count = 0;
- int seqlen = sequence.length;
- while (j < fromCol && j < seqlen)
- {
- if (!Comparison.isGap(sequence[j]))
- {
- count++;
- }
- j++;
- }
-
- /*
- * find first and last residues between fromCol and toCol
- */
- int firstPos = 0;
- int lastPos = 0;
- boolean foundFirst = false;
-
- while (j <= toCol && j < seqlen)
- {
- if (!Comparison.isGap(sequence[j]))
- {
- count++;
- if (!foundFirst)
- {
- firstPos = count;
- foundFirst = true;
- }
- lastPos = count;
- }
- j++;
- }
-
- if (firstPos == 0)
- {
- /*
- * no residues in this range
- */
- return null;
- }
-
- /*
- * adjust for sequence start coordinate
- */
- firstPos += start - 1;
- lastPos += start - 1;
-
- return new Range(firstPos, lastPos);
- }
-
- /**
- * Returns the range of sequence positions included in the given alignment
- * position range. If no positions are included (the range is entirely gaps),
- * then returns null. The cursor parameter may provide a starting position in
- * the neighbourhood of the search (which may be left of, right of, or
- * overlapping the search region).
- *
- * @param fromCol
- * start column of region (0..)
- * @param toCol
- * end column of region (0..)
- * @param curs
- * @return
- */
- protected Range findPositions(int fromCol, int toCol, SequenceCursor curs)
- {
- if (!isValidCursor(curs))
- {
- /*
- * wrong or invalidated cursor, compute de novo
- */
- return findPositions(fromCol, toCol);
- }
-
- /*
- * keep this simple...first step from cursor to fromCol...
- */
- final int seqlen = sequence.length;
- int resNo = curs.residuePosition;
- int col = curs.columnPosition - 1; // from base 1 to base 0
- if (col != fromCol)
- {
- int delta = col > fromCol ? -1 : 1;
- while (col != fromCol && col >= 0 && col < seqlen)
- {
- if (!Comparison.isGap(sequence[col]))
- {
- resNo += delta;
- }
- col += delta;
- }
- }
-
- if (col < fromCol || col == seqlen)
- {
- /*
- * sequence lies to the left of the target region
- */
- return null;
- }
-
- /*
- * resNo is now the residue at fromCol (if not gapped), else the one
- * before it (if delta == 1), else the one after (if delta == -1);
- * we want the residue before fromCol
- */
- if (!Comparison.isGap(sequence[fromCol]))
- {
- resNo--;
- }
- else if (curs.columnPosition > fromCol)
- {
- resNo -= 2;
- }
-
- /*
- * now first and last residues between fromCol and toCol
- */
- int firstPos = 0;
- int lastPos = 0;
- boolean foundFirst = false;
-
- while (col <= toCol && col < seqlen)
- {
- if (!Comparison.isGap(sequence[col]))
- {
- resNo++;
- if (!foundFirst)
- {
- firstPos = resNo;
- foundFirst = true;
- }
- lastPos = resNo;
- }
- col++;
- }
-
- if (firstPos == 0)
- {
- /*
- * no residues in this range
- */
- return null;
- }
-
- return new Range(firstPos, lastPos);
- }
-
- /**
* Returns an int array where indices correspond to each residue in the
* sequence and the element value gives its position in the alignment
*
* {@inheritDoc}
*/
@Override
- public List<SequenceFeature> findFeatures(int from, int to,
+ public List<SequenceFeature> findFeatures(int fromColumn, int toColumn,
String... types)
{
+ int startPos = findPosition(fromColumn - 1); // convert base 1 to base 0
+ int endPos = findPosition(toColumn - 1);
+
+ List<SequenceFeature> result = new ArrayList<>();
if (datasetSequence != null)
{
- return datasetSequence.findFeatures(from, to, types);
+ result = datasetSequence.getFeatures().findFeatures(startPos, endPos,
+ types);
}
- return sequenceFeatureStore.findFeatures(from, to, types);
+ else
+ {
+ result = sequenceFeatureStore.findFeatures(startPos, endPos, types);
+ }
+
+ /*
+ * if the start or end column is gapped, startPos or endPos may be to the
+ * left or right, and we may have included adjacent or enclosing features;
+ * remove any that are not enclosing, non-contact features
+ */
+ if (endPos > this.end || Comparison.isGap(sequence[fromColumn - 1])
+ || Comparison.isGap(sequence[toColumn - 1]))
+ {
+ ListIterator<SequenceFeature> it = result.listIterator();
+ while (it.hasNext())
+ {
+ SequenceFeature sf = it.next();
+ int featureStartColumn = findIndex(sf.getBegin());
+ int featureEndColumn = findIndex(sf.getEnd());
+ boolean noOverlap = featureStartColumn > toColumn
+ || featureEndColumn < fromColumn;
+
+ /*
+ * reject an 'enclosing' feature if it is actually a contact feature
+ */
+ if (sf.isContactFeature() && featureStartColumn < fromColumn
+ && featureEndColumn > toColumn)
+ {
+ noOverlap = true;
+ }
+ if (noOverlap)
+ {
+ it.remove();
+ }
+ }
+ }
+
+ return result;
}
/**