Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,664,030	G→A	57.0%	L398L (CTG→TTG)	guaA ←	GMP synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,664,030	0	G	A	57.0%	21.2 / 66.4	86	L398L (CTG→TTG)	guaA	GMP synthetase
Reads supporting (aligned to +/- strand):  ref base G (16/21);  new base A (11/38);  total (27/59)
Fisher's exact test for biased strand distribution p-value = 5.98e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

ACTCTTTCTTCACTTCACCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTC  >  minE/1663964‑1664094                                                                   |                                                                 aCTCTTTCTTCACTTCACCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCa                                                                <  1:917442/69‑1 (MQ=255) aCTCTTTCTTCACTTCACCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCa                                                                <  1:1047196/69‑1 (MQ=255) aCTCTTTCTTCACTTCACCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCa                                                                <  1:1258924/69‑1 (MQ=255) aCTCTTTCTTCACTTCACCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCa                                                                <  1:1286350/69‑1 (MQ=255) aCTCTTTCTTCACTTCACCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCa                                                                <  1:3251110/69‑1 (MQ=255) aCTCTTTCTTCACTTCACCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCa                                                                <  1:3072808/69‑1 (MQ=255) aCTCTTTCTTCACTTCACCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCa                                                                <  1:3025479/69‑1 (MQ=255) aCTCTTTCTTCACTTCACCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCa                                                                <  1:217614/69‑1 (MQ=255) aCTCTTTCTTCACTTCACCAAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCa                                                                <  1:1971810/69‑1 (MQ=255)  ctctTTCTTCACTTCTCCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCa                                                                <  1:2596148/68‑1 (MQ=255)                   ccAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGccc                                               <  1:446647/68‑1 (MQ=255)                   ccAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGccc                                               <  1:1215090/68‑1 (MQ=255)                   ccAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGccc                                               <  1:2466639/68‑1 (MQ=255)                       aaCACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTc                                          <  1:3081904/69‑1 (MQ=255)                       aaCACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTc                                          <  1:1673375/69‑1 (MQ=255)                       aaCACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTc                                          <  1:2920766/69‑1 (MQ=255)                       aaCACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTc                                          <  1:2139310/69‑1 (MQ=255)                       aaCACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTc                                          <  1:2188998/69‑1 (MQ=255)                       aaCACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTc                                          <  1:2276346/69‑1 (MQ=255)                       aaCACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTc                                          <  1:243841/69‑1 (MQ=255)                       aaCACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTc                                          <  1:1090982/69‑1 (MQ=255)                       aaCACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTc                                          <  1:1513509/69‑1 (MQ=255)                       aaCACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTc                                          <  1:1504727/69‑1 (MQ=255)                       aaCACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTc                                          <  1:1367041/69‑1 (MQ=255)                       aaCACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTc                                          <  1:3185094/69‑1 (MQ=255)                       aaCACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTc                                          <  1:14793/69‑1 (MQ=255)                       aaCACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTc                                          <  1:3085974/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:2819081/1‑69 (MQ=255)                         cacGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCa                                        >  1:431136/1‑69 (MQ=255)                         cacGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCa                                        >  1:2337563/1‑69 (MQ=255)                         cacGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCa                                        >  1:2639289/1‑69 (MQ=255)                         cacGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCa                                        >  1:2846999/1‑69 (MQ=255)                         cacGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCa                                        >  1:1885355/1‑69 (MQ=255)                             aaCGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGaa                                     >  1:1235681/1‑67 (MQ=255)                             aaCGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGGcc                                    >  1:3122633/1‑69 (MQ=255)                             aaCGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGAcc                                    >  1:3184735/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:2459164/1‑69 (MQ=255)                             aaCGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGAcc                                    >  1:179438/1‑69 (MQ=255)                             aaCGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGAcc                                    >  1:1917976/1‑69 (MQ=255)                               cGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTCCAGACCaa                                  >  1:2426778/1‑69 (MQ=255)                               cGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTCCAGACCaa                                  >  1:1309051/1‑69 (MQ=255)                               cGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTCCAGACCaa                                  >  1:1286442/1‑69 (MQ=255)                               cGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTCCAGACCaa                                  >  1:3160164/1‑69 (MQ=255)                               cGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTCCAGACCaa                                  >  1:1638860/1‑69 (MQ=255)                               cGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTCCAGACCaa                                  >  1:1254743/1‑69 (MQ=255)                               cGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTCCAGACCaa                                  >  1:18754/1‑69 (MQ=255)                               cGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTCCAGACCaa                                  >  1:626340/1‑69 (MQ=255)                               cGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTCCAGACCaa                                  >  1:203428/1‑69 (MQ=255)                               cGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTCCAGACCaa                                  >  1:806339/1‑69 (MQ=255)                               cGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTCCAGACCaa                                  >  1:814228/1‑69 (MQ=255)                                              ccGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGccc                                               <  1:2997472/41‑1 (MQ=255)                                                 ggAACGGGTGACGGTACAACATGTCATACGGTAGGCCCAGCTc                                          <  1:1613144/43‑1 (MQ=255)                                                        gTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAg         <  1:3083738/69‑1 (MQ=255)                                                        gTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAg         <  1:1201965/69‑1 (MQ=255)                                                         tGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAg         <  1:823626/68‑1 (MQ=255)                                                         tGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGc        <  1:1926311/69‑1 (MQ=255)                                                         tGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGc        <  1:2289271/69‑1 (MQ=255)                                                         tGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGc        <  1:18336/69‑1 (MQ=255)                                                         tGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGc        <  1:1924579/69‑1 (MQ=255)                                                               tACAACATGTCATACGGTAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTc  <  1:74780/69‑1 (MQ=255)                                                               tACAACATGTCATACGGTAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTc  <  1:429232/69‑1 (MQ=255)                                                               tACAACATGTCATACGGTAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTc  <  1:744288/69‑1 (MQ=255)                                                               tACAACATGTCATACGGTAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTc  <  1:882541/69‑1 (MQ=255)                                                               tACAACATGTCATACGGTAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTc  <  1:3392861/69‑1 (MQ=255)                                                               tACAACATGTCATACGGTAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTc  <  1:1093061/69‑1 (MQ=255)                                                               tACAACATGTCATACGGTAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTc  <  1:995702/69‑1 (MQ=255)                                                               tACAACATGTCATACGGTAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTc  <  1:1379565/69‑1 (MQ=255)                                                               tACAACATGTCATACGGTAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTc  <  1:1423014/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:1963477/69‑1 (MQ=255)                                                               tACAACATGTCATACGGTAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTc  <  1:2798081/69‑1 (MQ=255)                                                               tACAACATGTCATACGGTAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTc  <  1:2116762/69‑1 (MQ=255)                                                               tACAACATGTCATACGGTAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTc  <  1:2698047/69‑1 (MQ=255)                                                               tACAACATGTCATACGGTAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTc  <  1:2277340/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)                                                                   |                                                                 ACTCTTTCTTCACTTCACCCAGAACACGAACGCCAAGGCCTGGTCCCGGGAACGGGTGACGGTACAGCATGTCGTACGGCAGGCCCAGCTCCAGACCAATCTTACGCACTTCGTCTTTGAACAGCTCTTTC  >  minE/1663964‑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: ‑