Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,036,664 T→C 38.0% A685A (GCT→GCC rsxC → fused predicted 4Fe‑4S ferredoxin‑type protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,036,6640TC38.0% 39.0 / 30.1 50A685A (GCT→GCCrsxCfused predicted 4Fe‑4S ferredoxin‑type protein
Reads supporting (aligned to +/- strand):  ref base T (17/14);  new base C (7/12);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 2.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

CGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCT  >  minE/1036600‑1036721
                                                                |                                                         
cgcgcAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCg                                                       >  1:304309/1‑69 (MQ=255)
          ggAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGc                                             <  1:2869454/69‑1 (MQ=255)
                     aggCTAATGCGGAGCCAGAAGAACAGATTGATCCGCGCAAAGCCGCCGTCGAAGCGGCTATTGCCcgt                                   >  1:2802934/2‑67 (MQ=11)
                     aggCTAATGCGGAGCCAGAAGAACAGATTGATCCGCGCAAAGCCGCCGTCGAAGCGGCTATTGCCcgt                                   >  1:762483/2‑67 (MQ=11)
                     aggCTAATGCGGAGCCAGAAGAACAGATTGATCCGCGCAAAGCCGCCGTCGAAGCGGCTATTGCCcgt                                   >  1:2943849/2‑67 (MQ=11)
                     aggCTAATGCGGAGCCAGAAGAACAGATTGATCCGCGCAAAGCCGCCGTCGAAGCGGCTATTGCCcgt                                   >  1:2712753/2‑67 (MQ=11)
                     aggCTAATGCGGAGCCAGAAGAACAGATTGATCCGCGCAAAGCCGCCGTCGAAGCGGCTATTGCCcgt                                   >  1:2375354/2‑67 (MQ=11)
                     aggCTAATGCGGAGCCAGAAGAACAGATTGATCCGCGCAAAGCCGCCGTCGAAGCGGCTATTGCCcgt                                   >  1:2998915/2‑67 (MQ=11)
                     aggCTAATGCGGAGCCAGAAGAACAGATTGATCCGCGCAAAGCCGCCGTCGAAGCGGCTATTGCCcgt                                   >  1:2280364/2‑67 (MQ=11)
                      ggCTAATGCGGAGCCAGAAGAACAGATTGATCCGCGCAAAGCCGCTGTCGAAGCGGCTATTGCCcgtga                                 <  1:2286288/69‑4 (MQ=11)
                      ggCTAATGCGGAGCCAGAAGAACAGATTGATCCGCGCAAAGCCGCTGTCGAAGCGGCTATTGCCcgtga                                 <  1:2800153/69‑4 (MQ=11)
                      ggCTAATGCGGAGCCAGAAGAACAGATTGATCCGCGCAAAGCCGCTGTCGAAGCGGCTATTGCCcgtga                                 <  1:2975563/69‑4 (MQ=11)
                      ggCTAATGCGGAGCCAGAAGAACAGATTGATCCGCGCAAAGCCGCTGTCGAAGCGGCTATTGCCcgtga                                 <  1:2421909/69‑4 (MQ=11)
                      ggCTAATGCGGAGCCAGAAGAACAGATTGATCCGCGCAAAGCCGCTGTCGAAGCGGCTATTGCCcgtga                                 <  1:2296654/69‑4 (MQ=11)
                      ggCTAATGCGGAGCCAGAAGAACAGATTGATCCGCGCAAAGCCGCTGTCGAAGCGGCTATTGCCcgtga                                 <  1:3291422/69‑4 (MQ=11)
                      ggCTAATGCGGAGCCAGAAGAACAGATTGATCCGCGCAAAGCCGCTGTCGAAGCGGCTATTGCCcgtga                                 <  1:2202809/69‑4 (MQ=11)
                      ggCTAATGCGGAGCCAGAAGAACAGATTGATCCGCGCAAAGCCGCTGTCGAAGCGGCTATTGCCcgtga                                 <  1:427945/69‑4 (MQ=11)
                      ggCTAATGCGGAGCCAGAAGAACAGATTGATCCGCGCAAAGCCGCTGTCGAAGCGGCTATTGCCcgtga                                 <  1:1227361/69‑4 (MQ=11)
                      ggCTAATGCGGAGCCAGAAGAACAGATTGATCCGCGCAAAGCCGCTGTCGAAGCGGCTATTGCCcgtga                                 <  1:857537/69‑4 (MQ=11)
                      ggCTAATGCGGAGCCAGAAGAACAGATTGATCCGCGCAAAGCCGCTGTCGAAGCGGCTATTGCCcgtga                                 <  1:955797/69‑4 (MQ=11)
                            tGCGGAGCCAGAAGAACAGATTGATCCGCGCAAAGCCGCTGTCGAAGCGGCTATTGCCcgtga                                 <  1:752547/63‑4 (MQ=11)
                              cGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCa                          >  1:2889386/1‑68 (MQ=255)
                              cGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCa                          >  1:3334150/1‑68 (MQ=255)
                              cGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCa                          >  1:748598/1‑68 (MQ=255)
                              cGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGAGGCTATTGCACGCGCCAAAGCa                          >  1:308637/1‑68 (MQ=255)
                               ggAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGc                                             <  1:2009271/48‑1 (MQ=255)
                                gAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGc                                             <  1:2862165/47‑1 (MQ=255)
                                   ccAGAAGAACAGATTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCGCGCAAg                    <  1:604434/69‑1 (MQ=255)
                                   ccAGAAGAACAGATTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCGCGCAAg                    <  1:378281/69‑1 (MQ=255)
                                   ccAGAAGAACAGATTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCGCGCAAg                    <  1:2541943/69‑1 (MQ=255)
                                   ccAGAAGAACAGATTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCGCGCAAg                    <  1:2735994/69‑1 (MQ=255)
                                   ccAGAAGAACAGATTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCGCGCAAg                    <  1:2927672/69‑1 (MQ=255)
                                       aagaaCAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTgg                <  1:3222271/69‑1 (MQ=255)
                                       aagaaCAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTgg                <  1:1840138/69‑1 (MQ=255)
                                        agaaCAGGTTGATCCGCGCAAAGCTGCTGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTgg                >  1:1491741/1‑68 (MQ=255)
                                        agaaCAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGc                                             <  1:2267947/39‑1 (MQ=255)
                                        agaaCAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCAcgcg                                  >  1:3272631/1‑50 (MQ=255)
                                        agaaCAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTgg                >  1:1220818/1‑68 (MQ=255)
                                        agaaCAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTgg                >  1:769768/1‑68 (MQ=255)
                                        agaaCAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTgg                >  1:1364305/1‑68 (MQ=255)
                                        agaaCAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTgg                >  1:693893/1‑68 (MQ=255)
                                        agaaCAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTgg                >  1:1046894/1‑68 (MQ=255)
                                        agaaCAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTgg                >  1:1980714/1‑68 (MQ=255)
                                        agaaCAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTgg                >  1:380630/1‑68 (MQ=255)
                                        agaaCAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTgg                >  1:2311160/1‑68 (MQ=255)
                                        agaaCAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTgg                >  1:231427/1‑68 (MQ=255)
                                        agaaCAGGTTGATACGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTgg                >  1:3113409/1‑68 (MQ=255)
                                               gTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCaa        <  1:3003422/69‑1 (MQ=255)
                                                     ccGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCt  <  1:863412/69‑1 (MQ=255)
                                                             aGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTgg                <  1:1172886/47‑1 (MQ=255)
                                                                |                                                         
CGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCT  >  minE/1036600‑1036721

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: