Ignore:
Timestamp:
Jan 26, 2013, 7:53:56 AM (10 years ago)
Author:
Peter
Message:

refs #745. Implement a workaround when bam_nt16_rev_table is not
available in -lbam.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/0.10-stable/yat/omic/BamRead.cc

    r2969 r2972  
    149149  std::string BamRead::sequence(void) const
    150150  {
    151     // FIXME bam_nt16_rev_table not available in old lbam
    152151    std::string result(sequence_length(), ' ');
     152#ifdef HAVE_BAM_NT16_REV_TABLE
    153153    for (size_t i=0; i<result.size(); ++i)
    154154      result[i] = bam_nt16_rev_table[sequence(i)];
     155#else
     156    std::string table = "=ACMGRSVTWYHKDBN";
     157    for (size_t i=0; i<result.size(); ++i)
     158      result[i] = table[sequence(i)];
     159#endif
    155160    return result;
    156161  }
Note: See TracChangeset for help on using the changeset viewer.