src/api/BamConstants.h

   1 // ***************************************************************************
   2 // BamConstants.h (c) 2011 Derek Barnett
   3 // Marth Lab, Department of Biology, Boston College
   4 // All rights reserved.
   5 // ---------------------------------------------------------------------------
   6 // Last modified: 5 April 2011 (DB)
   7 // ---------------------------------------------------------------------------
   8 // Provides basic constants for handling BAM files.
   9 // ***************************************************************************
  10
  11 #ifndef BAM_CONSTANTS_H
  12 #define BAM_CONSTANTS_H
  13
  14 #include <string>
  15
  16 /*! \namespace BamTools::Constants
  17     \brief Provides basic constants for handling BAM files.
  18 */
  19
  20 namespace BamTools {
  21 namespace Constants {
  22
  23 const int BAM_SIZEOF_INT = 4;
  24
  25 // header magic number
  26 const char* const  BAM_HEADER_MAGIC = "BAM\1";
  27 const unsigned int BAM_HEADER_MAGIC_LENGTH = 4;
  28
  29 // BAM alignment core size
  30 const int BAM_CORE_SIZE = 32;
  31 const int BAM_CORE_BUFFER_SIZE = 8;
  32
  33 // BAM alignment flags
  34 const int BAM_ALIGNMENT_PAIRED              = 1;
  35 const int BAM_ALIGNMENT_PROPER_PAIR         = 2;
  36 const int BAM_ALIGNMENT_UNMAPPED            = 4;
  37 const int BAM_ALIGNMENT_MATE_UNMAPPED       = 8;
  38 const int BAM_ALIGNMENT_REVERSE_STRAND      = 16;
  39 const int BAM_ALIGNMENT_MATE_REVERSE_STRAND = 32;
  40 const int BAM_ALIGNMENT_READ_1              = 64;
  41 const int BAM_ALIGNMENT_READ_2              = 128;
  42 const int BAM_ALIGNMENT_SECONDARY           = 256;
  43 const int BAM_ALIGNMENT_QC_FAILED           = 512;
  44 const int BAM_ALIGNMENT_DUPLICATE           = 1024;
  45
  46 // CIGAR constants
  47 const char* const BAM_CIGAR_LOOKUP = "MIDNSHP";
  48 const int BAM_CIGAR_MATCH    = 0;
  49 const int BAM_CIGAR_INS      = 1;
  50 const int BAM_CIGAR_DEL      = 2;
  51 const int BAM_CIGAR_REFSKIP  = 3;
  52 const int BAM_CIGAR_SOFTCLIP = 4;
  53 const int BAM_CIGAR_HARDCLIP = 5;
  54 const int BAM_CIGAR_PAD      = 6;
  55
  56 // TODO: Add support for 'X' and '=' ops
  57
  58 const char BAM_CIGAR_MATCH_CHAR    = 'M';
  59 const char BAM_CIGAR_INS_CHAR      = 'I';
  60 const char BAM_CIGAR_DEL_CHAR      = 'D';
  61 const char BAM_CIGAR_REFSKIP_CHAR  = 'N';
  62 const char BAM_CIGAR_SOFTCLIP_CHAR = 'S';
  63 const char BAM_CIGAR_HARDCLIP_CHAR = 'H';
  64 const char BAM_CIGAR_PAD_CHAR      = 'P';
  65
  66 const int BAM_CIGAR_SHIFT    = 4;
  67 const int BAM_CIGAR_MASK     = ((1 << BAM_CIGAR_SHIFT) - 1);
  68
  69 // BAM tag types
  70 const char BAM_TAG_TYPE_ASCII  = 'A';
  71 const char BAM_TAG_TYPE_UINT8  = 'c';
  72 const char BAM_TAG_TYPE_INT8   = 'C';
  73 const char BAM_TAG_TYPE_UINT16 = 's';
  74 const char BAM_TAG_TYPE_INT16  = 'S';
  75 const char BAM_TAG_TYPE_UINT32 = 'i';
  76 const char BAM_TAG_TYPE_INT32  = 'I';
  77 const char BAM_TAG_TYPE_FLOAT  = 'f';
  78 const char BAM_TAG_TYPE_STRING = 'Z';
  79 const char BAM_TAG_TYPE_HEX    = 'H';
  80
  81 const size_t BAM_TAG_TAGSIZE  = 2;
  82 const size_t BAM_TAG_TYPESIZE = 1;
  83
  84 // DNA bases
  85 const char* const BAM_DNA_LOOKUP = "=ACMGRSVTWYHKDBN";
  86 const unsigned char BAM_BASECODE_EQUAL = 0;
  87 const unsigned char BAM_BASECODE_A     = 1;
  88 const unsigned char BAM_BASECODE_C     = 2;
  89 const unsigned char BAM_BASECODE_G     = 4;
  90 const unsigned char BAM_BASECODE_T     = 8;
  91 const unsigned char BAM_BASECODE_N     = 15;
  92
  93 const char BAM_DNA_EQUAL   = '=';
  94 const char BAM_DNA_A       = 'A';
  95 const char BAM_DNA_C       = 'C';
  96 const char BAM_DNA_G       = 'G';
  97 const char BAM_DNA_T       = 'T';
  98 const char BAM_DNA_N       = 'N';
  99 const char BAM_DNA_DEL     = '-';
 100 const char BAM_DNA_PAD     = '*';
 101
 102 // zlib constants
 103 const int GZIP_ID1   = 31;
 104 const int GZIP_ID2   = 139;
 105 const int CM_DEFLATE = 8;
 106 const int FLG_FEXTRA = 4;
 107 const int OS_UNKNOWN = 255;
 108 const int BGZF_XLEN  = 6;
 109 const int BGZF_ID1   = 66;
 110 const int BGZF_ID2   = 67;
 111 const int BGZF_LEN   = 2;
 112 const int GZIP_WINDOW_BITS    = -15;
 113 const int Z_DEFAULT_MEM_LEVEL = 8;
 114
 115 // BZGF constants
 116 const int BGZF_BLOCK_HEADER_LENGTH = 18;
 117 const int BGZF_BLOCK_FOOTER_LENGTH = 8;
 118 const int BGZF_MAX_BLOCK_SIZE      = 65536;
 119 const int BGZF_DEFAULT_BLOCK_SIZE  = 65536;
 120
 121 } // namespace Constants
 122 } // namespace BamTools
 123
 124 #endif // BAM_CONSTANTS_H