Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,664,043	C→T	64.5%	L393L (CTG→CTA)	guaA ←	GMP synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,664,043	0	C	T	64.5%	46.4 / 44.7	76	L393L (CTG→CTA)	guaA	GMP synthetase
Reads supporting (aligned to +/- strand):  ref base C (16/11);  new base T (11/38);  total (27/49)
Fisher's exact test for biased strand distribution p-value = 2.35e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

CCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTC  >  minE/1663982‑1664094                                                              |                                                    ccAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGccc                                               <  1:1215090/68‑1 (MQ=255) ccAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGccc                                               <  1:446647/68‑1 (MQ=255) ccAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGccc                                               <  1:2466639/68‑1 (MQ=255)     aaCACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTc                                          <  1:1090982/69‑1 (MQ=255)     aaCACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTc                                          <  1:3085974/69‑1 (MQ=255)     aaCACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTc                                          <  1:3081904/69‑1 (MQ=255)     aaCACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTc                                          <  1:3185094/69‑1 (MQ=255)     aaCACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTc                                          <  1:2920766/69‑1 (MQ=255)     aaCACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTc                                          <  1:243841/69‑1 (MQ=255)     aaCACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTc                                          <  1:2276346/69‑1 (MQ=255)     aaCACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTc                                          <  1:2188998/69‑1 (MQ=255)     aaCACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTc                                          <  1:2139310/69‑1 (MQ=255)     aaCACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTc                                          <  1:1673375/69‑1 (MQ=255)     aaCACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTc                                          <  1:1513509/69‑1 (MQ=255)     aaCACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTc                                          <  1:1504727/69‑1 (MQ=255)     aaCACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTc                                          <  1:14793/69‑1 (MQ=255)     aaCACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTc                                          <  1:1367041/69‑1 (MQ=255)     aaCACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTc                                          <  1:1431082/69‑1 (MQ=255)      acacGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTc                                          <  1:1630358/68‑1 (MQ=255)       gacGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCa                                        >  1:2709860/2‑69 (MQ=255)       cacGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCa                                        >  1:2337563/1‑69 (MQ=255)       cacGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCa                                        >  1:1885355/1‑69 (MQ=255)       cacGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCa                                        >  1:2639289/1‑69 (MQ=255)       cacGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCa                                        >  1:431136/1‑69 (MQ=255)       cacGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCa                                        >  1:2819081/1‑69 (MQ=255)       cacGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCa                                        >  1:2846999/1‑69 (MQ=255)           aaCGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGaa                                     >  1:1235681/1‑67 (MQ=255)           aaCGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGGcc                                    >  1:3122633/1‑69 (MQ=255)           aaCGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGAcc                                    >  1:2459164/1‑69 (MQ=255)           aaCGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGAcc                                    >  1:283012/1‑69 (MQ=255)           aaCGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGAcc                                    >  1:3017798/1‑69 (MQ=255)           aaCGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGAcc                                    >  1:3208853/1‑69 (MQ=255)           aaCGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGAcc                                    >  1:1917976/1‑69 (MQ=255)           aaCGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGAcc                                    >  1:3184735/1‑69 (MQ=255)           aaCGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGAcc                                    >  1:179438/1‑69 (MQ=255)             cGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTCCAGACCaa                                  >  1:1254743/1‑69 (MQ=255)             cGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTCCAGACCaa                                  >  1:1286442/1‑69 (MQ=255)             cGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTCCAGACCaa                                  >  1:626340/1‑69 (MQ=255)             cGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTCCAGACCaa                                  >  1:3160164/1‑69 (MQ=255)             cGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTCCAGACCaa                                  >  1:2426778/1‑69 (MQ=255)             cGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTCCAGACCaa                                  >  1:1309051/1‑69 (MQ=255)             cGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTCCAGACCaa                                  >  1:1638860/1‑69 (MQ=255)             cGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTCCAGACCaa                                  >  1:18754/1‑69 (MQ=255)             cGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTCCAGACCaa                                  >  1:814228/1‑69 (MQ=255)             cGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTCCAGACCaa                                  >  1:203428/1‑69 (MQ=255)             cGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTCCAGACCaa                                  >  1:806339/1‑69 (MQ=255)                            ccGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGccc                                               <  1:2997472/41‑1 (MQ=255)                               ggAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTc                                          <  1:1613144/43‑1 (MQ=255)                                      gTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAg         <  1:1201965/69‑1 (MQ=255)                                      gTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAg         <  1:3083738/69‑1 (MQ=255)                                       tGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAg         <  1:823626/68‑1 (MQ=255)                                       tGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGc        <  1:1926311/69‑1 (MQ=255)                                       tGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGc        <  1:1924579/69‑1 (MQ=255)                                       tGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGc        <  1:18336/69‑1 (MQ=255)                                       tGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGc        <  1:2289271/69‑1 (MQ=255)                                             tACAACATGTCATACGGTAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTc  <  1:429232/69‑1 (MQ=255)                                             tACAACATGTCATACGGTAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTc  <  1:3392861/69‑1 (MQ=255)                                             tACAACATGTCATACGGTAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTc  <  1:882541/69‑1 (MQ=255)                                             tACAACATGTCATACGGTAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTc  <  1:74780/69‑1 (MQ=255)                                             tACAACATGTCATACGGTAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTc  <  1:995702/69‑1 (MQ=255)                                             tACAACATGTCATACGGTAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTc  <  1:744288/69‑1 (MQ=255)                                             tACAACATGTCATACGGTAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTc  <  1:3111593/69‑1 (MQ=255)                                             tACAACATGTCATACGGTAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTc  <  1:3055034/69‑1 (MQ=255)                                             tACAACATGTCATACGGTAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTc  <  1:3052728/69‑1 (MQ=255)                                             tACAACATGTCATACGGTAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTc  <  1:2979555/69‑1 (MQ=255)                                             tACAACATGTCATACGGTAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTc  <  1:2798081/69‑1 (MQ=255)                                             tACAACATGTCATACGGTAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTc  <  1:2698047/69‑1 (MQ=255)                                             tACAACATGTCATACGGTAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTc  <  1:2543452/69‑1 (MQ=255)                                             tACAACATGTCATACGGTAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTc  <  1:2539140/69‑1 (MQ=255)                                             tACAACATGTCATACGGTAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTc  <  1:2337841/69‑1 (MQ=255)                                             tACAACATGTCATACGGTAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTc  <  1:2277340/69‑1 (MQ=255)                                             tACAACATGTCATACGGTAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTc  <  1:2116762/69‑1 (MQ=255)                                             tACAACATGTCATACGGTAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTc  <  1:1963477/69‑1 (MQ=255)                                             tACAACATGTCATACGGTAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTc  <  1:1423014/69‑1 (MQ=255)                                             tACAACATGTCATACGGTAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTc  <  1:1379565/69‑1 (MQ=255)                                             tACAACATGTCATACGGTAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTc  <  1:1093061/69‑1 (MQ=255)                                                              |                                                    CCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTC  >  minE/1663982‑1664094

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