Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,317,540 T→C 57.2% F36F (TTT→TTC fpr → ferredoxin‑NADP reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,317,5400TC57.2% 10.5 / 28.6 42F36F (TTT→TTCfprferredoxin‑NADP reductase
Reads supporting (aligned to +/- strand):  ref base T (11/7);  new base C (5/19);  total (16/26)
Fisher's exact test for biased strand distribution p-value = 1.14e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

CTGGACCGACGCCCTGTTTAGTCTCACCGTTCACGCCCCCGTGCTTCCGTTTACCGCCGGGCAATTTACCAAGCTTGGCCTTGAAATCGACGGCGAACGCGTCCAGCGCGCCTA  >  minE/2317474‑2317587
                                                                  |                                               
cTGGACCGACGCCCTGTTTAGTCTCACCGTTCACGCCCCCGTGCTTCCGTTTACCGCCGGGCAATTTa                                                >  1:2341710/1‑68 (MQ=255)
cTGGACCGACGCCCTGTTTAGTCTCACCGTTCACGCCCCCGTGCTTCCGTTTACCGCCGGGCAATTTa                                                >  1:1326274/1‑68 (MQ=255)
cTGGACCGACGCCCTGTTTAGTCTCACCGTTCACGCCCCCGTGCTTCCGTTTACCGCCGGGCAATTTa                                                >  1:2296795/1‑68 (MQ=255)
cTGGACCGACGCCCTGTTTAGTCTCACCGTTCACGCCCCCGTGCTTCCGTTTACCGCCGGGCAATTTa                                                >  1:2913653/1‑68 (MQ=255)
                  tAGTCTCACCGTTCACGCCCCCGTGCTTCCGTTTACCGCCGGGCAATTTACCAAGCTTGGCCTTGaa                               <  1:1043442/67‑1 (MQ=255)
                  tAGTCTCACCGTTCACGCCCCCGTGCTTCCGTTTACCGCCGGGCAATTTACCAAGCTTGGCCTTGaa                               <  1:143865/67‑1 (MQ=255)
                  tAGTCTCACCGTTCACGCCCCCGTGCTTCCGTTTACCGCCGGGCAATTTACCAAGCTTGGCCTTGaa                               <  1:2597219/67‑1 (MQ=255)
                  tAGTCTCACCGTTCACGCCCCCGTGCTTCCGTTTACCGCCGGGCAATTTACCAAGCTTGGCCTTGaa                               <  1:2156255/67‑1 (MQ=255)
                  tAGTCTCACCGTTCACGCCCCCGTGCTTCCGTTTACCGCCGGGCAATTTACCAAGCTTGGCCTTGa                                <  1:970639/66‑1 (MQ=255)
                  tAGTCTCACCGTTCACGCCCCCGTGCTTCCGTTTACCACCGGGCAATTTACCAAGCTTGGCCTTGa                                <  1:899101/66‑1 (MQ=255)
                       tcACCGTTCACGCCCCCGTGCTTCCGTTTACCGCCGGGCAATTTACCAAGCTTgg                                      >  1:267459/1‑55 (MQ=255)
                           cGTTCACGCTCCCGTGCTTCCGTTTACCGCCGTGCAATTCACCaaa                                           <  1:525294/46‑2 (MQ=37)
                           cGTTCACGCTCCCGTGCTTCCGTTTACCGCCGGGCAATTCACCaaa                                           <  1:3372703/46‑2 (MQ=38)
                           cGTTCACGCTCCCGTGCTTCCGTTTACCGCCGGGCAATTCACCaaa                                           <  1:3366819/46‑2 (MQ=38)
                           cGTTCACGCTCCCGTGCTTCCGTTTACCGCCGGGCAATTCACCaaa                                           <  1:3211916/46‑2 (MQ=38)
                           cGTTCACGCTCCCGTGCTTCCGTTTACCGCCGGGCAATTCACCaaa                                           <  1:3063995/46‑2 (MQ=38)
                           cGTTCACGCTCCCGTGCTTCCGTTTACCGCCGGGCAATTCACCaaa                                           <  1:466535/46‑2 (MQ=38)
                           cGTTCACGCTCCCGTGCTTCCGTTTACCGCCGGGCAATTCACCaaa                                           <  1:828616/46‑2 (MQ=38)
                           cGTTCACGCTCCCGTGCTTCCGTTTACCGCCGGGCAATTCACCaaa                                           <  1:947502/46‑2 (MQ=38)
                           cGTTCACGCTCCCGTGCTTCCGTTTACCGCCGGGCAATTCACCaaa                                           <  1:2444518/46‑2 (MQ=38)
                           cGTTCACGCTCCCGTGCTTCCGTTTACCGCCGGGCAATTCACCaaa                                           <  1:1218662/46‑2 (MQ=38)
                           cGTTCACGCTCCCGTGCTTCCGTTTACCGCCGGGCAATTCACCaaa                                           <  1:1503073/46‑2 (MQ=38)
                           cGTTCACGCTCCCGTGCTTCCGTTTACCGCCGGGCAATTCACCaaa                                           <  1:1947659/46‑2 (MQ=38)
                           cGTTCACGCTCCCGTGCTTCCGTTTACCGCCGGGCAATTCACCaaa                                           <  1:1939657/46‑2 (MQ=38)
                           cGTTCACGCTCCCGTGCTTCCGTTTACCGCCGGGCAATTCACCaaa                                           <  1:1667374/46‑2 (MQ=38)
                           cGTTCACGCTCCCGTGCTTCCGTTTACAGCCGGGCAATTCACCaaa                                           <  1:2083129/46‑2 (MQ=255)
                            gTTCACGCTCCCGTGCTTCCGTTTACCGCCGGGCAATTCACCaaa                                           <  1:619692/45‑2 (MQ=38)
                            gTTCACGCTCCCGTGCTTCCGTTTACCGCCGGGCAATTCACCaaa                                           <  1:336115/45‑2 (MQ=38)
                            gTTCACGCTCCCGTGCTTCCGTTTACCGCCGGGCAATTCACCaaa                                           <  1:1768504/45‑2 (MQ=38)
                                 cGCTCCCGTGCTTCCGTTTACCGCCGGGCAATTCACCaaa                                           <  1:3160765/40‑2 (MQ=37)
                                     cccGTGCTTCCGTTTACCGCCGGGCAATTTACCAAGCTTGGCCTTGAAATCGACGGCGAACGCGTCCAg          <  1:2407012/69‑1 (MQ=255)
                                          gCTTCCGTTTACCGCCGGGCAATTTACCAAGCTTGGCCTTGAAATCGACGGCGAACGCGTcc            >  1:3187053/1‑62 (MQ=255)
                                          gCTTCCGTTTACCGCCGGGCAATTTACCAAGCTTGGCCTTGAAATCGACGGCGAACGCGTCCAgcgcg      >  1:1044414/1‑68 (MQ=255)
                                          gCTTCCGTTTACCGCCGGGCAATTTACCAAGCTTGGCCTTGAAATCGACGGCGAACGCGTCCAgcgcg      >  1:1495831/1‑68 (MQ=255)
                                          gCTTCCGTTTACCGCCGGGCAATTTACCAAGCTTGGCCTTGAAATCGACGGCGAACGCGTCCAgcgcg      >  1:1560991/1‑68 (MQ=255)
                                          gCTTCCGTTTACCGCCGGGCAATTTACCAAGCTTGGCCTTGAAATCGACGGCGAACGCGTCCAgcgcg      >  1:2843219/1‑68 (MQ=255)
                                          gCTTCCGTTTACCGCCGGGCAATTTACCAAGCTTGGCCTTGAAATCGACGGCGAACGCGTCCAgcgcg      >  1:2171867/1‑68 (MQ=255)
                                             tCCGTTTACCGCCGGGCAATTCACCAAACTTGGCCTTGAAATTGATGGCGACCGCGTTCAGCGCGCCt   >  1:2343843/1‑68 (MQ=255)
                                             tCCGTTTACCGCCGGGCAATTCACCAAACTTGGCCTTGAAATTGATGGCGAACGCGTTCAGCGCGCCTa  >  1:1691265/1‑69 (MQ=255)
                                             tCCGTTTACCGCCGGGCAATTCACCAAACTTGGCCTTGAAATTGATGGCGAACGCGTTCAGCGCGCCTa  >  1:3122918/1‑69 (MQ=255)
                                             tCCGTTTACCGCCGGGCAATTCACCAAACTTGGCCTTGAAATTGATGGCGAACGCGTTCAGCGCGCCTa  >  1:78127/1‑69 (MQ=255)
                                             tCCGTTTACCGCCGGGCAATTCACCAAACTTGGCCTTGAAATTGATGGCGAACGCGTTCAGCGCGCCTa  >  1:2210337/1‑69 (MQ=255)
                                                                  |                                               
CTGGACCGACGCCCTGTTTAGTCTCACCGTTCACGCCCCCGTGCTTCCGTTTACCGCCGGGCAATTTACCAAGCTTGGCCTTGAAATCGACGGCGAACGCGTCCAGCGCGCCTA  >  minE/2317474‑2317587

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: