Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,409,168 A→C 12.1% intergenic (+73/‑107) rnr → / → rlmB exoribonuclease R, RNase R/23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,409,1680AC12.1% 104.9 / 7.6 33intergenic (+73/‑107)rnr/rlmBexoribonuclease R, RNase R/23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (17/12);  new base C (2/2);  total (19/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.40e-01

GCCATCGGCGAAAACGCAGAAAATAGCTGCAGCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCA  >  NC_000913/4409020‑4409309
                                                                                                                                                    |                                                                                                                                             
gCCATCGGCGAAAACGCAGAAAATAGCTGCAGCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGa                                                                                                                                               >  4:538752/1‑149 (MQ=255)
 ccATCGGCGAAAACGCAGAAAATAGCTGCAGCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCACGCCGCATCCGGCATGaa                                                                                                                                              <  2:396179/149‑1 (MQ=255)
    tCGGCGAAAACGCAGAAAATAGCTGCAGCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGc                                                                                                                                           >  1:567527/1‑149 (MQ=255)
    tCGGCGAAAACGCAGAAAATAGCTGCAGCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGc                                                                                                                                           >  2:135587/1‑149 (MQ=255)
                          cTGCAGCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCAGCGTGAACGCCTTATCCGGCCTACATAATCAc                                                                                                                     >  2:518660/1‑149 (MQ=255)
                           tGCAGCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACg                                                                                                                    <  3:538752/149‑1 (MQ=255)
                                cGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAAt                                                                                                               <  1:135587/149‑1 (MQ=255)
                                             gcgTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCa                                                                                                  >  4:135998/1‑149 (MQ=255)
                                                    gcgAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCa                                                                                           >  1:130398/1‑149 (MQ=255)
                                                    gcgAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCa                                                                                           <  2:130398/149‑1 (MQ=255)
                                                       aagaaAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGGTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGt                                                                                        >  4:135483/1‑149 (MQ=255)
                                                         gaaAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAg                                                                                      >  2:413596/1‑149 (MQ=255)
                                                          aaaaaaGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTg                                                                                                    <  4:498802/134‑1 (MQ=255)
                                                          aaaaaaGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTg                                                                                                    >  3:498802/1‑134 (MQ=255)
                                                            aaaaGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGcc                                                                                   >  1:124957/1‑149 (MQ=255)
                                                                        gTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAg                                                                       <  2:567527/149‑1 (MQ=255)
                                                                         tGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGt                                                                      <  2:124957/149‑1 (MQ=255)
                                                                                   ccTCTTTAAAAGAATAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAAc                                                                           >  1:520032/1‑134 (MQ=255)
                                                                                   ccTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAAc                                                                           <  2:520032/134‑1 (MQ=255)
                                                                                    ctctTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCaa                                                           >  4:269825/1‑149 (MQ=255)
                                                                                    ctctTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTCAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCaa                                                           >  4:472601/1‑149 (MQ=255)
                                                                                     tctTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCaaa                                                          <  1:413596/149‑1 (MQ=255)
                                                                                           aaaGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCACGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTaa                                                    >  4:53584/1‑149 (MQ=255)
                                                                                              gaagaGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAaaca                                                 >  2:1094/1‑149 (MQ=255)
                                                                                                      ttAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTAc                                         <  3:135998/149‑1 (MQ=255)
                                                                                                                      aTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATg                         <  3:269825/149‑1 (MQ=255)
                                                                                                                       tGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGa                        >  2:682938/1‑149 (MQ=255)
                                                                                                                       tGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTCATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAaacaac                                               >  4:652901/1‑126 (MQ=255)
                                                                                                                       tGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCGTATCCCGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAaacaac                                               <  3:652901/126‑1 (MQ=255)
                                                                                                                         cgcTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGAtt                      <  4:457586/149‑1 (MQ=255)
                                                                                                                                    ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCa                                                                <  3:641044/95‑1 (MQ=17)
                                                                                                                                    ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCa                                                                >  4:641044/2‑96 (MQ=17)
                                                                                                                                     gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAg                                                                       <  1:460167/88‑1 (MQ=17)
                                                                                                                                     gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAg                                                                       >  2:460167/1‑88 (MQ=17)
                                                                                                                                             cggcgTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCa  >  1:151147/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
GCCATCGGCGAAAACGCAGAAAATAGCTGCAGCGACCAAAGCGAAGCGTGCGGCGAAGAAAAAAGTGGCAGAGTGATCAATACCCTCTTTAAAAGAAGAGGGTTAGATTGCTGACAAAATGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCA  >  NC_000913/4409020‑4409309

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