Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,206,243	G→A	29.7%	A356A (GCC→GCT)	tldD ←	predicted peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,206,243	0	G	A	29.7%	64.8 / 28.0	67	A356A (GCC→GCT)	tldD	predicted peptidase
Reads supporting (aligned to +/- strand):  ref base G (19/28);  new base A (8/12);  total (27/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTA  >  minE/2206192‑2206306                                                    |                                                                cGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCg                                                >  1:2517363/1‑69 (MQ=255) cGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCg                                                >  1:40986/1‑69 (MQ=255)      tACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGt                                           >  1:3163526/1‑69 (MQ=255)               gcaTATAGGTGTTGGTCATACGCGGCATGGGCAGATGAGCGTAGGATTCACGGCGACCGTTGCCAGTCg                                  <  1:2121541/69‑1 (MQ=255)               gcaTATAGGTGTTGGTCATACGCGGCATGGGCAGATGAGCGTAGGATTCACGGCGACCGTTGCCAGTCg                                  <  1:98271/69‑1 (MQ=255)               gcaTATAGGTGTTGGTCATACGCGGCATGGGCAGATGAGCGTAGGATTCACGGCGACCGTTGCCAGTCg                                  <  1:2217372/69‑1 (MQ=255)               gcaTATAGGTGTTGGTCATACGCGGCATGGGCAGATGAGCGTAGGATTCACGGCGACCGTTGCCAGTCg                                  <  1:671520/69‑1 (MQ=255)               gcaTATAGGTGTTGGTCATACGCGGCATGGGCAGATGAGCGTAGGATTCACGGCGACCGTTGCCAGTCg                                  <  1:2268518/69‑1 (MQ=255)               gcaTATAGGTGTTGGTCATACGCGGCATGGGCAGATGAGCGTAGGATTCACGGCGACCGTTGCCAGTCg                                  <  1:1272583/69‑1 (MQ=255)               gcaTATAGGTGTTGGTCATACGCGGCATGGGCAGATGAGCGTAGGATTCACGGCGACCGTTGCCAGTCg                                  <  1:2325959/69‑1 (MQ=255)               gcaTATAGGTGTTGGTCATACGCGGCATGGGCAGATGAGCGTAGGATTCACGGCGACCGTTGCCAGTCg                                  <  1:3345063/69‑1 (MQ=255)               gcaTATAGGTGTTGGTCATACGCGGCATGGGCAGATGAGCGTAGGATTCACGGCGACCGTTGCCAGTCg                                  <  1:158487/69‑1 (MQ=255)               gcaTATAGGTGTTGGTCATACGCGGCATGGGCAGATGAGCGTAGGATTCACGGCGACCGTTGCCAGTCg                                  <  1:3189700/69‑1 (MQ=255)                    taGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTc                             >  1:1014578/1‑69 (MQ=255)                    taGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTc                             >  1:2062239/1‑69 (MQ=255)                    taGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTc                             >  1:2262656/1‑69 (MQ=255)                    taGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTc                             >  1:1390504/1‑69 (MQ=255)                    taGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTc                             >  1:47965/1‑69 (MQ=255)                      ggTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCAt                           >  1:1506784/1‑69 (MQ=255)                      ggTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCAt                           >  1:1208116/1‑69 (MQ=255)                       gtgtTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATc                          >  1:2222028/1‑69 (MQ=255)                       gtgtTGGTCATACGCGGCATGGGCAGATGAGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCAt                           >  1:920454/1‑68 (MQ=255)                       gtgtTGGTCATACGCGGCATGGGCAGATGAGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATc                          >  1:265837/1‑69 (MQ=255)                       gtgtTGGTCATACGCGGCATGGGCAGATGAGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATc                          >  1:50391/1‑69 (MQ=255)                       gtgtTGGTCATACGCGGCATGGGCAGATGAGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATc                          >  1:996903/1‑69 (MQ=255)                       gtgtTGGTCATACGCGGCATGGGCAGATGAGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATc                          >  1:1078891/1‑69 (MQ=255)                       gtgtTGGTCATACGCGGCATGGGCAGATGAGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATc                          >  1:1113896/1‑69 (MQ=255)                       gtgtTGGTCATACGCGGCATGGGCAGATGAGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATc                          >  1:146762/1‑69 (MQ=255)                       gtgtTGGTCATACGCGGCATGGGCAGATGAGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATc                          >  1:176923/1‑69 (MQ=255)                         gtTGGTCATACGCGGCATGGGCAGATGAGCGTAGGATTCACGGCGACCGTTGCCAGTCg                                  <  1:868724/59‑1 (MQ=255)                            ggTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGt                                    >  1:2976138/1‑54 (MQ=255)                            ggTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCa                      >  1:2216854/1‑68 (MQ=255)                            ggTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCa                      >  1:1434250/1‑68 (MQ=255)                            ggTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCa                      >  1:350606/1‑68 (MQ=255)                            ggTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCa                      >  1:1863347/1‑68 (MQ=255)                             gTCATACGCGGCATGGGCAGATGAGCGTAGGATTCACGGCGACCGTTGACAGTCg                                  <  1:4024/55‑1 (MQ=255)                               cATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCa                   <  1:1170885/68‑1 (MQ=255)                               cATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCa                   <  1:670739/68‑1 (MQ=255)                                   cgcgGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAAcgc              >  1:3020837/1‑69 (MQ=255)                                   cgcgGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAAcgc              >  1:3147366/1‑69 (MQ=255)                                   cgcgGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAAcgc              >  1:3092887/1‑69 (MQ=255)                                     cgGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAAcgcg             <  1:1281421/68‑1 (MQ=255)                                     cgGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAAcgcg             <  1:3027885/68‑1 (MQ=255)                                      ggCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAAcgcg             <  1:3170264/67‑1 (MQ=255)                                         aTGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCTCGTTg        <  1:1156255/69‑1 (MQ=255)                                         aTGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTg        <  1:648619/69‑1 (MQ=255)                                         aTGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTg        <  1:2676997/69‑1 (MQ=255)                                              cAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGtt    <  1:796525/68‑1 (MQ=255)                                              cAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGtt    <  1:869795/68‑1 (MQ=255)                                              cAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGtt    <  1:903694/68‑1 (MQ=255)                                              cAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGtt    <  1:907701/68‑1 (MQ=255)                                              cAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGtt    <  1:479789/68‑1 (MQ=255)                                              cAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGtt    <  1:1114782/68‑1 (MQ=255)                                              cAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGtt    <  1:1932434/68‑1 (MQ=255)                                              cAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGtt    <  1:1864310/68‑1 (MQ=255)                                              cAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGtt    <  1:1530142/68‑1 (MQ=255)                                              cAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGtt    <  1:1577318/68‑1 (MQ=255)                                              cAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGtt    <  1:3092133/68‑1 (MQ=255)                                              cAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGtt    <  1:2437329/68‑1 (MQ=255)                                              cAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGtt    <  1:1658758/68‑1 (MQ=255)                                               aGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGtt    <  1:483139/67‑1 (MQ=255)                                               aGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTa  <  1:3391732/69‑1 (MQ=255)                                               aGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTa  <  1:1965654/69‑1 (MQ=255)                                               aGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTa  <  1:2743691/69‑1 (MQ=255)                                               aGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTa  <  1:1151488/69‑1 (MQ=255)                                               aGATGGGCATAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTa  <  1:1655302/69‑1 (MQ=255)                                               aGATGGGCATAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATAAAACGCGCGTTGAGTTTa  <  1:2035714/69‑1 (MQ=255)                                                    |                                                                CGATTTACCCGGCAGCATATAGGTGTTGGTCATACGCGGCATGGGCAGATGGGCGTAGGATTCACGGCGACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTA  >  minE/2206192‑2206306

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