Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,317,828 T→A 36.3% D132E (GAT→GAA fpr → ferredoxin‑NADP reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,317,8280TA36.3% 53.9 / 52.1 66D132E (GAT→GAAfprferredoxin‑NADP reductase
Reads supporting (aligned to +/- strand):  ref base T (30/12);  new base A (16/8);  total (46/20)
Fisher's exact test for biased strand distribution p-value = 7.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGCTGGCAACCGGTACAGCGATTGGCCCTTATTTATCGATTCTGCAACTAGGTAAAGATTTAGATCGCTTCAAAAATCTGGTCCTGGTGCACGCCGCACGTTATGCCGCCGACTTAAGCTAT  >  minE/2317762‑2317885
                                                                  |                                                         
gATGCTGGCAACCGGTACAGCGACTGGCCCTTATTTATCGATTCTGCAACTAGGTAAAGATTTAGATCg                                                         >  1:73442/1‑69 (MQ=255)
         aaCCGGTACAGCGATTGGCCCTTATTTATCGATTCTGCAACTAGGTAAAGATTTAGATCGCTTCaaaa                                                 <  1:1140596/68‑1 (MQ=255)
          aCCGGTACAGCGATTGGCCCTTATTTATCGATTCTGCAACTAGGTAAAGATTTAGATCGCTTCaaaa                                                 <  1:989171/67‑1 (MQ=255)
           ccGGTACAGCGATTGGCCCTTATTTATCGATTCTGCAATTGGGCAAAGATTTAGAACGCTTCaaaaaca                                              >  1:380373/1‑67 (MQ=255)
           ccGGTACAGCGATTGGCCCTTATTTATCGATTCTGCAATTGGGCAAAGATTTAGAACGCTTCaaaaaca                                              >  1:224569/1‑67 (MQ=255)
           ccGGTACAGCGATTGGCCCTTATTTATCGATTCTGCAATTGGGCAAAGATTTAGAACGCTTCaaaaaca                                              >  1:2300130/1‑67 (MQ=255)
           ccGGTACAGCGATTGGCCCTTATTTATCGATTCTGCAATTGGGCAAAGATTTAGAACGCTTCaaaaaca                                              >  1:2735227/1‑67 (MQ=255)
           ccGGTACAGCGATTGGCCCTTATTTATCGATTCTGCAATTGGGCAAAGATTTAGAACGCTTCaaaaaca                                              >  1:2773687/1‑67 (MQ=255)
           ccGGTACAGCGATTGGCCCTTATTTATCGATTCTGCAATTGGGCAAAGATTTAGAACGCTTCaaaaaca                                              >  1:2800931/1‑67 (MQ=255)
           ccGGTACAGCGATTGGCCCTTATTTATCGATTCTGCAATTGGGCAAAGATTTAGAACGCTTCaaaaaca                                              >  1:3111335/1‑67 (MQ=255)
           ccGGTACAGCGATTGGCCCTTATTTATCGATTCTGCAATTGGGCAAAGATTTAGAACGCTTCaaaaaca                                              >  1:3222868/1‑67 (MQ=255)
           ccGGTACAGCGATTGGCCCTTATTTATCGATTCTGCAATTGGGCAAAGATTTAGAACGCTTCaaaaaca                                              >  1:1775845/1‑67 (MQ=255)
           ccGGTACAGCGATTGGCCCTTATTTATCGATTCTGCAATTGGGCAAAGATTTAGAACGCTTCaaaaaca                                              >  1:1749680/1‑67 (MQ=255)
           ccGGTACAGCGATTGGCCCTTATTTATCGATTCTGCAATTGGGCAAAGATTTAGAACGCTTCaaaaaca                                              >  1:1597167/1‑67 (MQ=255)
           ccGGTACAGCGATTGGCCCTTATTTATCGATTCTGCAATTGGGCAAAGATTTAGAACGCTTCaaaaaca                                              >  1:3243082/1‑67 (MQ=255)
           ccGGTACAGCGATTGGCCCTTATTTATCGATTCTGCAATTGGGCAAAGATTTAGAACGCTTCaaaaaca                                              >  1:1375030/1‑67 (MQ=255)
           ccGGTACAGCGATTGGCCCTTATTTATCGATTCTGCAATTGGGCAAAGATTTAGAACGCTTCaaaaaca                                              >  1:671474/1‑67 (MQ=255)
           ccGGTACAGCGATTGGCCCTTATTTATCGATTCTGCAATTGGGCAAAGATTTAGAACGCTTCaaaaaca                                              >  1:733390/1‑67 (MQ=255)
           ccGGTACAGCGATTGGCCCTTATTTATCGATTCTGCAATTGGGCAAAGATTTAGAACGCTTCaaaaaca                                              >  1:994642/1‑67 (MQ=255)
             ggTACAGCGATTGGCCCTTATTTATCGATTCTGCAACTAGGTAAAGATTTAGATCGCTTCAAAAATCTg                                            >  1:27183/1‑69 (MQ=255)
             ggTACAGCGATTGGCCCTTATTTATCGATTCTGCAACTAGGTAAAGATTTAGATCGCTTCAAAAATCTg                                            >  1:427503/1‑69 (MQ=255)
             ggTACAGCGATTGGCCCTTATTTATCGATTCTGCAACTAGGTAAAGATTTAGATCGCTTCAAAAATCTg                                            >  1:592615/1‑69 (MQ=255)
             ggTACAGCGATTGGCCCTTATTTATCGATTCTGCAACTAGGTAAAGATTTAGATCGCTTCAAAAATCTg                                            >  1:966535/1‑69 (MQ=255)
             ggTACAGCGATTGGCCCTTATTTATCGATTCTGCAACTAGGTAAAGATTTAGATCGCTTCAAAAATCTg                                            >  1:1372444/1‑69 (MQ=255)
             ggTACAGCGATTGGCCCTTATTTATCGATTCTGCAACTAGGAAAAGATTTAGATAGCTTCAAAAATCTg                                            >  1:2395396/1‑69 (MQ=255)
                   gCGATTGGCCCTTATTTATCTATTCTGCAACTAGGTAAAGATTTAGATCGCTTCAAAAATCTGGTCCTg                                      >  1:1209520/1‑69 (MQ=255)
                   gCGATTGGCCCTTATTTATCGATTCTGCAACTAGGTAAAGATTTAGATCGCTTCAAAAATCTGGTCCTg                                      >  1:2948257/1‑69 (MQ=255)
                   gCGATTGGCCCTTATTTATCGATTCTGCAACTAGGTAAAGATTTAGATCGCTTCAAAAATCTGGTCCTg                                      >  1:531863/1‑69 (MQ=255)
                   gCGATTGGCCCTTATTTATCGATTCTGCAACTAGGTAAAGATTTAGATCGCTTCAAAAATCTGGTCCTg                                      >  1:1005162/1‑69 (MQ=255)
                    cGATTGGCCCTTATTTATCGATTCTGCAACTAGGTAAAGATTTAGATCGCTTCAAAAATCTGGTCCTgg                                     >  1:3077145/1‑69 (MQ=255)
                    cGATTGGCCCTTATTTATCGATTCTGCAACTAGGTAAAGATTTAGATCGCTTCAAAAATCTGGTCCTgg                                     >  1:1992299/1‑69 (MQ=255)
                    cGATTGGCCCTTATTTATCGATTCTGCAACTAGGTAAAGATTTAGATCGCTTCAAAAATCTGGTCCTgg                                     >  1:3093184/1‑69 (MQ=255)
                                    atCGATTCTGCAACTAGGTAAAGATTTAGATCGCTTCAAAAATCTGGTCCTGGTGCACGCCGCACGTTa                     >  1:2914300/1‑69 (MQ=255)
                                    atCGATTCTGCAACTAGGTAAAGATTTAGATCGCTTCAAAAATCTGGTCCTGGTGCACGCCGCACGTTa                     >  1:563574/1‑69 (MQ=255)
                                    atCGATTCTGCAACTAGGTAAAGATTTAGATCGCTTCAAAAATCTGGTCCTGGTGCACGCCGCACGTTa                     >  1:12587/1‑69 (MQ=255)
                                    atCGATTCTGCAACTAGGTAAAGATTTAGATCGCTTCAAAAATCTGGTCCTGGTGCACGCCGCACGTTa                     >  1:2342095/1‑69 (MQ=255)
                                    atCGATTCTGCAACTAGGTAAAGATTTAGATCGCTTCAAAAATCTGGTCCTGGTGCACGCCGCACGTTa                     >  1:760025/1‑69 (MQ=255)
                                    atCGATTCTGCAACTAGGTAAAGATTTAGATCGCTTCAAAAATCTGGTCCTGGTGCACGCCGCACGTTa                     >  1:1572085/1‑69 (MQ=255)
                                    atCGATTCTGCAACTAGGTAAAGATTTAGATCGCTTCAAAAATCTGGTCCTGGTGCACGCCGCACGTTa                     >  1:3277940/1‑69 (MQ=255)
                                       gATTCTGCAACTAGGTAAAGATTTAGATCGCTTCAAAAATCTGGTCCTGGTGCACGCCGCACGTTAt                    <  1:1580328/67‑1 (MQ=255)
                                       gATTCTGCAACTAGGTAAAGATTTAGATCGCTTCAAAAATCTGGGCCTGGTGCACGCCGCACGTTAt                    <  1:1074884/67‑1 (MQ=255)
                                                   aGGTAAAGATTTAGATCGCTTCAAAAATCTGGTCCTGGTGCACGCCGCACGTTATGCCGCCGACTTaa       >  1:1838501/1‑68 (MQ=255)
                                                       aaaGATTTAGATCGCTTCAAAAATCTGGTCCTGGTGCACGCCGCACGTTATGCCGCCGACTTAAGCTAt  >  1:3215020/1‑69 (MQ=255)
                                                        aaGATTTAGATCGCTTCAAAAATCTGGTCCTGGTGCACGCCGCACGt                       >  1:1493684/1‑47 (MQ=255)
                                                        aaGATTTAGATCGCTTCAAAAATCTGGTCCTGGTGCACGCCGCACGt                       >  1:1314148/1‑47 (MQ=255)
                                                        aaGATTTAGATCGCTTCAAAAATCTGGTCCTGGTGCACGCCGCACGt                       <  1:93731/47‑1 (MQ=255)
                                                        aaGATTTAGATCGCTTCAAAAATCTGGTCCTGGTGCACGCCGCACGt                       <  1:1200940/47‑1 (MQ=255)
                                                        aaGATTTAGATCGCTTCAAAAATCTGGTCCTGGTGCACGCCGCACGt                       <  1:1875895/47‑1 (MQ=255)
                                                        aaGATTTAGATCGCTTCAAAAATCTGGTCCTGGTGCACGCCGCACGt                       >  1:117061/1‑47 (MQ=255)
                                                        aaGATTTAGATCGCTTCAAAAATCTGGTCCTGGTGCACGCCGCACGt                       >  1:1912651/1‑47 (MQ=255)
                                                        aaGATTTAGATCGCTTCAAAAATCTGGTCCTGGTGCACGCCGCACGt                       >  1:2742182/1‑47 (MQ=255)
                                                        aaGATTTAGATCGCTTCAAAAATCTGGTCCTGGTGCACGCCGCACGt                       <  1:269514/47‑1 (MQ=255)
                                                        aaGATTTAGATCGCTTCAAAAATCTGGTCCTGGTGCACGCCGCACGt                       <  1:2606258/47‑1 (MQ=255)
                                                        aaGATTTAGATCGCTTCAAAAATCTGGTCCTGGTGCACGCCGCACGt                       >  1:2479190/1‑47 (MQ=255)
                                                        aaGATTTAGATCGCTTCAAAAATCTGGTCCTGGTGCACGCCGCACGt                       <  1:2370228/47‑1 (MQ=255)
                                                        aaGATTTAGATCGCTTCAAAAATCTGGTCCTGGTGCACGCCGCACGt                       <  1:2294237/47‑1 (MQ=255)
                                                        aaGATTTAGATCGCTTCAAAAATCTGGTCCTGGTGCACGCCGCACGt                       >  1:2122417/1‑47 (MQ=255)
                                                         aGATTTAGATCGCTTCAAAAATCTGGTCCTGGTGCACGCCGCACGt                       <  1:745568/46‑1 (MQ=255)
                                                         aGATTTAGAACGCTTCAAAAACATGGTGCTGGTTCACGCCGCACGATATGCCGCCGACTTaa       <  1:3292693/62‑1 (MQ=255)
                                                         aGATTTAGAACGCTTCAAAAACATGGTGCTGGTTCACGCCGCACGATATGCCGCCGACTTaa       <  1:3292918/62‑1 (MQ=255)
                                                         aGATTTAGAACGCTTCAAAAACATGGTGCTGGTTCACGCCGCACGATATGCCGCCGACTTaa       <  1:3411581/62‑1 (MQ=255)
                                                         aGATTTAGAACGCTTCAAAAACATGGTGCTGGTTCACGCCGCACGATATGCCGCCGACTTaa       <  1:192325/62‑1 (MQ=255)
                                                         aGATTTAGAACGCTTCAAAAACATGGTGCTGGTTCACGCCGCACGATATGCCGCCGACTTaa       <  1:2009800/62‑1 (MQ=255)
                                                         aGAATTAGAACGCTTCAAAAACATGGTGCTGGTTCACCCCGCACGATATGCCGCCGACTTaa       <  1:2086308/62‑1 (MQ=255)
                                                          gATTTAGAACGCTTCAAAAACATGGTGCTGGTTCACGCCGCACGATATGCCGCCGACTTaa       <  1:3162947/61‑1 (MQ=255)
                                                           aTTTAGAACGCTTCAAAAACATGGTGCTGGTTCACGCCGCACGATATGCCGCCGACTTaa       <  1:2780711/60‑1 (MQ=255)
                                                                  |                                                         
GATGCTGGCAACCGGTACAGCGATTGGCCCTTATTTATCGATTCTGCAACTAGGTAAAGATTTAGATCGCTTCAAAAATCTGGTCCTGGTGCACGCCGCACGTTATGCCGCCGACTTAAGCTAT  >  minE/2317762‑2317885

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: