Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,233,560 C→T 6.9% intergenic (‑327/‑198) lysC ← / → pgi lysine‑sensitive aspartokinase 3/glucosephosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,233,5600CT 6.9% 179.3 / 3.3 58intergenic (‑327/‑198)lysC/pgilysine‑sensitive aspartokinase 3/glucosephosphate isomerase
Reads supporting (aligned to +/- strand):  ref base C (27/27);  new base T (2/2);  total (29/29)
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.13e-01

TGAACGTGGCCAGCCGTTCAATCACCTCGGCGATGCACCCCCTCAGGTGTTATCACAGGACTGGCTCCTCCAACACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGG  >  NC_000913/4233414‑4233707
                                                                                                                                                  |                                                                                                                                                   
tgaACGTGGCCAGCCGTTCAATCACCTCGGCGATGCACCCCCTCAGGTGTTATCACAGGACTGGCTCCTCCAACACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCgg                                                                                                                                                   >  1:356024/1‑149 (MQ=255)
 gaACGTGGCCAGCCGTTCAATCACCTCGGCGATGCACCCCCTCAGGTGTTATCACAGGACTGGCTCCTCCAACACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCggg                                                                                                                                                  <  3:192641/149‑1 (MQ=255)
  aaCGTGGCCAGCCGTTCAATCACCTCGGCGATGCACCCCCTCAGGTGTTATCACAGGACTGGCTCCTCCAACACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCgggg                                                                                                                                                 >  2:530113/1‑149 (MQ=255)
   aCGTGGCCAGCCGTTCAATCACCTCGGCGATGCACCCCCTCAGGTGTTATCACAGGACTGGCTCCTCCAACACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGc                                                                                                                                                <  4:324860/149‑1 (MQ=255)
               gTTCAATCACCTCGGCGATGCACCCCCTCAGGTGTTATCACAGGACTGGCTCCTCCAACACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGa                                                                                                                                    <  2:88570/149‑1 (MQ=255)
                  cAATCACCTCGGCGATGCACCCCCTCAGGTGTTATCACAGGACTGGCTCCTCCAACCCCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATgg                                                                                                                                 <  2:356024/149‑1 (MQ=255)
                     tCACCTCGGCGATGCACCCCCTCAGGTGTTATCACAGGACTGGCTCCTCCAACACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGt                                                                                                                              >  4:647300/1‑149 (MQ=255)
                     tCACCTCGGCGATGCACCCCCTCAGGTGTTATCACAGGACTGGCTCCTCCAACACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGt                                                                                                                              >  4:150656/1‑149 (MQ=255)
                      cACCTCGGCGATGCACCCCCTCAGGTGTTATCACAGGACTGGCTCCTCCAACACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTc                                                                                                                             <  3:247592/149‑1 (MQ=255)
                          tCGGCGATGCACCCCCTCAGGTGTTATCACAGGACTGGCTCCTCCAACACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGt                                                                                                                                          <  2:7024/132‑1 (MQ=255)
                          tCGGCGATGCACCCCCTCAGGTGTTATCACAGGACTGGCTCCTCCAACACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGt                                                                                                                                          >  1:7024/1‑132 (MQ=255)
                          tCGGCGATGCACCCCCTCAGGTGTTATCACAGGACTGGCTCCTCCAACACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGc                                                                                                                         >  1:171965/1‑149 (MQ=255)
                               gATGCACCCCCTCAGGTGTTATCACAGGACTGGCTCCTCCAACACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGAtt                                                                                                                    >  3:123074/1‑149 (MQ=255)
                                 tGCACCCCCTCAGGTGTTATCACAGGACTGGCTCCTCCAACACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGAtttt                                                                                                                  <  2:171965/149‑1 (MQ=255)
                                   cACCCCCTCAGGTGTTATCACAGGACTGGCTCCTCCAACACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTc                                                                                                                >  4:321440/1‑149 (MQ=255)
                                    aCCCCCTCAGGTGTTATCACAGGACTGGCTCCTCCAACACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCa                                                                                                               <  1:139659/149‑1 (MQ=255)
                                           cAGGTGTTATCACAGGACTGGCTCCTCCAACACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTc                                                                                                        >  1:529309/1‑149 (MQ=255)
                                             ggTGTTATCACAGGACTGGCTCCTCCAACACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCt                                                                                                      >  2:2589/1‑149 (MQ=255)
                                              gtgtTATCACAGGACTGGCTCCTCCAACACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGTGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTg                                                                                                     >  2:228132/1‑149 (MQ=255)
                                              gtgtTATCACAGGACTGGCTCCTCCAACACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGTGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTg                                                                                                     >  4:133436/1‑149 (MQ=255)
                                                  tATCACAGGACTGGCTCCTCCAACACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCgg                                                                                                 <  4:478982/149‑1 (MQ=255)
                                                   aTCACAGGACTGGCTCCTCCAACACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGt                                                                                                <  3:150656/149‑1 (MQ=255)
                                                      acaGGACTGGCTCCTCCAACACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCag                                                                                             >  3:351002/1‑149 (MQ=255)
                                                      acaGGACTGGCTCCTCCAACACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCag                                                                                             >  2:529966/1‑149 (MQ=255)
                                                             tGGCTCCTCTAACACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTg                                                                                                                                           <  1:628791/96‑1 (MQ=255)
                                                             tGGCTCCTCCAACACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTg                                                                                                                                           >  2:628791/1‑96 (MQ=255)
                                                                ctcctcCAACACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAAc                                                                                                                                      <  1:79302/98‑1 (MQ=255)
                                                                ctcctcCAACACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAAc                                                                                                                                      >  2:79302/1‑98 (MQ=255)
                                                                    tccAACACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCtaata                                                                               >  1:439954/1‑149 (MQ=255)
                                                                       aaCACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAAt                                                                            <  2:224579/149‑1 (MQ=255)
                                                                       aaCACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAAt                                                                            <  2:439954/149‑1 (MQ=255)
                                                                       aaCACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAAt                                                                            <  2:397137/149‑1 (MQ=255)
                                                                                     ggCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCtaat                                                                                >  2:94802/1‑131 (MQ=255)
                                                                                     ggCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCtaat                                                                                <  1:94802/131‑1 (MQ=255)
                                                                                       cAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAg                                                            <  3:321440/149‑1 (MQ=255)
                                                                                       cAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAg                                                            <  1:530113/149‑1 (MQ=255)
                                                                                        aaCGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAgg                                                           >  4:120211/1‑149 (MQ=255)
                                                                                                tcttctGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTacac                                                   <  3:120211/149‑1 (MQ=255)
                                                                                                tcttctGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTacac                                                   <  3:647300/149‑1 (MQ=255)
                                                                                                  ttctGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGTGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTc                                                 <  3:133436/149‑1 (MQ=255)
                                                                                                  ttctGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGTGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTc                                                 <  1:228132/149‑1 (MQ=255)
                                                                                                      taccTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATg                                                                           <  3:287645/117‑1 (MQ=255)
                                                                                                      taccTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATg                                                                           >  4:287645/3‑119 (MQ=255)
                                                                                                             gcTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCt                                      <  4:351002/149‑1 (MQ=255)
                                                                                                               tAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTaa                                    >  4:239613/1‑149 (MQ=255)
                                                                                                                   gcAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGc                                >  1:5016/1‑149 (MQ=255)
                                                                                                                   gcAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGc                                >  3:318036/1‑149 (MQ=255)
                                                                                                                          gtaCATTTATAAATAAAGGGTGAGCGGGTCGTTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTGAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTaaaa                         >  4:424671/1‑149 (MQ=255)
                                                                                                                             cATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCa                      <  2:5016/149‑1 (MQ=255)
                                                                                                                              aTTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCAGTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCAt                     >  1:469230/1‑149 (MQ=255)
                                                                                                                               tttATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCa                               >  3:644755/1‑138 (MQ=255)
                                                                                                                               tttATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCa                               <  4:644755/138‑1 (MQ=255)
                                                                                                                               tttATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATc                    >  3:263836/1‑149 (MQ=255)
                                                                                                                                     aataaaGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACAttt              >  1:609464/1‑149 (MQ=255)
                                                                                                                                          aGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCtg         >  1:655281/1‑149 (MQ=255)
                                                                                                                                           gggTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCtgt        <  4:123074/149‑1 (MQ=255)
                                                                                                                                                 gCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTgg  <  4:318036/149‑1 (MQ=255)
                                                                                                                                                 gCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTgg  <  1:2589/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TGAACGTGGCCAGCCGTTCAATCACCTCGGCGATGCACCCCCTCAGGTGTTATCACAGGACTGGCTCCTCCAACACCGTTACTTGGGCAACGCGCCTCTTCTGGCCTGCGCTAGCGCAGGTAGTACATTTATAAATAAAGGGTGAGCGGGGCGGTTGTCAACGATGGGGTCATGCGGATTTTTCATCCACTCCTGGCGGTCAGTAGTTCAGCTAATAAATGCTTCACTGCGCTAAGGGTTTACACTCAACATTACGCTAACGGCACTAAAACCATCACATTTTTCTGTGACTGG  >  NC_000913/4233414‑4233707

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