Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,052,811	A→G	13.5%	R171R (CGT→CGC)	ydhF ←	predicted oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,052,811	0	A	G	13.5%	168.6 / 12.3	95	R171R (CGT→CGC)	ydhF	predicted oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (40/42);  new base G (9/4);  total (49/46)
Fisher's exact test for biased strand distribution p-value = 2.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

CCAGCAGTAACGGCTGATGCACCGGGGATATTTCCACCTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACGCACT  >  minE/1052747‑1052871                                                                 |                                                             ccAGCAGTAACGGCTGATGCACCGGGGATATTTCCACCTGATTAGTGGCAAGGGTAAACGGCAGACGt                                                           <  1:2647760/68‑1 (MQ=255) ccAGCAGTAACGGCTGATGCACCGGGGATATTTCCACCTGATTAGTGGCAAGGGTAAACGGCAGACGTg                                                          <  1:401625/69‑1 (MQ=255)          aCGGCTGATGCACCGGGGATATTTCCACCTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAAc                                                  >  1:1141082/1‑68 (MQ=255)          aCGGCTGATGCACCGGGGATATTTCCACCTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACa                                                 >  1:1616832/1‑69 (MQ=255)               tGATGCACCGGGGATATTTCCACCTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCa                                            <  1:1632699/69‑1 (MQ=255)                              atTTCCACCTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCTGGCGTa                             <  1:439069/69‑1 (MQ=255)                              atTTCCACCTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCg                               >  1:2692287/1‑67 (MQ=255)                              atTTCCACCTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCg                               >  1:2820344/1‑67 (MQ=255)                              atTTCCACCTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCg                               >  1:23642/1‑67 (MQ=255)                              atTTCCACCTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCg                               >  1:2257555/1‑67 (MQ=255)                              atTTCCACCTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCg                               >  1:3125056/1‑67 (MQ=255)                              atTTCCACCTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTa                             <  1:1239146/69‑1 (MQ=255)                               tttCCACCTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTaa                            >  1:2314480/1‑69 (MQ=255)                                ttCCACCTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTaaa                           >  1:2683968/1‑69 (MQ=255)                                ttCCACCTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTaaa                           >  1:3243115/1‑69 (MQ=255)                                ttCCACCTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTaaa                           >  1:2498355/1‑69 (MQ=255)                                 tCCACCTGATTAGTGGCAAGAGTAAATGGCAGGCGTGATTGCAACAGCGCAAATTGCGCAGGTGTAAAg                          <  1:2482729/69‑1 (MQ=255)                                 tCCACCTGATTAGTGGCAAGAGTAAATGGCAGGCGTGATTGCAACAGCGCAAATTGCGCAGGTGTAAAg                          <  1:1381641/69‑1 (MQ=255)                                 tCCACCTGATTAGTGGCAAGAGTAAATGGCAGGCGTGATTGCAACAGCGCAAATTGCGCAGGTGTAAAg                          <  1:704905/69‑1 (MQ=255)                                 tCCACCTGATTAGTGGCAAGAGTAAATGGCAGGCGTGATTGCAACAGCGCAAATTGCGCAGGTGTAAAg                          <  1:1916671/69‑1 (MQ=255)                                      cTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                     <  1:2911787/69‑1 (MQ=255)                                      cTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                     <  1:1302299/69‑1 (MQ=255)                                      cTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                     <  1:704843/69‑1 (MQ=255)                                      cTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                     <  1:1975040/69‑1 (MQ=255)                                      cTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                     <  1:2301494/69‑1 (MQ=255)                                      cTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                     <  1:1959746/69‑1 (MQ=255)                                      cTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                     <  1:3202855/69‑1 (MQ=255)                                      cTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                     <  1:3204867/69‑1 (MQ=255)                                      cTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                     <  1:2301983/69‑1 (MQ=255)                                      cTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                     <  1:264026/69‑1 (MQ=255)                                      cTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                     <  1:194686/69‑1 (MQ=255)                                      cTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                     <  1:1475339/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:1318748/1‑68 (MQ=255)                                            aGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGc                >  1:311767/1‑68 (MQ=255)                                            aGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGc                >  1:3016679/1‑68 (MQ=255)                                            aGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGc                >  1:3008682/1‑68 (MQ=255)                                            aGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGc                >  1:1548966/1‑68 (MQ=255)                                            aGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGc                >  1:1413519/1‑68 (MQ=255)                                            aGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGc                >  1:3137579/1‑68 (MQ=255)                                            aGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGc                >  1:35885/1‑68 (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:2905407/68‑1 (MQ=255)                                                  aaGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAAt          <  1:1973180/68‑1 (MQ=255)                                                  aaGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAAt          <  1:3296994/68‑1 (MQ=255)                                                  aaGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAAt          <  1:3354839/68‑1 (MQ=255)                                                  aaGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAAt          <  1:1655004/68‑1 (MQ=255)                                                  aaGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAAt          <  1:358961/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:2079972/3‑68 (MQ=255)                                                    gagTAAATGGCAGGCGTGATTGCAACAGCGCAAATTGCGCAGGTGTAAAGTTTGATACACCAAAATGa        >  1:1725659/3‑68 (MQ=255)                                                    gagTAAATGGCAGGCGTGATTGCAACAGCGCAAATTGCGCAGGTGTAAAGTTTGATACACCAAAATGa        >  1:226143/3‑68 (MQ=255)                                                    gagTAAATGGCAGGCGTGATTGCAACAGCGCAAATTGCGCAGGTGTAAAGTTTGATACACCAAAATGa        >  1:228634/3‑68 (MQ=255)                                                    gagTAAATGGCAGGCGTGATTGCAACAGCGCAAATTGCGCAGGTGTAAAGTTTGATACACCAAAATGa        >  1:2850419/3‑68 (MQ=255)                                                     ggTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACg      >  1:836800/1‑69 (MQ=255)                                                     ggTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACg      >  1:527407/1‑69 (MQ=255)                                                     ggTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACg      >  1:1510701/1‑69 (MQ=255)                                                     ggTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCAAAAATGAc       >  1:833407/1‑68 (MQ=255)                                                      gTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                     <  1:1157494/53‑1 (MQ=255)                                                      gTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                     <  1:1563327/53‑1 (MQ=255)                                                      gTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                     <  1:3363459/53‑1 (MQ=255)                                                      gTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                     <  1:1541985/53‑1 (MQ=255)                                                      gTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                     <  1:552441/53‑1 (MQ=255)                                                      gTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                     <  1:804291/53‑1 (MQ=255)                                                      gTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                     <  1:2824610/53‑1 (MQ=255)                                                      gTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                     <  1:653254/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:69901/1‑69 (MQ=255)                                                         aaCGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACGCACt  >  1:1416838/1‑69 (MQ=255)                                                         aaCGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACGCACt  >  1:741755/1‑69 (MQ=255)                                                         aaCGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACGCACt  >  1:1603071/1‑69 (MQ=255)                                                         aaCGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACGCACt  >  1:483056/1‑69 (MQ=255)                                                         aaCGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACGCACt  >  1:431285/1‑69 (MQ=255)                                                         aaCGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACGCACt  >  1:264728/1‑69 (MQ=255)                                                         aaCGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACGCACt  >  1:3307083/1‑69 (MQ=255)                                                         aaCGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACGCACt  >  1:3289600/1‑69 (MQ=255)                                                         aaCGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACGCACt  >  1:1856311/1‑69 (MQ=255)                                                         aaCGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACGCACt  >  1:1947906/1‑69 (MQ=255)                                                         aaCGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACGCACt  >  1:2031768/1‑69 (MQ=255)                                                         aaCGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACGCACt  >  1:2189387/1‑69 (MQ=255)                                                         aaCGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACGCACt  >  1:2372476/1‑69 (MQ=255)                                                              cAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGa                     <  1:2196366/45‑1 (MQ=255)                                                                 |                                                             CCAGCAGTAACGGCTGATGCACCGGGGATATTTCCACCTGATTAGTGGCAAGGGTAAACGGCAGACGTGATTGCAACAGGGCAAATTGCGCAGGCGTAAAGTTCGATACGCCAAAATGACGCACT  >  minE/1052747‑1052871

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