Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 466,148 G→A P389S (CCA→TCA)  ybaE ← uncharacterized protein YbaE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913466,1480GA97.1% 109.8 / ‑4.4 34P389S (CCA→TCA) ybaEuncharacterized protein YbaE
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (18/15);  minor base T (1/0);  total (19/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

ATTAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGATGCCGCTACTGACCTGGCTGACCCGTTGCAGCTCCTCCGGTTTGCCGATGGTGATTTGCACGGGATGCCGACAACTGGTTCCCAAATCTT  >  NC_000913/465922‑466390
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  
aTTAATCTGTCGCCCATCATGAGGTCTGCGTTTGCCTTTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTTTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGAt                                                                                                                                                                                                                            <  1:130541/251‑1 (MQ=255)
  tAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATc                                                                                                                                                                                                                           >  2:319795/1‑250 (MQ=255)
      cTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAACAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGtt                                                                                                                                                                                                                      >  1:298225/1‑251 (MQ=255)
        gTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGttt                                                                                                                                                                                                                     <  2:188183/250‑1 (MQ=255)
          cGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTtct                                                                                                                                                                                                                  >  1:352593/1‑251 (MQ=255)
              catcatGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGGATGAAGTTCTATCGGTAGGGGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGGGATCAGGGTTTCTCCgc                                                                                                                                                                                                              >  2:6241/1‑250 (MQ=255)
                     aGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACCGCCTTCAGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGGGATCCGGTTTTCTCCCACTTCTa                                                                                                                                                                                                        >  2:167420/1‑250 (MQ=255)
                                     gTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGAAATGCATGACTGGCGGTGATCAGGGTTTCTCCGACTTCTAAAGTTTGTAATAAAccc                                                                                                                                                                                       >  1:54802/1‑250 (MQ=255)
                                                           ttttGCGTTATGAAAAATAATTGTGAGTTCACTGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATc                                                                                                                                                                  <  2:352593/250‑1 (MQ=255)
                                                                           cctcttccgacctCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCa                                                                                                                                                 <  2:298834/239‑1 (MQ=255)
                                                                                      ttCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGGTCAGGGTTTCTCCGACTTCTAAAGTTTGTAATAAACCGGATTGATGAATAATTGAGATCAATTTTCGCGCCTGACAGAGGGAGAGTc                                                                                                                                      >  1:175334/1‑251 (MQ=255)
                                                                                        cacaGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCa                                                                                                                                                 >  1:298834/1‑238 (MQ=255)
                                                                                        cacaGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGGGAGTCg                                                                                                                                     >  2:186412/1‑250 (MQ=255)
                                                                                        cacaGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCg                                                                                                                                     <  1:186412/250‑1 (MQ=255)
                                                                                            gCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGAc                                                                                                                                >  1:121517/1‑251 (MQ=255)
                                                                                                   gCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATACACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGAATTCTAACGTTTGTAATCAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGGCGGGGACtttttcg                                                                                                                         >  1:64615/1‑248 (MQ=255)
                                                                                                                   gTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCGGCGTCAAGGTTTTCGGTAGGTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTGGCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCaa                                                                                                         <  1:140980/251‑1 (MQ=255)
                                                                                                                     aGGCGTTCTGCCATGGTATGAAGTAGTATCGGTAGGGGATAACCCAGCGTCAAGGTTTTCGGTAGTTATACTTCAGCCGGTCCCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACGGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGGAACGTCAAATAGCaaa                                                                                                        <  2:272736/250‑1 (MQ=255)
                                                                                                                           tCTGCCATGGTATGAAGTTCTATCGGTAGGGGATAAACCAGCGGCAAGGTTTTCGGTAGTTTGACTTAATCCGGTACCTGACACTGCGGAATAGTCCAGCCTGTCAGTACTGCATGACTGGCGGCGATCAGGGTTTCTCCGCCTTCTTATGTTTGTGCTTAACCGGGTTGATGGGTAATGGCGATCACTTTTCGCGCCTGCCAGAGGGGGAGTCGGGGGCTTTTgcgccccgtccacttggaaaaacctc                                                                                                  >  1:398889/1‑228 (MQ=255)
                                                                                                                              tcctgggTATGAAGTTCTATCGGTATTTGATAAACCATCGCTAAGTTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGGTTTCTCCGACTTCTAACGTTTGTAATAAACCTGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCTAGATGGAGAGTCGGGGACTTTTGCGCAACTTCAAATAGCAAAAACCTAAAc                                                                                               <  1:133819/245‑1 (MQ=255)
                                                                                                                               ccATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGaataa                                                                                                                                                                          <  2:195333/174‑1 (MQ=255)
                                                                                                                               ccATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGaataa                                                                                                                                                                          >  1:195333/1‑174 (MQ=255)
                                                                                                                               ccATGGTATGAAGTTCTATCGGTAGGGGATAAACCAGCGGCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGTAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGAGCACTTTGCGCGGCTGCCAGAGGGAGATTCGGGGACTTTGGCGCAACGTCAAATAGCAAATACCTAAACTg                                                                                             <  2:64615/251‑1 (MQ=255)
                                                                                                                                     tATGAAGTTCTATCGGTAGGGGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGGTTGATGAAAAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGGCTTTTGCGGAACGTCAAATAGCAAAAAACAAAACTGATGcgc                                                                                       >  2:341146/1‑249 (MQ=255)
                                                                                                                                          aGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGATGCCGCTACt                                                                                  <  1:319795/251‑1 (MQ=255)
                                                                                                                                          aGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCCGCCCGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCCAAAAACTAAACTGATGgcgccacc                                                                                  >  1:161618/1‑247 (MQ=255)
                                                                                                                                            ttCTATCGGTAGGCGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCCGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTTCCAGAGGGGGAGTCGGGGACTTTTTCGCCACGTTCAATACCCAAAACCTAAACTGGTGcccccactg                                                                                 >  2:134518/1‑243 (MQ=255)
                                                                                                                                                                          ggTTTTCGGTAGTTTGCCTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGATGCCGCTACTGACCTGGCTGACCCGTTGCAGCTCCTCCGGtt                                                  <  1:341146/251‑1 (MQ=255)
                                                                                                                                                                                             tcacccGGTCCCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGATGCCGCTACTGACCTGGCTGACCCGTTGCAGCTCCTCCGGTTTGCCGATGGTGATTTGCa                                <  2:212096/246‑1 (MQ=255)
                                                                                                                                                                                                       ccTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCCGAGGGAGAGTCGGGGGCTTTTGCGCAACGTCAAATAGCACAAACCTAAACTGATGCCGCTACTGCACTGGCTGACCCCGTGCAGCTCCTCCCGGTTGCCGCTGGGGGTTTTCACGGGAAGcccg                     >  2:139285/1‑249 (MQ=255)
                                                                                                                                                                                                              aTGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGGTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGGCTTTTTCGCAACGTCAAATAGCAAAAACCTAAACTGATTCCCCTACTGACCTGGCTTACCCGTTGCAGCTCCTCCGGTTTGCCGATGGGGATTTGCACGGGGTGCCGACAAccgg              >  1:138502/1‑248 (MQ=255)
                                                                                                                                                                                                                    aaTAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGATGCCGCTACTGACCTGGCTGACCCGTTGCAGCTCCTCCGGGTTGCCGATGGTGATTTGCACGGGATGCCGACAACTGGTTCCCa        <  1:95723/251‑1 (MQ=255)
                                                                                                                                                                                                                       aGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGGTTTTTCCGACTTCTAACGTTTGTAATAAAACGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAAACTAAACTGATGCCGCTACTGAACTGGCTGACCCGTTGCAGCTCCTCCGGTTTGCCGAGGGAGATTTGCACGGGATGCCGACAACTGGTTCCCAAAt     >  1:103432/1‑251 (MQ=255)
                                                                                                                                                                                                                          ccAGCCTGACAGTAATGCGTGTCTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCTGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCTGGGACTTTTGCGCAACGGCAAATAGCAAAACCCTAAACTTATGCCGCTACTGACCTGGCTGACCCGTTGCAGCTCCTCCGGTTTGCCGATGGTGATTTGCTCGGTATGCCGACCACTGGTTCCCAAAGCtt  <  2:1742/251‑1 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                  
ATTAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGATGCCGCTACTGACCTGGCTGACCCGTTGCAGCTCCTCCGGTTTGCCGATGGTGATTTGCACGGGATGCCGACAACTGGTTCCCAAATCTT  >  NC_000913/465922‑466390

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