test/jalview/datamodel/CigarArrayTest.java

   1 /*
   2  * Jalview - A Sequence Alignment Editor and Viewer ($$Version-Rel$$)
   3  * Copyright (C) $$Year-Rel$$ The Jalview Authors
   4  *
   5  * This file is part of Jalview.
   6  *
   7  * Jalview is free software: you can redistribute it and/or
   8  * modify it under the terms of the GNU General Public License
   9  * as published by the Free Software Foundation, either version 3
  10  * of the License, or (at your option) any later version.
  11  *
  12  * Jalview is distributed in the hope that it will be useful, but
  13  * WITHOUT ANY WARRANTY; without even the implied warranty
  14  * of MERCHANTABILITY or FITNESS FOR A PARTICULAR
  15  * PURPOSE.  See the GNU General Public License for more details.
  16  *
  17  * You should have received a copy of the GNU General Public License
  18  * along with Jalview.  If not, see <http://www.gnu.org/licenses/>.
  19  * The Jalview Authors are detailed in the 'AUTHORS' file.
  20  */
  21 package jalview.datamodel;
  22
  23 import static org.testng.Assert.assertEquals;
  24
  25 import jalview.gui.JvOptionPane;
  26
  27 import org.testng.annotations.BeforeClass;
  28 import org.testng.annotations.Test;
  29
  30 public class CigarArrayTest
  31 {
  32   @BeforeClass(alwaysRun = true)
  33   public void setUpJvOptionPane()
  34   {
  35     JvOptionPane.setInteractiveMode(false);
  36     JvOptionPane.setMockResponse(JvOptionPane.CANCEL_OPTION);
  37   }
  38
  39   @Test(groups = "Functional")
  40   public void testConstructor()
  41   {
  42     SequenceI seq1 = new Sequence("sq1",
  43             "ASFDDABACBACBACBACBACBACBABCABCBACBABCAB");
  44     Sequence seq2 = new Sequence("sq2",
  45             "TTTTTTACBCBABCABCABCABCBACBACBABCABCABCBA");
  46
  47     // construct alignment
  48     AlignmentI al = new Alignment(new SequenceI[] { seq1, seq2 });
  49
  50     // hide columns
  51     HiddenColumns hc = new HiddenColumns();
  52     hc.hideColumns(3, 6);
  53     hc.hideColumns(16, 20);
  54
  55     // select group
  56     SequenceGroup sg1 = new SequenceGroup();
  57     sg1.addSequence(seq1, false);
  58     sg1.setStartRes(2);
  59     sg1.setEndRes(23);
  60
  61     // Cigar array meanings:
  62     // M = match
  63     // D = deletion
  64     // I = insertion
  65     // number preceding M/D/I is the number of residues which
  66     // match/are deleted/are inserted
  67     // In the CigarArray constructor only matches or deletions are created, as
  68     // we are comparing a sequence to its own subsequence (the group) + hidden
  69     // columns.
  70
  71     CigarArray cig = new CigarArray(al, hc, sg1);
  72     String result = cig.getCigarstring();
  73     assertEquals(result, "1M4D9M5D3M");
  74
  75     // group starts at hidden cols
  76     sg1.setStartRes(3);
  77     cig = new CigarArray(al, hc, sg1);
  78     result = cig.getCigarstring();
  79     assertEquals(result, "4D9M5D3M");
  80
  81     // group starts at last but 1 hidden col
  82     sg1.setStartRes(5);
  83     cig = new CigarArray(al, hc, sg1);
  84     result = cig.getCigarstring();
  85     assertEquals(result, "2D9M5D3M");
  86
  87     // group starts at last hidden col
  88     sg1.setStartRes(6);
  89     cig = new CigarArray(al, hc, sg1);
  90     result = cig.getCigarstring();
  91     assertEquals(result, "1D9M5D3M");
  92
  93     // group starts just after hidden region
  94     sg1.setStartRes(7);
  95     cig = new CigarArray(al, hc, sg1);
  96     result = cig.getCigarstring();
  97     assertEquals(result, "9M5D3M");
  98
  99     // group ends just before start of hidden region
 100     sg1.setStartRes(5);
 101     sg1.setEndRes(15);
 102     cig = new CigarArray(al, hc, sg1);
 103     result = cig.getCigarstring();
 104     assertEquals(result, "2D9M");
 105
 106     // group ends at start of hidden region
 107     sg1.setEndRes(16);
 108     cig = new CigarArray(al, hc, sg1);
 109     result = cig.getCigarstring();
 110     assertEquals(result, "2D9M1D");
 111
 112     // group ends 1 after start of hidden region
 113     sg1.setEndRes(17);
 114     cig = new CigarArray(al, hc, sg1);
 115     result = cig.getCigarstring();
 116     assertEquals(result, "2D9M2D");
 117
 118     // group ends at end of hidden region
 119     sg1.setEndRes(20);
 120     cig = new CigarArray(al, hc, sg1);
 121     result = cig.getCigarstring();
 122     assertEquals(result, "2D9M5D");
 123
 124     // group ends just after end of hidden region
 125     sg1.setEndRes(21);
 126     cig = new CigarArray(al, hc, sg1);
 127     result = cig.getCigarstring();
 128     assertEquals(result, "2D9M5D1M");
 129
 130     // group ends 2 after end of hidden region
 131     sg1.setEndRes(22);
 132     cig = new CigarArray(al, hc, sg1);
 133     result = cig.getCigarstring();
 134     assertEquals(result, "2D9M5D2M");
 135   }
 136 }