Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 461,006 T→G 100% intergenic (‑83/‑626) gltA ← / → sdhC citrate synthase/succinate dehydrogenase, membrane subunit, binds cytochrome b556

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE461,0060TG88.7% 75.4 / 1.1 35intergenic (‑83/‑626)gltA/sdhCcitrate synthase/succinate dehydrogenase, membrane subunit, binds cytochrome b556
Reads supporting (aligned to +/- strand):  ref base T (1/3);  new base G (31/0);  total (32/3)
Fisher's exact test for biased strand distribution p-value = 6.11e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCTCGCCGCTCTGTGAAAGAGGGGAAAACCTGGGTACAGAGCTCTGGGCGCTTGCAGG  >  minE/460940‑461071
                                                                  |                                                                 
gcCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGtttttt                                                                 <  1:1057568/69‑1 (MQ=255)
gcCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGtttttt                                                                 <  1:2560284/69‑1 (MQ=255)
gcCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGtttttt                                                                 <  1:1800697/69‑1 (MQ=255)
 ccTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGtttttt                                                                 <  1:3103989/68‑1 (MQ=255)
           tAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCTcgc                                                      >  1:1774753/1‑69 (MQ=255)
                       aaCTTAAATTTGCCTTCGCACATCAACCTGTCTTTACCCGTTTTTTATTTGGCTCGCCGCTCTGTGaaa                                          <  1:1138715/69‑1 (MQ=255)
                        aCTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTGTTATTTGGCTCGCCGCTCTGTGAAag                                         >  1:932892/1‑69 (MQ=255)
                        aCTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTGTTATTTGGCTCGCCGCTCTGTGAAag                                         >  1:535493/1‑69 (MQ=255)
                        aCTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTGTTATTTGGCTCGCCGCTCTGTGAAag                                         >  1:3339949/1‑69 (MQ=255)
                        aCTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTGTTATTTGGCTCGCCGCTCTGTGAAag                                         >  1:2596151/1‑69 (MQ=255)
                        aCTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTGTTATTTGGCTCGCCGCTCTGTGAAag                                         >  1:2457237/1‑69 (MQ=255)
                        aCTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTGTTATTTGGCTCGCCGCTCTGTGAAag                                         >  1:2119353/1‑69 (MQ=255)
                        aCTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTGTTATTTGGCTCGCCGCTCTGTGAAag                                         >  1:2044344/1‑69 (MQ=255)
                        aCTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTGTTATTTGGCTCGCCGCTCTGTGAAag                                         >  1:1705211/1‑69 (MQ=255)
                        aCTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTGTTATTTGGCTCGCCGCTCTGTGAAag                                         >  1:1561326/1‑69 (MQ=255)
                        aCTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTGTTATTTGGCTCGCCGCTCTGTGAAag                                         >  1:1400806/1‑69 (MQ=255)
                        aCTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTGTTATTTGGCTCGCCGCTCTGTGAAag                                         >  1:1226270/1‑69 (MQ=255)
                        aCTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTGTTATTTGGCTCGACGCTCTGTGAAag                                         >  1:1469693/1‑69 (MQ=255)
                                         cacaTCAACCTGGCTTTACCCGTTTTTTATTTGGCTCGCCGCTCTGTGAAAGAGGGGAAAACCTggg                          <  1:2253258/67‑1 (MQ=255)
                                                       tttACCCGTTTTTTATTTGGCTCGCCGCTCTGTGAAAGAGGGGAAAACCTGGGTACAGAGCTCTGGgcg          <  1:877105/69‑1 (MQ=255)
                                                          aCCCGTTTGTTATTTGGCTCGCCGCTCTGTGAAagag                                       >  1:2486850/1‑37 (MQ=255)
                                                          aCCCGTTTGTTATTTGGCTCGCCGCTCTGTGAAAGAGGGGAAAACCTGGGTACAGAGCTCTggg            >  1:3059433/1‑64 (MQ=255)
                                                          aCCCGTTTGTTATTTGGCTCGCCGCTCTGTGAAAGAGGGGAAAACCTGGGTACAGAGCTCTGGgcgc         >  1:2921928/1‑67 (MQ=255)
                                                          aCCCGTTTGTTATTTGGCTCGCCGCTCTGTGAAAGAGGGGAAAACCTGGGTACAGAGCTCTGGgcgc         >  1:1646170/1‑67 (MQ=255)
                                                          aCCCGTTTGTTATTTGGCTCGCCGCTCTGTGAAAGAGGGGAAAACCTGGGTACAGAGCTCTGGgcgc         >  1:1233648/1‑67 (MQ=255)
                                                          aCCCGTTTGTTATTTGGCTCGCCGCTCTGTGAAAGAGGGGAAAACCTGGGTACAGAGCTCTGGgcgc         >  1:3386838/1‑67 (MQ=255)
                                                          aCCCGTTTGTTATTTGGCTCGCCGCTCTGTGAAAGAGGGGAAAACCTGGGTACAGAGCTCTGGgcgc         >  1:486599/1‑67 (MQ=255)
                                                          aCCCGTTTGTTATTTGGCTCGCCGCTCTGTGAAAGAGGGGAAAACCTGGGTACAGAGCTCTGGgcgc         >  1:1214975/1‑67 (MQ=255)
                                                          aCCCGTTTGTTATTTGGCTCGCCGCTCTGTGAAAGAGGGGAAAACCTGGGTACAGAGCTCTGGgcgc         >  1:859327/1‑67 (MQ=255)
                                                          aCCCGTTTGTTATTTGGCTCGCCGCTCTGTGAAAGAGGGGAAAACCTGGGTACAGAGCTCTGGgcgc         >  1:882058/1‑67 (MQ=255)
                                                          aCCCGTTTGTTATTTGGCTCGCCGCTCTGTGAAAGAGGGGAAAACCTGGGTACAGAGCTCTGGgcgc         >  1:2689225/1‑67 (MQ=255)
                                                          aCCCGTTTGTTATTTGGCTCGCCGCTCTGTGAAAGAGGGGAAAACCTGGGTACAGAGCTCTGGgcgc         >  1:2324409/1‑67 (MQ=255)
                                                          aCCCGTTTGTTATTTGGCTCGCCGCTCTGTGAAAGAGGGGAAAACCTGGGTACAGAGCTCTGGgcgc         >  1:2185015/1‑67 (MQ=255)
                                                          aCCCGTTTGTTATTTGGCTCGCCGCTCTGTGAAAGAGGGGAAAACCTGGGTACAGAGCTCTGGgcgc         >  1:1546912/1‑67 (MQ=255)
                                                          aCCCGTTTGTTATTTGGCTCGCCGCTCTGTGAAAGAGGGGAAAACCTGGGTACAGAGCTCTGGgcgc         >  1:2094867/1‑67 (MQ=255)
                                                          aCCCGTTTGTTATTTGGCTCGCCGCTCTGTGAAAGAGGGGAAAACCTGGGTACAGAGCTCTGGgcg          >  1:378302/1‑66 (MQ=255)
                                                          aCCCGTTTGTTATTTGGCTCGCCGCTCTGTGAAAGAGGGGAAAACCTGGGTACAGAGCTCTGGGCGCt        >  1:3124423/1‑68 (MQ=255)
                                                          aCCCGTTTGTTATTTGGCTCGCCGCCCTGTGAAAGAGGGGAAAACCTGGGTACAGAGCTCTGGgcgc         >  1:2883218/1‑67 (MQ=255)
                                                          aCCCGTTTGTTATTTGGCTCGCAGCTCTGTGAAAGAGGGGAAAACCTGGGTACAGAGCTCTGGgcgc         >  1:925362/1‑67 (MQ=255)
                                                                 ttttATTTGGCTCGCCGCTCTGTTAAAGAGGGGAAAACCTGGGTACAGAGCTCTGGGCGCTTGCAgg  >  1:1011352/1‑67 (MQ=255)
                                                                 ttttATTTGGCTCGCCGCTCTGTGAAAGAGGGGAAAACCTGGGTACAGAGCTCTGGGCGCTTGCAgg  >  1:2688903/1‑67 (MQ=255)
                                                                 ttttATTTGGCTCGCCGCTCTGTGAAAGAGGGGAAAACCTGGGTACAGAGCTCTGGGCGCTTGCAgg  >  1:2635269/1‑67 (MQ=255)
                                                                 ttttATTTGGCTCGCCGCTCTGTGAAAGAGGGGAAAACCTGGGTACAGAGCTCTGGGCGCTTGCAgg  >  1:25063/1‑67 (MQ=255)
                                                                 ttttATTTGGCTCGCCGCTCTGTGAAAGAGGGGAAAACCTGGGTACAGAGCTCTGGGCGCTTGCAgg  >  1:2457807/1‑67 (MQ=255)
                                                                 ttttATTTGGCTCGCCGCTCTGTGAAAGAGGGGAAAACCTGGGTACAGAGCTCTGGGCGCTTGCAgg  >  1:1501512/1‑67 (MQ=255)
                                                                 ttttATTTGGCTCGCCGCTCTGTGAAAGAGGGGAAAACCTGGGTACAGAGCTCTGGGCGCTTGCAgg  >  1:648079/1‑67 (MQ=255)
                                                                  |                                                                 
GCCTTATTGCGTAAGACTGCCGGAACTTAAATTTGCCTTCGCACATCAACCTGGCTTTACCCGTTTTTTATTTGGCTCGCCGCTCTGTGAAAGAGGGGAAAACCTGGGTACAGAGCTCTGGGCGCTTGCAGG  >  minE/460940‑461071

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: