Changeset 3026
- Timestamp:
- Apr 6, 2013, 9:31:02 AM (10 years ago)
- Location:
- trunk
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/test/bam.cc
r2987 r3026 129 129 << " expected " << 8 << "\n"; 130 130 } 131 uint8_t q = bam.qual(24); 132 if (q != 72) { 133 suite.add(false); 134 suite.err() << "error: quality: " << static_cast<int>(q) 135 << " expected " << 72 << "\n"; 136 } 137 uint8_t new_q = 60; 138 bam.qual(24, new_q); 139 q = bam.qual(24); 140 if (q != new_q) { 141 suite.add(false); 142 suite.err() << "error: quality: " << static_cast<int>(q) 143 << " expected " << new_q << "\n"; 144 } 131 145 test_cigar(suite, bam, in.header()); 132 146 } -
trunk/yat/omic/BamRead.cc
r2993 r3026 166 166 167 167 168 void BamRead::qual(size_t i, uint8_t q) 169 { *(bam1_qual(bam_)+i)=q; } 170 171 168 172 std::string BamRead::sequence(void) const 169 173 { -
trunk/yat/omic/BamRead.h
r2993 r3026 172 172 173 173 /** 174 \brief set quality of a base 175 176 \param i base to modify 177 \param q new quality 178 179 \since New in yat 0.11 180 */ 181 void qual(size_t i, uint8_t q); 182 183 /** 174 184 \return query name 175 185 … … 228 238 229 239 \see bam_nt16_table 240 241 \since New in yat 0.11 230 242 */ 231 243 void sequence(size_t i, uint8_t x);
Note: See TracChangeset
for help on using the changeset viewer.