Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,732,510 C→T 51.7% C54C (TGC→TGT yfhL → predicted 4Fe‑4S cluster‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,732,5100CT51.7% 1.3 / 106.0 116C54C (TGC→TGTyfhLpredicted 4Fe‑4S cluster‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (11/45);  new base T (22/38);  total (33/83)
Fisher's exact test for biased strand distribution p-value = 6.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 121 total aligned reads displayed.

GCGATAAGTGTACCGAATGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGCCCGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATAAATT  >  minE/1732452‑1732578
                                                          |                                                                    
gagATAAGTGTACCGAATGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGCCCGATcccc                                                             >  1:1424779/3‑68 (MQ=255)
gCGATAAGTGTACCGAATGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGCCCGATcccc                                                             >  1:2432647/1‑68 (MQ=255)
gCGATAAGTGTACCGAATGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGCCCGATcccc                                                             >  1:457114/1‑68 (MQ=255)
gCGATAAGTGTACCGAATGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGCCCGATcccc                                                             >  1:1042531/1‑68 (MQ=255)
 cGATAAGTGTACCGAATGCGTAGGGTACTACGAGACACCAACCTGCCAGAAGGTGTGCCCGATCCCCaa                                                           <  1:1331282/69‑1 (MQ=255)
 cGATAAGTGTACCGAATGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGCCCGATCCCCaa                                                           <  1:3097363/69‑1 (MQ=255)
 cGATAAGTGTACCGAATGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGCCCGATCCCCaa                                                           <  1:3069406/69‑1 (MQ=255)
 cGATAAGTGTACCGAATGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGCCCGATCCCCaa                                                           <  1:2836444/69‑1 (MQ=255)
 cGATAAGTGTACCGAATGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGCCCGATCCCCaa                                                           <  1:2724722/69‑1 (MQ=255)
 cGATAAGTGTACCGAATGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGCCCGATCCCCaa                                                           <  1:1014846/69‑1 (MQ=255)
 cGATAAGTGTACCGAATGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGCCCGATCCCCaa                                                           <  1:2320971/69‑1 (MQ=255)
 cGATAAGTGTACCGAATGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGCCCGATCCCCaa                                                           <  1:2286505/69‑1 (MQ=255)
 cGATAAGTGTACCGAATGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGCCCGATCCCCaa                                                           <  1:2229719/69‑1 (MQ=255)
 cGATAAGTGTACCGAATGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGCCCGATCCCCaa                                                           <  1:2180247/69‑1 (MQ=255)
 cGATAAGTGTACCGAATGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGCCCGATCCCCaa                                                           <  1:2165420/69‑1 (MQ=255)
 cGATAAGTGTACCGAATGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGCCCGATCCCCaa                                                           <  1:1682945/69‑1 (MQ=255)
 cGATAAGTGTACCGAATGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGCCCGATCCCCaa                                                           <  1:3404612/69‑1 (MQ=255)
 cGATAAGTGTACCGAATGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGCCCGATCCCCaa                                                           <  1:1042786/69‑1 (MQ=255)
 cGATAAGTGTACCGAATGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGCCCGATCCCCaa                                                           <  1:140544/69‑1 (MQ=255)
 cGATAAGTGTACCGAATGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGCCCGATCCCCaa                                                           <  1:787710/69‑1 (MQ=255)
  gATAAGTGTACCGAATGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGCCCGATCCCCaa                                                           <  1:780989/68‑1 (MQ=255)
      aGTGTACCGAATGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGCCCGATCCCCAATACTa                                                      <  1:619965/69‑1 (MQ=255)
      aGTGTACCGAATGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGCCCGATCCCCAATACTa                                                      <  1:3235099/69‑1 (MQ=255)
      aGTGTACCGAATGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGCCCGATCCCCAATACTa                                                      <  1:2985807/69‑1 (MQ=255)
      aGTGTACCGAATGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGCCCGATCCCCAATACTa                                                      <  1:1997262/69‑1 (MQ=255)
      aGTGTACCGAATGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGCCCGATCCCCAATACTa                                                      <  1:706779/69‑1 (MQ=255)
      aGTGTACCGAATGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGCCCGATCCCCAATACTa                                                      <  1:753005/69‑1 (MQ=255)
            ccGAATGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGCCCGATCCCCAATACTATtgt                                                  >  1:2557906/1‑67 (MQ=255)
            ccGAATGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGCCCGATCCCCAATACTATtgt                                                  >  1:2115005/1‑67 (MQ=255)
                 tGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGa                                            >  1:1781979/1‑68 (MQ=255)
                 tGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGa                                            >  1:34738/1‑68 (MQ=255)
                 tGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGAt                                           >  1:3052527/1‑69 (MQ=255)
                 tGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGAt                                           >  1:2669218/1‑69 (MQ=255)
                 tGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGAt                                           >  1:2609018/1‑69 (MQ=255)
                 tGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGAt                                           >  1:1149410/1‑69 (MQ=255)
                 tGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGAt                                           >  1:1113076/1‑69 (MQ=255)
                 tGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGAt                                           >  1:2840716/1‑69 (MQ=255)
                 tGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGAt                                           >  1:381207/1‑69 (MQ=255)
                 tGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGAt                                           >  1:285441/1‑69 (MQ=255)
                 tGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGAt                                           >  1:2903667/1‑69 (MQ=255)
                 tGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGAt                                           >  1:1573211/1‑69 (MQ=255)
                 tGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGAt                                           >  1:1236066/1‑69 (MQ=255)
                 tGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGAt                                           >  1:1499753/1‑69 (MQ=255)
                 tGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGAt                                           >  1:2121251/1‑69 (MQ=255)
                 tGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGAt                                           >  1:1990594/1‑69 (MQ=255)
                 tGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACGATTGTGAAAGa                                            >  1:2724172/1‑68 (MQ=255)
                 tGCGTAGGGCACTACGAGACACAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATc                                          >  1:2128009/1‑69 (MQ=255)
                     tAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGCCCGATCCCCAATACTATTGTGAAAGATCCgg                                       <  1:876990/69‑1 (MQ=255)
                     tAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGCCCGATCCCCAATACTATTGTGAAAGATCCgg                                       <  1:6412/69‑1 (MQ=255)
                      aGGGCACTACGAGACACCAACCTGCCAGAAGGGGTGCCCGATCCCCAATACTATTGTGAAAGATCCgg                                       <  1:1066995/68‑1 (MQ=255)
                                 agaCACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCga                             <  1:2525163/67‑1 (MQ=255)
                                 agaCACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCga                             <  1:3170353/67‑1 (MQ=255)
                                 agaCACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCga                             <  1:3390002/67‑1 (MQ=255)
                                 agaCACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCga                             <  1:3338156/67‑1 (MQ=255)
                                 agaCACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCga                             <  1:272366/67‑1 (MQ=255)
                                 agaCACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCga                             <  1:2693052/67‑1 (MQ=255)
                                 agaCACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCga                             <  1:1473501/67‑1 (MQ=255)
                                 agaCACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCga                             <  1:62986/67‑1 (MQ=255)
                                 agaCACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCga                             <  1:1827618/67‑1 (MQ=255)
                                 agaCACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCga                             <  1:646311/67‑1 (MQ=255)
                                 agaCACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCga                             <  1:1767545/67‑1 (MQ=255)
                                 agaCACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCga                             <  1:2053522/67‑1 (MQ=255)
                                 agaCACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCga                             <  1:2090569/67‑1 (MQ=255)
                                 agaCACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCga                             <  1:218038/67‑1 (MQ=255)
                                 agaCACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCga                             <  1:2186609/67‑1 (MQ=255)
                                 agaCACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCga                             <  1:266273/67‑1 (MQ=255)
                                 agaCACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCga                             <  1:227457/67‑1 (MQ=255)
                                 agaCACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCga                             <  1:814290/67‑1 (MQ=255)
                                 agaCACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCga                             <  1:1411510/67‑1 (MQ=255)
                                 agaCACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCga                             <  1:92277/67‑1 (MQ=255)
                                 agaCACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCga                             <  1:2441018/67‑1 (MQ=255)
                                 agaCACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCga                             <  1:2456072/67‑1 (MQ=255)
                                 agaCACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCga                             <  1:1805933/67‑1 (MQ=255)
                                 agaCACCAACCTGCCAGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCga                             <  1:1152455/67‑1 (MQ=255)
                                    cacCAACCTGCCAGAAGGTGTGCCCGATCCCCAATACTATTGTGAAAGATCCGGCGCa                                   <  1:762690/58‑1 (MQ=255)
                                    cacCAACCTGCCAGAAGGTGTGCCCGATCCCCAATACTATTGTGAAAGATCCGGCGCa                                   <  1:970849/58‑1 (MQ=255)
                                    cacCAACCTGCCAGAAGGTGTGCCCGATCCCCAATACTATTGTGAAAGATCCGGCGCa                                   <  1:417109/58‑1 (MQ=255)
                                    cacCAACCTGCCAGAAGGTGTGCCCGATCCCCAATACTATTGTGAAAGATCCGGCGCa                                   <  1:1358539/58‑1 (MQ=255)
                                    cacCAACCTGCCAGAAGGTGTGCCCGATCCCCAATACTATTGTGAAAGATCCGGCGCa                                   <  1:3180204/58‑1 (MQ=255)
                                    cacCAACCTGCCAGAAGGTGTGCCCGATCCCCAATACTATTGTGAAAGATCCGGCGCa                                   <  1:2844884/58‑1 (MQ=255)
                                        aaCCTGCCAGAAGGTGTGCCCGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACagaag                     >  1:2100144/1‑68 (MQ=255)
                                         aCCTGCCAGAAGGTGTGCCCGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACAgaagaa                   <  1:3338989/69‑1 (MQ=255)
                                         aCCTGCCAGAAGGTGTGCCCGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACAgaagaa                   <  1:1701575/69‑1 (MQ=255)
                                                aGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCGAGACGGAAGAACAGtt              <  1:2611525/67‑1 (MQ=255)
                                                aGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCGAGACGGAAGAACAGtt              <  1:1014894/67‑1 (MQ=255)
                                                aGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCGAGACGGAAGAACAGtt              <  1:115549/67‑1 (MQ=255)
                                                aGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCGAGACGGAAGAACAGtt              <  1:319714/67‑1 (MQ=255)
                                                aGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCGAGACGGAAGAACAGtt              <  1:1828428/67‑1 (MQ=255)
                                                aGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCGAGACGGAAGAACAGtt              <  1:1952982/67‑1 (MQ=255)
                                                aGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCGAGACGGAAGAACAGtt              <  1:3101962/67‑1 (MQ=255)
                                                aGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCGAGACGGAAGAACAGtt              <  1:2393860/67‑1 (MQ=255)
                                                aGAAGGTGTGTCCGATCCCCAATACTATTGTGAAAGATCCGCAACATGTCGAGACGGAAGAACAGtt              <  1:2629372/67‑1 (MQ=255)
                                                  aaGGTGTGCCCGATCCCCAATACTATTGTGAAAGATCCgg                                       <  1:168728/40‑1 (MQ=255)
                                                      tgtgCCCGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATa      <  1:2264571/69‑1 (MQ=255)
                                                      tgtgCCCGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATa      <  1:3154449/69‑1 (MQ=255)
                                                      tgtgCCCGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATa      <  1:122750/69‑1 (MQ=255)
                                                       gtgCCCGATCCCCAATACTATTGTTAAAGATCCGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATa      <  1:896779/68‑1 (MQ=255)
                                                          cccGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATAAAtt  >  1:1296139/1‑69 (MQ=255)
                                                          cccGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATAAAtt  >  1:1190798/1‑69 (MQ=255)
                                                          cccGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATAAAtt  >  1:502907/1‑69 (MQ=255)
                                                          |                                                                    
GCGATAAGTGTACCGAATGCGTAGGGCACTACGAGACACCAACCTGCCAGAAGGTGTGCCCGATCCCCAATACTATTGTGAAAGATCCGGCGCATGTCGAGACAGAAGAACAGTTGTGGGATAAATT  >  minE/1732452‑1732578

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: