Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,983,953	G→C	12.7%	A93A (GCG→GCC)	yggV →	dITP/XTP pyrophosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,983,953	0	G	C	12.7%	97.3 / 4.9	39	A93A (GCG→GCC)	yggV	dITP/XTP pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base G (26/8);  new base C (3/2);  total (29/10)
Fisher's exact test for biased strand distribution p-value = 5.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTACCGGACGACCAACGTCAGG  >  minE/1983884‑1984020                                                                      |                                                                    tggtctggCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGAcgcg                                                                     >  1:3118601/1‑70 (MQ=255) tggtctggCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGa                                                                    >  1:2571640/1‑71 (MQ=255)   gtctggCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGAcgcg                                                                     <  1:546813/68‑1 (MQ=255)                 gtgCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCt                                                    >  1:3100210/1‑71 (MQ=255)                 gtgCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCt                                                    >  1:2734528/1‑71 (MQ=255)                    cTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGc                                                  >  1:2279745/1‑70 (MQ=255)                            cgcgcCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCaaa                                               >  1:3123415/1‑65 (MQ=255)                            cgcgcCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAActgc                                         >  1:563328/1‑71 (MQ=255)                            cgcgcCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAActgc                                         >  1:544077/1‑71 (MQ=255)                             gcgcCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAActgct                                        >  1:1930579/1‑71 (MQ=255)                             gcgcCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAActgct                                        >  1:3234372/1‑71 (MQ=255)                                ccGGGCTTTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGa                                     <  1:992213/71‑1 (MQ=255)                                ccGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTgg                                      >  1:1728362/1‑70 (MQ=255)                                   ggCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAAc                                  <  1:506241/71‑1 (MQ=255)                                    gCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAAc                                  >  1:1301806/1‑70 (MQ=255)                                    gCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAAc                                  >  1:169327/1‑70 (MQ=255)                                         tACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCaaaaa                                             >  1:3286338/1‑54 (MQ=255)                                         tACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCaaaaa                                             >  1:496449/1‑54 (MQ=255)                                         tACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATg                              >  1:413127/1‑69 (MQ=255)                                         tACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGa                             >  1:726259/1‑70 (MQ=255)                                         tACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGa                             >  1:4084183/1‑70 (MQ=255)                                         tACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGa                             <  1:368056/70‑1 (MQ=255)                                         tACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGa                             >  1:2807957/1‑70 (MQ=255)                                                cgcgTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTa                     <  1:670343/71‑1 (MQ=255)                                                   gTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTACCg                  <  1:254329/71‑1 (MQ=255)                                                    ttATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTACCg                  <  1:726214/70‑1 (MQ=255)                                                     tATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTTCTGGAAACAATGAAAGACGTACCGga                >  1:3121199/1‑71 (MQ=255)                                                     tATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTACCGga                >  1:2152344/1‑71 (MQ=255)                                                     tATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTACCGga                >  1:4221463/1‑71 (MQ=255)                                                      aTTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTACCGgac               >  1:1648346/1‑71 (MQ=255)                                                      aTTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTACCGga                >  1:1970104/1‑70 (MQ=255)                                                       ttCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTACCGgac               >  1:2622069/1‑70 (MQ=255)                                                             tGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACa                                 >  1:2774104/1‑46 (MQ=255)                                                             tGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACa                                 >  1:2765040/1‑46 (MQ=255)                                                             agAAGACGCCACCGATCTAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTACCGGACGACCAAc        >  1:2845900/2‑71 (MQ=255)                                                             agAAGACGCCACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTACCGGACGACCAAc        >  1:3279317/2‑71 (MQ=255)                                                             agAAGACGCCACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTACCGGACGACCAAc        >  1:3140890/2‑71 (MQ=255)                                                               aaGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAAc                                  <  1:1840003/43‑1 (MQ=255)                                                               aaGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAAc                                  <  1:615717/43‑1 (MQ=255)                                                                aGACGCCACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTACCGGACGACCAACGt      <  1:3825807/70‑1 (MQ=255)                                                                aGACGCCACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTACCGGACGACCAACGt      <  1:3135216/70‑1 (MQ=255)                                                                   cgcgACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTACCGGACGACCaa         >  1:1296789/1‑64 (MQ=255)                                                                   cgcgACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTACCGGACGACCAACGTCAgg  <  1:1906376/71‑1 (MQ=255)                                                                   cgcgACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTACCGGACGACCAACGTCAgg  <  1:2865173/71‑1 (MQ=255)                                                                      |                                                                    TGGTCTGGCGGTAGATGTGCTTGGCGGCGCGCCGGGCATTTACTCCGCGCGTTATTCCGGTGAAGACGCGACCGATCAAAAGAATCTGCAAAAACTGCTGGAAACAATGAAAGACGTACCGGACGACCAACGTCAGG  >  minE/1983884‑1984020

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑