Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,036,505 G→A 100% A632A (GCG→GCA 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,5050GA98.2% 121.3 / ‑5.5 48A632A (GCG→GCArsxCfused predicted 4Fe‑4S ferredoxin‑type protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (29/18);  minor base C (0/1);  total (29/19)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCG  >  minE/1036438‑1036572
                                                                   |                                                                   
cGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcacg                                                                   >  1:2997194/1‑70 (MQ=25)
cGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcacg                                                                   <  1:3306334/70‑1 (MQ=25)
     ccAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCt                                                             >  1:626171/1‑71 (MQ=21)
      cAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTg                                                            >  1:1691315/1‑71 (MQ=21)
         aagaaCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcacgc                                                                  >  1:1957299/1‑62 (MQ=25)
         aagaaCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCaa                                                                >  1:1563277/1‑64 (MQ=25)
          agaaCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCAGTGCCAAAGCACGCAAGCTg                                                            >  1:3191490/1‑67 (MQ=255)
           gaaCAGATCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGaa                                                         >  1:2760775/1‑69 (MQ=255)
                           gcAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGcc                                       <  1:1744277/71‑1 (MQ=14)
                              aaGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAg                                    >  1:3762797/1‑71 (MQ=25)
                              aaGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAg                                    >  1:425037/1‑71 (MQ=25)
                              aaGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAg                                    >  1:1286372/1‑71 (MQ=25)
                                   gccgTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAg                                >  1:3145908/1‑70 (MQ=255)
                                    ccgTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGcc                                       <  1:4229799/62‑1 (MQ=14)
                                     cgTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaac                             <  1:1240979/71‑1 (MQ=255)
                                         gAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaacaac                          <  1:4159477/70‑1 (MQ=255)
                                         gAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAacaaca                         <  1:2915994/71‑1 (MQ=255)
                                              cGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCg                    <  1:3642724/71‑1 (MQ=255)
                                               gCTATTGCCCGTGCCAAAGCCCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGaac                             <  1:3274149/61‑1 (MQ=255)
                                               gCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCg                    <  1:664645/70‑1 (MQ=255)
                                               gCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCg                    <  1:2785947/70‑1 (MQ=255)
                                               gCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCg                    <  1:3968991/70‑1 (MQ=255)
                                                 tATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGAt                  >  1:305658/1‑70 (MQ=255)
                                                    tGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCcgc              >  1:1531241/1‑71 (MQ=255)
                                                    tGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCcgc              >  1:1788886/1‑71 (MQ=255)
                                                     gCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGa                               >  1:2604101/1‑53 (MQ=255)
                                                     gCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCcgcg             >  1:1667186/1‑71 (MQ=255)
                                                     gCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCcgcg             >  1:1227596/1‑71 (MQ=255)
                                                      cccGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGa                   <  1:3525823/64‑1 (MQ=255)
                                                      cccGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGa                   <  1:917563/64‑1 (MQ=255)
                                                      cccGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGa                   <  1:2316622/64‑1 (MQ=255)
                                                      cccGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCcgcg             <  1:2613876/70‑1 (MQ=255)
                                                        cGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATc                 >  1:2107620/1‑64 (MQ=255)
                                                        cGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATc                 >  1:1248955/1‑64 (MQ=255)
                                                        cGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCaa          <  1:1738119/71‑1 (MQ=255)
                                                         gTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCTCaaa         <  1:2964298/71‑1 (MQ=255)
                                                         gTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCaaa         <  1:425204/71‑1 (MQ=255)
                                                          tGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCaaa         <  1:4242266/70‑1 (MQ=255)
                                                            ccAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGa                               >  1:3118580/1‑46 (MQ=255)
                                                            ccAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGa                               >  1:2600145/1‑46 (MQ=255)
                                                            ccAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGa                               >  1:2723222/1‑46 (MQ=255)
                                                            ccAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGa                               >  1:2304922/1‑46 (MQ=255)
                                                            ccAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGa                               >  1:3749038/1‑46 (MQ=255)
                                                            ccAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGa                               >  1:3688866/1‑46 (MQ=255)
                                                            ccAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGa                               >  1:3633999/1‑46 (MQ=255)
                                                            ccAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGa                               >  1:1823360/1‑46 (MQ=255)
                                                             cAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCaa          <  1:1683490/66‑1 (MQ=255)
                                                              aaaGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAgccgc    >  1:1704835/1‑71 (MQ=255)
                                                               aaGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAg        >  1:653030/1‑66 (MQ=255)
                                                               aaGCACACAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAg        >  1:571755/1‑66 (MQ=255)
                                                                aGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGccgccg  >  1:2915859/1‑71 (MQ=255)
                                                                   |                                                                   
CGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCG  >  minE/1036438‑1036572

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: