Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,087,373 T→C 7.3% intergenic (+186/+732) PP_1865 → / ← aroH transposase/3‑deoxy‑7‑phosphoheptulonate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,087,3730TC 7.3% 157.5 / 2.5 56intergenic (+186/+732)PP_1865/aroHtransposase/3‑deoxy‑7‑phosphoheptulonate synthase
Reads supporting (aligned to +/- strand):  ref base T (27/24);  new base C (2/2);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01

GGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGT  >  NC_002947/2087232‑2087520
                                                                                                                                               |                                                                                                                                                   
ggAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATctccctc                                                                                                                                              >  5:161755/1‑149 (MQ=255)
                 gCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCg                                                                                                                             <  5:71880/149‑1 (MQ=255)
                 gCTCCCATACTCCGGCCCCTACGCTGCGGTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGATCAAGGGGTGCATGCCCCAACCTATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTc                                                                                                                              <  1:110861/149‑1 (MQ=17)
                       atcatcCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTTCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTGATCTACCTCCCCAATGCCTCCACACGGcctcc                                                                                                                       >  2:293603/1‑149 (MQ=17)
                                  ccTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCGCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGc                                                                                                                                      >  3:101028/1‑123 (MQ=255)
                                  ccTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCCCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGc                                                                                                                                      <  4:101028/123‑1 (MQ=255)
                                      cGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCt                                                                                                         <  8:197688/148‑1 (MQ=255)
                                                cGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTg                                                                                               <  1:153461/148‑1 (MQ=255)
                                                cGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCACGCCCCATACAAGGGGCCTCAGCCCTTGACGGGCATCCATGCCCGTCATCACCCAACGCAATGCCTCCACTCGGCCTCCTCAAGTCGCAATCTGCGTCGCCTg                                                                                               >  4:86807/1‑148 (MQ=17)
                                                 ggggTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCAACCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCAACTCCCACCGCAATACCTCCACTCGGCCTCCAGAAGTCGCAATCAGCGTCGCCTggc                                                                                             >  4:272272/1‑148 (MQ=14)
                                                       cccTCACTCCGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTTGGGGCCTCAGCGCTTGACGTGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGt                                                                                                     <  5:226453/136‑1 (MQ=11)
                                                        ccTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACcgc                                                                                      >  1:191662/1‑149 (MQ=255)
                                                          tCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACcgcgc                                                                                    >  2:80880/1‑149 (MQ=255)
                                                              tCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCt                                                                                 >  8:174113/1‑148 (MQ=255)
                                                                cGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCt                                                                                 >  3:256089/1‑146 (MQ=255)
                                                                    cTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAg                                                                           >  5:335850/1‑148 (MQ=255)
                                                                          cgcGGGAGCACGCGCTGGAGCCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGAGTACCGCGCGCTTAGAAGAAAGATc                                                                    <  2:245398/147‑1 (MQ=255)
                                                                          cccGCGAGCACGCGCTGCAGGCCCCAACCATTCGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAAAGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCGGGGTCTACCGCGCGCGTAGAAGCAAGATc                                                                    <  2:245395/149‑1 (MQ=17)
                                                                           ccGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATc                                                                    <  3:272272/148‑1 (MQ=255)
                                                                              ggAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCaaaa                                                                <  7:174113/149‑1 (MQ=255)
                                                                                 gCACGGGCTGCAGGCGCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATGTCGCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCg                                                             <  7:312987/149‑1 (MQ=255)
                                                                                   aCGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGAt                                                                     >  6:96144/1‑139 (MQ=255)
                                                                                   aCGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGAt                                                                     <  5:96144/139‑1 (MQ=255)
                                                                                    ccgcgCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGa                                                                                             <  1:191656/112‑1 (MQ=1)
                                                                                    ccgcgCTGCAGGCCCCATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGa                                                                                             >  2:191656/2‑113 (MQ=1)
                                                                                             aGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAAtt                                                  <  3:49770/149‑1 (MQ=255)
                                                                                              ggCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGATTACCGCGCGCTTAGAAGc                                                                          >  3:63916/1‑124 (MQ=255)
                                                                                              ggCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGATTACCGCGCGCTTAGAAGc                                                                          <  4:63916/124‑1 (MQ=255)
                                                                                                ccccATGGATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCgg                                               >  4:146469/1‑149 (MQ=255)
                                                                                                ccccATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCgg                                               >  1:224137/1‑149 (MQ=255)
                                                                                                 cccATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGt                                              >  4:291671/1‑149 (MQ=255)
                                                                                                 cccATCCTGGGGCCCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGa                                            >  1:61462/1‑149 (MQ=255)
                                                                                                  ccatccatGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGt                                              >  2:225541/1‑148 (MQ=255)
                                                                                                    atcaatGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTCGACTACCGCGCGCGTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTCAt                                           >  6:34518/1‑149 (MQ=255)
                                                                                                              gCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCa                                  >  8:195084/1‑149 (MQ=255)
                                                                                                                cTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCt                                                                                                <  7:240146/85‑1 (MQ=255)
                                                                                                                cTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCt                                                                                                >  8:240146/1‑85 (MQ=255)
                                                                                                                  cAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCa                                  >  5:172470/1‑145 (MQ=255)
                                                                                                                  cAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCa                                  <  6:172470/145‑1 (MQ=255)
                                                                                                                   aGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATcggccg                                     >  3:171549/1‑141 (MQ=255)
                                                                                                                   aGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATcggccg                                     <  4:171549/141‑1 (MQ=255)
                                                                                                                     cgcTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGa                           <  6:161755/149‑1 (MQ=255)
                                                                                                                      gcTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGAt                          <  1:293603/149‑1 (MQ=255)
                                                                                                                       cTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATg                         <  4:256089/149‑1 (MQ=255)
                                                                                                                          gACGGGCATCCATGCCCGTCATCTCCATCCGCAATGCCTCCAATAGGCCTCCTGAAGTCTCAATCTGCGTCGCCTGGACTACCGCGCGCTTACAAGCAAGATCAAAAGCGAGAAGCGAAATCGTTGATCGGCCGGCACCGCCGATGCa                       >  7:150117/1‑148 (MQ=25)
                                                                                                                               gCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGc                 >  3:187091/1‑149 (MQ=255)
                                                                                                                                   ccATGCCCGTCACCTCCCTCCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAgg             <  2:166238/149‑1 (MQ=18)
                                                                                                                                     aTGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGa               <  7:167705/145‑1 (MQ=255)
                                                                                                                                     aTGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGa               >  8:167705/1‑145 (MQ=255)
                                                                                                                                        cccGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCt         <  1:80880/148‑1 (MQ=255)
                                                                                                                                         ccGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGc       >  3:158198/1‑149 (MQ=255)
                                                                                                                                         ccGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGc       >  2:28105/1‑149 (MQ=255)
                                                                                                                                         ccGTCACCTCCCTCCGCAAGCCCTGCGCTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGc       >  1:348206/1‑149 (MQ=21)
                                                                                                                                          cGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTTGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGATGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGcc      <  5:34518/149‑1 (MQ=255)
                                                                                                                                          cGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGAAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGcc      >  4:25447/1‑149 (MQ=255)
                                                                                                                                              aTCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGt  <  4:158198/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
GGAACCAAAAACCGTCGGCTCCCATCATCCGGCCCCTACGCTGCGCTCCGGGGTCCCCTCGCTCCGGCCTTGCTCCCGGGAGCACGCGCTGCAGGCCCCATCCATGGGGCCTCAGCGCTTGACGGGCATCCATGCCCGTCATCTCCCTCCGCAATGCCTCCACTCGGCCTCCTGAAGTCGCAATCTGCGTCGCCTGGACTACCGCGCGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGT  >  NC_002947/2087232‑2087520

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