Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,855,562 T→C 51.3% A56A (GCA→GCG yjeS ← predicted Fe‑S electron transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,855,5620TC51.3% 2.8 / 71.4 82A56A (GCA→GCGyjeSpredicted Fe‑S electron transport protein
Reads supporting (aligned to +/- strand):  ref base T (22/18);  new base C (27/15);  total (49/33)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGATCACGCGCAGCGTACCGGGCAATAACTCATGAGGGCGAGCGCGCAGCATACCGTGACGTGCCATCCAGTCCATTTCGCCGTGGTATTGTTTGTCCAGCCATGCTTGCAGTTTGGGCTCGGACTCGCT  >  minE/2855500‑2855630
                                                              |                                                                    
cTGATCACGCGCAGCGTACCGGGCAATAATTCGTGAGGACGAGCACGCAGCATACCGTGACGCGCCa                                                                  >  1:448005/1‑67 (MQ=255)
cTGATCACGCGCAGCGTACCGGGCAATAATTCGTGAGGACGAGCACGCAGCATACCGTGACGCGCCa                                                                  >  1:1834429/1‑67 (MQ=255)
cTGATCACGCGCAGCGTACCGGGCAATAATTCGTGAGGACGAGCACGCAGCATACCGTGACGCGCCa                                                                  >  1:3157092/1‑67 (MQ=255)
cTGATCACGCGCAGCGTACCGGGCAATAATTCGTGAGGACGAGCACGCAGCATACCGTGACGCGCCa                                                                  >  1:320041/1‑67 (MQ=255)
cTGATCACGCGCAGCGTACCGGGCAATAATTCGTGAGGACGAGCACGCAGCATACCGTGACGCGCCa                                                                  >  1:2956647/1‑67 (MQ=255)
cTGATCACGCGCAGCGTACCGGGCAATAATTCGTGAGGACGAGCACGCAGCATACCGTGACGCGCCa                                                                  >  1:2913820/1‑67 (MQ=255)
cTGATCACGCGCAGCGTACCGGGCAATAATTCGTGAGGACGAGCACGCAGCATACCGTGACGCGCCa                                                                  >  1:1981507/1‑67 (MQ=255)
cTGATCACGCGCAGCGTACCGGGCAATAATTCGTGAGGACGAGCACGCAGCATACCGTGACGCGCCa                                                                  >  1:3323941/1‑67 (MQ=255)
cTGATCACGCGCAGCGTACCGGGCAATAATTCGTGAGGACGAGCACGCAGCATACCGTGACGCGCCa                                                                  >  1:2875527/1‑67 (MQ=255)
cTGATCACGCGCAGCGTACCGGGCAATAATTCGTGAGGACGAGCACGCAGCATACCGTGACGCGCCa                                                                  >  1:2868119/1‑67 (MQ=255)
cTGATCACGCGCAGCGTACCGGGCAATAATTCGTGAGGACGAGCACGCAGCATACCGTGACGCGCCa                                                                  >  1:1524430/1‑67 (MQ=255)
cTGATCACGCGCAGCGTACCGGGCAATAATTCGTGAGGACGAGCACGCAGCATACCGTGACGCGCCa                                                                  >  1:3356041/1‑67 (MQ=255)
cTGATCACGCGCAGCGTACCGGGCAATAATTCGTGAGGACGAGCACGCAGCATACCGTGACGCGCCa                                                                  >  1:351204/1‑67 (MQ=255)
cTGATCACGCGCAGCGTACCGGGCAATAATTCGTGAGGACGAGCACGCAGCATACCGTGACGCGCCa                                                                  >  1:2114643/1‑67 (MQ=255)
cTGATCACGCGCAGCGTACCGGGCAATAATTCGTGAGGACGAGCACGCAGCATACCGTGACGCGCCa                                                                  >  1:212050/1‑67 (MQ=255)
cTGATCACGCGCAGCGTACCGGGCAATAATTCGTGAGGACGAGCACGCAGCATACCGTGACGCGCCa                                                                  >  1:634352/1‑67 (MQ=255)
cTGATCACGCGCAGCGTACCGGGCAATAATTCGTGAGGACGAGCACGCAGCATACCGTGACGCGCCa                                                                  >  1:1235635/1‑67 (MQ=255)
cTGATCACGCGCAGCGTACCGGGCAATAATTCGTGAGGACGAGCACGCAGCATACCGTGACGCGCCa                                                                  >  1:1223239/1‑67 (MQ=255)
cTGATCACGCGCAGCGTACCGGGCAATAATTCGTGAGGACGAGCACGCAGCATACCGTGACGCGCCa                                                                  >  1:814945/1‑67 (MQ=255)
cTGATCACGCGCAGCGTACCGGGCAATAATTCGTGAGGACGAGCACGCAGCATACCGTGACGCGCCa                                                                  >  1:817203/1‑67 (MQ=255)
cTGATCACGCGCAGCGTACCGGGCAATAATTCGTGAGGACGAGCACGCAGCATACCGTGACGCGCCa                                                                  >  1:272140/1‑67 (MQ=255)
cTGATCACGCGCAGCGTACCGGGCAATAATTCGTGAGGACGAGCACGCAGCATACCGTGACGCGCCa                                                                  >  1:1105908/1‑67 (MQ=255)
cTGATCACGCGCAGCGTACCGGGCAATAATTCGTGAGGACGAGCACGCAGCATACCGTGACGCGCCa                                                                  >  1:2191906/1‑67 (MQ=255)
cTGATCACGCGCAGCGTACCGGGCAATAATTCGTGAGGACGAGCACGCAGCATACCGTGACGCGCCa                                                                  >  1:2289477/1‑67 (MQ=255)
cTGATCACGCGCAGCGTACCGGGCAATAATTCGTGAGGACGAGCACGCAGCATACCGTGACGCGCCa                                                                  >  1:1029644/1‑67 (MQ=255)
cTGATCACGCGCAGCGTACCGGGCAATAACTCATGAGGGCGAGCGCGCAGCATACCGTGACGTGCCa                                                                  >  1:987121/1‑67 (MQ=255)
cTGATCACGCGCAGCGTACCGGCCAATAATTCGTGAGGACGAGCACGCAGCATACCGTGACGCGCCa                                                                  >  1:2425370/1‑67 (MQ=255)
cTGATCACGCGCAGCGTACAGGGCAATAATTCGTGAGGACGAGCACGCAGCATACCGTGACGCGCCa                                                                  >  1:1552137/1‑67 (MQ=255)
     cACGCGCAGCGTACCGGGCAATAACTCATGAGGGCGAGCGCGCAGCATACCGTGACGTGCCATCCAGTc                                                           >  1:3022753/1‑69 (MQ=255)
     cACGCGCAGCGTACCGGGCAATAACTCATGAGGGCGAGCGCGCAGCATACCGTGACGTGCCATCCAGTc                                                           >  1:1756249/1‑69 (MQ=255)
     cACGCGCAGCGTACCGGGCAATAACTCATGAGGGCGAGCGCGCAGCATACCGTGACGTGCCATCCAGTc                                                           >  1:3326531/1‑69 (MQ=255)
     cACGCGCAGCGTACCGGGCAATAACTCATGAGGGCGAGCGCGCAGCATACCGTGACGTGCCATCCAGTc                                                           >  1:605929/1‑69 (MQ=255)
     cACGCGCAGCGTACCGGGCAATAACTCATGAGGGCGAGCGCGCAGCATACCGTGACGTGCCATCCAGTc                                                           >  1:2550294/1‑69 (MQ=255)
     cACGCGCAGCGTACCGGGCAATAACTCATGAGGGCGAGCGCGCAGCATACCGTGACGTGCCATCCAGTc                                                           >  1:9102/1‑69 (MQ=255)
                       cAATAACTCATGAGGGCGAGCGCGCAGCATACCGTGACGTGCCATCCAGTCCATTTCGCCGTTGTAttg                                         <  1:1835434/69‑1 (MQ=255)
                       cAATAACTCATGAGGGCGAGCGCGCAGCATACCGTGACGTGCCATCCAGTCCATTTCGCCGTGGTAttg                                         <  1:301540/69‑1 (MQ=255)
                       cAATAACTCATGAGGGCGAGCGCGCAGCATACCGTGACGTGCCATCCAGTCCATTTCGCCGTGGTAttg                                         <  1:1777711/69‑1 (MQ=255)
                       cAATAACTCATGAGGGCGAGCGCGCAGCATACCGTGACGTGCCATCCAGTCCATTTCGCCGTGGTAttg                                         <  1:1774300/69‑1 (MQ=255)
                       cAATAACTCATGAGGGCGAGCGCGCAGCATACCGTGACGTGCCATCCAGTCCATTTCGCCGTGGTAttg                                         <  1:611475/69‑1 (MQ=255)
                       cAATAACTCATGAGGGCGAGCGCGCAGCATACCGTGACGTGCCATCCAGTCCATTTCGCCGTGGTAttg                                         <  1:1123981/69‑1 (MQ=255)
                       cAATAACTCATGAGGGCGAGCGCGCAGCATACCGTGACGTGCCATCCAGTCCATTTCGCCGTGGTAttg                                         <  1:949609/69‑1 (MQ=255)
                       cAATAACTCATGAGGGCGAGCGCGCAGCATACCGTGACGTGCAATCCAGTCCATTTCGCCGTGGTAttg                                         <  1:2281555/69‑1 (MQ=255)
                           aaCTCATGAGGGCGAGCGCGCAGCATACCGTGACGTGCCATCCAGTCCATTTCGCCGTg                                               >  1:1637824/1‑59 (MQ=255)
                           aaCTCATGAGGGCGAGCGCGCAGCATACCGTGACGTGCCATCCAGTCCATTTCGCCGTGGTAttgtttg                                     >  1:2269928/1‑69 (MQ=255)
                           aaCTCATGAGGGCGAGCGCGCAGCATACCGTGACGTGCCATCCAGTCCATTTCGCCGTGGTAttgtttg                                     >  1:1076416/1‑69 (MQ=255)
                           aaCTCATGAGGGCGAGCGCGCAGCATACCGTGACGTGCCATCCAGTCCATTTCGCCGTGGTAttgtttg                                     >  1:1168397/1‑69 (MQ=255)
                           aaCTCATGAGGGCGAGCGCGCAGCATACCGTGACGTGCCATCCAGTCCATTTCGCCGTGGTAttgtttg                                     >  1:3198122/1‑69 (MQ=255)
                           aaCTCATGAGGGCGAGCGCGCAGCATACCGTGACGTGCCATCCAGTCCATTTCGCCGTGGTAttgtttg                                     >  1:3186475/1‑69 (MQ=255)
                           aaCTCATGAGGGCGAGCGCGCAGCATACCGTGACGTGCCATCCAGTCCATTTCGCCGTGGTAttgtttg                                     >  1:1885459/1‑69 (MQ=255)
                           aaCTCATGAGGGCGAGCGCGCAGCATACCGTGACGTGCCATCCAGTCCATTTCGCCGTGGTAttgtttg                                     >  1:209813/1‑69 (MQ=255)
                           aaCTCATGAGGGCGAGCGCGCAGCATACCGTGACGTGCCATCCAGTCCATTTCGCCGTGGTAttgtttg                                     >  1:2226114/1‑69 (MQ=255)
                             cTCATGAGGGCGAGCGCGCAGCATACCGTGACGTGCCATCCAGTCCATTTCGCCGTGGt                                             >  1:2237174/1‑59 (MQ=255)
                             cTCATGAGGGCGAGCGCGCAGCATACCGTGACGTGCCATCCAGTCCATTTCGCCGTGGt                                             >  1:2813615/1‑59 (MQ=255)
                             cTCATGAGGGCGAGCGCGCAGCATACCGTGACGTGCCATCCAGTCCATTTCGCCGTGGt                                             >  1:3026730/1‑59 (MQ=255)
                             cTCATGAGGGCGAGCGCGCAGCATACCGTGACGTGCCATCCAGTCCATTTCGCCGTGGt                                             >  1:221195/1‑59 (MQ=255)
                                  gAGGGCGAGCGCGCAGCATACCGTGACGTGCCATCCAGTCCATTTCGCCGTGGTATTGTTTGTCCAGcc                              >  1:649360/1‑69 (MQ=255)
                                          gcgcgcAGCATACCGTGACGTGCCATCCAGTCCATTTCGCCGTGGTATTGTTTGTCCAGCCATGCTTGc                      <  1:1801211/69‑1 (MQ=255)
                                           cgcgcAGCATACCGTGACGTGCCATCCAGTCCATTTCGCCGTGGTATTGTTTGTCCAGCCATGCTTGCa                     <  1:1495813/69‑1 (MQ=255)
                                           cgcgcAGCATACCGTGACGTGCCATCCAGTCCATTTCGCCGTGGTATTGTTTGTCCAGCCATGCTTGCa                     <  1:1103648/69‑1 (MQ=255)
                                           cgcgcAGCATACCGTGACGTGCCATCCAGTCCATTTCGCCGTGGTATTGTTTGTCCAGCCATGCTTGCa                     <  1:888653/69‑1 (MQ=255)
                                           cgcgcAGCATACCGTGACGTGCCATCCAGTCCATTTCGCCGTGGTATTGTTTGTCCAGCCATGCTTGCa                     <  1:847136/69‑1 (MQ=255)
                                           cgcgcAGCATACCGTGACGTGCCATCCAGTCCATTTCGCCGTGGTATTGTTTGTCCAGCCATGCTTGCa                     <  1:608478/69‑1 (MQ=255)
                                           cgcgcAGCATACCGTGACGTGCCATCCAGTCCATTTCGCCGTGGTATTGTTTGTCCAGCCATGCTTGCa                     <  1:3249761/69‑1 (MQ=255)
                                           cgcgcAGCATACCGTGACGTGCCATCCAGTCCATTTCGCCGTGGTATTGTTTGTCCAGCCATGCTTGCa                     <  1:2176169/69‑1 (MQ=255)
                                                    tACCGTGACGTGCCATCCAGTCCATTTCGCCGTGGTATTGTTTGTCCAGCCATGCTTGCAGTTTGGGCt            >  1:983713/1‑69 (MQ=255)
                                                    tACCGTGACGCGCCATCCAGTCCATTTCGCCGTGGTATTGTTTGTCCAGCCATGCTTGCAGTTTGGGCt            <  1:1794150/69‑1 (MQ=255)
                                                    tACCGTGACGCGCCATCCAGTCCATTTCGCCGTGGTATTGTTTGTCCAGCCATGCTTGCAGTTTGGGCt            <  1:2113178/69‑1 (MQ=255)
                                                    tACCGTGACGCGCCATCCAGTCCATTTCGCCGTGGTATTGTTTGTCCAGCCATGCTTGCAGTTTGGGCt            <  1:2136284/69‑1 (MQ=255)
                                                    tACCGTGACGCGCCATCCAGTCCATTTCGCCGTGGTATTGTTTGTCCAGCCATGCTTGCAGTTTGGGCt            <  1:2905078/69‑1 (MQ=255)
                                                    tACCGTGACGCGCCATCCAGTCCATTTCGCCGTGGTATTGTTTGTCCAGCCATGCTTGCAGTTTGGGCt            <  1:1962012/69‑1 (MQ=255)
                                                    tACCGTGACGCGCCATCCAGTCCATTTCGCCGTGGTATTGTTTGTCCAGCCATGCTTGCAGTTTGGGCt            <  1:1143397/69‑1 (MQ=255)
                                                    tACCGTGACGCGCCATCCAGTCCATTTCGCCGTGGTATTGTTTGTCCAGCCATGCTTGCAGTTTGGGCt            <  1:1958103/69‑1 (MQ=255)
                                                    tACCGTGACGCGCCATCCAGTCCATTTCGCCGTGGTATTGTTTGTCCAGCCATGCTTGCAGTTTGGGCt            <  1:73458/69‑1 (MQ=255)
                                                    tACCGTGACGCGCCATCCAGTCCATTTCGCCGTGGTATTGTTTGTCCAGCCATGCTTGCAGTTTGGGCt            <  1:3010836/69‑1 (MQ=255)
                                                    tACCGTGACGCGCCATCCAGTCCATTTCGCCGTGGTATTGTTTGTCCAGCCATGCTTGCAGTTTGGGCt            <  1:1397191/69‑1 (MQ=255)
                                                    tACCGTGACGCGCCATCCAGTCCATTTCGCCGTGGTATTGTTTGTCCAGCCATGCTTGCAGTTTGGGCt            <  1:50802/69‑1 (MQ=255)
                                                    tACCGTGACGCGCCATCCAGTCCATTTCGCCGTGGTATTGTTTGTCCAGCCATGCTTGCAGTTTGGGCt            <  1:2214738/69‑1 (MQ=255)
                                                    tACCGTGACGCGCCATCCAGTCCATTTCGCCGTGGTATTGTTTGTCCAGCCATGCTTGCAGTTTGGGCt            <  1:3167763/69‑1 (MQ=255)
                                                    tACCGTGACGCGCCATCCAGTCCATTTCGCCGTGGTATTGTTTGTCCAGCCATGCTTGCAGTTTGGGCt            <  1:3241118/69‑1 (MQ=255)
                                                     aCCGTGACGCGCCATCCAGTCCATTTCGCCGTGGTATTGTTTGTCCAGCCATGCTTGCAGTTTGGGCt            <  1:195257/68‑1 (MQ=255)
                                                          gACGTGCCATCCAGTCCATTTCGCCGTGGTATTGTTTGTCCAGCCATGCTTGCAGTTTGGGCTCGGACt      <  1:1331621/69‑1 (MQ=255)
                                                           aCGTGCCATCCAGTCCATTTCGCCGTGGTATTGTTTGTCCAGCCATGCTTGCAGTTTGGGCTCGGACt      <  1:1447746/68‑1 (MQ=255)
                                                              tGCCATCCAGTCCATTTCGCCGTGGTATTGTTTGTCCAGCCATTCTTGCAGTTTGGGCTCGGACTCGCt  <  1:1632606/69‑1 (MQ=255)
                                                              tGCCATCCAGTCCATTTCGCCGTGGTATTGTTTGTCCAGCCATGCTTGCAGTTTGGGCTCGGACTCGCt  <  1:1327879/69‑1 (MQ=255)
                                                              tGCCATCCAGTCCATTTCGCCGTGGTATTGTTTGTCCAGCCATGCTTGCAGTTTGGGCTCGGACTCGCt  <  1:691276/69‑1 (MQ=255)
                                                              tGCCATCCAGTCCATTTCGCCGTGGTATTGTTTGTCCAGCCATGCTTGCAGTTTGGGCTCGGACTCGCt  <  1:1438826/69‑1 (MQ=255)
                                                              tGCCATCCAGTCCATTTCGCCGTGGTATTGTTTGTCCAGCCATGCTTGCAGTTTGGGCTCGGACTCGCt  <  1:1094243/69‑1 (MQ=255)
                                                              tGCCATCCAGTCCATTTCGCCGTGGTATTGTTTGTCCAGCCATGCTTGCAGTTTGGGCTCGGACTCGCt  <  1:2301068/69‑1 (MQ=255)
                                                              |                                                                    
CTGATCACGCGCAGCGTACCGGGCAATAACTCATGAGGGCGAGCGCGCAGCATACCGTGACGTGCCATCCAGTCCATTTCGCCGTGGTATTGTTTGTCCAGCCATGCTTGCAGTTTGGGCTCGGACTCGCT  >  minE/2855500‑2855630

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: