Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,569,075 T→G 5.2% T40T (ACA→ACC mtlA ← fused mannitol‑specific PTS enzyme IIABC components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,569,0750TG 5.2% 204.3 / 4.0 77T40T (ACA→ACCmtlAfused mannitol‑specific PTS enzyme IIABC components
Reads supporting (aligned to +/- strand):  ref base T (37/36);  new base G (2/2);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

CAGCAGCGGCAGGAGATAAGTGATCATCGGCCCGACCAGCTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGTGATGATACCCCACGCGATAAACGCGCCGATATTTGGCATCACCATGTTG  >  minE/2569006‑2569144
                                                                     |                                                                     
cagcagCGGCAGGAGATAAGTGATCATCGGCCCGACCAGCTTCGCCAGCGTCTCGTTCGGTAACCACCCTg                                                                      >  1:4559175/1‑71 (MQ=255)
  gcagcGGCAGGAGATAAGTGATCATCGGCCCGACCAGCTTCGCCAGCGTCTCGTTCGGTAACCACCCTGtt                                                                    <  1:3973611/71‑1 (MQ=255)
        gCAGGAGATAAGTGTTCATCGGCCCGACCAGCTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGaat                                                               >  1:5135964/1‑70 (MQ=255)
         cAGGAGATAAGTGATCATCGGCCCGACCAGCTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGaataa                                                             <  1:3059095/71‑1 (MQ=255)
         cAGGAGATAAGTGATCATCGGCCCGACCAGCTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGa                                                                 >  1:3751525/1‑67 (MQ=255)
          aGGAGATAAGTGATCATCGGCCCGACCAGCTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGaataaa                                                            >  1:2261688/1‑71 (MQ=255)
           ggAGATAAGTGATCATCTGCCCGACCAGCTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGaataaat                                                           <  1:4093894/71‑1 (MQ=255)
            gagaTAAGTGATCATCGGCCCGACCAGCTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGaataaata                                                          <  1:776289/71‑1 (MQ=255)
              gaTAAGTGATCATCGGCCCGACCAGCTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGAATAAATAAc                                                        <  1:5445046/71‑1 (MQ=255)
               aTAAGTGATCATCGGCCCGACCAGCTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGAATAAATAAcg                                                       >  1:1841671/1‑71 (MQ=255)
                tAAGTGATCATCGGCCCGACCAGCTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGaataaat                                                           <  1:4488618/66‑1 (MQ=255)
                 aaGTGATCATCGGCCCGACCAGCTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGAataaataa                                                         >  1:174590/1‑67 (MQ=255)
                     gATCATCGGCCCGACCAGCTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGtga                                                 >  1:5623203/1‑71 (MQ=255)
                      atcatcGGCCCGACCAGCTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTg                                                                   <  1:940132/52‑1 (MQ=255)
                      atcatcGGCCCGACCAGCTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTg                                                                   <  1:2120624/52‑1 (MQ=255)
                      atcatcGGCCCGACCAGCTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGtgat                                                >  1:2336609/1‑71 (MQ=255)
                       tcatcGGCCCGACCAGCTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGtgatg                                               <  1:3854699/71‑1 (MQ=255)
                         atcGGCCCGACCAGCTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGtgat                                                >  1:3869856/1‑68 (MQ=255)
                          tcGGCCCGACCAGCTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGaataaat                                                           <  1:544750/56‑1 (MQ=255)
                          tcGGCCCGACCAGCTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGaataaat                                                           <  1:2686657/56‑1 (MQ=255)
                          tcGGCCCGACCAGCTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGtgatga                                              <  1:2582751/69‑1 (MQ=255)
                          tcGGCCCGACCAGCTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGTGATGATa                                            <  1:2056632/71‑1 (MQ=255)
                           cGGCCCGACCAGCTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGaataaat                                                           <  1:5166023/55‑1 (MQ=255)
                                 gACCAGCTTCGCGAGCGTCTCGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGtgat                                                <  1:4067389/60‑1 (MQ=255)
                                 gACCAGCTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGtgat                                                <  1:2157232/60‑1 (MQ=255)
                                 gACCAGCTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGtgat                                                <  1:5670246/60‑1 (MQ=255)
                                  aCCAGCTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGaataa                                                             >  1:1691802/1‑46 (MQ=255)
                                  aCCAGCTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGAATAAATAAcg                                                       >  1:316555/1‑52 (MQ=255)
                                  aCCAGCTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGTGATGATACCCCAcgc                                    <  1:3845802/71‑1 (MQ=255)
                                  aCCAGCTTCCCCAGCGTCTCGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGTGATGATACCCCAcgc                                    <  1:1166793/71‑1 (MQ=255)
                                   ccAGCTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGTGATGATACCCCAcgc                                    >  1:4679329/1‑70 (MQ=255)
                                    cAGCTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGtga                                                 >  1:4446481/1‑56 (MQ=255)
                                    cAGCTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGTGATGATACCCCACGCGa                                  >  1:67731/1‑71 (MQ=255)
                                    aaGCTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGtga                                                 >  1:70154/2‑56 (MQ=255)
                                     aGCTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGTGATGATACCCCACGCGAt                                 <  1:5698130/71‑1 (MQ=255)
                                       cTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGaataaat                                                           <  1:1243928/43‑1 (MQ=255)
                                       cTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGaataaat                                                           <  1:4327753/43‑1 (MQ=255)
                                       cTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGaat                                                               <  1:5641014/39‑1 (MQ=255)
                                       cTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGtgatga                                              >  1:2670525/1‑56 (MQ=255)
                                       cTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGTGATGATACCCCACGCGATa                                >  1:2675980/1‑70 (MQ=255)
                                        ttCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGTGATGATAccc                                         >  1:4633505/1‑60 (MQ=255)
                                        ttCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGTGATGATAccc                                         >  1:648603/1‑60 (MQ=255)
                                        ttCGACAGCGTCTCGTTCGGTAACCACCCTGTTGGAATAAATAAcg                                                       >  1:2311176/1‑46 (MQ=255)
                                         tCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGAataaataa                                                         >  1:2835131/1‑43 (MQ=255)
                                          cGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGTGATGATACCCCACGCGATAAAcg                            >  1:2617703/1‑71 (MQ=255)
                                          cGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGAAAAAATAACGCTGTGATGATACGCCACGCGATa                                <  1:29603/67‑1 (MQ=255)
                                           gCCAGCGTCTCGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGTGATGATACCCCACGCGATAAAcgc                           <  1:2425969/71‑1 (MQ=255)
                                              aGCGTCTCGTTCGTTAACCACCCTGTTGGAATAAATAACGCGGtgat                                                <  1:1450676/47‑1 (MQ=255)
                                              aGCGTCTCGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGTGATGATACCCCACGCGATAAACGCGcc                        >  1:4570591/1‑71 (MQ=255)
                                              aGCGTCTCGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGTGATGATACCCCACGCGATAAACGCGcc                        >  1:5613050/1‑71 (MQ=255)
                                                  tctcGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGTGATGATACCCCACGCGATAAACGCGCCGata                    <  1:5634626/71‑1 (MQ=255)
                                                   ttcGTTAGGCAGCCAACCGGTGGGAATAAATAACGCAGTAATAAAACCCCAGGCAATAAAAGCACCAata                    >  1:2891923/2‑70 (MQ=255)
                                                    tcGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGtgatga                                              >  1:2236692/1‑43 (MQ=255)
                                                    tcGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGtgatga                                              >  1:3661003/1‑43 (MQ=255)
                                                    tcGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGtgatg                                               >  1:4925902/1‑42 (MQ=255)
                                                    tcGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGtgatg                                               >  1:912037/1‑42 (MQ=255)
                                                    tcGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGtgat                                                >  1:1030171/1‑41 (MQ=255)
                                                    tcGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGt                                                   <  1:4253430/38‑1 (MQ=255)
                                                    tcGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGt                                                   <  1:2760253/38‑1 (MQ=255)
                                                    tcGTTCGGTAAACACCCTGTTGGAATAAATAACGCGGtgatg                                               >  1:1065264/1‑42 (MQ=255)
                                                     cGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGTGATGATACCCCACGCGATAAACGCGCCGATAttt                 >  1:4618268/1‑71 (MQ=255)
                                                     cGTTAGGCAGCCAACCGGTGGGAATAAATAACGCAGTAATAAAACCCCAGGCAATAAAAGCACCAATAtt                  <  1:2015882/70‑1 (MQ=255)
                                                      gTTCGGTAACCACCCTGTTGGAATAAATAACGCGGTGATGATACCCCACGCGATAAACGCGCCGATAttt                 <  1:517938/70‑1 (MQ=255)
                                                       ttCGGTAACCACCCTGTTGGAATAAATAACGCGGtgat                                                <  1:1339191/38‑1 (MQ=255)
                                                       ttCGGTAACCACCCTGTTGGAATAAATAACGCGGtgat                                                <  1:4801529/38‑1 (MQ=255)
                                                       ttCGGTAACCACCCTGTTGGAATAAATAACGCGGtgat                                                <  1:821719/38‑1 (MQ=255)
                                                       ttCGGTAACCACCCTGTTGGAATAAATAACGCGGTGATGATACCCCACGCGATAAACGCGCCGatat                   <  1:4230686/67‑1 (MQ=255)
                                                         cGGTAACCACCCTGTTGGAATAAATAACGCGGtga                                                 >  1:3114775/1‑35 (MQ=255)
                                                         cGGTAACCACCCTGTTGGAATAAATAACGCGGTGATGATACCCCACGCGATAAACGCGCCGATATTTGGc              <  1:4888212/70‑1 (MQ=255)
                                                              gccaacCGGTGGGAATAAATAACGCAGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGAtg           <  1:2863261/67‑2 (MQ=255)
                                                               ccaacCGGTGGGAATAAATAACGCAGTAATAAAACCCCAGGCAATAAAAGCACCAATATTGGGGATGACCa       >  1:4707540/1‑71 (MQ=255)
                                                                caccCTGTTGGAATAAATAACGCGGTGATGATACCCCACGCGATAAACGCGCCGata                    >  1:255166/1‑57 (MQ=255)
                                                                caccCTGTTGGAATAAATAACGCGGTGATGATACCCCACGCGATAAACGCGCCGATATTTGGCATCACCAt      >  1:2490011/1‑71 (MQ=255)
                                                                caccCTGTTGGAATAAATAACGCGGTGATGATACCCCACGCGATAAACGCGCCGATATTTGGCATCACCAt      >  1:51359/1‑71 (MQ=255)
                                                                caccCTGTTGGAATAAATAACGCGGTGATGATACCCCACGCGATAAACGCGCCGATATTTGGCATCACCAt      >  1:5082293/1‑71 (MQ=255)
                                                                    cTGTTGGAATAAATAACGCGGTGATGATACCCCAcgc                                    <  1:2227238/37‑1 (MQ=255)
                                                                    cTGTTGGAATAAATAACGCGGTGATGATACCCCACGCGATAAACGCGCCGATATTTGGCATCACCATGTTg  >  1:4581898/1‑71 (MQ=255)
                                                                     |                                                                     
CAGCAGCGGCAGGAGATAAGTGATCATCGGCCCGACCAGCTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGTGATGATACCCCACGCGATAAACGCGCCGATATTTGGCATCACCATGTTG  >  minE/2569006‑2569144

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: