Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,249,123 G→A 12.8% intergenic (‑170/+150) treA ← / ← dhaH periplasmic trehalase/fused predicted dihydroxyacetone‑specific PTS enzyme HPr component and EI component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,249,1230GA12.8% 103.8 / 6.9 39intergenic (‑170/+150)treA/dhaHperiplasmic trehalase/fused predicted dihydroxyacetone‑specific PTS enzyme HPr component and EI component
Reads supporting (aligned to +/- strand):  ref base G (17/17);  new base A (2/3);  total (19/20)
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.34e-01

CTGTCAGGCATTTTAGAAAAGGGAATTAAAGGGATAAATAACGGCGAAATAGAGTGATCAAATAACGTCTGTTTATTAGCCAGCCAGCGGCGTTGGGTTCGTGATTAGCGAAGAAAAATCATCATAAATTGGTGAAA  >  W3110S.gb/1249057‑1249193
                                                                  |                                                                      
cTGTCAGGCATTTTAGAAAAGGGAATTAAAGGGATAAATAACGGCGAAATAGAGTGATCAAATAACGTCTg                                                                    >  1:17504/1‑71 (MQ=255)
   tCAGGCATTTTAGAAAAGGGAATTAAAGGGATAAATAACGGCGAAATAGAGTGATCAAATAACGTCTGtt                                                                  <  1:3911698/70‑1 (MQ=255)
   tCAGGCATTTTAGAAAAGGGAATTAAAGGGATAAATAACGGCGAAATAGAGTGATCAAATAACGTCTGtt                                                                  <  1:2861429/70‑1 (MQ=255)
   tCAGGCATTTTAGAAAAGGGAATTAAAGGGATAAATAACGGCGAAATAGAGTGATCAAATAACGTCTGtt                                                                  <  1:1801994/70‑1 (MQ=255)
   tCAGGCATTTTAGAAAAGGGAATTAAAGGGATAAATAACGGCGAAATAGAGTGATCAAATAACGTCTGtt                                                                  <  1:3347477/70‑1 (MQ=255)
   tCAGGCATTTTAGAAAAGGGAATTAAAGGGATAAATAACGGCGAAATAGAGTGATCAAATAACGTCTGtt                                                                  <  1:3860987/70‑1 (MQ=255)
     aGGCATTTTAGAAAAGGGAATTAAAGGGATAAATAACGGCGAAATAGAGTGATCAAATAACGTCTGTttat                                                               <  1:3954480/71‑1 (MQ=255)
     aGGCATTTTAGAAAAGGGAATTAAAGGGATAAATAACGGCGAAATAGAGTGATCAAATAACGTCTGTttat                                                               <  1:727645/71‑1 (MQ=255)
         aTTTTAGAAAAGGGAATTAAAGGGATAAATAACGGCGAAATAGAGTGATCAAATAACGTCTGTttat                                                               >  1:3028843/1‑67 (MQ=255)
         aTTTTAGAAAAGGGAATTAAAGGGATAAATAACGGCGAAATAGAGTGATCAAATAACGTCTGTTTATTagc                                                           >  1:1869194/1‑71 (MQ=255)
          ttttAGAAAAGGGAATTAAAGGGATAAATAACGGCGAAATAGAGTGATCAAATAACGTCTGtt                                                                  <  1:1937889/63‑1 (MQ=255)
          ttttAGAAAAGGGAATTAAAGGGATAAATAACGGCGAAATAGAGTGATCAAATAACGTCTGtt                                                                  <  1:761212/63‑1 (MQ=255)
            ttAGAAAAGGGAATTAAAGGGATAAATAACGGCGAAATAGAGTGATCAAATAACGTCTGTTTATTagccag                                                        >  1:3873723/1‑71 (MQ=255)
            ttAGAAAAGGGAATTAAAGGGATAAATAACGGCGAAATAGAGTGATCAAATAACGTCTGTTTATTagccag                                                        >  1:2205903/1‑71 (MQ=255)
            ttAGAAAAGGGAATTAAAGGGATAAATAACGGCGAAATAGAGTGATCAAATAACGTCTGTTTATTagccag                                                        >  1:3272731/1‑71 (MQ=255)
            ttAGAAAAGGGAATTAAAGGGATAAATAACGGCGAAATAGAGTGATCAAATAACATCTGTTTATTagccag                                                        >  1:466114/1‑71 (MQ=255)
               gAAAAGGGAATTAAAGGGATAAATAACGGCGAAATAGAGTGATCAAATAACGTCTGTTTATTagccagc                                                       <  1:2748781/69‑1 (MQ=255)
                    gggAATTAAAGGGATAAATAACGGCGAAATAGAGTGATCAAATAACGTCTGTTTATTAGCCAGCCAgcggc                                                >  1:2027805/1‑71 (MQ=255)
                                  taaataaCGGCGAAATAGAGTGATCAAATAACGTCTGttt                                                                 <  1:4048527/40‑1 (MQ=255)
                                  taaataaCGGCGAAATAGAGTGATCAAATAACGTCTGttt                                                                 <  1:30331/40‑1 (MQ=255)
                                  taaataaCGGCGAAATAGAGTGATCAAATAACGTCTGttt                                                                 <  1:1320077/40‑1 (MQ=255)
                                  taaataaCGGCGAAATAGAGTGATCAAATAACGTCTGTTTATTAGCCAGCCAGCGGCGtt                                             <  1:4068920/60‑1 (MQ=255)
                                  taaataaCGGCGAAATAGAGTGATCAAATAACGTCTGTTTATTAGCCAGCCAGCGGCGtt                                             <  1:1139253/60‑1 (MQ=255)
                                  taaataaCGGCGAAATAGAGTGATCAAATAACGTCTGTTTATTAGCCAGCCAGCGGCGTTGGGTTCGTGAt                                  >  1:667391/1‑71 (MQ=255)
                                   aaataaCGGCGAAATAGAGTGATCAAATAACGTCTGTTTATTAGCCAGCCAGCGGCGTTGGGTTCGTGAt                                  >  1:3056207/1‑70 (MQ=255)
                                    aataaCGGCGAAATAGAGTGATCAAATAACGTCTGTTTATTAGCCAGCCAGCGGCGTTGGGTTCGTGAt                                  >  1:245341/1‑69 (MQ=255)
                                    aataaCGGCGAAATAGAGTGATCAAATAACATCTGTTTATTAGCCAGCCAGCGGCGTTGGGTTCGTg                                    >  1:3406808/1‑67 (MQ=255)
                                           gCGAAATAGAGTGATCAAATAACGTCTGTTTATTAGCCAGCCAGCGGCGTTGGGTTCGTGa                                   >  1:3597083/1‑61 (MQ=255)
                                             gAAATAGAGTGATCAAATAACGTCTGTTTATTAGCCAGCCAGCGGCGTTGGGTTCGTGATTAGcgacgat                        >  1:2420044/1‑69 (MQ=255)
                                              aaaTAGAGTGATCAAATAACGTCTGTTTATTAGCCAGCCAGCGGCGTTGGGTTCGTGATTAGCGAAGaaaa                      <  1:2056480/71‑1 (MQ=255)
                                                   gagTGATCAAATAACGTCTGTTTATTAGCCAGCCAGCGGCGTTGGGTTCGTGAt                                  >  1:1257360/1‑54 (MQ=255)
                                                   gagTGATCAAATAACGTCTGTTTATTAGCCAGCCAGCGGCGTTGGGTTCGTGAt                                  >  1:397613/1‑54 (MQ=255)
                                                   gagTGATCAAATAACGTCTGTTTATTAGCCAGCCAGCGGCGTTGGGTTCGTGAt                                  >  1:2667619/1‑54 (MQ=255)
                                                        aTCAAATAACATCTGTTTATTAGCCAGCCAGCGGCGt                                              <  1:3910682/37‑1 (MQ=255)
                                                        aTCAAATAACATCTGTTTATTAGCCAGCCAGCGGCGt                                              <  1:550882/37‑1 (MQ=255)
                                                        aTCAAATAACATCTGTTTATTAGCCAGCCAGCGGCGt                                              <  1:254510/37‑1 (MQ=255)
                                                          cAAATAACGTCTGTTTATTAGCCAGCCAGCGGCGTTGGGTTCGTGATTAGCGAAGAAAAATCATCATAAAt          >  1:981739/1‑71 (MQ=255)
                                                           aaaTAACGTCTGTTTATTAGCCAGCCAGCGGCGTTGGGTTCGTGATTAGCGAAGAAAAATCATCATAAAtt         >  1:1755322/1‑71 (MQ=255)
                                                            aaTAACGTCTGTTTATTAGCCAGCCAGCGGCGTTGGGTTCGTGATTAGCGAAGAAAAATCATCATAAATTg        <  1:167303/71‑1 (MQ=255)
                                                                  gTCTGTTTATTAGCCAGCCAGCGGCGTTGGGTTCGTGATTAGCGAAGAAAAATCATCATAAATTGGTGaaa  <  1:937795/71‑1 (MQ=255)
                                                                  |                                                                      
CTGTCAGGCATTTTAGAAAAGGGAATTAAAGGGATAAATAACGGCGAAATAGAGTGATCAAATAACGTCTGTTTATTAGCCAGCCAGCGGCGTTGGGTTCGTGATTAGCGAAGAAAAATCATCATAAATTGGTGAAA  >  W3110S.gb/1249057‑1249193

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