Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,508,541 T→A N587K (AAT→AAA rlhA → 23S rRNA 5‑hydroxycytidine C2501 synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,508,5410TA100.0% 67.6 / NA 20N587K (AAT→AAArlhA23S rRNA 5‑hydroxycytidine C2501 synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (12/8);  total (12/8)

GGCGAAGTCCCGGTGATGATCACCAAGCATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCG  >  NC_000913/1508464‑1508615
                                                                             |                                                                          
ggCGAAGTCCCGGTGATGATCACCAAGCATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAAAATc                                                                         <  2:283900/81‑1 (MQ=255)
      gTCCCGGTGATGATCACCAAGCATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAAAATCAAAAGc                                                                   >  2:283930/1‑81 (MQ=255)
       tCCCGGTGATGATCACCAAGCATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAAAATCAAAAGCt                                                                  <  2:283901/81‑1 (MQ=255)
        cccGGTGATGATCACCAAGCATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAAAATCAAAAGCTg                                                                 >  2:283929/1‑81 (MQ=255)
                     aCCAAGCATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAAAATCAAAAGCTGGAAGGCGACGcc                                                     <  2:283903/80‑1 (MQ=255)
                     aCCAAGCATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAAAATCAAAAGCTGGAAGGCGACGcc                                                     <  1:283902/80‑1 (MQ=255)
                            aTTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAAAATCAAAAGCTGGAAGGCGACGCCAATGCAAc                                             >  2:283937/1‑81 (MQ=255)
                                       tttGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAAAATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGc                                  >  2:283942/1‑81 (MQ=255)
                                           ccTTTAATCTGTGCCCGAAACAGGCGAAAGGCAAAATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATg                              <  2:283904/81‑1 (MQ=255)
                                              ttAATCTGTGCCCGAAACAGGCGAAAGGCAAAATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAg                           >  1:283926/1‑81 (MQ=255)
                                               tAATCTGTGCCCGAAACAGGCGAAAGGCAAAATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAg                           >  1:283939/1‑80 (MQ=255)
                                                           gAAACAGGCGAAAGGCAAAATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTa                >  1:283945/1‑79 (MQ=255)
                                                               cAGGCGAAAGGCAAAATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGttt          <  1:283905/81‑1 (MQ=255)
                                                               cAGGCGAAAGGCAAAATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGttt          <  1:283906/81‑1 (MQ=255)
                                                                   cGAAAGGCAAAATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGAtt      >  2:283947/1‑81 (MQ=255)
                                                                   cGAAAGGCAAAATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGAtt      <  2:283907/81‑1 (MQ=255)
                                                                    gAAAGGCAAAATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTg     >  2:283933/1‑81 (MQ=255)
                                                                    gAAAGGCAAAATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTg     >  2:283908/1‑81 (MQ=255)
                                                                      aaGGCAAAATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTgcc   >  2:283944/1‑81 (MQ=255)
                                                                       aGGCAAAATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTgccg  >  2:283938/1‑81 (MQ=255)
                                                                             |                                                                          
GGCGAAGTCCCGGTGATGATCACCAAGCATTGTCTGCGCTTTGCCTTTAATCTGTGCCCGAAACAGGCGAAAGGCAATATCAAAAGCTGGAAGGCGACGCCAATGCAACTGGTTAACGGCGATGAAGTATTAACGCTAAAGTTTGATTGCCG  >  NC_000913/1508464‑1508615

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: