Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,121,861	A→G	17.0%	A122A (GCA→GCG)	yraP →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,121,861	0	A	G	17.0%	99.4 / 10.8	64	A122A (GCA→GCG)	yraP	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (16/37);  new base G (8/3);  total (24/40)
Fisher's exact test for biased strand distribution p-value = 1.44e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GGCGTAGACGGTGCCAACGAAGTGTATAACGAGATTCGTCAGGGCCAGCCGATTGGTCTGGGCGAAGCATCTAACGATACGTGGATCACCACCAAAGTGCGTTCGCAGCTCTTAACCAGCGACCTGGT  >  minE/2121793‑2121920                                                                     |                                                            ggCGTAGACGGTGCCAACGAAGTGTATAACGAGATTCGTCAGGGCCAGCCGATTGGTCTGGGCGAAGCa                                                             >  1:1530438/1‑69 (MQ=255)   cGTAGACGGTGCCAACGAAGTGTATAACGAGATTCGTCAGGGCCAGCCGATTGGTCTGGGCGAAGCGTc                                                           <  1:2069514/69‑1 (MQ=255)   cGTAGACGGTGCCAACGAAGTGTATAACGAGATTCGTCAGGGCCAGCCGATTGGTCTGGGCGAAGCGTc                                                           <  1:3123429/69‑1 (MQ=255)   cGTAGACGGTGCCAACGAAGTGTATAACGAGATTCGTCAGGGCCAGCCGATTGGTCTGGGCGAAGCGTc                                                           <  1:2716242/69‑1 (MQ=255)          ggTGCCAACGAAGTGTATAACGAGATTCGTCAGGGCCAGCCGATTGGTCTGGGCGAAGCATCTAACGAt                                                    >  1:1745397/1‑69 (MQ=255)             gCCAACTAAGTGTATAACGAGATTCGTCAGGGCCAGCCGATTGGTCTTGGCGAAGCATCTAACGATACg                                                 <  1:623017/69‑1 (MQ=255)             gCCAACGAAGTGTATAACGAGATTCGTCAGGGCCAGCCGATTGGTCTGGGCGAAGCATCTAAc                                                       <  1:1933884/63‑1 (MQ=255)             gCCAACGAAGTGTATAACGAGATTCGTCAGGGCCAGCCGATTGGTCTGGGCGAAGCATCTAAc                                                       <  1:2306142/63‑1 (MQ=255)             gCCAACGAAGTGTATAACGAGATTCGTCAGGGCCAGCCGATTGGTCTGGGCGAAGCATCTAAc                                                       <  1:375641/63‑1 (MQ=255)             gCCAACGAAGTGTATAACGAGATTCGTCAGGGCCAGCCGATTGGTCTGGGCGAAGCATCTAAc                                                       <  1:3179294/63‑1 (MQ=255)             gCCAACGAAGTGTATAACGAGATTCGTCAGGGCCAGCCGATTGGTCTGGGCGAAGCATCTAACGATACg                                                 <  1:1774263/69‑1 (MQ=255)             gCCAACGAAGTGTATAACGAGATTCGTCAGGGCCAGCCGATTGGTCTGGGCGAAGCATCTAACGATACg                                                 <  1:1013772/69‑1 (MQ=255)                aaCGAAGTGTATAACGAGATTCGTCAGGGCCAGCCGATTGGTCTGGGCGAAGCATCTAACGATACGTgg                                              <  1:1056208/69‑1 (MQ=255)                             aCGAGATTCGTCAGGGCCAGCCGATTGGTCTGGGCGAAGCATCTAACGATACGTGGATCACCACCaaa                                  >  1:2885107/1‑68 (MQ=255)                             aCGAGATTCGTCAGGGCCAGCCGATTGGTCTGGGCGAAGCATCTAACGATACGTGGATCACCACCAAAg                                 >  1:3248337/1‑69 (MQ=255)                             aCGAGATTCGTCAGGGCCAGCCGATTGGTCTGGGCGAAGCATCTAACGATACGTGGATCACCACCAAAg                                 >  1:2733398/1‑69 (MQ=255)                             aCGAGATTCGTCAGGGCCAGCCGATTGGTCTGGGCGAAGCATCTAACGATACGTGGATCACCACCAAAg                                 >  1:2935839/1‑69 (MQ=255)                             aCGAGATTCGTCAGGGCCAGCCGATTGGTCTGGGCGAAGCATCTAACGATACGTGGATCACCACCAAAg                                 >  1:2899211/1‑69 (MQ=255)                             aCGAGATTCGTCAGGGCCAGCCGATTGGTCTGGGCGAAGCATCTAACGATACGTGGATCACCACCAAAg                                 >  1:762015/1‑69 (MQ=255)                                    tCGTCAGGGCCAGCCGATTGGTCTGGGCGAAGCATCTAACGATACGTGGATCACCACCaa                                   <  1:2347952/60‑1 (MQ=255)                                    tCGTCAGGGCCAGCCGATTGGTCTGGGCGAAGCATCTAACGATACGTGGATCACCACCaa                                   <  1:2409742/60‑1 (MQ=255)                                    tCGTCAGGGCCAGCCGATTGGTCTGGGCGAAGCATCTAACGATACGTGGATCACCACCaa                                   <  1:2009898/60‑1 (MQ=255)                                    tCGTCAGGGCCAGCCGATTGGTCTGGGCGAAGCATCTAACGATACGTGGATCACCACCaa                                   <  1:3059145/60‑1 (MQ=255)                                    tCGTCAGGGCCAGCCGATTGGTCTGGGCGAAGCATCTAACGATACGTGGATCACCACCaa                                   <  1:331792/60‑1 (MQ=255)                                    tCGTCAGGGCCAGCCGATTGGTCTGGGCGAAGCATCTAACGATACGTGGATCACCACCaa                                   <  1:680300/60‑1 (MQ=255)                                     cGTCAGGGCCAGCCGATTGGTCTGGGCGAAGCATCTAACGATACGTGGATCACCACCaa                                   <  1:2465670/59‑1 (MQ=255)                                        cAGGGCCAGCCGATTGGTCTGGGCGAAGCATCTAAc                                                       <  1:3128623/36‑1 (MQ=255)                                        cAGGGCCAGCCGATTGGTCTGGGCGAAGCATCTAAc                                                       <  1:2542054/36‑1 (MQ=255)                                        cAGGGCCAGCCGATTGGTCTGGGCGAAGCATCTAAc                                                       <  1:453107/36‑1 (MQ=255)                                        cAGGGCCAGCCGATTGGTCTGGGCGAAGCATCTAAc                                                       <  1:247171/36‑1 (MQ=255)                                        cAGGGCCAGCCGATTGGTCTGGGCGAAGCATCTAAc                                                       <  1:1236759/36‑1 (MQ=255)                                        cAGGGCCAGCCGATTGGTCTGGGCGAAGCATCTAAc                                                       <  1:1137832/36‑1 (MQ=255)                                        cAGGGCCAGCCGATTGGTCTGGGCGAAGCATCTAAc                                                       <  1:1083675/36‑1 (MQ=255)                                        cAGGGCCAGCCGATTGGTCTGGGCGAAGCATCTAAc                                                       <  1:2108363/36‑1 (MQ=255)                                        cAGGGCCAGCCGATTGGTCTGGGCGAAGCATCTAAc                                                       <  1:1760570/36‑1 (MQ=255)                                        cAGGGCCAGCCGATTGGTCTGGGCGAAGCATCTAAc                                                       <  1:1841672/36‑1 (MQ=255)                                        cAGGGCCAGCCGATTGGGCTGGGCGAAGCATCTAAc                                                       <  1:2229204/36‑1 (MQ=255)                                            gCCAGCCGATTGGTCTGGGCGAAGCATCTAACGATACGTGGATCACCACCAAAGTGCGTTCGCAGctc                   >  1:377039/1‑68 (MQ=255)                                            gCCAGCCGATTGGTCTGGGCGAAGCATCTAACGATACGTGGATCACCACCAAAGTGCGTTCGCAGctc                   >  1:387615/1‑68 (MQ=255)                                            gCCAGCCGATTGGTCTGGGCGAAGCATCTAACGATACGTGGATCACCACCAAAGTGCGTTCGCAGctc                   >  1:3135110/1‑68 (MQ=255)                                            gCCAGCCGATTGGTCTGGGCGAAGCATCTAACGATACGTGGATCACCACCAAAGTGCGTTCGCAGctc                   >  1:1871324/1‑68 (MQ=255)                                            gCCAGCCGATTGGTCTGGGCGAAGCATCTAACGATACGTGGATCACCACCAAAGTGCGTTCGCAGctc                   >  1:2048578/1‑68 (MQ=255)                                            gCCAGCCGATTGGTCTGGGCGAAGCATCTAACGATACGTGGATCACCACCAAAGTGCGTTCGCAGctc                   >  1:1217734/1‑68 (MQ=255)                                            gCCAGCCGATTGGTCTGGGCGAAGCATCTAACGATACGTGGATCACCACCAAAGTGCGTTCGCAGc                     >  1:2569017/1‑66 (MQ=255)                                              cAGCCGATTGGTCTGGGCGAAGCATCTAACGATACGTGGATCACCACCAAAGTGCGTTCGCAGCTCTTa                <  1:450943/69‑1 (MQ=255)                                              cAGCCGATTGGTCTGGGCGAAGCATCTAACGATACGTGGATCACCACCAAAGTGCGTTCGCAGCTCTTa                <  1:220799/69‑1 (MQ=255)                                              cAGCCGATTGGTCTGGGCGAAGCATCTAACGATACGTGGATCACCACCAAAGTGCGTTCGCAGCTCTTa                <  1:2316400/69‑1 (MQ=255)                                              cAGCCGATTGGTCTGGGCGAAGCATCTAACGATACGTGGATCACCACCAAAGTGCGTTCGCAGCTCTTa                <  1:3320044/69‑1 (MQ=255)                                              cAGCCGATTGGTCTGGGCGAAGCATCTAACGATACGTGGATCACCACCAAAGTGCGTTCGCAGCTCTTa                <  1:385340/69‑1 (MQ=255)                                              cAGCCGATTGGTCTGGGCGAAGCATCTAACGATACGTGGATCACCACCAAAGTGCGTTCGCAGCTCTTa                <  1:375879/69‑1 (MQ=255)                                                     ttGGTCTGGGCGAAGCATCTAACGATACGTGGATCACCACCAAAGTGCGTTCGCAGCTCTTAACCAGCg         <  1:172369/69‑1 (MQ=255)                                                     ttGGTCTGGGCGAAGCATCTAACGATACGTGGATCACCACCAAAGTGCGTTCGCAGCTCTTAACCAGCg         <  1:2645393/69‑1 (MQ=255)                                                     ttGGTCTGGGCGAAGCATCTAACGATACGTGGATCACCACCAAAGTGCGTTCGCAGCTCTTAACCAGCg         <  1:3205126/69‑1 (MQ=255)                                                     ttGGTCTGGGCGAAGCATCTAACGATACGTGGATCACCACCAAAGTGCGTTCGCAGCTCTTAACCAGCg         <  1:2270257/69‑1 (MQ=255)                                                      tGGTCTGGGCGAAGCATCTAACGATACGTGGATCACCACCAAAGTGCGTTCGCAGCTCTTAACCAGCg         <  1:490565/68‑1 (MQ=255)                                                             ggCGAAGCATCTAACGATACGTGGATCACCACCAAAGTGCGTTCGCAGCTCTTAACCAGCGACCTGGt  >  1:563449/1‑68 (MQ=255)                                                             ggCGAAGCATCTAACGATACGTGGATCACCACCAAAGTGCGTTCGCAGCTCTTAACCAGCGACCTGGt  >  1:1336908/1‑68 (MQ=255)                                                                gAAGCGTCTAACGATACGTGGATCACCACCAAAGTACGTTCGCAGCTGTTAAcc             >  1:1083042/1‑54 (MQ=255)                                                                gAAGCGTCTAACGATACGTGGATCACCACCAAAGTACGTTCGCAGCTGTTAAcc             >  1:80461/1‑54 (MQ=255)                                                                gAAGCGTCTAACGATACGTGGATCACCACCAAAGTACGTTCGCAGCTGTTAAcc             >  1:1484680/1‑54 (MQ=255)                                                                gAAGCGTCTAACGATACGTGGATCACCACCAAAGTACGTTCGCAGCTGTTAAcc             >  1:359526/1‑54 (MQ=255)                                                                gAAGCGTCTAACGATACGTGGATCACCACCAAAGTACGTTCGCAGCTGTTAAcc             >  1:3361011/1‑54 (MQ=255)                                                                gAAGCGTCTAACGATACGTGGATCACCACCAAAGTACGTTCGCAGCTGTTAAcc             >  1:1709332/1‑54 (MQ=255)                                                                gAAGCGTCTAACGATACGTGGATCACCACCAAAGTACGTTCGCAGCTGTTAAcc             >  1:1816079/1‑54 (MQ=255)                                                                gAAGCGTCTAACGATACGTGGATCACCACCAAAGTACGTTCGCAGCTGTTAAc              >  1:2468647/1‑53 (MQ=255)                                                                     |                                                            GGCGTAGACGGTGCCAACGAAGTGTATAACGAGATTCGTCAGGGCCAGCCGATTGGTCTGGGCGAAGCATCTAACGATACGTGGATCACCACCAAAGTGCGTTCGCAGCTCTTAACCAGCGACCTGGT  >  minE/2121793‑2121920

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