Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,176,188	G→A	48.4%	A60A (GCC→GCT)	ydjA ←	predicted oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,176,188	0	G	A	48.4%	2.1 / 56.0	66	A60A (GCC→GCT)	ydjA	predicted oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (16/18);  new base A (22/10);  total (38/28)
Fisher's exact test for biased strand distribution p-value = 8.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

ACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAAAAATGCCACGGTTGCATGGACT  >  minE/1176128‑1176249                                                             |                                                              aCGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGaa                                                       >  1:1387626/1‑69 (MQ=255) aCGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGaa                                                       >  1:2852747/1‑69 (MQ=255) aCGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGa                                                        >  1:1754343/1‑68 (MQ=255)      cTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCt                                                  <  1:2555223/69‑1 (MQ=255)      cTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCt                                                  <  1:1946521/69‑1 (MQ=255)           tCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTcgc                                               <  1:959702/67‑1 (MQ=255)            tgATTGCTTTGTTATCACTACCGGCAGCAAT‑‑CCCCCTGCTCCAGTACAGCGCTGAAACGCTCGCGccc                                           <  1:1844612/67‑1 (MQ=255)             gATTGCTTTGTTATCACTACCGGCAGCAATCGCCCCCTGCTCCAGTACAGCGCTGAAACGCTCGCGccc                                           <  1:2452670/69‑1 (MQ=255)             gATTGCTTTGTTATCACTACCGGCAGCAATCGCCCCCTGCTCCAGTACAGCGCTGAAACGCTCGCGccc                                           <  1:1020997/69‑1 (MQ=255)             gATTGCTTTGTTATCACTACCGGCAGCAATCGCCCCCTGCTCCAGTACAGCGCTGAAACGCTCGCGccc                                           <  1:2126429/69‑1 (MQ=255)             gATTGCTTTGTTATCACTACCGGCAGCAATCGCCCCCTGCTCCAGTACAGCGCTGAAACGCTCGCGccc                                           <  1:3210936/69‑1 (MQ=255)             gATTGCTTTGTTATCACTACCGGCAGCAATCGCCCCCTGCTCCAGTACAGCGCTGAAACGCTCGCGccc                                           <  1:911365/69‑1 (MQ=255)             gATTGCTTTGTTATCACTACCGGCAGCAATCGCCCCCTGCTCCAGTACAGCGCTGAAACGCTCGCGccc                                           <  1:1347373/69‑1 (MQ=255)             gATTGCTTTGTTATCACTACCGGCAGCAATCGCCCCCTGCTCCAGTACAGCGCTGAAACGCTCGCGccc                                           <  1:1426500/69‑1 (MQ=255)             gATTGCTTTGTTATCACTACCGGCAGCAATCGCCCCCTGCTCCAGTACAGCGCTGAAACGCTCGCGccc                                           <  1:64313/69‑1 (MQ=255)              attgCTTTGTTATCACTACCGGCAGCAATCGCCCCCTGCTCCAGTACAGCGCTGAAACGCTCGCGccc                                           <  1:1429487/65‑1 (MQ=255)                aGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACg                                                    <  1:1604937/57‑1 (MQ=255)                aGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACg                                                    <  1:1431903/57‑1 (MQ=255)                aGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACg                                                    <  1:1949125/57‑1 (MQ=255)                aGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACg                                                    <  1:19605/57‑1 (MQ=255)                aGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACg                                                    <  1:772449/57‑1 (MQ=255)                aGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGcccttc                                        <  1:2922418/69‑1 (MQ=255)                aGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGcccttc                                        >  1:279822/1‑69 (MQ=255)                aGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGcccttc                                        >  1:2637417/1‑69 (MQ=255)                aGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGcccttc                                        <  1:2635034/69‑1 (MQ=255)                           caCTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATc                              >  1:3394549/1‑68 (MQ=255)                            aCTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATca                             >  1:2483942/1‑68 (MQ=255)                            aCTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATca                             >  1:557378/1‑68 (MQ=255)                            aCTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATca                             >  1:397380/1‑68 (MQ=255)                            aCTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATca                             >  1:47138/1‑68 (MQ=255)                            aCTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATca                             >  1:1499862/1‑68 (MQ=255)                            aCTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATca                             >  1:1630988/1‑68 (MQ=255)                            aCTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATca                             >  1:2219429/1‑68 (MQ=255)                            aCTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATca                             >  1:170304/1‑68 (MQ=255)                                           cGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAAAAATGCCACgg             >  1:63946/1‑69 (MQ=255)                                             cccccTGCTCCAGTACAGCGCTGAACCGCTCGCGcccttcccc                                     >  1:127799/1‑43 (MQ=255)                                             cccccTGCTCCAGTACAGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAGAAATGCCACGGt            >  1:1698093/1‑68 (MQ=255)                                             cccccTGCTCCAGTACAGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAGAAATGCCACGGt            >  1:1360762/1‑68 (MQ=255)                                             cccccTGCTCCAGTACAGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAGAAATGCCACGGt            >  1:809910/1‑68 (MQ=255)                                             cccccTGCTCCAGTACAGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAGAAATGCCACGGt            >  1:136660/1‑68 (MQ=255)                                             cccccTGCTCCAGTACAGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAGAAATGCCACGGt            >  1:3256613/1‑68 (MQ=255)                                             cccccTGCTCCAGTACAGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAGAAATGCCACGGt            >  1:2372947/1‑68 (MQ=255)                                             cccccTGCTCCAGTACAGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAGAAATGCCACGGt            >  1:3055327/1‑68 (MQ=255)                                             cccccTGCTCCAGTACAGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAGAAATGCCACGGt            >  1:3044568/1‑68 (MQ=255)                                             cccccTGCTCCAGTACAGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAGAAATGCCACGGt            >  1:2007523/1‑68 (MQ=255)                                             cccccTGCTCCAGTACAGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAGAAATGCCACGGt            >  1:2601045/1‑68 (MQ=255)                                             cccccTGCTCCAGTACAGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAGAAATGCCACGGt            >  1:2573458/1‑68 (MQ=255)                                             cccccTGCTCCAGTACAGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAGAAATGCCACGGt            >  1:2228874/1‑68 (MQ=255)                                             cccccTGCTCCAGTACAGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCA‑AAAGAAATGCCACGGTTg          >  1:1984937/1‑69 (MQ=255)                                              ccccTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAAAAATGCCACGGTTg          <  1:3103224/69‑1 (MQ=255)                                              ccccTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAAAAATGCCACGGTTg          <  1:3209356/69‑1 (MQ=255)                                              ccccTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAAAAATGCCACGGTTg          <  1:1688072/69‑1 (MQ=255)                                              ccccTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAAAAATGCCACGGTTg          <  1:165196/69‑1 (MQ=255)                                              ccccTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAAAAATGCCACGGTTg          <  1:2610976/69‑1 (MQ=255)                                               cccTGTTCCAGTACGGCGCTGAAACGCTCGCGGCCTTCCCCTTCAATCACAAAAAAATGCCACGGTTg          <  1:1696742/68‑1 (MQ=255)                                                  tgctCCAGTACAGCGCTGAAACGCTCTCGCCCTTCCCCTTCAATCACAAAGAAATGCCACGGTTg          >  1:2289450/4‑65 (MQ=255)                                                  tgctCCAGTACAGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAGAAATGCCACGGTTg          >  1:887348/4‑65 (MQ=255)                                                  tgctCCAGTACAGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAGAAATGCCACGGTTg          >  1:818123/4‑65 (MQ=255)                                                  tgctCCAGTACAGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAGAAATGCCACGGTTg          >  1:786586/4‑65 (MQ=255)                                                  tgctCCAGTACAGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAGAAATGCCACGGTTg          >  1:2483704/4‑65 (MQ=255)                                                  tgctCCAGTACAGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAGAAATGCCACGGTTg          >  1:2872903/4‑65 (MQ=255)                                                  tgctCCAGTACAGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAGAAATGCCACGGTTg          >  1:3327862/4‑65 (MQ=255)                                                  tgctCCAGTACAGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAGAAATGCCACGGTTg          >  1:1822777/4‑65 (MQ=255)                                                  tGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAAAAATGCCACGGTTGCATg      >  1:2465879/1‑69 (MQ=255)                                                      ccAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAAAAATGCCACGGTTGCATGGACt  <  1:723115/69‑1 (MQ=255)                                                      ccAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAAAAATGCCACGGTTGCATGGACt  <  1:623000/69‑1 (MQ=255)                                                             |                                                              ACGGGCTTTGTCGATAGCTTTGTCATCACTACCGGCAGCAATCGCCCCCTGTTCCAGTACGGCGCTGAAACGCTCGCGCCCTTCCCCTTCAATCACAAAAAAATGCCACGGTTGCATGGACT  >  minE/1176128‑1176249

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