Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,703,672 G→A 33.3% A355T (GCC→ACC)  aroB → 3‑dehydroquinate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,703,6720GA33.3% 70.1 / 50.1 57A355T (GCC→ACC) aroB3‑dehydroquinate synthase
Reads supporting (aligned to +/- strand):  ref base G (18/20);  new base A (8/11);  total (26/31)
Fisher's exact test for biased strand distribution p-value = 7.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

ATTCTTCCGTTGGCAATTGGTAAGAGTGAAGTTCGCAGCGGCGTTTCGCACGAGCTTGTTCTTAACGCCATTGCCGATTGTCAATCAGCGTAACAACAAGAAAGGTCAGGCCGCTTATCAAGCGGTCTATTAGC  >  minE/2703606‑2703739
                                                                  |                                                                   
aTTCTTCCGTTGGCAATTGGTAAGAGTGAAGTTCGCAGCGGCGTTTCGCACGAGCTTGTTCTTAACGCCAt                                                                 >  1:1010129/1‑71 (MQ=255)
aTTCTTCCGTTGGCAATTGGTAAGAGTGAAGTTCGCAGCGGCGTTTCGCACGAGCTTGTTCTTAACGCCAt                                                                 >  1:1079622/1‑71 (MQ=255)
aTTCTTCCGTTGGCAATTGGTAAGAGTGAAGTTCGCAGCGGCGTTTCGCACGAGCTTGTTCTTAACACCAt                                                                 <  1:1395816/71‑1 (MQ=255)
   cttcCGTTGGCAATTGGTAAGAGTGAAGTTCGCAGCGGCGTTTCGCACGAGCTTGTTCTTAACGCCATTGc                                                              >  1:1291147/1‑71 (MQ=255)
   cttcCGTTGGCAATTGGTAAGAGTGAAGTTCGCAGCGGCGTTTCGCACGAGCTTGTTCTTAACACCATTGc                                                              >  1:1167882/1‑71 (MQ=255)
      ccGTTGGCAATTGGTAAGAGTGAAGTTCGCAGCGGCGTTTCGCACGAGCTTGTTCTTAACACCATTGCCg                                                            <  1:1158156/70‑1 (MQ=255)
             cAATTGGTAAGAGTGAAGTTCGCAGCGGCGTTTCGCACGAGCTTGTTCTTAACGCCATTGCCGATTGTCaa                                                    <  1:1538323/71‑1 (MQ=255)
              aaTTGGTAAGAGTGAAGTTCGCAGCGGCGTTTCGCACGAGCTTGTTCTTAACACCATTGCCGATTGTCAAt                                                   <  1:694697/71‑1 (MQ=255)
                  ggTAAGAGTGAAGTTCGCAGCGGCGTTTCGCACGAGCTTGTTCTTAACGCCATTGCCGATTGTCAATCAg                                                >  1:1454045/1‑70 (MQ=255)
                  ggTAAGAGTGAAGTTCGCAGCGGCGTTTCGCACGAGCTTGTTCTTAACACCATTGCCGATTGTCAATCAGc                                               >  1:865339/1‑71 (MQ=255)
                      agagTGAAGTTCGCAGCGGCGTTTCGCACGAGCTTGTTCTTAACGCCATTGCCGATTGTCAATCAGCGTa                                            <  1:1943158/70‑1 (MQ=255)
                      agagTGAAGTTCGCAGCGGCGTTTCGCACGAGCTTGTTCTTAACACCATTGCCGATTGTCAATCAGCGTaa                                           >  1:766469/1‑71 (MQ=255)
                       gagTGAAGTTCGCAGCGGCGTTTCGCACGAGCTTGTTCTTAACGCCATTGCCGATTGTCAATCAGCGTaac                                          <  1:1816978/71‑1 (MQ=255)
                         gTGAAGTTCGCAGCGGCGTTTCGCACGAGCTTGTTCTTAACACCATTGCCGATTGTCAATCAGCGTaacaa                                        >  1:1714029/1‑71 (MQ=255)
                         gTGAAGTTCGCAGCGGCGTTTCGCACGAGCTTGATCTTAACGCCATTGCCGATTGTCaa                                                    <  1:1897727/59‑1 (MQ=255)
                          tGAAGTTCGCAGCGGCGTTTCGCACGAGCTTGTTCTTAACGc                                                                    <  1:2152325/42‑1 (MQ=255)
                          tGAAGTTCGCAGCGGCGTTTCGCACGAGCTTGTTCTTAACGc                                                                    <  1:1805253/42‑1 (MQ=255)
                          tGAAGTTCGCAGCGGCGTTTCGCACGAGCTTGTTCTTAACGCCATTGCCGATTGTCAATc                                                  >  1:1075578/1‑60 (MQ=255)
                          tGAAGTTCGCAGCGGCGTTTCGCACGAGCTTGTTCTTAACGCCATTGCCGATTGTCAATc                                                  >  1:1099415/1‑60 (MQ=255)
                              gTTCGCAGCGGCGTTTCGCACGAGCTTGTTCTTAACGCCATTGCCTATTGTCAATCAGCGTAACAACAAGa                                   <  1:475335/71‑1 (MQ=255)
                              gTTCGCAGCGGCGTTTCGCACGAGCTTGTTCTTAACGCCATTGCCGATTGTCAATCAGCGTAACAACAAGa                                   <  1:1114196/71‑1 (MQ=255)
                               ttCGCAGCGGCGTTTCGCACGAGCTTGTTCTTAACGCCATTGCCGATTGTCAATCAGCGTAAcaacaa                                     <  1:653270/68‑1 (MQ=255)
                               ttCGCAGCGGCGTTTCGCACGAGCTTGTTCTTAACGCCATTGCCGATTGTCAATCAGCGTAACAACAAGa                                   <  1:854709/70‑1 (MQ=255)
                                tCGCAGCGGCGTTTCGCACGAGCTTGTTCTTAACACCATTGCCGATTGTCAATCAGCGTaac                                          >  1:1398266/1‑62 (MQ=255)
                                 cGCAGCGGCGTTTCGCACGAGCTTGTTCTTAACGCCATTGCCGATTGTCAATCAGCGTAacaaca                                      >  1:1873565/1‑65 (MQ=255)
                                    aGCGGCGTTTCGCACGAGCTTGTTCTTAACGCCATTGCCGATTGTCAATCAGCGTAACAACAAGAAAGGTc                             <  1:916186/71‑1 (MQ=255)
                                    aGCGGCGTTTCGCACGAGCTTGTTCTTAACACCATTGCCGATTGTCAATCAGCGTAACAACAAGAAAGGt                              >  1:1146193/1‑70 (MQ=255)
                                     gcggcgTTTCGCACGAGCTTGTTCTTAACGCCATTGCCGATTGTCAATCAGCGTAACAACAAGAAAGGt                              >  1:2018258/1‑69 (MQ=255)
                                     gcggcgTTTCGCACGAGCTTGTTCTTAACGCCATTGCCGATTGTCAATCAGCGTAACAACAAGAAAGGGCa                            <  1:1211077/71‑1 (MQ=255)
                                        gcgTTTCGCACGAGCTTGTTCTTAACGCCATTGCCGATTGTCAATCAGCGTAACAACAAGAAAGGTCAGGc                         >  1:1556049/1‑71 (MQ=255)
                                         cgTTTCGCACGAGCTTGTTCTTAACACCATTGCCGATTGTCAATCAGCGTAACAACAAGAAAGGTCAGGcc                        <  1:260424/71‑1 (MQ=255)
                                             tCGCACGAGCTTGTTCTTAACGCCATTGCCGATTGTCAATCAGCGTAACAACAAGAAAGGTCAGGCCGCt                     >  1:1936661/1‑70 (MQ=255)
                                             tCGCACGAGCTTGTTCTTAACACCATTGCCGATTGTCAATCAGCGTAACAACAAGAAAGGTCAGGCCGCt                     >  1:430234/1‑70 (MQ=255)
                                              cGCACGAGCTTGTTCTTAACACCATTGCCTATTGTCAATCAGCGTAACAACAAGAAAGGTCAGGCCGCtt                    <  1:1733340/70‑1 (MQ=255)
                                                  cGAGCTTGTTCTTAACACCATTGCCGATTGTCAATCAGCGTAACAACAAGa                                   >  1:442354/1‑51 (MQ=255)
                                                  cGAGCTTGTTCTTAACACCATTGCCGATTGTCAATCAGCGTAACAACAAGAAAGGTCAGGCCGCTTATCaa               <  1:348672/71‑1 (MQ=255)
                                                   gAGCTTGTTCTTAACGCCATTGCCGATTGTCAATCAGCGTAACAACAAGAAAGGTCAGGCCGCTTATCAAg              <  1:73579/71‑1 (MQ=255)
                                                   gAGCTTGTTCTTAACGCCATTGCCGATTGTCAATCAGCGTAACAACAAGAAAGGTCAGGCCGCTTATCAAg              <  1:1800314/71‑1 (MQ=255)
                                                     gCTTGTTCTTAACGCCATTGCCGATTGTCAATCAGCGTAACAACAAGAAAGGTCAGGCCGCTTATCAAGCg            <  1:273851/71‑1 (MQ=255)
                                                         gTTCTTAACGCCATTGCCGATTGTCAATCAGCGTAACAACAAGAAAGGTCAGGCCGCTTATCa                <  1:1719619/63‑1 (MQ=255)
                                                         gTTCTTAACGCCATTGCCGATTGTCAATCAGCGTAACAACAAGAAAGGTCAGGCCGCTTATCAAGCGGTc         >  1:1327563/1‑70 (MQ=255)
                                                         gTTCTTAACACCATTGCCGATTGTCAATCAGCGTAACAACAAGAAAGGTCAGGCCGCTTATCAAGCGGt          <  1:81594/69‑1 (MQ=255)
                                                         gTTCTTAACACCATTGCCGATTGTCAATCAGCGTAACAACAAGAAAGGTCAGGCCGCTTATCAAGCGGTc         <  1:1239074/70‑1 (MQ=255)
                                                         gTTCTTAACACCATTGCCGATTGTCAATCAGCGTAACAACAAGAAAGGTCAGGCCGCTTATCAAGCGGTc         <  1:1000676/70‑1 (MQ=255)
                                                         gTTCTTAACACCATTGCCGATTGTCAATCAGCGTAACAACAAGAAAGGTCAGGCCGCTTATCAAGCGGTc         <  1:1590160/70‑1 (MQ=255)
                                                          ttCTTAACGCCATTGCCGATTGTCAATCAGCGTAACAACAAGAAAGGt                              <  1:326920/48‑1 (MQ=255)
                                                          ttCTTAACGCCATTGCCGATTGTCAATCAGCGTAACAACAAGAAAGGt                              <  1:1917889/48‑1 (MQ=255)
                                                          ttCTTAACACCATTGCCGATTGTCAATCAGCGTAACAACAAGAAAgg                               <  1:1774502/47‑1 (MQ=255)
                                                            cTTAACGCCATTGCCGATTGTCAATCAGCGTAACAACAAGAAAGGTCAGGCCGCTTATCAAGCGGTCTa       >  1:2191021/1‑69 (MQ=255)
                                                            cTTAACGCCATTGCCGATTGTCAATCAGCGTAACAACAAGAAAGGTCAGGCCGCTTATCAAGCGGTCTAt      <  1:271293/70‑1 (MQ=255)
                                                            cTTAACGACATTGCCGATTGTCAATCAGCGTAACAACAAGAAAGGTCAg                           >  1:151664/1‑49 (MQ=255)
                                                             ttAACGCCATTGCCGATTGTCAATCAGCGTAACAACAAGAAAGGTCAGGCCGCTTa                   >  1:1603232/1‑56 (MQ=255)
                                                               aaCGCCATTGCCTATTGTCAATCAGCGTAACAACAAGAAAGGTCAGGCCGCTTATCAAGCGGTCTAtt     <  1:213269/68‑1 (MQ=255)
                                                               aaCGCCATTGCCGATTGTCAATCAGCGTAACAACAAGAAAg                                >  1:1931198/1‑41 (MQ=255)
                                                               aaCGCCATTGCCGATTGTCAATCAGCGTAACAACAAGAAAg                                >  1:917216/1‑41 (MQ=255)
                                                               aaCGCCATTGCCGATTGTCAATCAGCGTAACAACAAGAAAGGTCAGGCCGCTTATCAAGCGGTCTATTAGc  >  1:328447/1‑71 (MQ=255)
                                                               aaCGCCATTGCCGATTGTCAATCAGCGTAACAACAAGAAAGGTCAGGCCGCTTATCAAGCGGTCTATTAGc  >  1:1142225/1‑71 (MQ=255)
                                                                  |                                                                   
ATTCTTCCGTTGGCAATTGGTAAGAGTGAAGTTCGCAGCGGCGTTTCGCACGAGCTTGTTCTTAACGCCATTGCCGATTGTCAATCAGCGTAACAACAAGAAAGGTCAGGCCGCTTATCAAGCGGTCTATTAGC  >  minE/2703606‑2703739

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