Package org.moltools.apps.probemaker.seq

Examples of org.moltools.apps.probemaker.seq.Probe


    ProbeDesignTask task = pm.getProbeDesignTask(false, true, null, new DefaultNamingScheme(), new AllAcceptor(),
        new NoSelector(), new DefaultProbeDesigner(new PairwiseTagAllocator()), new DefaultTSSPairDesigner(new DefaultTSSConstructor()), new StreamErrorHandler(System.out,true));

    task.run();

    Probe p = proj.getProbes().getSequenceAt(0);
    Probe p2 = proj.getProbes().getSequenceAt(1);

    //Move dummy to arm stage
    proj.getSettings().getCandidateModulesStage1().clear();
    proj.getSettings().addCandidateModuleStage2(dummy);
    task = pm.getProbeDesignTask(false, true, null, new DefaultNamingScheme(), new AllAcceptor(),
        new NoSelector(), new DefaultProbeDesigner(new PairwiseTagAllocator()), new DefaultTSSPairDesigner(new DefaultTSSConstructor()), new StreamErrorHandler(System.out,true));
    task.run();

    Probe p3 = proj.getProbes().getSequenceAt(0);
    Probe p4 = proj.getProbes().getSequenceAt(1);

//  Check that both probes of fair rank  
    assertTrue("Probe 1 (stage 1) expected to be of fair quality",ProbeMakerPropertyUtils.getRank(p)==ProbeMakerConstants.OK_QUALITY); //$NON-NLS-1$
    assertTrue("Probe 2 (stage 1) expected to be of fair quality",ProbeMakerPropertyUtils.getRank(p2)==ProbeMakerConstants.OK_QUALITY); //$NON-NLS-1$
View Full Code Here


    for (Iterator<Probe> pi = probes.iterator(); pi.hasNext(); ) {
      if (isAborted()) {
        return;
      }
      try {
        Probe p = pi.next();
        setStatus(CoreMessages.getString("ProbeAnalysisTask.MESSAGE_ANALYZING") + p.getName(),0); //$NON-NLS-1$
        reAnalyzeProbe(p);
        increaseProgress();
      }
      catch (Exception ex) {
        //Print the error
        handleError(ex);
      }
    }
    //Recalculate the probe qualities and finish up the work
    for (Iterator<Probe> i = probes.iterator(); i.hasNext(); ) {
      Probe p = i.next();
      Analyzer.calculateAndSetQuality(p);
    }
   
    elapsedTime = new Date().getTime() - startTime;
    setCompleted();
View Full Code Here

  }

  protected Element getProbeListElement(String name, DefiniteSequenceDB<Probe> probes, TagLibrary[] libraries) {
    Element probeList = new Element(name,ns);
    for (Iterator<Probe> i = probes.iterator(); i.hasNext(); ) {
      Probe pr = i.next();
      Element probeE = getProbeElement(PROBE,pr,libraries);
      probeList.addContent(probeE);
    }
    return probeList;
  }
View Full Code Here

          NucleotideSequence tag = getTag(subEs[pos]);
          tags.put(tag,null);
        }
      }     

      Probe p = new DefaultProbe(tssP,upstream,downstream,type);
      p.setID(id);
      p.setName(probeName);
      ProbeMakerPropertyUtils.setRank(p,rank);
      if (pg!=null)
        DesignUtils.addGroup(p,pg);

      //Go through and add tags
      int pos = 0;
      for (Iterator ti = tags.keySet().iterator();ti.hasNext();pos++) {
        NucleotideSequence t = (NucleotideSequence) ti.next();
        p.addTagAt(t,pos);       
        project.getTagAllocationTable().use(t,p,pos);       
      }
      tags.clear();

      addMessagesFromElement(p,probeE.getChild("messageList",ns));       //$NON-NLS-1$
View Full Code Here

        for (int i = 1; i < l.size() - 1; i++) {
          tags.add(new SimpleChangeableNucleotideSequence(id+"|Tag " + i,  //$NON-NLS-1$
              l.get(i), NucleotideSequence.DNA,  "Tag " + i)); //$NON-NLS-1$
        }
        PropertyAcceptorNucleotideSequence tss3 = ProbeMakerSequenceFactory.createTSS(new SimpleChangeableNucleotideSequence("3' TSS", l.get(l.size() - 1), NucleotideSequence.DNA,"3' TSS"),null); //$NON-NLS-1$
        Probe p = new DefaultProbe(ProbeMakerSequenceFactory.createTSSPair(tss5, tss3, null),0,0,NucleotideSequence.DNA);
        for (int i = 0;i<tags.size();i++) {
          p.addTagAt(tags.get(i),i)
        }
       
        p.setName(name);
        p.setID(id);
        p.setTagsAllocated(true);
        return p;
      }
    };
  }
View Full Code Here

    PropertyAcceptorNucleotideSequence tss5 = ProbeMakerSequenceFactory.createTSS(seq,target.getID());   
    ProbeMakerPropertyUtils.addMessage(tss5, new Message("Test",new MessageType("test_type","Test type"),Message.ERROR));     //$NON-NLS-1$ //$NON-NLS-2$ //$NON-NLS-3$
    TSSPair pair = new TSSPair(tss5,null,target);
    ProbeMakerPropertyUtils.addMessage(pair, new Message("Test",new MessageType("test_type","Test type"),Message.ERROR)); //$NON-NLS-1$ //$NON-NLS-2$ //$NON-NLS-3$

    Probe probe = new DefaultProbe(pair,0,0,NucleotideSequence.DNA);

    ProbeMakerPropertyUtils.setRank(probe, ProbeMakerConstants.BAD_QUALITY);

    assertTrue("Incorrect rank", ProbeMakerConstants.BAD_QUALITY == ProbeMakerPropertyUtils.getRank(probe)); //$NON-NLS-1$

    assertTrue("Incorrect target ID", target.getID().equals(DesignUtils.getTargetID(probe))); //$NON-NLS-1$
    assertTrue("Incorrect target", target.equals(pair.getTarget())); //$NON-NLS-1$
    assertTrue("Incorrect target", target.equals(probe.getTarget())); //$NON-NLS-1$
    assertTrue("Incorrect message count", 2 == ProbeMakerPropertyUtils.getAllMessages(probe).size()); //$NON-NLS-1$
    assertTrue("Wrong tss pair",pair == probe.getTSSPair()); //$NON-NLS-1$

    ProbeMakerPropertyUtils.addMessage(probe, new Message("Test",new MessageType("test_type","Test type"),Message.ERROR)); //$NON-NLS-1$ //$NON-NLS-2$ //$NON-NLS-3$

    assertTrue("Incorrect message count", 3 == ProbeMakerPropertyUtils.getAllMessages(probe).size()); //$NON-NLS-1$
View Full Code Here

  }

  public void testGroup() {
   
    Probe p1 = new DefaultProbe(NucleotideSequence.DNA);
    p1.setID("P1");
   
    TargetGroup tg = new TargetGroup("tg");
    ProbeGroup pg = new ProbeGroup(tg, "PG");
    assertEquals(0, pg.getMembers().size());
    pg.addMember(p1);
    assertEquals(1, pg.getMembers().size());
    p1.setName("Probe 1");   
    pg.removeMember(p1);
    assertEquals(0, pg.getMembers().size());
 
View Full Code Here

    setStatus("Writing probes",0);       //$NON-NLS-1$

    //Make a list of probe groups
    List<ProbeGroup> p_groups = new ArrayList<ProbeGroup>();
    for (Iterator<Probe> i = proj.getProbes().iterator(); i.hasNext(); ) {
      Probe pr = i.next();
      ProbeGroup pg = (ProbeGroup) DesignUtils.getGroupWithType(pr,ProbeGroup.GROUP_TYPE);
      if (pg != null) {
        if (!p_groups.contains(pg)) {
          p_groups.add(pg);
        }
      }
    }

    //Print the probe groups - targetGroup;groupTagLibrary;groupTagNumber;
    p.println(TITLE_DESCRIPTOR + PROBE_GROUPS);
    for (Iterator<ProbeGroup> i = p_groups.iterator(); i.hasNext(); ) {
      ProbeGroup pg = i.next();
      TargetGroup tg = pg.getTargetGroup();
      int tg_no = tg == null ? -1 : groups.indexOf(tg);
      int lib_no = -1;
      int tag_no = -1;
      String tagID = proj.getTagAllocationTable().getGroupTagID(pg);
      if (tagID != null) {
        for (int l = 0; l < libraries.length; l++) {
          try {
            NucleotideSequence s = libraries[l].getSequenceByID(tagID);
            lib_no = l;
            tag_no = libraries[l].indexOfSequence(s);
          }
          catch (IDNotFoundException inx) {
            //do nothing
          }
        }
      }
      p.println(tio.getOutString(new String[] { String.valueOf(tg_no),String.valueOf(lib_no),String.valueOf(tag_no) }));
    }
    p.println(TITLE_DESCRIPTOR + LISTEND);

    //Print the probes - description;ID;type;upStream;downStream;target;group number;quality;messageCountArms;messageCountProbe;;5'seq;Tm;3'seq;Tm;tag 1 lib;tag 1 pos;...;...;tag n lib;tag n pos;
    p.println(TITLE_DESCRIPTOR + PROBES);
    DefiniteSequenceDB<Probe> probes = proj.getProbes();
    for (Iterator<Probe> i = probes.iterator(); i.hasNext(); ) {
      Probe pr = i.next();
      PropertyAcceptorNucleotideSequence fivePrime;
      PropertyAcceptorNucleotideSequence threePrime;
      try {
        fivePrime = (PropertyAcceptorNucleotideSequence) pr.getTSSPair().getSequence(TSSPair.KEY_FIVE_PRIME);
      }
      catch (ClusterException e) {
        fivePrime = null;
      }
      try {
        threePrime = (PropertyAcceptorNucleotideSequence) pr.getTSSPair().getSequence(TSSPair.KEY_THREE_PRIME);
      }
      catch (ClusterException e) {
        threePrime = null;
      }

      int targetNo = -1;
      try {
        targetNo = pr.getTarget() == null ? -1 : proj.getTargets().indexOfSequence(proj.getTargets().getSequenceByID(pr.getTarget().getID()));
      }
      catch (IDNotFoundException ix) {
        //let targetNo be -1
      }
      Group<?> pg = DesignUtils.getGroupWithType(pr,ProbeGroup.GROUP_TYPE);

      String outStr = tio.getOutString(new String[] {
          pr.getName().replace(VALUE_SEPARATOR, VALUE_SEPARATOR_REPLACEMENT),
          pr.getID() != null ? pr.getID().replace(VALUE_SEPARATOR, VALUE_SEPARATOR_REPLACEMENT) : "No ID",                            //$NON-NLS-1$
              String.valueOf(pr.getType()),
              String.valueOf(pr.getUpstreamBlockCount()),
              String.valueOf(pr.getDownstreamBlockCount()),
              String.valueOf(targetNo),
              String.valueOf(pg == null ? -1 : p_groups.indexOf(pg)),
              String.valueOf(ProbeMakerPropertyUtils.getRank(pr)),
              String.valueOf(fivePrime == null ? 0 : ProbeMakerPropertyUtils.getAllMessages(fivePrime).size()),
              String.valueOf(threePrime == null ? 0 : ProbeMakerPropertyUtils.getAllMessages(threePrime).size()),
              String.valueOf(ProbeMakerPropertyUtils.getOwnMessages(pr.getTSSPair()).size()),
              String.valueOf(ProbeMakerPropertyUtils.getOwnMessages(pr).size()),
              fivePrime == null ? "" : fivePrime.seqString(), //$NON-NLS-1$
                  String.valueOf(fivePrime == null ? 0 : ProbeMakerPropertyUtils.getHybridizationTemp(fivePrime)),
                  threePrime == null ? "" : threePrime.seqString(), //$NON-NLS-1$
                      String.valueOf(threePrime == null ? 0 : ProbeMakerPropertyUtils.getHybridizationTemp(threePrime))
      });
      for (int t = 0; t < pr.getTags().size(); t++) {
        outStr += VALUE_SEPARATOR;
        NucleotideSequence tag = pr.getTagAt(t);
        //This isn't perfect but will work (if noone calls their tags 'Spacer')
        if (tag.getID().startsWith("Spacer")) { //$NON-NLS-1$
          outStr += "Spacer"; //$NON-NLS-1$
          outStr += VALUE_SEPARATOR;
          outStr += tag.seqString();
        }
        else {
          int lib_no = -1;
          int tag_no = -1;
          for (int l = 0; l < libraries.length; l++) {
            try {
              NucleotideSequence s = libraries[l].getSequenceByID(tag.getID());
              if (s!=null) {
                lib_no = l;
                tag_no = libraries[l].indexOfSequence(s);
              }
            }
            catch (IDNotFoundException inx) {
              //do nothing
            }
          }
          outStr += String.valueOf(lib_no);
          outStr += VALUE_SEPARATOR;
          outStr += String.valueOf(tag_no);
        }
      }
      p.println(outStr);

      if (fivePrime != null) {
        for (Iterator<Message> mi = ProbeMakerPropertyUtils.getAllMessages(fivePrime).iterator();mi.hasNext();) {
          Message m = mi.next();
          MessageType mt = m.getType();
          p.println(tio.getOutString(new String[] {mt.getType(),mt.getDescription(),m.getMessage(),String.valueOf(m.getSeverity())}));
        }
      }
      if (threePrime != null) {
        for (Iterator<Message> mi = ProbeMakerPropertyUtils.getAllMessages(threePrime).iterator();mi.hasNext();) {
          Message m = mi.next();
          MessageType mt = m.getType();
          p.println(tio.getOutString(new String[] {mt.getType(),mt.getDescription(),m.getMessage(),String.valueOf(m.getSeverity())}));
        }
      }
      for (Iterator<Message> mi = ProbeMakerPropertyUtils.getOwnMessages(pr.getTSSPair()).iterator();mi.hasNext();) {
        Message m = mi.next();
        MessageType mt = m.getType();
        p.println(tio.getOutString(new String[] {mt.getType(),mt.getDescription(),m.getMessage(),String.valueOf(m.getSeverity())}));
      }
      for (Iterator<Message> mi = ProbeMakerPropertyUtils.getOwnMessages(pr).iterator();mi.hasNext();) {
View Full Code Here

        return false;
      }

      setStatus("Selecting probe"); //$NON-NLS-1$
      htmllog.printLine("Selecting probe"); //$NON-NLS-1$
      Probe selected = analyzer.selectProbe(tempCandidates);
      if (selected == null) {
        //Nothing
      }
      else {
        //We also have to set the selected probe's group to be the correct one.
View Full Code Here

    m.setData(DINAMeltSecStructModule.KEY_NA_CONC,new Double(0.01));
    m.setData(DINAMeltSecStructModule.KEY_DG_LIMIT,new Double(-7.5));
    m.setData(DINAMeltSecStructModule.KEY_DG_SPAN,new Double(2));
   
    new ProbeAnalysisTask(project,null).run();
    Probe p1 = project.getProbes().getSequenceAt(0);
    assertEquals(ProbeMakerConstants.OK_QUALITY,ProbeMakerPropertyUtils.getRank(p1));
    m.setData(DINAMeltSecStructModule.KEY_NA_CONC,new Double(1));
    new ProbeAnalysisTask(project,null).run();
    p1 = project.getProbes().getSequenceAt(0);
    assertEquals(ProbeMakerConstants.BAD_QUALITY,ProbeMakerPropertyUtils.getRank(p1));
View Full Code Here

TOP

Related Classes of org.moltools.apps.probemaker.seq.Probe

Copyright © 2018 www.massapicom. All rights reserved.
All source code are property of their respective owners. Java is a trademark of Sun Microsystems, Inc and owned by ORACLE Inc. Contact coftware#gmail.com.