Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,153,446 A→G 10.5% intergenic (‑349/+249) ppc ← / ← argE phosphoenolpyruvate carboxylase/acetylornithine deacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,153,4460AG10.5% 106.1 / 3.2 53intergenic (‑349/+249)ppc/argEphosphoenolpyruvate carboxylase/acetylornithine deacetylase
Reads supporting (aligned to +/- strand):  ref base A (19/26);  new base G (5/2);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 2.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTA  >  NC_000913/4153303‑4153594
                                                                                                                                               |                                                                                                                                                    
cTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTgg                                                                                                                                                 <  1:143220/149‑1 (MQ=255)
 ttATGGAAATGTTAAAAAATCGCCCCAAGTACCCCCAAAGGTGTAGGTCGGATAAGATGCGCAAGTCTCGCATCCGACATTACTGCTTCATTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTggg                                                                                                                                                <  1:163733/149‑1 (MQ=18)
 ttATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTggg                                                                                                                                                >  2:579914/1‑149 (MQ=255)
 ttATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTggg                                                                                                                                                >  2:684740/1‑149 (MQ=255)
 ttATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTggg                                                                                                                                                >  2:494740/1‑149 (MQ=255)
    gggAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCAt                                                                                                                                             <  2:652414/148‑1 (MQ=255)
      ggaaTGTTAAAAGCTCGCCCCAAGTAACACCAGAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCg                                                                                                                                           <  2:420950/147‑1 (MQ=255)
      gAACTGTTAAAAAACCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCg                                                                                                                                           <  2:265405/149‑1 (MQ=255)
           gTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGa                                                                                                                                      <  2:128414/149‑1 (MQ=255)
            ttAAAAAATCGCCCCAAGTAACACCTAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAg                                                                                                                                     >  1:347794/1‑149 (MQ=255)
            ttAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAg                                                                                                                                     >  1:623099/1‑149 (MQ=255)
             tAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCGTCGTTTGAGt                                                                                                                                    >  1:15323/1‑149 (MQ=255)
             tAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCAGCCTACAGTTGGGCATCGTTTGAGt                                                                                                                                    >  1:425909/1‑149 (MQ=255)
             tAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGGATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGt                                                                                                                                    >  1:655646/1‑149 (MQ=255)
                          cAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGTCATCGTTTGAGTCACCGtcggccgg                                                                                                                       >  1:531578/1‑149 (MQ=255)
                           aaGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGTATCCGACATTATTGCGGCACGGGAGAGTGGCGACAGTGCCGGATGCGGCGCGGGCGCCTTATCCGGCCTACGGTGTGGCATCGTTTGAGTAACGGccggtcggg                                                                                                                      >  2:659840/1‑148 (MQ=14)
                              tAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATaa                                                                                                                   >  1:258533/1‑149 (MQ=255)
                                acacCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGa                                                                                                                 <  1:619039/149‑1 (MQ=255)
                                 cacCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGAt                                                                                                                <  1:187099/149‑1 (MQ=12)
                                 cacCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTAGGTCGGATAAGAt                                                                                                                <  2:530569/149‑1 (MQ=12)
                                    cAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATgcg                                                                                                             <  1:643065/149‑1 (MQ=255)
                                       gggTGGAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCaa                                                                                                          <  1:211313/148‑1 (MQ=255)
                                         gtgtAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGt                                                                                                        >  1:171826/1‑149 (MQ=39)
                                            gaGGTCGGATAAGGTGGGCAAGTATCGCAGCCGACACTATTGCGGCACTGGAGTTTGGCAGCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATc                                                                                                     <  2:679834/148‑1 (MQ=11)
                                             aGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCg                                                                                                    <  2:596783/149‑1 (MQ=40)
                                             aGGCCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCg                                                                                                    <  1:246733/149‑1 (MQ=255)
                                                      aaGATGTGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACa                                                                                           <  1:397577/149‑1 (MQ=38)
                                                        gATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGACGTCGGTTGAGTCACTGTCGGTCGAATAAGATGCGCAAGCCTCGCGTCCGAccta                                                                                         >  2:632427/1‑146 (MQ=11)
                                                           gcgcAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCCCTGTCGGTCGGAGAAGATGCGCAAGTATCGCATCCGGCattatt                                                                                      >  1:189825/1‑149 (MQ=255)
                                                             gcAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCCGATAAGATCCGCAAGCATCGCGTCCGACATTATTGc                                                                                    >  2:145784/1‑149 (MQ=255)
                                                             gcAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGACGGATAAGATGCGCCAGTATCGCATCCGACATTATTGc                                                                                    >  2:113487/1‑149 (MQ=255)
                                                             gcAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGACTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGc                                                                                    <  1:494740/149‑1 (MQ=255)
                                                               aaGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCACGTATCGCATCCGACattatagagc                                                                                  >  2:323093/1‑144 (MQ=12)
                                                                      gCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGa                                                                           <  1:573326/149‑1 (MQ=38)
                                                                            gACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTgg                                                                     <  1:391714/149‑1 (MQ=38)
                                                                            gACATTATTGCGGCACAGGAGTTTGGCAACAGTGCCGGATGGGGCGCGGGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCGCTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTgg                                                                     <  2:225155/149‑1 (MQ=255)
                                                                               attattGCCGCACTGGAGTTTCGCAACGGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCaa                                                                  <  2:35373/149‑1 (MQ=255)
                                                                                 tattGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACa                                                                <  2:347794/149‑1 (MQ=38)
                                                                                 tattGCGGCACAGGAGTTTGGCGACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGCGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACa                                                                <  2:567851/149‑1 (MQ=255)
                                                                                         cACTGGGGTTTGGCAACAGTGCCGGATGCGGCGCGGGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGa                                                        <  2:65471/149‑1 (MQ=255)
                                                                                              gAGTTTGGCAACGGTGCAGGATGTGGCGCGGGCTCCTTATCCGGCATACAGGTGGGCATCGTTTGCGTCGCTGTCGGTCGGATAAGATGCGCTAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATgcgg                                                   <  2:434730/149‑1 (MQ=14)
                                                                                                   tGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCAATGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGGGGTTGGCAACAGTGCCGGATGCGGCGCGa                                              >  1:311402/1‑149 (MQ=255)
                                                                                                    cgCAACAGTGCGGGATGGGGTGGGATGACCATATCTGACCAACACTCGCGCATCGATTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAg                                             <  2:218496/148‑1 (MQ=14)
                                                                                                      cAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTg                                                                      >  1:554973/1‑122 (MQ=39)
                                                                                                          aGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCtt                                       >  1:474602/1‑149 (MQ=255)
                                                                                                              ccGGATACGGGGGGGGCGCCGTATCCGGCCTACCGGTGGCTATTGTTTGGGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATcc                                   <  2:14607/149‑1 (MQ=11)
                                                                                                                          ccaGCGCCTTATCCGGCCAACAGTTGGCCATCGTTTGAGTCATTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGAGGCGGCGCGAGCGCCTTATCCGGCTTACGGTTg                       <  1:289201/147‑1 (MQ=14)
                                                                                                                            aGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTggg                     >  2:341635/1‑149 (MQ=255)
                                                                                                                              cgcCTTATCCGGACTACAGTTGGTCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCa                   <  1:153891/149‑1 (MQ=255)
                                                                                                                                 cTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCGTCg                >  2:28874/1‑149 (MQ=255)
                                                                                                                                      ccGGCCTACACTTGCCCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGa           <  1:186151/149‑1 (MQ=255)
                                                                                                                                       cGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAg          >  1:18023/1‑149 (MQ=255)
                                                                                                                                       cGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGGTTGGCACCAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGCTGAGCGGCGTTTGAg          >  2:38885/1‑149 (MQ=14)
                                                                                                                                             ataGTTGGGCATGGTTTGGGTCGCTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCCACAGTGCCGGATGCGGCGCGAGCGTCTTTTCCGGTCTATGGTTGTGCATCGTTTGAGTCACTg    <  2:63437/147‑1 (MQ=14)
                                                                                                                                               aGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCGCTGTa  >  1:535376/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
CTTATGGAAATGTTAAAAAATCGCCCCAAGTAACACCAAAGGTGTAGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACAGTTGGGCATCGTTTGAGTCACTGTCGGTCGGATAAGATGCGCAAGTATCGCATCCGACATTATTGCGGCACTGGAGTTTGGCAACAGTGCCGGATGCGGCGCGAGCGCCTTATCCGGCCTACGGTTGGGCATCGTTTGAGTCACTGTA  >  NC_000913/4153303‑4153594

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