Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,052,826	G→C	13.6%	A166A (GCC→GCG)	ydhF ←	predicted oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,052,826	0	G	C	13.6%	207.3 / 19.9	95	A166A (GCC→GCG)	ydhF	predicted oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (42/40);  new base C (9/4);  total (51/44)
Fisher's exact test for biased strand distribution p-value = 2.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGATGCACCGGGGATATTTCCACCTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACGCACTTTGCCGCTCTGA  >  minE/1052761‑1052883                                                                  |                                                          tGATGCACCGGGGATATTTCCACCTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCa                                                        <  1:1632699/69‑1 (MQ=255)                atTTCCACCTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCTGGCGTa                                         <  1:439069/69‑1 (MQ=255)                atTTCCACCTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCg                                           >  1:3125056/1‑67 (MQ=255)                atTTCCACCTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCg                                           >  1:23642/1‑67 (MQ=255)                atTTCCACCTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCg                                           >  1:2820344/1‑67 (MQ=255)                atTTCCACCTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCg                                           >  1:2257555/1‑67 (MQ=255)                atTTCCACCTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCg                                           >  1:2692287/1‑67 (MQ=255)                atTTCCACCTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTa                                         <  1:1239146/69‑1 (MQ=255)                 tttCCACCTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTaa                                        >  1:2314480/1‑69 (MQ=255)                  ttCCACCTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTaaa                                       >  1:3243115/1‑69 (MQ=255)                  ttCCACCTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTaaa                                       >  1:2683968/1‑69 (MQ=255)                  ttCCACCTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTaaa                                       >  1:2498355/1‑69 (MQ=255)                   tCCACCTGATTAGTGGCAAGAGTAAATGGCAGGCGTGATTGCAACAGCGCAAATTGCGCAGGTGTAAAg                                      <  1:1381641/69‑1 (MQ=255)                   tCCACCTGATTAGTGGCAAGAGTAAATGGCAGGCGTGATTGCAACAGCGCAAATTGCGCAGGTGTAAAg                                      <  1:704905/69‑1 (MQ=255)                   tCCACCTGATTAGTGGCAAGAGTAAATGGCAGGCGTGATTGCAACAGCGCAAATTGCGCAGGTGTAAAg                                      <  1:2482729/69‑1 (MQ=255)                   tCCACCTGATTAGTGGCAAGAGTAAATGGCAGGCGTGATTGCAACAGCGCAAATTGCGCAGGTGTAAAg                                      <  1:1916671/69‑1 (MQ=255)                        cTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                                 <  1:2301983/69‑1 (MQ=255)                        cTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                                 <  1:2301494/69‑1 (MQ=255)                        cTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                                 <  1:2911787/69‑1 (MQ=255)                        cTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                                 <  1:264026/69‑1 (MQ=255)                        cTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                                 <  1:1975040/69‑1 (MQ=255)                        cTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                                 <  1:1959746/69‑1 (MQ=255)                        cTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                                 <  1:3202855/69‑1 (MQ=255)                        cTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                                 <  1:194686/69‑1 (MQ=255)                        cTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                                 <  1:3204867/69‑1 (MQ=255)                        cTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                                 <  1:1475339/69‑1 (MQ=255)                        cTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                                 <  1:1302299/69‑1 (MQ=255)                        cTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                                 <  1:704843/69‑1 (MQ=255)                        cTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                                 <  1:493521/69‑1 (MQ=255)                        cTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                                 <  1:1338886/69‑1 (MQ=255)                         tGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                                 <  1:1533344/68‑1 (MQ=255)                              aggggCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACg                             >  1:727089/4‑67 (MQ=255)                              aGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACg                             >  1:1509399/1‑67 (MQ=255)                              aGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGc                            >  1:311767/1‑68 (MQ=255)                              aGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGc                            >  1:1413519/1‑68 (MQ=255)                              aGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGc                            >  1:3008682/1‑68 (MQ=255)                              aGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGc                            >  1:3016679/1‑68 (MQ=255)                              aGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGc                            >  1:1548966/1‑68 (MQ=255)                              aGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGc                            >  1:1318748/1‑68 (MQ=255)                              aGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGc                            >  1:35885/1‑68 (MQ=255)                              aGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGc                            >  1:3137579/1‑68 (MQ=255)                                    aaGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAAt                      <  1:2905407/68‑1 (MQ=255)                                    aaGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAAt                      <  1:798268/68‑1 (MQ=255)                                    aaGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAAt                      <  1:1417111/68‑1 (MQ=255)                                    aaGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAAt                      <  1:2704772/68‑1 (MQ=255)                                    aaGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAAt                      <  1:358961/68‑1 (MQ=255)                                    aaGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAAt                      <  1:1973180/68‑1 (MQ=255)                                    aaGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAAt                      <  1:1655004/68‑1 (MQ=255)                                    aaGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAAt                      <  1:3354839/68‑1 (MQ=255)                                    aaGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAAt                      <  1:3296994/68‑1 (MQ=255)                                    aaGGGTAAACGGCAGACGTGATTCCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAAt                      <  1:859448/68‑1 (MQ=255)                                      gagTAAATGGCAGGCGTGATTGCAAGAGCGCAAATTGCGCAGGTGTAAAGTTTGATACACCAAAATGa                    >  1:1770159/3‑68 (MQ=255)                                      gagTAAATGGCAGGCGTGATTGCAACAGCGCAAATTGCGCAGGTGTAAAGTTTGATACACGAAAATGa                    >  1:1832556/3‑68 (MQ=255)                                      gagTAAATGGCAGGCGTGATTGCAACAGCGCAAATTGCGCAGGTGTAAAGTTTGATACACCAAAGTGa                    >  1:195276/3‑68 (MQ=255)                                      gagTAAATGGCAGGCGTGATTGCAACAGCGCAAATTGCGCAGGTGTAAAGTTTGATACACCAAAATg                     >  1:1164629/3‑67 (MQ=255)                                      gagTAAATGGCAGGCGTGATTGCAACAGCGCAAATTGCGCAGGTGTAAAGTTTGATACACCAAAATGa                    >  1:228634/3‑68 (MQ=255)                                      gagTAAATGGCAGGCGTGATTGCAACAGCGCAAATTGCGCAGGTGTAAAGTTTGATACACCAAAATGa                    >  1:226143/3‑68 (MQ=255)                                      gagTAAATGGCAGGCGTGATTGCAACAGCGCAAATTGCGCAGGTGTAAAGTTTGATACACCAAAATGa                    >  1:1725659/3‑68 (MQ=255)                                      gagTAAATGGCAGGCGTGATTGCAACAGCGCAAATTGCGCAGGTGTAAAGTTTGATACACCAAAATGa                    >  1:2850419/3‑68 (MQ=255)                                      gagTAAATGGCAGGCGTGATTGCAACAGCGCAAATTGCGCAGGTGTAAAGTTTGATACACCAAAATGa                    >  1:2079972/3‑68 (MQ=255)                                       ggTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACg                  >  1:836800/1‑69 (MQ=255)                                       ggTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACg                  >  1:1510701/1‑69 (MQ=255)                                       ggTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACg                  >  1:527407/1‑69 (MQ=255)                                       ggTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCAAAAATGAc                   >  1:833407/1‑68 (MQ=255)                                        gTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                                 <  1:3363459/53‑1 (MQ=255)                                        gTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                                 <  1:1563327/53‑1 (MQ=255)                                        gTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                                 <  1:1541985/53‑1 (MQ=255)                                        gTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                                 <  1:552441/53‑1 (MQ=255)                                        gTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                                 <  1:653254/53‑1 (MQ=255)                                        gTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                                 <  1:1157494/53‑1 (MQ=255)                                        gTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                                 <  1:2824610/53‑1 (MQ=255)                                        gTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                                 <  1:804291/53‑1 (MQ=255)                                         tAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                                 <  1:1361329/52‑1 (MQ=255)                                         tAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACGCa                <  1:1994274/69‑1 (MQ=255)                                           aaCGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                                 <  1:1079648/50‑1 (MQ=255)                                           aaCGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACGCAc               >  1:760744/1‑68 (MQ=255)                                           aaCGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACGCACt              >  1:1416838/1‑69 (MQ=255)                                           aaCGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACGCACt              >  1:2372476/1‑69 (MQ=255)                                           aaCGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACGCACt              >  1:264728/1‑69 (MQ=255)                                           aaCGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACGCACt              >  1:2189387/1‑69 (MQ=255)                                           aaCGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACGCACt              >  1:741755/1‑69 (MQ=255)                                           aaCGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACGCACt              >  1:1603071/1‑69 (MQ=255)                                           aaCGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACGCACt              >  1:69901/1‑69 (MQ=255)                                           aaCGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACGCACt              >  1:2031768/1‑69 (MQ=255)                                           aaCGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACGCACt              >  1:1947906/1‑69 (MQ=255)                                           aaCGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACGCACt              >  1:1856311/1‑69 (MQ=255)                                           aaCGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACGCACt              >  1:483056/1‑69 (MQ=255)                                           aaCGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACGCACt              >  1:3289600/1‑69 (MQ=255)                                           aaCGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACGCACt              >  1:431285/1‑69 (MQ=255)                                           aaCGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACGCACt              >  1:3307083/1‑69 (MQ=255)                                                cAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                                 <  1:2196366/45‑1 (MQ=255)                                                           gCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACGCACTTTGCCGCTCTGa  >  1:2019340/1‑65 (MQ=255)                                                           gCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACGCACTTTGCCGCTCTGa  >  1:3009415/1‑65 (MQ=255)                                                           gCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACGCACTTTGCCGCTCTGa  >  1:2075841/1‑65 (MQ=255)                                                           gCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACGCACTTTGCCGCTCTGa  >  1:2302806/1‑65 (MQ=255)                                                                  |                                                          TGATGCACCGGGGATATTTCCACCTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACGCACTTTGCCGCTCTGA  >  minE/1052761‑1052883

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: ‑