Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,732,733 G→A 47.8% L27L (CTG→CTA yheN → predicted intracellular sulfur oxidation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,732,7330GA47.8% 6.8 / 82.1 94L27L (CTG→CTAyheNpredicted intracellular sulfur oxidation protein
Reads supporting (aligned to +/- strand):  ref base G (24/25);  new base A (33/12);  total (57/37)
Fisher's exact test for biased strand distribution p-value = 2.04e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGGTGACCGGGCCAGCATACGGTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATCGGGAAGGGGTCTATAACGCTAACCA  >  minE/2732668‑2732798
                                                                 |                                                                 
gtggtgACCGGGCCAGCATACGGTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGAt                                                                 >  1:2195531/1‑68 (MQ=255)
gtggtgACCGGGCCAGCATACGGTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGAt                                                                 >  1:1809468/1‑68 (MQ=255)
gtggtgACCGGGCCAGCATACGGTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGAt                                                                 >  1:1445162/1‑68 (MQ=255)
  ggtgACCGGGCCAGCATACGGTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGAt                                                                 >  1:2151690/1‑66 (MQ=255)
   gtgACCGGGCCAGCATACGGTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCa                                                             <  1:2151657/69‑1 (MQ=255)
   gtgACCGGGCCAGCATACGGTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCa                                                             <  1:249304/69‑1 (MQ=255)
   gtgACCGGGCCAGCATACGGTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCa                                                             <  1:2496214/69‑1 (MQ=255)
   gtgACCGGGCCAGCATACGGTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCa                                                             <  1:682520/69‑1 (MQ=255)
   gtgACCGGGCCAGCATACGGTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCa                                                             <  1:538259/69‑1 (MQ=255)
   gtgACCGGGCCAGCATACGGTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCa                                                             <  1:2952207/69‑1 (MQ=255)
   gtgACCGGGCCAGCATACGGTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCa                                                             <  1:1527725/69‑1 (MQ=255)
   gtgACCGGGCCAGCATACGGTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCa                                                             <  1:3237577/69‑1 (MQ=255)
     gACCGGGCCAGCATACGGTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCa                                                             <  1:2655137/67‑1 (MQ=255)
         gggCCAGCATACGGTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGc                                                       <  1:1679358/69‑1 (MQ=255)
         gggCCAGCATACGGTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGc                                                       <  1:1764119/69‑1 (MQ=255)
         gggCCAGCATACGGTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGc                                                       <  1:443551/69‑1 (MQ=255)
         gggCCAGCATACGGTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGc                                                       <  1:678747/69‑1 (MQ=255)
          ggCCAGCATACGGTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGc                                                       <  1:1788753/68‑1 (MQ=255)
                 aTACGGTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGa                                                  >  1:2459832/1‑66 (MQ=255)
                 aTACGGTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGa                                                  >  1:1812264/1‑66 (MQ=255)
                 aTACGGTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAg                                                 >  1:2398075/1‑67 (MQ=255)
                 aTACGGTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAg                                                 >  1:3247148/1‑67 (MQ=255)
                      gTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAagca                                          <  1:1827375/69‑1 (MQ=255)
                      gTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAagca                                          <  1:2100599/69‑1 (MQ=255)
                      gTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAagca                                          <  1:2470664/69‑1 (MQ=255)
                      gTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAagca                                          <  1:2074620/69‑1 (MQ=255)
                      gTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAagca                                          <  1:1147932/69‑1 (MQ=255)
                      gTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAagca                                          <  1:2186523/69‑1 (MQ=255)
                      gTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAagca                                          <  1:2450051/69‑1 (MQ=255)
                      gTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAagca                                          <  1:2693239/69‑1 (MQ=255)
                      gTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAagca                                          <  1:2775411/69‑1 (MQ=255)
                      gTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAagca                                          <  1:586331/69‑1 (MQ=255)
                      gTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAagca                                          <  1:1733794/69‑1 (MQ=255)
                      gTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAagca                                          <  1:2955477/69‑1 (MQ=255)
                               aGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCAt                                                    >  1:3408190/1‑50 (MQ=255)
                               aGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAagca                                          >  1:96400/1‑60 (MQ=255)
                               aGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGtttttt                                 >  1:349938/1‑69 (MQ=255)
                               aGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGtttttt                                 >  1:689858/1‑69 (MQ=255)
                               aGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGtttttt                                 >  1:2444399/1‑69 (MQ=255)
                               aGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGtttttt                                 >  1:847621/1‑69 (MQ=255)
                               aGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGtttttt                                 >  1:889350/1‑69 (MQ=255)
                               aGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGtttttt                                 >  1:1128421/1‑69 (MQ=255)
                               aGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGtttttt                                 >  1:2332599/1‑69 (MQ=255)
                               aGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGtttttt                                 >  1:2271545/1‑69 (MQ=255)
                               aGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGtttttt                                 >  1:1056466/1‑69 (MQ=255)
                               aGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGtttttt                                 >  1:2156328/1‑69 (MQ=255)
                               aGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGtttttt                                 >  1:1662117/1‑69 (MQ=255)
                               aGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGtttttt                                 >  1:2150394/1‑69 (MQ=255)
                               aGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGtttttt                                 >  1:1937723/1‑69 (MQ=255)
                               aGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGtttt                                   >  1:2469163/1‑67 (MQ=255)
                               aGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTAATAGAAGAAGGCCATGAGTTAAGCAGCGtttttt                                 >  1:896153/1‑69 (MQ=255)
                               aGGCGAGTAGTGCTATTCAGTTTGCGCAGGCGCTAAt                                                                 >  1:128446/1‑37 (MQ=38)
                                  cGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCt                              <  1:738306/69‑1 (MQ=255)
                                  cGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCt                              <  1:1562301/69‑1 (MQ=255)
                                  cGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCt                              <  1:669122/69‑1 (MQ=255)
                                  cGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCt                              <  1:1412962/69‑1 (MQ=255)
                                  cGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCt                              <  1:825245/69‑1 (MQ=255)
                                  cGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCt                              <  1:2059184/69‑1 (MQ=255)
                                  cGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCt                              <  1:2214592/69‑1 (MQ=255)
                                   gAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCt                              <  1:1738073/68‑1 (MQ=255)
                                      tagtGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATc                           >  1:83268/1‑68 (MQ=255)
                                      tagtGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATc                           >  1:1444800/1‑68 (MQ=255)
                                      tagtGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATc                           >  1:3361876/1‑68 (MQ=255)
                                        gtGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATc                           >  1:1533597/1‑66 (MQ=255)
                                        gtGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATc                           >  1:736403/1‑66 (MQ=255)
                                        gtGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATc                           >  1:604878/1‑66 (MQ=255)
                                        gtGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATc                           >  1:896478/1‑66 (MQ=255)
                                        gtGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATc                           >  1:1494815/1‑66 (MQ=255)
                                        gtGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATc                           >  1:473088/1‑66 (MQ=255)
                                        gtGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATc                           >  1:2961415/1‑66 (MQ=255)
                                        gtGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATc                           >  1:3345374/1‑66 (MQ=255)
                                          gCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATCGGGaa                      >  1:1998609/1‑69 (MQ=255)
                                           cTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATCGGGAAg                     >  1:1503190/1‑69 (MQ=255)
                                           cTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATCGGGAAg                     >  1:3002264/1‑69 (MQ=255)
                                                  gTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATCGGGAAGGGGTc                >  1:975685/1‑67 (MQ=255)
                                                  gTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATCGGGAAGGGGTCt               >  1:1521455/1‑68 (MQ=255)
                                                   tttGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATCGGGAAGGGGTCtat             <  1:2814800/69‑1 (MQ=255)
                                                   tttGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTCTCTATCGGGAAGGGGTCtat             <  1:320116/69‑1 (MQ=255)
                                                     tGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGttt                                    >  1:406741/1‑42 (MQ=255)
                                                     tGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGttt                                    >  1:2288623/1‑42 (MQ=255)
                                                     tGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGttt                                    >  1:921910/1‑42 (MQ=255)
                                                     tGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGttt                                    >  1:2156090/1‑42 (MQ=255)
                                                     tGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGttt                                    >  1:209644/1‑42 (MQ=255)
                                                     tGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGttt                                    >  1:1165384/1‑42 (MQ=255)
                                                     tGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGttt                                    >  1:1246866/1‑42 (MQ=255)
                                                     tGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGttt                                    >  1:1926366/1‑42 (MQ=255)
                                                     tGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGttt                                    >  1:698278/1‑42 (MQ=255)
                                                     tGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGttt                                    >  1:1664600/1‑42 (MQ=255)
                                                     tGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGttt                                    >  1:130830/1‑42 (MQ=255)
                                                     tGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGttt                                    >  1:449442/1‑42 (MQ=255)
                                                     tGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGttt                                    >  1:2993392/1‑42 (MQ=255)
                                                     tGCGCAGGCGCTAATAGCAGAAGGCCATGAGTTAAGCAGCGttt                                    >  1:3286169/1‑42 (MQ=255)
                                                     tGCGCAGGCGCTAATAACAGAAGGCCATGAGTTAAGCAGCGttt                                    >  1:1999057/1‑42 (MQ=255)
                                                              gcTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATCGGGAAGGGGTCTATAACGCTAACCa  <  1:903156/69‑1 (MQ=255)
                                                                 gATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATCGGGAAGGGGTCTATaa           >  1:2427476/1‑57 (MQ=255)
                                                                 |                                                                 
GTGGTGACCGGGCCAGCATACGGTACGCAACAGGCGAGTAGTGCTTTTCAGTTTGCGCAGGCGCTGATAGCAGATGGCCATGAGTTAAGCAGCGTCTTTTTCTATCGGGAAGGGGTCTATAACGCTAACCA  >  minE/2732668‑2732798

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