cn3d_blast.hpp
上传用户:yhdzpy8989
上传日期:2007-06-13
资源大小:13604k
文件大小:4k
源码类别:

生物技术

开发平台:

C/C++

  1. /*
  2.  * ===========================================================================
  3.  * PRODUCTION $Log: cn3d_blast.hpp,v $
  4.  * PRODUCTION Revision 1000.0  2003/10/29 18:30:40  gouriano
  5.  * PRODUCTION PRODUCTION: IMPORTED [ORIGINAL] Dev-tree R1.9
  6.  * PRODUCTION
  7.  * ===========================================================================
  8.  */
  9. /*  $Id: cn3d_blast.hpp,v 1000.0 2003/10/29 18:30:40 gouriano Exp $
  10. * ===========================================================================
  11. *
  12. *                            PUBLIC DOMAIN NOTICE
  13. *               National Center for Biotechnology Information
  14. *
  15. *  This software/database is a "United States Government Work" under the
  16. *  terms of the United States Copyright Act.  It was written as part of
  17. *  the author's official duties as a United States Government employee and
  18. *  thus cannot be copyrighted.  This software/database is freely available
  19. *  to the public for use. The National Library of Medicine and the U.S.
  20. *  Government have not placed any restriction on its use or reproduction.
  21. *
  22. *  Although all reasonable efforts have been taken to ensure the accuracy
  23. *  and reliability of the software and data, the NLM and the U.S.
  24. *  Government do not and cannot warrant the performance or results that
  25. *  may be obtained by using this software or data. The NLM and the U.S.
  26. *  Government disclaim all warranties, express or implied, including
  27. *  warranties of performance, merchantability or fitness for any particular
  28. *  purpose.
  29. *
  30. *  Please cite the author in any work or product based on this material.
  31. *
  32. * ===========================================================================
  33. *
  34. * Authors:  Paul Thiessen
  35. *
  36. * File Description:
  37. *      module for aligning with BLAST and related algorithms
  38. *
  39. * ===========================================================================
  40. */
  41. #ifndef CN3D_BLAST__HPP
  42. #define CN3D_BLAST__HPP
  43. #include <corelib/ncbistl.hpp>
  44. #include <list>
  45. BEGIN_SCOPE(Cn3D)
  46. class Sequence;
  47. class BlockMultipleAlignment;
  48. class BLASTer
  49. {
  50. public:
  51.     BLASTer(void) { }
  52.     ~BLASTer(void) { }
  53.     static const std::string BLASTResidues;
  54.     typedef std::list < BlockMultipleAlignment * > AlignmentList;
  55.     // creates new pairwise alignments (as two-row BlockMultipleAlignments), each of which has
  56.     // the given master and one of the given sequences. If the alignment algorithm fails to
  57.     // align the new sequence, it will include a null-alignment for that sequence. If usePSSM is
  58.     // true, will do BLAST/PSSM, otherwise, BlastTwoSequences. In both cases, alignment is restricted
  59.     // to locations based on the multiple and alignTo/From values.
  60.     void CreateNewPairwiseAlignmentsByBlast(const BlockMultipleAlignment *multiple,
  61.         const AlignmentList& toRealign, AlignmentList *newAlignments, bool usePSSM);
  62.     // calculate the self-hits for the given alignment - e.g., the score of each row when
  63.     // aligned with the multiple using BLAST/PSSM. The scores will be placed in each row, -1.0
  64.     // if no significant alignment is found.
  65.     void CalculateSelfHitScores(const BlockMultipleAlignment *multiple);
  66. };
  67. // utility functions
  68. extern int LookupBLASTResidueNumberFromCharacter(char r);
  69. extern int LookupBLASTResidueNumberFromThreaderResidueNumber(char r);
  70. END_SCOPE(Cn3D)
  71. #endif // CN3D_BLAST__HPP
  72. /*
  73. * ---------------------------------------------------------------------------
  74. * $Log: cn3d_blast.hpp,v $
  75. * Revision 1000.0  2003/10/29 18:30:40  gouriano
  76. * PRODUCTION: IMPORTED [ORIGINAL] Dev-tree R1.9
  77. *
  78. * Revision 1.9  2003/02/03 19:20:02  thiessen
  79. * format changes: move CVS Log to bottom of file, remove std:: from .cpp files, and use new diagnostic macros
  80. *
  81. * Revision 1.8  2003/01/22 14:47:30  thiessen
  82. * cache PSSM in BlockMultipleAlignment
  83. *
  84. * Revision 1.7  2002/08/01 12:51:36  thiessen
  85. * add E-value display to block aligner
  86. *
  87. * Revision 1.6  2002/06/13 13:32:39  thiessen
  88. * add self-hit calculation
  89. *
  90. * Revision 1.5  2002/05/22 17:17:09  thiessen
  91. * progress on BLAST interface ; change custom spin ctrl implementation
  92. *
  93. * Revision 1.4  2002/05/17 19:10:27  thiessen
  94. * preliminary range restriction for BLAST/PSSM
  95. *
  96. * Revision 1.3  2002/03/28 14:06:02  thiessen
  97. * preliminary BLAST/PSSM ; new CD startup style
  98. *
  99. * Revision 1.2  2001/09/27 15:36:47  thiessen
  100. * decouple sequence import and BLAST
  101. *
  102. * Revision 1.1  2001/09/18 03:09:38  thiessen
  103. * add preliminary sequence import pipeline
  104. *
  105. */