Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,036,647 G→A 35.7% V680I (GTT→ATT)  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,6470GA35.7% 60.5 / 37.1 67V680I (GTT→ATT) rsxCfused predicted 4Fe‑4S ferredoxin‑type protein
Reads supporting (aligned to +/- strand):  ref base G (27/16);  new base A (7/17);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 1.12e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

CTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAA  >  minE/1036582‑1036715
                                                                 |                                                                    
cTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTg                                                                   <  1:2441299/69‑1 (MQ=255)
     gCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATcc                                                               >  1:2800490/1‑68 (MQ=255)
     gCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCcg                                                              >  1:1135643/1‑69 (MQ=255)
     gCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCcg                                                              >  1:2550010/1‑69 (MQ=255)
     gCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCcg                                                              >  1:2207789/1‑69 (MQ=255)
     gCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCcg                                                              >  1:343848/1‑69 (MQ=255)
     gCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCcg                                                              >  1:2133436/1‑69 (MQ=255)
     gCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCcg                                                              >  1:343971/1‑69 (MQ=255)
     gCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCcg                                                              >  1:47651/1‑69 (MQ=255)
     gCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCcg                                                              >  1:1571583/1‑69 (MQ=255)
     gCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCcg                                                              >  1:153770/1‑69 (MQ=255)
          aGCCAAAGCGCGCAATCGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCaa                                                         <  1:1557112/69‑1 (MQ=255)
          aGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCaa                                                         <  1:1640959/69‑1 (MQ=255)
          aGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCaa                                                         <  1:853290/69‑1 (MQ=255)
          aGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCaa                                                         <  1:788752/69‑1 (MQ=255)
          aGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCaa                                                         <  1:216916/69‑1 (MQ=255)
          aGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCaa                                                         <  1:2573398/69‑1 (MQ=255)
          aGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCaa                                                         <  1:2109993/69‑1 (MQ=255)
          aGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCaa                                                         <  1:1978570/69‑1 (MQ=255)
          aGCCAAAGCGCGCAAACGGGAACAGCAACCGGATAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCaa                                                         <  1:1414262/69‑1 (MQ=255)
                  cgcgcAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCg                                                 >  1:304309/1‑69 (MQ=255)
                            ggAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGc                                       <  1:2869454/69‑1 (MQ=255)
                                       aggCTAATGCGGAGCCAGAAGAACAGATTGATCCGCGCAAAGCCGCCGTCGAAGCGGCTATTGCCcgt                             >  1:2712753/2‑67 (MQ=11)
                                       aggCTAATGCGGAGCCAGAAGAACAGATTGATCCGCGCAAAGCCGCCGTCGAAGCGGCTATTGCCcgt                             >  1:2802934/2‑67 (MQ=11)
                                       aggCTAATGCGGAGCCAGAAGAACAGATTGATCCGCGCAAAGCCGCCGTCGAAGCGGCTATTGCCcgt                             >  1:2280364/2‑67 (MQ=11)
                                       aggCTAATGCGGAGCCAGAAGAACAGATTGATCCGCGCAAAGCCGCCGTCGAAGCGGCTATTGCCcgt                             >  1:2375354/2‑67 (MQ=11)
                                       aggCTAATGCGGAGCCAGAAGAACAGATTGATCCGCGCAAAGCCGCCGTCGAAGCGGCTATTGCCcgt                             >  1:2943849/2‑67 (MQ=11)
                                       aggCTAATGCGGAGCCAGAAGAACAGATTGATCCGCGCAAAGCCGCCGTCGAAGCGGCTATTGCCcgt                             >  1:2998915/2‑67 (MQ=11)
                                       aggCTAATGCGGAGCCAGAAGAACAGATTGATCCGCGCAAAGCCGCCGTCGAAGCGGCTATTGCCcgt                             >  1:762483/2‑67 (MQ=11)
                                        ggCTAATGCGGAGCCAGAAGAACAGATTGATCCGCGCAAAGCCGCTGTCGAAGCGGCTATTGCCcgtga                           <  1:2975563/69‑4 (MQ=11)
                                        ggCTAATGCGGAGCCAGAAGAACAGATTGATCCGCGCAAAGCCGCTGTCGAAGCGGCTATTGCCcgtga                           <  1:3291422/69‑4 (MQ=11)
                                        ggCTAATGCGGAGCCAGAAGAACAGATTGATCCGCGCAAAGCCGCTGTCGAAGCGGCTATTGCCcgtga                           <  1:857537/69‑4 (MQ=11)
                                        ggCTAATGCGGAGCCAGAAGAACAGATTGATCCGCGCAAAGCCGCTGTCGAAGCGGCTATTGCCcgtga                           <  1:427945/69‑4 (MQ=11)
                                        ggCTAATGCGGAGCCAGAAGAACAGATTGATCCGCGCAAAGCCGCTGTCGAAGCGGCTATTGCCcgtga                           <  1:955797/69‑4 (MQ=11)
                                        ggCTAATGCGGAGCCAGAAGAACAGATTGATCCGCGCAAAGCCGCTGTCGAAGCGGCTATTGCCcgtga                           <  1:2800153/69‑4 (MQ=11)
                                        ggCTAATGCGGAGCCAGAAGAACAGATTGATCCGCGCAAAGCCGCTGTCGAAGCGGCTATTGCCcgtga                           <  1:1227361/69‑4 (MQ=11)
                                        ggCTAATGCGGAGCCAGAAGAACAGATTGATCCGCGCAAAGCCGCTGTCGAAGCGGCTATTGCCcgtga                           <  1:2421909/69‑4 (MQ=11)
                                        ggCTAATGCGGAGCCAGAAGAACAGATTGATCCGCGCAAAGCCGCTGTCGAAGCGGCTATTGCCcgtga                           <  1:2296654/69‑4 (MQ=11)
                                        ggCTAATGCGGAGCCAGAAGAACAGATTGATCCGCGCAAAGCCGCTGTCGAAGCGGCTATTGCCcgtga                           <  1:2286288/69‑4 (MQ=11)
                                        ggCTAATGCGGAGCCAGAAGAACAGATTGATCCGCGCAAAGCCGCTGTCGAAGCGGCTATTGCCcgtga                           <  1:2202809/69‑4 (MQ=11)
                                              tGCGGAGCCAGAAGAACAGATTGATCCGCGCAAAGCCGCTGTCGAAGCGGCTATTGCCcgtga                           <  1:752547/63‑4 (MQ=11)
                                                cGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCa                    >  1:2889386/1‑68 (MQ=255)
                                                cGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCa                    >  1:748598/1‑68 (MQ=255)
                                                cGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCa                    >  1:3334150/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:2735994/69‑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: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:1046894/1‑68 (MQ=255)
                                                          agaaCAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTgg          >  1:1220818/1‑68 (MQ=255)
                                                          agaaCAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTgg          >  1:1364305/1‑68 (MQ=255)
                                                          agaaCAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTgg          >  1:769768/1‑68 (MQ=255)
                                                          agaaCAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTgg          >  1:231427/1‑68 (MQ=255)
                                                          agaaCAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTgg          >  1:693893/1‑68 (MQ=255)
                                                          agaaCAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTgg          >  1:2311160/1‑68 (MQ=255)
                                                          agaaCAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTgg          >  1:1980714/1‑68 (MQ=255)
                                                          agaaCAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTgg          >  1:380630/1‑68 (MQ=255)
                                                          agaaCAGGTTGATACGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTgg          >  1:3113409/1‑68 (MQ=255)
                                                                 gTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCaa  <  1:3003422/69‑1 (MQ=255)
                                                                 |                                                                    
CTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTGCACGCGCCAAAGCACGCAAGCTGGAACAGCAA  >  minE/1036582‑1036715

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: