Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,029,739 (A)8→7 coding (231/1191 nt) rlmI ← 23S rRNA m(5)C1962 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,029,7320A.94.4% 200.6 / 5.4 54coding (238/1191 nt)rlmI23S rRNA m(5)C1962 methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base . (26/25);  total (27/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCT  >  NC_000913/1029608‑1029867
                                                                                                                            |                                                                                                                                       
ccGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCg                                                                                                                         >  2:277043/1‑140 (MQ=255)
                  aTACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCa                                                                                                       >  2:857563/1‑140 (MQ=255)
                   tACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCaa                                                                                                      <  2:903923/140‑1 (MQ=255)
                   tACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCaa                                                                                                      >  2:756903/1‑140 (MQ=255)
                     cccGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAAc                                                                                                    >  2:210311/1‑140 (MQ=255)
                     cccGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAAc                                                                                                    >  1:532322/1‑140 (MQ=255)
                       cGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGGAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGt                                                                                                  >  2:739588/1‑140 (MQ=255)
                       cGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGt                                                                                                  >  2:356723/1‑140 (MQ=255)
                        ggCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTc                                                                                                        <  2:429482/133‑1 (MQ=255)
                        ggCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTc                                                                                                        >  1:429482/1‑133 (MQ=255)
                        ggCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTc                                                                                                 >  1:82653/1‑140 (MQ=255)
                        ggCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTc                                                                                                 >  1:17750/1‑140 (MQ=255)
                        ggCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGGCAAACGTc                                                                                                 >  1:463687/1‑140 (MQ=255)
                         gCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTcc                                                                                                >  1:800573/1‑140 (MQ=255)
                              ccATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGAcg                                                                                           <  2:463687/140‑1 (MQ=255)
                                    gATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGcc                                                                                     <  2:800573/140‑1 (MQ=255)
                                       tctcCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCgg                                                                                  <  1:739588/140‑1 (MQ=255)
                                                aTTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGaa                                                                         >  2:50591/1‑140 (MQ=255)
                                                     aCGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTgg                                                                    <  1:756903/140‑1 (MQ=255)
                                                          aGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAAt                                                               >  2:109610/1‑140 (MQ=255)
                                                              gTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAgc                                                           >  1:588113/1‑140 (MQ=255)
                                                               tCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGaa                                                                         >  2:520579/1‑125 (MQ=255)
                                                               tCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAgcg                                                          >  1:440174/1‑140 (MQ=255)
                                                               tCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAgcg                                                          >  2:728854/1‑140 (MQ=255)
                                                               tCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGAAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGaa                                                                         <  1:520579/125‑1 (MQ=255)
                                                                 gAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAACAATGTCGATAGACTCAGACGGGGCAAACGTCCAGACGCGCGCCCGGGTTTTCGGAGCTGGGGAATAAAcgcc                                                        >  2:948618/1‑140 (MQ=255)
                                                                 gAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGcgcc                                                        <  1:815898/140‑1 (MQ=255)
                                                                  aGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGcgccg                                                       <  2:17750/140‑1 (MQ=255)
                                                                  aGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGcgccg                                                       <  2:371014/140‑1 (MQ=255)
                                                                      cATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTg                                                   >  2:518327/1‑140 (MQ=255)
                                                                        tCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCt                                                 <  2:440174/140‑1 (MQ=255)
                                                                         cTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTa                                                <  1:356723/140‑1 (MQ=255)
                                                                                    agccagTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTccc                                     >  2:354865/1‑140 (MQ=255)
                                                                                    agccagTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTccc                                     <  1:50591/140‑1 (MQ=255)
                                                                                      ccagTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCtg                                   <  2:932250/140‑1 (MQ=255)
                                                                                      ccagTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCtg                                   <  2:532322/140‑1 (MQ=255)
                                                                                         gTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCtgatg                                <  1:109610/140‑1 (MQ=255)
                                                                                         gTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCtgatg                                <  1:277043/140‑1 (MQ=255)
                                                                                         gTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCtgatg                                <  1:70676/140‑1 (MQ=255)
                                                                                              gCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATcaa                           <  1:857563/140‑1 (MQ=255)
                                                                                               ccATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATcaac                          >  2:17141/1‑140 (MQ=255)
                                                                                                cATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCtgatg                                >  2:308472/1‑133 (MQ=255)
                                                                                                cATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCtgatg                                <  1:308472/133‑1 (MQ=255)
                                                                                                   ttttGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTccc                                     >  2:447584/1‑125 (MQ=255)
                                                                                                   ttttGTGCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTccc                                     <  1:447584/125‑1 (MQ=255)
                                                                                                      tgtgCTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCg                   <  1:539185/140‑1 (MQ=255)
                                                                                                          cTTGTTGCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATgg               >  1:38559/1‑140 (MQ=255)
                                                                                                               tgCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCa          >  1:262591/1‑140 (MQ=255)
                                                                                                                gCAAACGGCGGGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCAc         <  1:122128/140‑1 (MQ=255)
                                                                                                                     cggcggGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAgg    >  2:56264/1‑140 (MQ=255)
                                                                                                                      ggcggGAAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAgg    <  2:395529/140‑1 (MQ=255)
                                                                                                                      ggcggGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAggc   <  1:420440/140‑1 (MQ=255)
                                                                                                                      ggcggGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAggc   <  1:210311/140‑1 (MQ=255)
                                                                                                                       gcggGAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAggct  <  1:56264/140‑1 (MQ=255)
                                                                                                                            |                                                                                                                                       
CCGAAACGATCGATAGTAATACCCGGCAGGCCATCAGATTCTCCGGCGATTAAACGATAGCTGTCGAGGCCATCTTTTTGCGCCAGCCAGTCACGCCATTTTTGTGCTTGTTGCAAACGGCGGGAAAAAAAAGCAATGTCGATAGACTCAGACGGGTCAAACGTCCAGACGCGCGCCCGGATTTGCGAAGCTGGCGAATAAGCGCCGCGTGCTAACCATTTTCCCTGATGATCAACAATATCGATGGTTTCACCGAGGCT  >  NC_000913/1029608‑1029867

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: