();
if (refs == null || entry == null)
{
return rfs;
}
for (int i = 0; i < refs.length; i++)
{
if (comparator.matches(entry, refs[i]))
{
rfs.add(refs[i]);
}
}
return rfs;
}
interface DbRefComp
{
public boolean matches(DBRefEntry refa, DBRefEntry refb);
}
/**
* match on all non-null fields in refa
*/
// TODO unused - remove?
public static DbRefComp matchNonNullonA = new DbRefComp()
{
@Override
public boolean matches(DBRefEntry refa, DBRefEntry refb)
{
if (refa.getSource() == null
|| DBRefUtils.getCanonicalName(refb.getSource()).equals(
DBRefUtils.getCanonicalName(refa.getSource())))
{
if (refa.getVersion() == null
|| refb.getVersion().equals(refa.getVersion()))
{
if (refa.getAccessionId() == null
|| refb.getAccessionId().equals(refa.getAccessionId()))
{
if (refa.getMap() == null
|| (refb.getMap() != null && refb.getMap().equals(
refa.getMap())))
{
return true;
}
}
}
}
return false;
}
};
/**
* either field is null or field matches for all of source, version, accession
* id and map.
*/
// TODO unused - remove?
public static DbRefComp matchEitherNonNull = new DbRefComp()
{
@Override
public boolean matches(DBRefEntry refa, DBRefEntry refb)
{
if (nullOrEqualSource(refa.getSource(), refb.getSource())
&& nullOrEqual(refa.getVersion(), refb.getVersion())
&& nullOrEqual(refa.getAccessionId(), refb.getAccessionId())
&& nullOrEqual(refa.getMap(), refb.getMap()))
{
return true;
}
return false;
}
};
/**
* accession ID and DB must be identical. Version is ignored. Map is either
* not defined or is a match (or is compatible?)
*/
// TODO unused - remove?
public static DbRefComp matchDbAndIdAndEitherMap = new DbRefComp()
{
@Override
public boolean matches(DBRefEntry refa, DBRefEntry refb)
{
if (refa.getSource() != null && refb.getSource() != null
&& DBRefUtils.getCanonicalName(refb.getSource()).equals(
DBRefUtils.getCanonicalName(refa.getSource())))
{
// We dont care about version
if (refa.getAccessionId() != null && refb.getAccessionId() != null
// FIXME should be && not || here?
|| refb.getAccessionId().equals(refa.getAccessionId()))
{
if ((refa.getMap() == null || refb.getMap() == null)
|| (refa.getMap() != null && refb.getMap() != null && refb
.getMap().equals(refa.getMap())))
{
return true;
}
}
}
return false;
}
};
/**
* accession ID and DB must be identical. Version is ignored. No map on either
* or map but no maplist on either or maplist of map on a is the complement of
* maplist of map on b.
*/
// TODO unused - remove?
public static DbRefComp matchDbAndIdAndComplementaryMapList = new DbRefComp()
{
@Override
public boolean matches(DBRefEntry refa, DBRefEntry refb)
{
if (refa.getSource() != null && refb.getSource() != null
&& DBRefUtils.getCanonicalName(refb.getSource()).equals(
DBRefUtils.getCanonicalName(refa.getSource())))
{
// We dont care about version
if (refa.getAccessionId() != null && refb.getAccessionId() != null
|| refb.getAccessionId().equals(refa.getAccessionId()))
{
if ((refa.getMap() == null && refb.getMap() == null)
|| (refa.getMap() != null && refb.getMap() != null))
{
if ((refb.getMap().getMap() == null && refa.getMap().getMap() == null)
|| (refb.getMap().getMap() != null
&& refa.getMap().getMap() != null && refb
.getMap().getMap().getInverse()
.equals(refa.getMap().getMap())))
{
return true;
}
}
}
}
return false;
}
};
/**
* accession ID and DB must be identical. Version is ignored. No map on both
* or or map but no maplist on either or maplist of map on a is equivalent to
* the maplist of map on b.
*/
// TODO unused - remove?
public static DbRefComp matchDbAndIdAndEquivalentMapList = new DbRefComp()
{
@Override
public boolean matches(DBRefEntry refa, DBRefEntry refb)
{
if (refa.getSource() != null && refb.getSource() != null
&& DBRefUtils.getCanonicalName(refb.getSource()).equals(
DBRefUtils.getCanonicalName(refa.getSource())))
{
// We dont care about version
// if ((refa.getVersion()==null || refb.getVersion()==null)
// || refb.getVersion().equals(refa.getVersion()))
// {
if (refa.getAccessionId() != null && refb.getAccessionId() != null
|| refb.getAccessionId().equals(refa.getAccessionId()))
{
if (refa.getMap() == null && refb.getMap() == null)
{
return true;
}
if (refa.getMap() != null
&& refb.getMap() != null
&& ((refb.getMap().getMap() == null && refa.getMap()
.getMap() == null) || (refb.getMap().getMap() != null
&& refa.getMap().getMap() != null && refb
.getMap().getMap().equals(refa.getMap().getMap()))))
{
return true;
}
}
}
return false;
}
};
/**
* accession ID and DB must be identical, or null on a. Version is ignored. No
* map on either or map but no maplist on either or maplist of map on a is
* equivalent to the maplist of map on b.
*/
public static DbRefComp matchDbAndIdAndEitherMapOrEquivalentMapList = new DbRefComp()
{
@Override
public boolean matches(DBRefEntry refa, DBRefEntry refb)
{
if (refa.getSource() != null && refb.getSource() != null
&& DBRefUtils.getCanonicalName(refb.getSource()).equals(
DBRefUtils.getCanonicalName(refa.getSource())))
{
// We dont care about version
if (refa.getAccessionId() == null
|| refa.getAccessionId().equals(refb.getAccessionId()))
{
if (refa.getMap() == null || refb.getMap() == null)
{
return true;
}
if ((refa.getMap() != null && refb.getMap() != null)
&& (refb.getMap().getMap() == null && refa.getMap()
.getMap() == null)
|| (refb.getMap().getMap() != null
&& refa.getMap().getMap() != null && (refb
.getMap().getMap().equals(refa.getMap().getMap()))))
{
return true;
}
}
}
return false;
}
};
/**
* accession ID only must be identical.
*/
public static DbRefComp matchId = new DbRefComp()
{
@Override
public boolean matches(DBRefEntry refa, DBRefEntry refb)
{
if (refa.getAccessionId() != null && refb.getAccessionId() != null
&& refb.getAccessionId().equals(refa.getAccessionId()))
{
return true;
}
return false;
}
};
/**
* Parses a DBRefEntry and adds it to the sequence, also a PDBEntry if the
* database is PDB.
*
* Used by file parsers to generate DBRefs from annotation within file (eg
* Stockholm)
*
* @param dbname
* @param version
* @param acn
* @param seq
* where to annotate with reference
* @return parsed version of entry that was added to seq (if any)
*/
public static DBRefEntry parseToDbRef(SequenceI seq, String dbname,
String version, String acn)
{
DBRefEntry ref = null;
if (dbname != null)
{
String locsrc = DBRefUtils.getCanonicalName(dbname);
if (locsrc.equals(DBRefSource.PDB))
{
/*
* Check for PFAM style stockhom PDB accession id citation e.g.
* "1WRI A; 7-80;"
*/
Regex r = new com.stevesoft.pat.Regex(
"([0-9][0-9A-Za-z]{3})\\s*(.?)\\s*;\\s*([0-9]+)-([0-9]+)");
if (r.search(acn.trim()))
{
String pdbid = r.stringMatched(1);
String chaincode = r.stringMatched(2);
if (chaincode == null)
{
chaincode = " ";
}
// String mapstart = r.stringMatched(3);
// String mapend = r.stringMatched(4);
if (chaincode.equals(" "))
{
chaincode = "_";
}
// construct pdb ref.
ref = new DBRefEntry(locsrc, version, pdbid + chaincode);
PDBEntry pdbr = new PDBEntry();
pdbr.setId(pdbid);
pdbr.setType(PDBEntry.Type.PDB);
pdbr.setProperty(new Hashtable());
pdbr.setChainCode(chaincode);
// pdbr.getProperty().put("CHAIN", chaincode);
seq.addPDBId(pdbr);
}
else
{
System.err.println("Malformed PDB DR line:" + acn);
}
}
else
{
// default:
ref = new DBRefEntry(locsrc, version, acn);
}
}
if (ref != null)
{
seq.addDBRef(ref);
}
return ref;
}
/**
* Returns true if either object is null, or they are equal
*
* @param o1
* @param o2
* @return
*/
public static boolean nullOrEqual(Object o1, Object o2)
{
if (o1 == null || o2 == null)
{
return true;
}
return o1.equals(o2);
}
/**
* canonicalise source string before comparing. null is always wildcard
*
* @param o1
* - null or source string to compare
* @param o2
* - null or source string to compare
* @return true if either o1 or o2 are null, or o1 equals o2 under
* DBRefUtils.getCanonicalName
* (o1).equals(DBRefUtils.getCanonicalName(o2))
*/
public static boolean nullOrEqualSource(String o1, String o2)
{
if (o1 == null || o2 == null)
{
return true;
}
return DBRefUtils.getCanonicalName(o1).equals(
DBRefUtils.getCanonicalName(o2));
}
/**
* Selects just the DNA or protein references from a set of references
*
* @param selectDna
* if true, select references to 'standard' DNA databases, else to
* 'standard' peptide databases
* @param refs
* a set of references to select from
* @return
*/
public static DBRefEntry[] selectDbRefs(boolean selectDna,
DBRefEntry[] refs)
{
return selectRefs(refs, selectDna ? DBRefSource.DNACODINGDBS
: DBRefSource.PROTEINDBS);
// could attempt to find other cross
// refs here - ie PDB xrefs
// (not dna, not protein seq)
}
/**
* Returns the (possibly empty) list of those supplied dbrefs which have the
* specified source database, with a case-insensitive match of source name
*
* @param dbRefs
* @param source
* @return
*/
public static List searchRefsForSource(DBRefEntry[] dbRefs,
String source)
{
List matches = new ArrayList();
if (dbRefs != null && source != null)
{
for (DBRefEntry dbref : dbRefs)
{
if (source.equalsIgnoreCase(dbref.getSource()))
{
matches.add(dbref);
}
}
}
return matches;
}
}