001/* 002 * BioJava development code 003 * 004 * This code may be freely distributed and modified under the 005 * terms of the GNU Lesser General Public Licence. This should 006 * be distributed with the code. If you do not have a copy, 007 * see: 008 * 009 * http://www.gnu.org/copyleft/lesser.html 010 * 011 * Copyright for this code is held jointly by the individual 012 * authors. These should be listed in @author doc comments. 013 * 014 * For more information on the BioJava project and its aims, 015 * or to join the biojava-l mailing list, visit the home page 016 * at: 017 * 018 * http://www.biojava.org/ 019 * 020 */ 021package demo; 022 023import org.biojava.nbio.core.alignment.matrices.SubstitutionMatrixHelper; 024import org.biojava.nbio.core.alignment.template.SubstitutionMatrix; 025import org.biojava.nbio.core.sequence.compound.AminoAcidCompound; 026 027public class DemoLoadSubstMax { 028 029 030 public static void main (String[] args){ 031 // that's the PAM250 matrix (named a bit unclear in AAindex...) 032 String max2="DAYM780301"; 033 SubstitutionMatrix<AminoAcidCompound> substMax2 = SubstitutionMatrixHelper.getMatrixFromAAINDEX(max2); 034 System.out.println("PAM250 matrix: "+ substMax2); 035 036 // and here BLOSUM62... 037 String max3="HENS920102"; 038 SubstitutionMatrix<AminoAcidCompound> substMax3 = SubstitutionMatrixHelper.getMatrixFromAAINDEX(max3); 039 System.out.printf("%s matrix: %s", max3, substMax3); 040 System.out.println(); 041 042 // This one I developed a while ago to be optimised for the alignment of distantly related sequences 043 String matrixName4 = "PRLA000101"; 044 SubstitutionMatrix<AminoAcidCompound> substMax4 = SubstitutionMatrixHelper.getMatrixFromAAINDEX(matrixName4); 045 System.out.printf("%s matrix: %s", matrixName4, substMax4); 046 System.out.println(); 047 } 048}