{
for (int j = 0; j < alAnot.length; j++)
{
- if (ds.getSequenceFeatures() != null)
+
+ String key = type2id(alAnot[j].label);
+ boolean isrna = alAnot[j].isValidStruc();
+
+ if (isrna)
{
- feature = ds.getSequenceFeatures()[0].type;
+ // hardwire to secondary structure if there is RNA secondary
+ // structure on the annotation
+ key = "SS";
}
- // ?bug - feature may still have previous loop value
- String key = type2id(feature);
-
if (key == null)
{
+
continue;
}
out.append(new Format("%-" + maxid + "s").form("#=GR "
+ printId(s[i], jvSuffix) + " " + key + " "));
ann = alAnot[j].annotations;
- boolean isrna = alAnot[j].isValidStruc();
String seq = "";
for (int k = 0; k < ann.length; k++)
{