Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 786,462 (A)5→4 coding (830/1053 nt) aroG → 3‑deoxy‑7‑phosphoheptulonate synthase, Phe‑sensitive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913786,4580A.97.1% 129.5 / 0.1 36coding (826/1053 nt)aroG3‑deoxy‑7‑phosphoheptulonate synthase, Phe‑sensitive
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base . (18/15);  total (19/17)
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.37e-01

CCATATCATTCTGCGCGGCGGTAAAGAGCCTAACTACAGCGCGAAGCACGTTGCTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGT  >  NC_000913/786319‑786595
                                                                                                                                           |                                                                                                                                         
ccATATCATTCTGCGCGGCGGTAAAGAGCCTAACTACAGCGCGAAGCGCGTGGCTGAAGTGAAAGAAGGGCCGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCaaaa                                                                                                                                        <  1:356505/143‑1 (MQ=255)
ccATATCATTCTGCGCGGCGGTAAAGAGCCTAACTACAGCGCGAAGCACGTTGCTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCaaaa                                                                                                                                        <  1:310664/143‑1 (MQ=255)
  agagcATTCTGCGCGGCGGTAAAGAGCCTAACTACAGCGCGAAGCACGTTGCTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCaaaagc                                                                                                                                      <  2:492188/139‑3 (MQ=255)
    atCATTCTGCGCGGCGGTAAAGAGCCTAACTACAGCGCGAAGCACGTTGCTGAAGTGAAAGAAGGGCTGAACAAAAAAGCACTTTAAGGGCAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCaaaagcag                                                                                                                                    >  1:570450/1‑139 (MQ=255)
      cATTCTGCGCGTCGGTAAGGAGCCTTACTACGGCACGAAGCACGTTGCTGAAGTGAAAGAAGGGTTGAACATAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAgat                                                                                                                                 <  2:572525/143‑1 (MQ=255)
        ttCTGCGCGGCGGTAAAGAGCCTAACTACAGCGCGAAGCACGTTGCTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAgatgg                                                                                                                               <  1:2984/143‑1 (MQ=255)
              gcggcgGTAAAGAGCCTAACTACAGCGCGAAGCACGTTGCTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGttt                                                                                                                         <  1:236587/143‑1 (MQ=255)
                  cggTAAAGAGCCTAACTACAGCGCGAAGCACGTTGCTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGttt                                                                                                                         <  1:400692/139‑1 (MQ=255)
                  cggTAAAGAGCCTAACTACAGCGCGAAGCACGTTGCTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGttt                                                                                                                         >  2:400692/1‑139 (MQ=255)
                      aaaGAGCCTAACTACAGCGCGAAGCACGTTGCTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAACAGATGGGTGTTTGTGCTgcc                                                                                                                 >  1:457058/1‑141 (MQ=255)
                        agagCCTAACTACAGCGCGAAGCACGTTGCTGAAGTGAAAGAAGGGCTGAACAAAGTAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACACCTCCAAAGCCGATGGATGCTTGTGCTGACGt                                                                                                               >  2:731819/1‑143 (MQ=255)
                           ggtgaaCGACAGCGCGAAGGACGTCGCTGGAGTGAAAGGAGTGCTGAACCAAGCAGGCCTGCCAGGACCGGCGATGGTCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTg                                                                                                            <  2:348051/139‑1 (MQ=255)
                            ccTAACTACAGCGCGAAGCACGTTGCTGAAGTGACAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACACGTCCAAACAATTCCAAACCAGAGGGATGTTTGTGCTGACGCTTGc                                                                                                           >  2:25514/1‑143 (MQ=255)
                               accTACAGCGCGAAGCACGTTGCTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCcag                                                                                                        <  1:758508/141‑1 (MQ=255)
                                  tACAGCGCGAAGCACGTTGCTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTCACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGCCGTTTGCcagcag                                                                                                     >  2:178150/1‑143 (MQ=255)
                                  ggcAGCGCGAAGCACGTGCCTGGACGGATAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCTACCAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCcagcag                                                                                                     <  2:228110/141‑1 (MQ=255)
                                    cAGCGCGAAGCACGTTGCTGAAGTGAAAGAAGGGCTGCACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCCAAAGCAGATGGACGTTTGTGCTGACGTTTCCccgcagct                                                                                                   >  2:171495/1‑141 (MQ=255)
                                          gAAGCACGTTGCTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCgg                                                                                             >  1:734282/1‑143 (MQ=255)
                                                           tGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCACAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCatta                                                                            >  2:361793/1‑143 (MQ=255)
                                                             aaaGAAGGGCTGAACAAAGCAGGCATGCAAGCACGGGGGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATGGCCGGGGGCGAAAAGGCCatgatt                                                                          >  1:403832/1‑143 (MQ=255)
                                                                  agGGCTGAACAAAGCAGGCCGGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGt                                                                     <  1:460054/143‑1 (MQ=255)
                                                                       tGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGAGGCGAAAAGGCCATTATTGGCGTGAtgg                                                                >  2:449501/1‑143 (MQ=255)
                                                                         aaCAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGAtggtg                                                              >  1:464238/1‑143 (MQ=255)
                                                                            aaaGCAGGCCTGCCAGCACAGGTGATGGTCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGaa                                                           <  1:153031/143‑1 (MQ=255)
                                                                              aGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACCCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGAGGGTGGaatg                                                         >  1:497855/1‑141 (MQ=255)
                                                                                   gCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAGGCCATc                                                    >  2:688512/1‑143 (MQ=255)
                                                                                   gCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATc                                                    <  1:401393/143‑1 (MQ=255)
                                                                                      tGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCAATTATTGGCGTGATGGTGGAACGCCCTCtgg                                                 >  1:335738/1‑143 (MQ=255)
                                                                                       gccacgacaGGTGATGATCGCTCCCAGCCATGCTAACTCGTGAAAACAGTTGAAAAGCAGATCGATGTTCGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCtggt                                                <  1:212133/137‑1 (MQ=255)
                                                                                                     tgatCGATTTCAGCCATGCTAACTCGTCCAAACAATTCGAAAGCAGATGGATGTTTGTGCTGACGTGTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATGATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCaga                                  <  2:759732/143‑1 (MQ=255)
                                                                                                        tCGATTTCAGCCATGCTAACTCGTCCAAAAAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGAAGGTGGAAAGCCATCTGGCGGACGGCAATCAGAGcc                               >  2:683902/1‑143 (MQ=255)
                                                                                                          gATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTc                             <  1:515816/143‑1 (MQ=255)
                                                                                                                aGCCATGCTAACTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCGGGTGGAAGGCAATCAGAGCCTAgagacc                       >  1:502461/1‑141 (MQ=255)
                                                                                                                      gCTAACTCGTCCAAACAATTCAATAGCAGATGGATGTTTGTGCTGACGTTTGCCAGGAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAg                 <  1:734696/143‑1 (MQ=255)
                                                                                                                         aaCTCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCg              <  1:43339/143‑1 (MQ=255)
                                                                                                                            tCGTGCAAGCACTTCAAAAGAAGATGGATGTTGGTGCTGACGTTTGCCAGCAGATGGCCAGTCGCTAAAGTTCCATTCTGGGCGTGATGGTGGGAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTg           <  2:192823/143‑1 (MQ=255)
                                                                                                                            tCGTCCAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTg           >  2:431895/1‑143 (MQ=255)
                                                                                                                               gccAAACAATTCTGGAGCAGATGGATGTGTGTGCTGATGTTTGCCACAAGATCGCCTGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGGAAGCCAGTTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGcc        <  2:88562/142‑1 (MQ=255)
                                                                                                                                ccAAACAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGACATTATTGGCGTGATGGTGGAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTa      >  1:267976/1‑143 (MQ=255)
                                                                                                                                     cAATTCAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGCAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGGGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCAGGCCTaagag  >  2:81542/1‑139 (MQ=255)
                                                                                                                                           |                                                                                                                                         
CCATATCATTCTGCGCGGCGGTAAAGAGCCTAACTACAGCGCGAAGCACGTTGCTGAAGTGAAAGAAGGGCTGAACAAAGCAGGCCTGCCAGCACAGGTGATGATCGATTTCAGCCATGCTAACTCGTCCAAACAATTCAAAAAGCAGATGGATGTTTGTGCTGACGTTTGCCAGCAGATTGCCGGTGGCGAAAAGGCCATTATTGGCGTGATGGTGGAAAGCCATCTGGTGGAAGGCAATCAGAGCCTCGAGAGCGGGGAGCCGCTGGCCTACGGT  >  NC_000913/786319‑786595

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