Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913843,1960TG33.9% 64.7 / 18.0 56F289V (TTT→GTT) rlmF23S rRNA m(6)A1618 methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (14/23);  new base G (18/1);  total (32/24)
Fisher's exact test for biased strand distribution p-value = 3.65e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTCGCGAAGCAGGTGATGTGGTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTG  >  NC_000913/843058‑843329
                                                                                                                                          |                                                                                                                                     
ttCGCGAAGCAGGTGATGTGGTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCgtt                                                                                                                                     >  2:1838814/1‑141 (MQ=255)
          aGGTGATGTGGTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGTTGTTGGCTGGa                                                                                                                           >  1:1535268/1‑141 (MQ=255)
          aGGTGATGTGGTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGTTATTggcgggg                                                                                                                           >  1:1535293/1‑135 (MQ=255)
          aGGTGATGTGGTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGTTATTTGCTggg                                                                                                                           >  1:1535281/1‑140 (MQ=255)
          aGGTGATGTGGTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGGTTTTGGCTggt                                                                                                                           >  1:1535232/1‑140 (MQ=255)
               aTGTGGTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGACTGGACCttt                                                                                                                      <  2:245416/141‑1 (MQ=255)
                  tgGTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATg                                                                                                                   <  2:47914/141‑1 (MQ=255)
                  tgGTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATg                                                                                                                   <  2:1849253/141‑1 (MQ=255)
                          aTCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGTTATTGCCTGGACCTTTATTAAcgacga                                                                                                           >  1:935924/1‑141 (MQ=255)
                           tCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTTCCTGGACCTTTATGAACGACGag                                                                                                          >  2:1865636/1‑141 (MQ=255)
                           tCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTTCCTGGACCATTATTAAAGACGag                                                                                                          >  2:1489885/1‑141 (MQ=255)
                           tCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTGTTTTGGCGGGTCCGTTATGAAAGACGag                                                                                                          >  2:220952/1‑141 (MQ=255)
                           tCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGTTATTGCCTGGACCTTTATGTAAGACGag                                                                                                          >  2:1489799/1‑141 (MQ=255)
                            cACTGGTAGCTCGTGGTGAAAACTTACAGCCGTTGTTTCGTGCCCTGACGGACGTGGGCGCGGTGAAGTTGGTTAAAAAAGAGATGGCCCAGGGGCAATAGCAGAGTCGGGTTATTGCCTGGACCTTTATGTACGACGagc                                                                                                         <  2:1024481/141‑1 (MQ=255)
                                ggTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGTTATTGGCTGGACCATTATGAAAGACGaccaacc                                                                                                     >  2:1700525/1‑138 (MQ=255)
                                  tATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGcgcc                                                                                                   <  1:248906/141‑1 (MQ=255)
                                  tATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGcgcc                                                                                                   <  1:249937/141‑1 (MQ=255)
                                  tATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGcgcc                                                                                                   <  1:249935/141‑1 (MQ=255)
                                  tATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGcgcc                                                                                                   <  1:249214/141‑1 (MQ=255)
                                  tATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTGTATTTGCGGGTCCATTATGTACGACGAGCAGcgcc                                                                                                   >  2:1419963/1‑141 (MQ=255)
                                  tATCTCGTGGTGAAAACTTACCGCCGTTGGATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGcgcc                                                                                                   <  1:602006/141‑1 (MQ=255)
                                       cGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGTTTTTTGCGGGACCTTTATTAAAGACGAGCAGCgccccccg                                                                                              >  2:1809630/1‑139 (MQ=255)
                                        gtggtgAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGTTTTTTGCTGGACCTTTATTAACGACGAGCAGCGCCCCCGtt                                                                                             >  1:770882/1‑141 (MQ=255)
                                        gtggtgAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGTTATTTCCGGGACCATTATGAAACACGAGCAGCGCCGCCGtt                                                                                             >  1:1438863/1‑141 (MQ=255)
                                        gtggtgAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGTTATTTCCATGTCCGTTATTTAAGACGATCAGCGCCGCCCtt                                                                                             >  1:2011651/1‑141 (MQ=255)
                                           gtgAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGTTATTGCCTGGGCCTTTTTTAACGACGAGCAGCGGCGCCGGTTTg                                                                                          >  2:1323893/1‑141 (MQ=255)
                                                    tACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGcc                                                                                 <  1:1419963/141‑1 (MQ=255)
                                                     aCCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCa                                                                                <  2:1690990/141‑1 (MQ=255)
                                                      ccgccgTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAg                                                                               <  2:1788937/141‑1 (MQ=255)
                                                                 tCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGt                                                                    <  2:1485382/141‑1 (MQ=255)
                                                                   gTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTTTTGCCTGGTCCTTTATGAACCACGAGCAGCGCCGCCGTTTTTTCAATCGCCCGCGTTTAATCGTTg                                                                  >  2:359486/1‑141 (MQ=255)
                                                                   gTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGGCTGGACCTTTTTGAACGACGAGCCGCGCCCCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTg                                                                  >  2:424003/1‑141 (MQ=255)
                                                                   gTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCGGGACCGTTATTAACGACGAGCCGCGCCGCCGGTTTGTCAATCGCCAGCGGTAAATCGTTg                                                                  >  1:949678/1‑141 (MQ=255)
                                                                   gTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGTTATTTGCTGGTCCGTTATTAAAGAAGAGCAGCGCCGCCGGTTTGTCAATAGCCCGCGTTaaatagttt                                                                  >  1:1895036/1‑135 (MQ=255)
                                                                   gTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGTTATTTCCTGGGCCTTTATTTAAGACCAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTg                                                                  >  1:949473/1‑141 (MQ=255)
                                                                            cggacgTGGGCGCGGTGAAGGTGGGTAAAACAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTGTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCcg                                                         <  2:292955/141‑1 (MQ=255)
                                                                                  tGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTa                                                   <  2:1676089/141‑1 (MQ=255)
                                                                                  tGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTa                                                   <  1:502492/141‑1 (MQ=255)
                                                                                  tGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCATTATTAACGACGAGCAGCGCCGCCGTTTTGTCAAAAGCCAGCGTTAAAGCGTTTGTTACCCCCGCGGTa                                                   >  2:485338/1‑141 (MQ=255)
                                                                                  tGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGTTTTTGGCTGGGCCTTTTTGAAACACGAGCAGCGCCCCCCTGTTTTGAAAAGCCACGGGTAAAGCGTTTGTTGCCCccgcggga                                                   >  1:180748/1‑139 (MQ=255)
                                                                                  tGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGTTATTGCCTGGACCATTATGAACGACGGGCCGCGCCCCCGTTTTGTAAAACCCCCGCGTTAAAACGTTGGTAACCCCGGCGGTa                                                   >  2:1153597/1‑141 (MQ=255)
                                                                                       gcgGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTGTATTTGCTGGGCCGTTATTAACCACGAGCCCCGCCGCCGGGTTGTGAATAGCCACCGGTAAATCGGTGGGTCCCCCCGCGGGAACggg                                              >  1:1684560/1‑140 (MQ=255)
                                                                                       gcgGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGGTATTTGCTGGTCCGTTATGAAAGGAGCGGCGCGCCCCCCTTTTTTGAATCGCCCCCGGTAAAACGGTGGGTACCCccgcggggaccgg                                              >  1:1684613/1‑134 (MQ=255)
                                                                                        cgGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGtt                                             <  1:1357520/141‑1 (MQ=255)
                                                                                        cgGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGtt                                             <  1:1357474/141‑1 (MQ=255)
                                                                                        cgGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGtt                                             <  1:1357445/141‑1 (MQ=255)
                                                                                                  ggtTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCGTTATTGCCGGGACCCTTATGAAACACGAGCAGCGCCCCCCGTTTTTCAAAACCCACCGGTAAAGCGGTTGTTCCCCCCCCGGTAACGGGTCAGCCGGGGt                                   >  1:1318898/1‑141 (MQ=255)
                                                                                                    ttAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTTCCTGGAACTTTAAGAAACACGAGCAGCGCCCCCGGTTTTTCAAACCCCAGCGTTAAAGCGTTGGTTACCCCCGCGGTAAAGGTTCAGCCGGGGTcg                                 >  2:586107/1‑141 (MQ=255)
                                                                                                      aaaaaaGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTcgcg                               <  1:1136373/141‑1 (MQ=255)
                                                                                                        aaaaGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTTCCTGGACCATTATTTACGACCAGCAGCGCCCCCGTGTTGTAAAACGCCAACGTTTAAGCGGTGGGTCCCCCCGCCGGAACGGTTCAGCCGGGGGCGcccc                             >  2:1998823/1‑141 (MQ=255)
                                                                                                                   cccAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAgcaggca                  <  1:701130/141‑1 (MQ=255)
                                                                                                                     cAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGTGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGgcagc                <  1:359486/141‑1 (MQ=255)
                                                                                                                            aaaaGCAGAGTCGCTTTATTGCCGGGACCATTATTAACGACGAGCAGCGCCGCCGGTTTGTCAATAGCCAGCGTTAAAGCGTTGGTTACCCCCGCGGTAACCGTTAAGCCGGGGTCGCGCCCGGCGCAGGCAGCACCGGAt         >  2:1770366/1‑141 (MQ=255)
                                                                                                                               aGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAg      <  2:700795/141‑1 (MQ=255)
                                                                                                                               aGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAg      <  2:770882/141‑1 (MQ=255)
                                                                                                                                  agagTCGCTTTATTGCCTGGGCCATTATGAACCACCATCAGCGCCCCCGTTTTTTCAATAGCCAGCGGTAAAGCGTTGGTTGCCCCGGCGGTGAACGTTCAGCCGGGGTGGCGCCCCGAGCCGGCAGCACCGGAttattat   >  2:1707418/1‑135 (MQ=255)
                                                                                                                                   gagTCGCTTTATTGCCGGGACCTTTATTTACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGGTGGTTCCCCCGGCGGTAACGGGTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTg  >  2:1396602/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                     
TTCGCGAAGCAGGTGATGTGGTTTACATCACTGGTATCTCGTGGTGAAAACTTACCGCCGTTGTATCGTGCCCTGACGGACGTGGGCGCGGTGAAGGTGGTTAAAAAAGAGATGGCCCAGGGGCAAAAGCAGAGTCGCTTTATTGCCTGGACCTTTATGAACGACGAGCAGCGCCGCCGTTTTGTCAATCGCCAGCGTTAAAGCGTTGGTTCCCCCGGCGGTAACGGTTCAGCCGGGGTCGCGCCCGGAGCAGGCAGCACCTGATAAGTCTG  >  NC_000913/843058‑843329

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑