Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,608,607 C→G 9.3% I141M (ATC→ATG)  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,6070CG 9.3% 181.1 / 9.6 54I141M (ATC→ATG) rimIribosomal‑protein‑S18‑alanine N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (24/25);  new base G (2/3);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.22e-01

GGAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTGATTACAGCGTCA  >  NC_000913/4608460‑4608742
                                                                                                                                                   |                                                                                                                                       
ggAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCaatcatca                                                                                                                                        <  2:221278/149‑1 (MQ=255)
    aaaCGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATGATGATgg                                                                                                                                     <  4:172246/148‑1 (MQ=255)
       cgcgGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATGATGATGGCGt                                                                                                                                  <  3:64121/148‑1 (MQ=255)
         cggcgTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGCTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATTGCGTTg                                                                                                                                <  4:254403/148‑1 (MQ=255)
                  gACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTa                                                                                                                      >  2:2829/1‑149 (MQ=255)
                  gACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTAAGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTa                                                                                                                      >  2:310159/1‑149 (MQ=255)
                   acacTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTAt                                                                                                                     <  1:207620/149‑1 (MQ=255)
                         tGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATa                                                                                                               <  2:297567/149‑1 (MQ=255)
                           gCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATAc                                                                                                              >  3:171165/1‑148 (MQ=255)
                            cTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTGGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACa                                                                                                             >  4:26521/1‑148 (MQ=255)
                            cTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACa                                                                                                             >  4:381713/1‑148 (MQ=255)
                              ggAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTa                                                                                                                      >  3:111257/1‑137 (MQ=255)
                              ggAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTa                                                                                                                      <  4:111257/137‑1 (MQ=255)
                               gAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCAAATCAGTATGTAATACAAGGt                                                                                                         <  2:162082/149‑1 (MQ=255)
                                  gTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAgg                                                                                                          <  3:268939/145‑1 (MQ=255)
                                  gTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAgg                                                                                                          >  4:268939/1‑145 (MQ=255)
                                  gTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGa                                                                                                      >  2:320912/1‑149 (MQ=255)
                                    ccGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTa                                                                                                                  <  3:279079/135‑1 (MQ=255)
                                    ccGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTa                                                                                                                  >  4:279079/1‑135 (MQ=255)
                                    ccGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGaa                                                                                                     >  3:191067/1‑148 (MQ=255)
                                          ttCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGaa                                                                                              >  2:145364/1‑149 (MQ=255)
                                             aaaCGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTg                                                                                           <  2:146452/149‑1 (MQ=255)
                                                     ccgccATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGa                                                                                   <  1:148758/149‑1 (MQ=255)
                                                             gCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCttt                                                                           <  1:371700/149‑1 (MQ=255)
                                                             gCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCtt                                                                            <  4:106361/148‑1 (MQ=255)
                                                                      gAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGTGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGAt                                                                  >  1:380/1‑149 (MQ=255)
                                                                      gAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGAt                                                                  >  2:131297/1‑149 (MQ=255)
                                                                       aaaGTTTCGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCTCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTACTACAAGGTGGAATACTGAAGTGGGACTGGATTTTCTTTGATGCCGAt                                                                  >  4:386369/1‑148 (MQ=255)
                                                                               ggCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATGATGATGGCGt                                                                                                                                  >  4:338448/1‑76 (MQ=255)
                                                                               ggCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATGATGATGGCGt                                                                                                                                  <  3:338448/76‑1 (MQ=255)
                                                                                gCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTg                                                         <  3:381713/148‑1 (MQ=255)
                                                                                 cTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGtt                                                       <  1:310159/149‑1 (MQ=255)
                                                                                 cTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGt                                                        <  4:192033/148‑1 (MQ=255)
                                                                                  tttAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGttt                                                      >  1:98692/1‑149 (MQ=255)
                                                                                     aaCGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTGGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTAc                                                    <  3:334560/148‑1 (MQ=255)
                                                                                      aCGAGACGACGACTCGCCGCAATTCCTACCCCACCACGGACGGTAGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTCGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTAcc                                                   <  4:245927/148‑1 (MQ=255)
                                                                                       cGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCtt                                                 <  1:145364/149‑1 (MQ=255)
                                                                                        gAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCttt                                                <  2:380/149‑1 (MQ=255)
                                                                                          ggCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTg                                               <  4:171165/148‑1 (MQ=255)
                                                                                                     gccgcAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCg                                   >  2:31260/1‑149 (MQ=255)
                                                                                                      ccgcAATTACTACCCCATCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGTTGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCg                                   >  4:271544/1‑148 (MQ=255)
                                                                                                         cAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGcc                                >  4:235156/1‑148 (MQ=255)
                                                                                                               ctacCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCg                         >  2:20765/1‑149 (MQ=255)
                                                                                                                 acCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCgg                        >  4:67062/1‑148 (MQ=255)
                                                                                                                  ccccACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTAcc                                                   <  2:334428/120‑1 (MQ=255)
                                                                                                                  ccccACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTAcc                                                   >  1:334428/1‑120 (MQ=255)
                                                                                                                            gacggTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTGa             <  3:400821/148‑1 (MQ=255)
                                                                                                                              cggTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTGATTa          <  2:98692/149‑1 (MQ=255)
                                                                                                                                gTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTGATTACa        <  1:2829/149‑1 (MQ=255)
                                                                                                                                 tCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCATTGATGCCGATGAAACGCTGTTTACCTTGGACTGATTCACCGGCCTGAAGCGGATGTTTCTTGATTACa        >  4:263747/1‑148 (MQ=255)
                                                                                                                                    cgAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTGATTACAGCg     >  4:384052/1‑148 (MQ=255)
                                                                                                                                     gAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCtt               >  1:275312/1‑137 (MQ=255)
                                                                                                                                     gAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCtt               <  2:275312/137‑1 (MQ=255)
                                                                                                                                      gcaaCGCCATGATGATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTGATTACAGCGTCa  >  1:321880/4‑149 (MQ=255)
                                                                                                                                       aGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTGATTACAGCGTCa  <  4:111584/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                       
GGAAAAACGCGGCGTGGCGACACTATGGCTGGAAGTCCGTGCTTCAAACGCTGCCGCCATTGCCCTGTACGAAAGTTTAGGCTTTAACGAGGCGACGATTCGCCGCAATTACTACCCCACCACGGACGGTCGCGAAGACGCCATCATCATGGCGTTGCCAATCAGTATGTAATACAAGGTGGAATAATGAAGTGGGACTGGATTTTCTTTGATGCCGATGAAACGCTGTTTACCTTTGACTCATTCACCGGCCTGCAGCGGATGTTTCTTGATTACAGCGTCA  >  NC_000913/4608460‑4608742

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