Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,569,078 T→G 5.3% P39P (CCA→CCC 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,0780TG 5.3% 201.2 / 4.1 75P39P (CCA→CCCmtlAfused mannitol‑specific PTS enzyme IIABC components
Reads supporting (aligned to +/- strand):  ref base T (36/35);  new base G (2/2);  total (38/37)
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.05e-01

GCAGCGGCAGGAGATAAGTGATCATCGGCCCGACCAGCTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGTGATGATACCCCACGCGATAAACGCGCCGATATTTGGCATCACCATGTTGCTG  >  minE/2569008‑2569147
                                                                      |                                                                     
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:5613050/1‑71 (MQ=255)
                                            aGCGTCTCGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGTGATGATACCCCACGCGATAAACGCGcc                           >  1:4570591/1‑71 (MQ=255)
                                                tctcGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGTGATGATACCCCACGCGATAAACGCGCCGata                       <  1:5634626/71‑1 (MQ=255)
                                                 ttcGTTAGGCAGCCAACCGGTGGGAATAAATAACGCAGTAATAAAACCCCAGGCAATAAAAGCACCAata                       >  1:2891923/2‑70 (MQ=255)
                                                  tcGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGtgatga                                                 >  1:3661003/1‑43 (MQ=255)
                                                  tcGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGtgatga                                                 >  1:2236692/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:51359/1‑71 (MQ=255)
                                                              caccCTGTTGGAATAAATAACGCGGTGATGATACCCCACGCGATAAACGCGCCGATATTTGGCATCACCAt         >  1:5082293/1‑71 (MQ=255)
                                                              caccCTGTTGGAATAAATAACGCGGTGATGATACCCCACGCGATAAACGCGCCGATATTTGGCATCACCAt         >  1:2490011/1‑71 (MQ=255)
                                                                  cTGTTGGAATAAATAACGCGGTGATGATACCCCAcgc                                       <  1:2227238/37‑1 (MQ=255)
                                                                  cTGTTGGAATAAATAACGCGGTGATGATACCCCACGCGATAAACGCGCCGATATTTGGCATCACCATGTTg     >  1:4581898/1‑71 (MQ=255)
                                                                     ttGGAATAAATAACGCGGTGATGATACCCCACGCGATAAACGCGCCGata                       <  1:4769913/50‑1 (MQ=255)
                                                                     ttGGAATAAATAACGCGGTGATGATACCCCACGCGATAAACGCGCCGATAttt                    <  1:292667/53‑1 (MQ=255)
                                                                     ttGGAATAAATAACGCGGTGATGATACCCCACGCGATAAACGCGCCGATAttt                    <  1:2644061/53‑1 (MQ=255)
                                                                     ttGGAATAAATAACGCGGTGATGATACCCCACGCGATAAACGCGCCGATAttt                    <  1:491361/53‑1 (MQ=255)
                                                                     ttGGAATAAATAACGCGGTGATGATACCCCACGCGATAAACGCGCCGATATTTGGCATCACCATGTTGCt   <  1:2837901/70‑1 (MQ=255)
                                                                     ttGGAATAAATAACGCGGTGATGATACCCCACGCGATAAACGCGCCGATATTTGGCATCACCATGTTGCt   <  1:24541/70‑1 (MQ=255)
                                                                     ttGGAATAAATAACGCGGTGATGATACCCCACGCGATAAACGCGCCGATATTTGGCATCACCATGTTGCt   <  1:1896580/70‑1 (MQ=255)
                                                                     ttGGAATAAATAACGCGGTGATGATACCCCACGCGATAAACGCGCCGATATTTGGCATCACCATGTTGCTg  <  1:618069/71‑1 (MQ=255)
                                                                     ttGGAATAAATAACGCGGTGATGATACCCCACGCGATAAACGCGCCGATATTTGGCATCACCATGTTGCTg  <  1:5613654/71‑1 (MQ=255)
                                                                     ttGGAATAAATAACGCGGTGATGATACCCCACGCGATAAACGCGCCGATATTTGGCATCACCATGTTGCTg  <  1:2391542/71‑1 (MQ=255)
                                                                      tGGAATAAATAACGCGGTGATGATACCCCACGCGATAAACGCGCCGATATTTGGCATCAc            <  1:3154253/60‑1 (MQ=255)
                                                                      |                                                                     
GCAGCGGCAGGAGATAAGTGATCATCGGCCCGACCAGCTTCGCCAGCGTCTCGTTCGGTAACCACCCTGTTGGAATAAATAACGCGGTGATGATACCCCACGCGATAAACGCGCCGATATTTGGCATCACCATGTTGCTG  >  minE/2569008‑2569147

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: