Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,036,495	T→C	34.1%	D209D (GAT→GAC)	ygiC →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,036,495	0	T	C	34.1%	57.7 / 45.3	82	D209D (GAT→GAC)	ygiC	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (20/34);  new base C (12/16);  total (32/50)
Fisher's exact test for biased strand distribution p-value = 6.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

GGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCGGGTTAGGTGAAAAAGGTCAGTTCACGGATTTACAGGATCAGGTAATTTCCAACCTGTTCAA  >  minE/2036427‑2036563                                                                     |                                                                     ggAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCgat                                                                      <  1:1065673/69‑1 (MQ=255) ggAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCgat                                                                      <  1:1734481/69‑1 (MQ=255) ggAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCgat                                                                      <  1:3310052/69‑1 (MQ=255) ggAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCgat                                                                      <  1:2290628/69‑1 (MQ=255) ggAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCgat                                                                      <  1:3091003/69‑1 (MQ=255) ggAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCgat                                                                      <  1:1536444/69‑1 (MQ=255) ggAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCgat                                                                      <  1:79056/69‑1 (MQ=255)       aTTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATc                                                                <  1:2096735/69‑1 (MQ=255)       aTTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATc                                                                <  1:1415372/69‑1 (MQ=255)         tCAGTTTTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg                                                              <  1:959009/69‑1 (MQ=255)         tCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg                                                              <  1:826034/69‑1 (MQ=255)         tCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg                                                              <  1:347117/69‑1 (MQ=255)         tCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg                                                              <  1:3413827/69‑1 (MQ=255)         tCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg                                                              <  1:1901502/69‑1 (MQ=255)         tCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg                                                              <  1:3325002/69‑1 (MQ=255)         tCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg                                                              <  1:433943/69‑1 (MQ=255)         tCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg                                                              <  1:1782722/69‑1 (MQ=255)         tCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg                                                              <  1:1649324/69‑1 (MQ=255)         tCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg                                                              <  1:1054349/69‑1 (MQ=255)         tCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg                                                              <  1:316250/69‑1 (MQ=255)         tCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg                                                              <  1:816666/69‑1 (MQ=255)         tCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg                                                              <  1:2791057/69‑1 (MQ=255)         tCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg                                                              <  1:3039850/69‑1 (MQ=255)         tCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg                                                              <  1:770255/69‑1 (MQ=255)         tCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg                                                              <  1:280680/69‑1 (MQ=255)         tCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTCCTACTGAGTTCCTCTACATCGATGATATCgg                                                              <  1:2887573/69‑1 (MQ=255)          cAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg                                                              >  1:682187/1‑68 (MQ=255)               tttGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCgat                                                                      >  1:457685/1‑55 (MQ=255)                 tGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCGGGTTAgg                                                        >  1:2881199/1‑67 (MQ=255)                 tGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCGGGTTAgg                                                        >  1:1337629/1‑67 (MQ=255)                 tGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCGGGTTAGGt                                                       >  1:2733605/1‑68 (MQ=255)                     ggACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCGGGTTAGGTGaaaa                                                  <  1:1577372/69‑1 (MQ=255)                     ggACTGCGCAACGGAAGCTGAAATTGCTACTGACTTCCTCTACATCGATGATATCGGGTTAGGTGaaaa                                                  <  1:2072818/69‑1 (MQ=255)                               aCGGAAGCTGAAATCGCTACCGAGTTCCTCTACATCGACGATATCGGGTTat                                                         <  1:375185/52‑2 (MQ=255)                               aCGGAAGCTGAAATCGCTACCGAGTTCCTCTACATCGACGATATCGGGTTAg                                                         <  1:2047522/52‑1 (MQ=255)                               aCGGAAGCTGAAATCGCTACCGAGTTCCTCTACATCGACGATATCGGGTTAg                                                         <  1:2699877/52‑1 (MQ=255)                               aCGGAAGCTGAAATCGCTACCGAGTTCCTCTACATCGACGATATCGGGTTAg                                                         <  1:2588552/52‑1 (MQ=255)                               aCGGAAGCTGAAATCGCTACCGAGTTCCTCTACATCGACGATATCGGGTTAg                                                         <  1:799926/52‑1 (MQ=255)                               aCGGAAGCTGAAATCGCTACCGAGTTCCTCTACATCGACGATATCGGGTTAg                                                         <  1:1197088/52‑1 (MQ=255)                               aCGGAAGCTGAAATCGCTACCGAGTTCCTCTACATCGACGATATCGGGTTAg                                                         <  1:2541532/52‑1 (MQ=255)                               aCGGAAGCTGAAATCGCTACCGAGTTCCTCTACATCGACGATATCGGGTTAg                                                         <  1:1470767/52‑1 (MQ=255)                               aCGGAAGCTGAAATCGCTACCGAGTTCCTCTACATCGACGATATCGGGTTAg                                                         <  1:695014/52‑1 (MQ=255)                               aCGGAAGCTGAAATCGCTACCGAGTTCCTCTACATCGACGATATCGGGTTAg                                                         <  1:237013/52‑1 (MQ=255)                               aCGGAAGCTGAAATCGCTACCGAGTTCCTCTACATCGACGATATCGGGTTAg                                                         <  1:506638/52‑1 (MQ=255)                               aCGGAAGCTGAAATCGCTACCGAGTTCCTCTACATCGACGATATCGGGTTAg                                                         <  1:2274175/52‑1 (MQ=255)                               aCGGAAGCTGAAATCGCTACCGAGTTCCTCTACATCGACGATATCGGGTTAg                                                         <  1:2253314/52‑1 (MQ=255)                               aCGGAAGCTGAAATCGCTACCGAGTTCCTCTACATCGACGATATCGGGTTAg                                                         <  1:3294844/52‑1 (MQ=255)                               aCGGAAGCTGAAATCGCTACCGAGTTCCTCTACATCGACGATATCGGGTTAg                                                         <  1:475340/52‑1 (MQ=255)                               aCGGAAGCTGAAATCGCTACCGAGTTCCTCTACATCGACGATATCGGGTTAg                                                         <  1:2186501/52‑1 (MQ=255)                                cGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg                                                              <  1:791378/46‑1 (MQ=255)                                cGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg                                                              <  1:1980001/46‑1 (MQ=255)                                      cTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCgg                                                              <  1:3263245/40‑1 (MQ=255)                                                    gAGTTCCTCTACATCGATGATATCGGGTTAGGTGAAAAAGGTCAGTTCACGGATTTACAGGATCAGGTa                   >  1:532608/1‑69 (MQ=255)                                                    gAGTTCCTCTACATCGATGATATCGGGTTAGGTGAAAAAGGTCAGTTCACGGATTTACAGGATCAGGTa                   <  1:2494280/69‑1 (MQ=255)                                                       ttCCTCTACATCGATGATATCGGGTTAGGTGAAAAAGGTCAGTTCACGGATTTACAGGATCAGGTaa                  <  1:762789/67‑1 (MQ=255)                                                       ttCCTCTACATCGATGATATCGGGTTAGGTGAAAAAGGTCAGTTCACGGATTTACAGGATCAGGTaa                  <  1:2784338/67‑1 (MQ=255)                                                       ttCCTCTACATCGATGATATCGGGTTAGGTGAAAAAGGTCAGTTCACGGATTTACAGGATCAGGTaa                  <  1:255346/67‑1 (MQ=255)                                                       ttCCTCTACATCGATGATATCGGGTTAGGTGAAAAAGGTCAGTTCACGGATTTACAGGATCAGGTaa                  <  1:1408257/67‑1 (MQ=255)                                                       ttCCTCTACATCGATGATATCGGGTTAGGTGAAAAAGGTCAGTTCACGGATTTACAGGATCAGGTaa                  <  1:2746299/67‑1 (MQ=255)                                                       ttCCTCTACATCGATGATATCGGGTTAGGTGAAAAAGGTCAGTTCACGGATTTACAGGATCAGGTaa                  <  1:2089372/67‑1 (MQ=255)                                                       ttCCTCTACATCGATGATATCGGGTTAGGTGAAAAAGGTCAGTTCACGGATTTACAGGATCAGGTaa                  <  1:1187732/67‑1 (MQ=255)                                                       ttCCTCTACATCGATGATATCGGGTTAGGTGAAAAAGGTCAGTTCACGGATTTACAGGATCAGGTaa                  <  1:1551292/67‑1 (MQ=255)                                                            ctACATCGACGATATCGGGTTAGGTGAAAAAGGTCAGTTc                                        >  1:1837746/1‑40 (MQ=255)                                                            ctACATCGACGATATCGGGTTAGGTGAAAAAGGTCAGTTc                                        >  1:1868789/1‑40 (MQ=255)                                                            ctACATCGACGATATCGGGTTAGGTGAAAAAGGTCAGTTc                                        >  1:1378328/1‑40 (MQ=255)                                                            ctACATCGACGATATCGGGTTAGGTGAAAAAGGTCAGTTc                                        >  1:1724977/1‑40 (MQ=255)                                                            ctACATCGACGATATCGGGTTAGGTGAAAAAGGTCAGTTc                                        >  1:2180394/1‑40 (MQ=255)                                                            ctACATCGACGATATCGGGTTAGGTGAAAAAGGTCAGTTc                                        >  1:2205972/1‑40 (MQ=255)                                                            ctACATCGACGATATCGGGTTAGGTGAAAAAGGTCAGTTc                                        >  1:2266124/1‑40 (MQ=255)                                                            ctACATCGACGATATCGGGTTAGGTGAAAAAGGTCAGTTc                                        >  1:605914/1‑40 (MQ=255)                                                            ctACATCGACGATATCGGGTTAGGTGAAAAAGGTCAGTTc                                        >  1:2839317/1‑40 (MQ=255)                                                            ctACATCGACGATATCGGGTTAGGTGAAAAAGGTCAGTTc                                        >  1:1545761/1‑40 (MQ=255)                                                            ctACATCGACGATATCGGGTTAGGTGAAAAAGGTCAGTTc                                        >  1:1071897/1‑40 (MQ=255)                                                            ctACATCGACGATATCGGGTTAGGTGAAAAAGGTCAGTTc                                        >  1:2568081/1‑40 (MQ=255)                                                             tACATCGATGATATCGGGTTAGGTGAAAAAGGTCAGTTCACGGATTTACAGGATCAGGTAATTTcc             >  1:110601/1‑66 (MQ=255)                                                             tACATCGATGATATCGGGTTAGGTGAAAAAGGTCAGTTCACGGATTTACAGGATCAGGTAATTTCCAAc          >  1:548337/1‑69 (MQ=255)                                                             tACATCGATGATATCGGGTTAGGTGAAAAAGGTCAGTTCACGGATTTACAGGATCAGGTAATTTCCAAc          >  1:1647769/1‑69 (MQ=255)                                                             tACATCGATGATATCGGGTTAGGTGAAAAAGGTCAGTTCACGGATTTACAGGATCAGGTAATTTCCAAc          >  1:618249/1‑69 (MQ=255)                                                             tACATCGATGATATCGGGTTAGGTGAAAAAGGTCAGTTCACGGATTTACAGGATCAGGTAATTTCCAAc          >  1:681965/1‑69 (MQ=255)                                                             tACATCGATGATATCGGGTTAGGTGAAAAAGGTCAGTTCACGGATTTACAGGATCAGGTAATTTCCAAc          >  1:771878/1‑69 (MQ=255)                                                             tACATCGATGATATCGGGTTAGGTGAAAAAGGTCAGTTCACGGATTTACAGGATCAGGTAATTTCCAAc          >  1:3400565/1‑69 (MQ=255)                                                             tACATCGATGATATCGGGTTAGGTGAAAAAGGTCAGTTCACGGATTTACAGGATCAGGTAATTTCCAAc          >  1:2337730/1‑69 (MQ=255)                                                             tACATCGATGATATCGGGTTAGGTGAAAAAGGTCAGTTCACGGATTTACAGGATCAGGTAATTTCCAAc          >  1:2393942/1‑69 (MQ=255)                                                             tACATCGATGATATCGGGTTAGGTGAAAAAGGTCAGTTCACGGATTTACAGGATCAGGTAATTTCCAAc          >  1:2455563/1‑69 (MQ=255)                                                             tACATCGATGATATCGGGTTAGGTGAAAAAGGTCAGTTCACGGATTTACAGGATCAGGTAATTTCCAAc          >  1:246937/1‑69 (MQ=255)                                                             tACATCGATGATATCGGGTTAGGTGAAAAAGGTCAGTTCACGGATTTACAGGATCAGGTAATTTCCAAc          >  1:2826985/1‑69 (MQ=255)                                                             tACATCGATGATATCGGGTTAGGTGAAAAAGGTCAGTTCACGGATTTACAGGATCAGGTAATTTCCAAc          >  1:2814572/1‑69 (MQ=255)                                                             tACATCGATGATATCGGGTTAGGTGAAAAAGGTCAGTTCACGGATTTACAGGATCAGGTAATTTCCAAc          >  1:2558619/1‑69 (MQ=255)                                                             tACATCGATGATATCGGGTTAGGTGAAAAAGGTCAGTTCACGGATTTACAGGATCAGGTAATTTCCAAc          >  1:2645927/1‑69 (MQ=255)                                                                 tCGATGATATCGGGTTAGGTGAAAAAGGTCAGTTCACGGATTTACAGGATCAGGTAATTTCCAACCTg       <  1:1140168/68‑1 (MQ=255)                                                                   gatgatATCGGGTTAGGTGAAAAAGGTCAGTTCACGGATTTACAGGATCAGGTAATTTcc             >  1:38415/1‑60 (MQ=255)                                                                   gatgatATCGGGTTAGGTGAAAAAGGTCAGTTCACGGATTTACAGGATCAGGTAATTTcc             >  1:1474469/1‑60 (MQ=255)                                                                   gatgatATCGGGTTAGGTGAAAAAGGTCAGTTCACGGATTTACAGGATCAGGTAATTTcc             >  1:1896080/1‑60 (MQ=255)                                                                   gatgatATCGGGTTAGGTGAAAAAGGTCAGTTCACGGATTTACAGGATCAGGTAATTTcc             >  1:2240811/1‑60 (MQ=255)                                                                   gatgatATCGGGTTAGGTGAAAAAGGTCAGTTCACGGATTTACAGGATCAGGTAATTTcc             >  1:2812147/1‑60 (MQ=255)                                                                     tgatATCGGGTTAGGTGAAAAAGGTCAGTTCACGGATTTACAGGATCAGGTAATTTCCAACCTGTTCaa  <  1:1467821/69‑1 (MQ=255)                                                                     tgatATCGGGTTAGGTGAAAAAGGTCAGTTCACGGATTTACAGGATCAGGTAATTTCCAACCTGTTCaa  <  1:1291936/69‑1 (MQ=255)                                                                     tgatATCGGGTTAGGTGAAAAAGGTCAGTTCACGGATTTACAGGATCAGGTAATTTCCAACCTGTTCaa  <  1:2139343/69‑1 (MQ=255)                                                                     tgatATCGGGTTAGGTGAAAAAGGTCAGTTCACGGATTTACAGGATCAGGTAATTTCCAACCTGTTCaa  <  1:2298732/69‑1 (MQ=255)                                                                     tgatATCGGGTTAGGTGAAAAAGGTCAGTTCACGGATTTACAGGATCAGGTAATTTCCAACCTGTTCaa  <  1:987200/69‑1 (MQ=255)                                                                     |                                                                     GGAACCATTCAGTATTTGCAGGACTGCGCAACGGAAGCTGAAATTGCTACTGAGTTCCTCTACATCGATGATATCGGGTTAGGTGAAAAAGGTCAGTTCACGGATTTACAGGATCAGGTAATTTCCAACCTGTTCAA  >  minE/2036427‑2036563

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