org.biojava.nbio.genome.homology.BlastHomologyHits Maven / Gradle / Ivy

Go to download

Show more of this group Show more artifacts with this name
Show all versions of biojava-genome Show documentation

There is a newer version: 7.1.3

/*
 *                    BioJava development code
 *
 * This code may be freely distributed and modified under the
 * terms of the GNU Lesser General Public Licence.  This should
 * be distributed with the code.  If you do not have a copy,
 * see:
 *
 *      http://www.gnu.org/copyleft/lesser.html
 *
 * Copyright for this code is held jointly by the individual
 * authors.  These should be listed in @author doc comments.
 *
 * For more information on the BioJava project and its aims,
 * or to join the biojava-l mailing list, visit the home page
 * at:
 *
 *      http://www.biojava.org/
 *
 */
package org.biojava.nbio.genome.homology;

import org.biojava.nbio.genome.query.BlastXMLQuery;

import java.io.File;
import java.util.ArrayList;
import java.util.LinkedHashMap;

/**
 *
 * @author Scooter Willis 
 */
public class BlastHomologyHits {

	static public LinkedHashMap> getMatches(File xmlBlastHits, double ecutoff) throws Exception {
		LinkedHashMap> homologyHits = new LinkedHashMap>();
		BlastXMLQuery blastXMLQuery = new BlastXMLQuery(xmlBlastHits.getAbsolutePath());
		LinkedHashMap> hits = blastXMLQuery.getHitsQueryDef(ecutoff);
		for (String accessionid : hits.keySet()) {
			String[] data = accessionid.split(" "); // deal with notes/comments in blast results
			String id = data[0];
			ArrayList uniprotProteinHits = hits.get(accessionid);
			homologyHits.put(id, uniprotProteinHits);

		}
		return homologyHits;
	}
}