Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,960,141 C→T 13.0% intergenic (‑71/+144) fbaA ← / ← pgk fructose‑bisphosphate aldolase, class II/phosphoglycerate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,960,1410CT13.0% 125.2 / 8.5 46intergenic (‑71/+144)fbaA/pgkfructose‑bisphosphate aldolase, class II/phosphoglycerate kinase
Reads supporting (aligned to +/- strand):  ref base C (22/18);  new base T (3/3);  total (25/21)
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.50e-01

GTCTCTTGTCCTGTATCGTCGGGCCTTATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCTCTTGTCGCCTTCCTGCAACTC  >  minE/1960071‑1960206
                                                                      |                                                                 
gTCTCTTGTCCTGTATCGTCGGGCCTTATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGc                                                                   >  1:1077152/1‑71 (MQ=255)
 tctcTTGTCCTGTATCGTCGGGCCTTATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCt                                                                  <  1:2318008/71‑1 (MQ=255)
 tctcTTGTCCTGTAGCGTCGGGCCTTATAATCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCt                                                                  <  1:276851/71‑1 (MQ=255)
  ctctTGTCCTGTATCGTCGGGCCTTATACTCGTTATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCt                                                                  >  1:3728762/1‑70 (MQ=255)
  ctctTGTCCTGTATCGTCGGGCCTTATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGc                                                                   >  1:1026831/1‑69 (MQ=255)
      tGTCCTGTATCGTCGGGCCTTATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCAcc                                                              >  1:3844460/1‑70 (MQ=255)
       gTCCTGTATCGTCGGGCCTTATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCAccc                                                             <  1:1725567/70‑1 (MQ=255)
       gTCCTGTATCGTCGGGCCTTATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCAc                                                               <  1:629279/68‑1 (MQ=255)
         ccTGTATCGTCGGGCCTTATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCAcc                                                              >  1:1938988/1‑67 (MQ=255)
            gTATCGTCGGGCCTTATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCa                                                       >  1:2949818/1‑71 (MQ=255)
            gTATCGTCGGGCCTTATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCa                                                       >  1:1306025/1‑71 (MQ=255)
            gTATCGTCGGGCCTTATACTCGGCATACTTCAAGTTGCATGTGCTGCGGCTGCGTTCGc                                                                   <  1:1544539/59‑1 (MQ=18)
             tATCGTCGGGCCTTATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTc                                                        <  1:2817802/69‑1 (MQ=255)
               tcgtcgGGCCTTATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCactt                                                    <  1:1626264/71‑1 (MQ=255)
                     ggCCTTATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCacttacat                                                >  1:219610/1‑67 (MQ=255)
                      gCCTTATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATg                                             >  1:3000152/1‑71 (MQ=40)
                      catttATACTCGTCATACTTCAAGTTTCATGTACTGCGACTGCGTTCGTTCACCCCAGTCACTTACTTAt                                              <  1:1508903/67‑1 (MQ=11)
                      catttATACTCGTCATACTTCAAGTTGCATGTACTGCGACTGCGTTCGTTCACCCCAGTCACTTACTTAt                                              <  1:3171209/67‑1 (MQ=11)
                      catttATACTCGTCATACTTCAAGTTGCATGTACTGCGACTGCGTTCGTTCACCCCAGTCACTTACTTATg                                             >  1:2087616/4‑71 (MQ=11)
                       ccTTATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGt                                            <  1:659104/71‑1 (MQ=39)
                        tttATACTCGTCATACTTCAAGTTGCATGTACTGCGACTGCGTTCGTTCa                                                                <  1:1398768/49‑1 (MQ=11)
                        tttATACTCGTCATACTTCAAGTTGCATGTACTGCGACTGCGTTCGTTCACCCCAGTCACTTACTTATGt                                            >  1:3319782/2‑70 (MQ=11)
                        cTTATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTa                                           <  1:2442243/71‑1 (MQ=38)
                         ttATACTCGTCATACTTCAAGTTGCATGTACTGCGACTGCGTTCGTTCACCCCAGTCa                                                       >  1:2032955/1‑58 (MQ=11)
                          tataCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAg                                         <  1:528091/71‑1 (MQ=34)
                          tataCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAg                                         >  1:1982469/1‑71 (MQ=34)
                           ataCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCactta                                                   >  1:724298/1‑60 (MQ=35)
                           ataCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCactta                                                   >  1:2991028/1‑60 (MQ=35)
                           ataCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTaa                                          >  1:325241/1‑69 (MQ=35)
                            taCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCAccc                                                             <  1:3448954/49‑1 (MQ=35)
                            taCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGc                                        <  1:993451/70‑1 (MQ=33)
                            taCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCt                                       <  1:182209/71‑1 (MQ=32)
                             aCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCa                                                           >  1:1160582/1‑50 (MQ=32)
                             aCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCa                                                           >  1:1284439/1‑50 (MQ=32)
                             aCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCa                                                           >  1:3700598/1‑50 (MQ=32)
                             aCTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTaa                                          <  1:983832/67‑1 (MQ=31)
                                             gTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCActctc                      <  1:812104/71‑1 (MQ=32)
                                             gTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCActctc                      >  1:2531184/1‑71 (MQ=32)
                                                gCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCTCtt                    <  1:3810529/70‑1 (MQ=35)
                                                  aTGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCTCTTGTc                 <  1:3420675/71‑1 (MQ=34)
                                                      gctgcGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCTCTTGTCGCCt             >  1:3168874/1‑71 (MQ=34)
                                                       ctgcgtctgcgtTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCTCTTGTCGCCtt            >  1:1036523/1‑71 (MQ=34)
                                                           gtctgcgtTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCTCTTGTCGCCt             >  1:1946594/1‑66 (MQ=35)
                                                             ctgcgtTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCTCTTGTCGCCTTCCTGc       >  1:2874660/1‑70 (MQ=35)
                                                              tgcgtTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCTCTTGTCGCCTTCCTGCaa     >  1:2932855/1‑71 (MQ=34)
                                                                cgtTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCTCTTGTCGCCTTCCTGCa      <  1:656424/68‑1 (MQ=35)
                                                                cgtTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCTCTTGTCGCCTTCCTGCAAc    <  1:3332329/70‑1 (MQ=35)
                                                                cgtTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCTCTTGTCGCCTTCCTGCAACt   >  1:3311305/1‑71 (MQ=34)
                                                                  ttCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCTCTTGTCGCCTTCCTGCAACTc  >  1:1636664/1‑70 (MQ=35)
                                                                      |                                                                 
GTCTCTTGTCCTGTATCGTCGGGCCTTATACTCGTCATACTTCAAGTTGCATGTGCTGCGTCTGCGTTCGCTCACCCCAGTCACTTACTTATGTAAGCTCCTGGGGATTCACTCTCTTGTCGCCTTCCTGCAACTC  >  minE/1960071‑1960206

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: