Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,608,413 C→T R77C (CGT→TGT)  rimI → ribosomal‑protein‑S18‑alanine N‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,608,4130CT98.0% 183.2 / ‑4.3 51R77C (CGT→TGT) rimIribosomal‑protein‑S18‑alanine N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (32/18);  total (33/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

AAAACGTTTGCCAGCAACCAGGGCGAGCGTTATCTCAACTTTCAGTTAACGCAAAACGGCAAAATGGCGGCGTTTGCGATTACGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGCGTCAGGGATTGGGAAGGGCGCTGCTGGAACATCTGATCGACGAACTGGAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGC  >  NC_000913/4608272‑4608552
                                                                                                                                             |                                                                                                                                           
aaaaCGTTTGCCAGCAACCAGGGCGAGCGTTATCTCAACTTTCAGTTAACGCAAAACGGCAAAATGGCGGCGTTTGCGATTACGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGAt                                                                                                                                    >  2:635013/1‑151 (MQ=255)
 aaaCGTTTGCCAGCAACCAGGGCGAGCGTTATCTCAACTTTCAGTTAACGCAAAACGGCAAAATGGCGGCGTTTGCGATTACGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGAtt                                                                                                                                   <  1:30692/151‑1 (MQ=255)
 aaaCGTTTGCCAGCAACCAGGGCGAGCGTTATCTCAACTTTCAGTTAACGCAAAACGGCAAAATGGCGGCGTTTGCGATTACGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGAtt                                                                                                                                   <  1:30746/151‑1 (MQ=255)
 aaaCGTTTGCCAGCAACCAGGGCGAGCGTTATCTCAACTTTCAGTTAACGCAAAACGGCAAAATGGCGGCGTTTGCGATTACGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGAtt                                                                                                                                   >  1:362382/1‑151 (MQ=255)
     gTTTGCCAGCAACCAGGGCGAGCGTTATCTCAACTTTCAGTTAACGCAAAACGGCAAAATGGCGGCGTTTGCGATTACGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGATTggg                                                                                                                                >  2:837314/1‑150 (MQ=255)
         gCCAGCAACCAGGGCGAGCGTTATCTCAACTTTCAGTTAACGCAAAACGGCAAAATGGCGGCGTTTGCGATTACGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGATTGGGAAggg                                                                                                                           >  1:33576/1‑151 (MQ=255)
          ccAGCAACCAGGGCGAGCGTTATCTCAACTTTCAGTTAACGCAAAACGGCAAAATGGCGGCGTTTGCGATTACGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAgg                                                                                                                                       >  2:387226/1‑138 (MQ=255)
          ccAGCAACCAGGGCGAGCGTTATCTCAACTTTCAGTTAACGCAAAACGGCAAAATGGCGGCGTTTGCGATTACGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAgg                                                                                                                                       <  1:387252/138‑1 (MQ=255)
             gCAACCAGGGCGAGCGTTATCTCAACTTTCAGTTAACGCAAAACGGCAAAATGGCGGCGTTTGCGATTACGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGATTGGGAAGGGCgct                                                                                                                       >  2:778391/1‑151 (MQ=255)
             gCAACCAGGGCGAGCGTTATCTCAACTTTCAGTTAACGCAAAACGGCAAAATGGCGGCGTTTGCGATTACGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGATTGGGAAGGGCgct                                                                                                                       <  2:462760/151‑1 (MQ=255)
              cAACCAGGGCGAGCGTTATCTCAACTTTCAGTTAACGCAAAACGGCAAAATGGCGGCGTTTGCGATTACGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGATTGGGAAGGGCgctg                                                                                                                      <  1:627736/151‑1 (MQ=255)
              cAACCAGGGCGAGCGTTATCTCAACTTTCAGTTAACGCAAAACGGCAAAATGGCGGCGTTTGCGATTACGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGATTGGGAAGGGCgctg                                                                                                                      >  2:327961/1‑151 (MQ=255)
                 ccAGGGCGAGCGTTATCTCAACTTTCAGTTAACGCAAAACGGCAAAATGGCGGCGTTTGCGATTACGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGATTGGGAAGGGCGctgctg                                                                                                                   >  2:794920/1‑151 (MQ=255)
                     ggCGAGCGTTATCTCAACTTTCAGTTAACGCAAAACGGCAAAATGGCGGCGTTTGCGATTACGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGATTGGGAAGGGCGCTGCTGGAAc                                                                                                               >  1:286917/1‑151 (MQ=255)
                               aTCTCAACTTTCAGTTAACGCAAAACGGCAAAATGGCGGCGTTTGCGATTACGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGATTGGGAAGGGCGCTGCTGGAACATCTGATcg                                                                                                      <  2:56617/150‑1 (MQ=255)
                               aTCTCAACTTTCAGTTAACGCAAAACGGCAAAATGGCGGCGTTTGCGATTACGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGATTGGGAAGGGCGCTGCTGGAACATCTGATcg                                                                                                      <  1:20951/150‑1 (MQ=255)
                                  tcAACTTTCAGTTAACGCAAAACGGCAAAATGGCGGCGTTTGCGATTACGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGATTGGGAAGGGCGCTGCTGGAACATCTGATcgacg                                                                                                   >  1:398939/1‑150 (MQ=255)
                                  tcAACTTTCAGTTAACGCAAAACGGCAAAATGGCGGCGTTTGCGATTACGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGATTGGGAAGGGCGCTGCTGGAACATCTGATcgacg                                                                                                   >  2:523544/1‑150 (MQ=255)
                                     aCTTTCAGTTAACGCAAAACGGCAAAATGGCGGCGTTTGCGATTACGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGATTGGGAAGGGCGCTGCTGGAACATCTGATCGACGAACt                                                                                               <  1:428416/151‑1 (MQ=255)
                                      cTTTCAGTTAACGCAAAACGGCAAAATGGCGGCGTTTGCGATTACGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGATTGGGAAGGGCGCTGCTGGAACATCTGATCGACGAACTg                                                                                              >  2:342742/1‑151 (MQ=255)
                                       tttCAGTTAACGCAAAACGGCAAAATGGCGGCGTTTGCGATTACGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGATTGGGAAGGGCGCTGCTGGAACATCTGATCGACGAACTgg                                                                                             >  2:453956/1‑151 (MQ=255)
                                        ttCAGTTAACGCAAAACGGCAAAATGGCGGCGTTTGCGATTACGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGATTGGGAAGGGCGCTGCTGGAACATCTGATCGACGAACTgg                                                                                             <  2:342262/150‑1 (MQ=255)
                                        ttCAGTTAACGCAAAACGGCAAAATGGCGGCGTTTGCGATTACGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGATTGGGAAGGGCGCTGCTGGAACATCTGATCGACGAACTgg                                                                                             >  1:798917/1‑150 (MQ=255)
                                         tCAGTTAACGCAAAACGGCAAAATGGCGGCGTTTGCGATTACGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGCGTCAGGGATTGGGAAGGGCGCTGCTGGAACATCTGATCGACGAACTGGa                                                                                            >  2:463892/1‑150 (MQ=255)
                                           aGTTAACGCAAAACGGCAAAATGGCGGCGTTTGCGATTACGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGATTGGGAAGGGCGCTGCTGGAACATCTGATCGACGAACTGGaaaa                                                                                         <  1:667103/151‑1 (MQ=255)
                                           aGTTAACGCAAAACGGCAAAATGGCGGCGTTTGCGATTACGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGATTGGGAAGGGCGCTGCTGGAACATCTGATCGACGAACTGGaaa                                                                                          >  1:495842/1‑150 (MQ=255)
                                            gTTAACGCAAAACGGCAAAATGGCGGCGTTTGCGATTACGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGATTGGGAAGGGCGCTGCTGGAACATCTGATCGACGAACTGGaaaaa                                                                                        <  1:828335/151‑1 (MQ=255)
                                                  gCAAAACGGCAAAATGGCGGCGTTTGCGATTACGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGATTGGGAAGGGCGCTGCTGGAACATCTGATCGACGAACTGGAAAAACgcggc                                                                                  >  1:117095/1‑151 (MQ=255)
                                                   cAAAACGGCAAAATGGCGGCGTTTGCGATTACGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGATTGGGAAGGGCGCTGCTGGAACATCTGATCGACGAACTGGAAAAACgcggcg                                                                                 <  1:457156/151‑1 (MQ=255)
                                                       aCGGCAAAATGGCGGCGTTTGCGATTACGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGATTGGGAAGGGCGCTGCTGGAACATCTGATCGACGAACTGGAAAAACGCGGCGTGGc                                                                             >  2:511979/1‑151 (MQ=255)
                                                           cAAAATGGCGGCGTTTGCGATTACGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGATTGGGAAGGGCGCTGCTGGAACATCTGATCGACGAACTGGAAAAACGCGGCGTGGCGac                                                                          >  2:653714/1‑150 (MQ=255)
                                                                       gTTTGCGATTACGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGATTGGGAAGGGCGCTGCTGGAACATCTGATCGACGAACTGGAAAAACGCGGCGTGGCGACACTATGGCTGGaa                                                             <  2:544684/151‑1 (MQ=255)
                                                                             gATTACGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGATTGGGAAGGGCGCTGCTGGAACATCTGATCGACGAACTGGAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGt                                                       >  1:294081/1‑151 (MQ=255)
                                                                               ttACGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGATTGGGAAGGGCGCTGCTGGAACATCTGATCGACGAACTGGAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGc                                                     >  2:165206/1‑151 (MQ=255)
                                                                               ttACGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGATTGGGAAGGGCGCTGCTGGAACATCTGATCGACGAACTGGAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGc                                                     >  2:793432/1‑151 (MQ=255)
                                                                                tACGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGATTGGGAAGGGCGCTGCTGGAACATCTGATCGACGAACTGGAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCt                                                    >  2:303041/1‑151 (MQ=255)
                                                                                 aCGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGATTGGGAAGGGCGCTGCTGGAACATCTGATCGACGAACTGGAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCtt                                                   <  2:540631/151‑1 (MQ=255)
                                                                                   gCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGATTGGGAAGGGCGCTGCTGGAACATCTGATCGACGAACTGGAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTc                                                  >  1:161830/1‑150 (MQ=255)
                                                                                    cAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGATTGGGAAGGGCGCTGCTGGAACATCTGATCGACGAACTGGAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCa                                                 >  2:752610/1‑150 (MQ=255)
                                                                                     aaGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGATTGGGAAGGGCGCTGCTGGAACATCTGATCGACGAACTGGAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCaa                                                >  2:877995/1‑150 (MQ=255)
                                                                                                gATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGATTGGGAAGGGCGCTGCTGGAACATCTGATCGACGAACTGGAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCa                                    <  1:773121/151‑1 (MQ=255)
                                                                                                    aaGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGTGTCAGGGATTGGGAAGAGCGCTGCTGGAACATCTGGTCGACGAACTGGAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGc                                >  2:138684/1‑151 (MQ=255)
                                                                                                                  aaTATTGCGGTCGATCCTGACTATCAGTGTCAGGGATTGGGAAGGGCGCTGCTGGAACATCTGATCGACGAACTGGAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGtt                  >  1:487090/1‑151 (MQ=255)
                                                                                                                     atTGCGGTCGATCCTGACTATCAGTGTCAGGGATTGGGAAGGGCGCTGCTGGAACATCTGATCGACGAACTGGAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTa                <  1:778442/150‑1 (MQ=255)
                                                                                                                      ttGCGGTCGATCCTGACTATCAGTGTCAGGGATTGGGAAGGGCGCTGCTGGAACATCTGATCGACGAACTGGAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAgg              >  1:393529/1‑151 (MQ=255)
                                                                                                                         cGGTCGATCCTGACTATCAGTGTCAGGGATTGGGAAGGGCGCTGCTGGAACATCTGATCGACGAACTGGAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCtt           <  1:334788/151‑1 (MQ=255)
                                                                                                                            tCGATCCTGACTATCAGTGTCAGGGATTGGGAAGGGCGCTGCTGGAACATCTGATCGACGAACTGGAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTa         >  1:525808/1‑150 (MQ=255)
                                                                                                                             cGATCCTGACTATCAGTGTCAGGGATTGGGAAGGGCGCTGCTGGAACATCTGATCGACGAACTGGAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTaa        <  1:102031/150‑1 (MQ=255)
                                                                                                                                tCCTGACTATCAGTGTCAGGGATTGGGAAGGGCGCTGCTGGAACATCTGATCGACGAACTGGAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACg      >  1:464018/1‑149 (MQ=255)
                                                                                                                                 ccTGACTATCAGTGTCAGGGATTGGGAAGGGCGCTGCTGGAACATCTGATCGACGAACTGGAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAgg   >  1:646399/1‑151 (MQ=255)
                                                                                                                                  cTGACTATCAGTGTCAGGGATTGGGAAGGGCGCTGCTGGAACATCTGATCGACGAACTGGAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGc  >  1:55207/1‑151 (MQ=255)
                                                                                                                                             |                                                                                                                                           
AAAACGTTTGCCAGCAACCAGGGCGAGCGTTATCTCAACTTTCAGTTAACGCAAAACGGCAAAATGGCGGCGTTTGCGATTACGCAAGTGGTGCTGGATGAAGCTACATTGTTCAATATTGCGGTCGATCCTGACTATCAGCGTCAGGGATTGGGAAGGGCGCTGCTGGAACATCTGATCGACGAACTGGAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGC  >  NC_000913/4608272‑4608552

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