src/api/BamReader.h

   1 // ***************************************************************************\r
   2 // BamReader.h (c) 2009 Derek Barnett, Michael Str�mberg\r
   3 // Marth Lab, Department of Biology, Boston College\r
   4 // All rights reserved.\r
   5 // ---------------------------------------------------------------------------\r
   6 // Last modified: 9 July 2010 (DB)\r
   7 // ---------------------------------------------------------------------------\r
   8 // Uses BGZF routines were adapted from the bgzf.c code developed at the Broad\r
   9 // Institute.\r
  10 // ---------------------------------------------------------------------------\r
  11 // Provides the basic functionality for reading BAM files\r
  12 // ***************************************************************************\r
  13 \r
  14 #ifndef BAMREADER_H\r
  15 #define BAMREADER_H\r
  16 \r
  17 // C++ includes\r
  18 #include <string>\r
  19 \r
  20 // BamTools includes\r
  21 #include "BamAux.h"\r
  22 \r
  23 namespace BamTools {\r
  24   \r
  25 class BamReader {\r
  26 \r
  27     // constructor / destructor\r
  28     public:\r
  29         BamReader(void);\r
  30         ~BamReader(void);\r
  31 \r
  32     // public interface\r
  33     public:\r
  34 \r
  35         // ----------------------\r
  36         // BAM file operations\r
  37         // ----------------------\r
  38 \r
  39         // close BAM file\r
  40         void Close(void);\r
  41         // returns whether reader is open for reading or not\r
  42         bool IsOpen(void) const;\r
  43         // performs random-access jump to reference, position\r
  44         bool Jump(int refID, int position = 0);\r
  45         // opens BAM file (and optional BAM index file, if provided)\r
  46         bool Open(const std::string& filename, const std::string& indexFilename = "");\r
  47         // returns file pointer to beginning of alignments\r
  48         bool Rewind(void);\r
  49         // sets a region of interest (with left & right bound reference/position)\r
  50         // attempts a Jump() to left bound as well\r
  51         // returns success/failure of Jump()\r
  52         bool SetRegion(const BamRegion& region);\r
  53         bool SetRegion(const int& leftRefID, const int& leftBound, const int& rightRefID, const int& rightBound);\r
  54 \r
  55         // ----------------------\r
  56         // access alignment data\r
  57         // ----------------------\r
  58 \r
  59         // retrieves next available alignment (returns success/fail)\r
  60         bool GetNextAlignment(BamAlignment& bAlignment);\r
  61         \r
  62         // retrieves next available alignment core data (returns success/fail)\r
  63         // ** DOES NOT parse any character data (read name, bases, qualities, tag data)\r
  64         //    these can be accessed, if necessary, from the supportData \r
  65         // useful for operations requiring ONLY positional or other alignment-related information\r
  66         bool GetNextAlignmentCore(BamAlignment& bAlignment);\r
  67 \r
  68         // ----------------------\r
  69         // access auxiliary data\r
  70         // ----------------------\r
  71 \r
  72         // returns SAM header text\r
  73         const std::string GetHeaderText(void) const;\r
  74         // returns number of reference sequences\r
  75         int GetReferenceCount(void) const;\r
  76         // returns vector of reference objects\r
  77         const BamTools::RefVector& GetReferenceData(void) const;\r
  78         // returns reference id (used for BamReader::Jump()) for the given reference name\r
  79         int GetReferenceID(const std::string& refName) const;\r
  80         // returns the name of the file associated with this BamReader\r
  81         const std::string GetFilename(void) const;\r
  82 \r
  83         // ----------------------\r
  84         // BAM index operations\r
  85         // ----------------------\r
  86 \r
  87         // creates index for BAM file, saves to file (default = bamFilename + ".bai")\r
  88         bool CreateIndex(bool useDefaultIndex = true);\r
  89         \r
  90     // private implementation\r
  91     private:\r
  92         struct BamReaderPrivate;\r
  93         BamReaderPrivate* d;\r
  94 };\r
  95 \r
  96 } // namespace BamTools\r
  97 \r
  98 #endif // BAMREADER_H\r