Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	115,773	T→A	41.0%	A76A (GCA→GCT)	yacH ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	115,773	0	T	A	41.0%	26.7 / 57.9	61	A76A (GCA→GCT)	yacH	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (13/23);  new base A (4/21);  total (17/44)
Fisher's exact test for biased strand distribution p-value = 1.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.49e-01

ATCAGCACCTGCGAAAGTAGGGCGTCGGGATACAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGGTGTTACGGGAGGTGCGGAAACCGTTGCAGGAAA  >  minE/115708‑115838                                                                  |                                                                  aTCAGCACCTGCGAAAGTAGGGCGTCGGGATACAGCGCGACGGGCGCGACCCATTGATCTATTTGTGc                                                                 <  1:2292634/68‑1 (MQ=255) aTCAGCACCTGCGAAAGTAGGGCGTCGGGATACAGCGCGACGGGCGCGACCCATTGATCTATTTGTGc                                                                 <  1:1927186/68‑1 (MQ=255) aTCAGCACCTGCGAAAGTAGGGCGTCGGGATACAGCGCGACGGGCGCGACCCATTGATCTATTTGTGc                                                                 <  1:4942/68‑1 (MQ=255)           gCGAAAGTAGGGCGTCGGGGTATATCTCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                                                      <  1:2668343/69‑1 (MQ=255)           gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCATTGCTGAATg                                                      <  1:3387176/69‑1 (MQ=255)           gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                                                      <  1:7026/69‑1 (MQ=255)           gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                                                      <  1:1945015/69‑1 (MQ=255)           gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                                                      <  1:1974296/69‑1 (MQ=255)           gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                                                      <  1:1840732/69‑1 (MQ=255)           gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                                                      <  1:1839019/69‑1 (MQ=255)           gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                                                      <  1:1766218/69‑1 (MQ=255)           gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                                                      <  1:2169835/69‑1 (MQ=255)           gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                                                      <  1:2178574/69‑1 (MQ=255)           gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                                                      <  1:2187023/69‑1 (MQ=255)           gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                                                      <  1:2282317/69‑1 (MQ=255)           gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                                                      <  1:1539974/69‑1 (MQ=255)           gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                                                      <  1:2819282/69‑1 (MQ=255)           gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                                                      <  1:288219/69‑1 (MQ=255)           gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                                                      <  1:3051442/69‑1 (MQ=255)           gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                                                      <  1:3081135/69‑1 (MQ=255)           gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                                                      <  1:3182131/69‑1 (MQ=255)           gCGAAAGTAGGGCGTCGGGGTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                                                      <  1:654405/69‑1 (MQ=255)             gAAAGTAGGGCGTCGGGGTATAGCGCGAGAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATg                                                      <  1:973161/67‑1 (MQ=255)                  tAGGGCGTCGGGATACAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGAttt                                               >  1:225097/1‑69 (MQ=255)                   aGGGCGTCGGGATACAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTg                                              >  1:2525818/1‑69 (MQ=255)                           ggggTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATGCGGATTTGACGAcggc                                      >  1:2646863/1‑69 (MQ=255)                           ggggTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATGCGGATTTGACGAcggc                                      >  1:2349546/1‑69 (MQ=255)                           ggggTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGTGCTGAATGCGGATTTGACGAcggc                                      >  1:1752683/1‑69 (MQ=255)                           ggggTATAGCGCGACAGGTGCAACCCATTGATCTATTTGAGCAGAGCTGAATGCGGATTTGACGAcggc                                      >  1:1924490/1‑69 (MQ=255)                           gggATACAGCGCGACGGGCGCGACCCATTGATCTATTtgtg                                                                  <  1:1161205/41‑1 (MQ=255)                           gggATACAGCGCGACGGGCGCGACCCATTGATCTATTtgtg                                                                  <  1:911778/41‑1 (MQ=255)                           gggATACAGCGCGACGGGCGCGACCCATTGATCTATTtgtg                                                                  <  1:2646368/41‑1 (MQ=255)                           gggATACAGCGCGA‑GGGCGCGACCCATTGATCTATTtgtg                                                                  <  1:3398690/40‑1 (MQ=39)                             gATACAGCGCGACGGGCGCGACCCATTGATCTATTtgtg                                                                  <  1:2672119/39‑1 (MQ=255)                                aCAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGgtg                                 <  1:144692/69‑1 (MQ=255)                                aCAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGgtg                                 <  1:14412/69‑1 (MQ=255)                                aCAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGgtg                                 <  1:1878231/69‑1 (MQ=255)                                aCAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGgtg                                 <  1:1437560/69‑1 (MQ=255)                                aCAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGgtg                                 <  1:2446344/69‑1 (MQ=255)                                aCAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGgtg                                 <  1:1349907/69‑1 (MQ=255)                                aCAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGgtg                                 <  1:2362157/69‑1 (MQ=255)                                aCAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGgtg                                 <  1:551754/69‑1 (MQ=255)                                aCAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGgtg                                 <  1:1914786/69‑1 (MQ=255)                                 cAGCGCGACGGGCGCGACCCATTGATCTATTtgtg                                                                  <  1:393848/35‑1 (MQ=255)                                 cAGCGCGACGGGCGCGACCCATTGATCTATTtgtg                                                                  <  1:3407067/35‑1 (MQ=255)                                 cAGCGCGACGGGCGCGACCCATTGATCTATTtgtg                                                                  <  1:1219932/35‑1 (MQ=255)                                 cAGCGCGACGGGCGCGACCCATTGATCTATTtgtg                                                                  <  1:1208716/35‑1 (MQ=255)                                 cAGCGCGACGGGCGCGACCCATTGATCTATTtgtg                                                                  <  1:1193323/35‑1 (MQ=255)                                 cAGCGCGACGGGCGCGACCCATTGATCTATTtgtg                                                                  <  1:1902935/35‑1 (MQ=255)                                 cAGCGCGACGGGCGCGACCCATTGATCTATTtgtg                                                                  <  1:2086310/35‑1 (MQ=255)                                 cAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGgtg                                 <  1:1400908/68‑1 (MQ=255)                                      cgACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGGTGTTACgg                           <  1:2037405/69‑1 (MQ=255)                                      cgACGGGCGAGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGgtg                                 <  1:2496546/63‑1 (MQ=255)                                                gACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGGTGTTACGGGAGGTGCGGa                 <  1:745353/69‑1 (MQ=255)                                                gACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGGTGTTACGGGAGGTGCGGa                 <  1:2079756/69‑1 (MQ=255)                                                gACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGGTGTTACGGGAGGTGCGGa                 <  1:2808791/69‑1 (MQ=255)                                                gACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGGTGTTACGGGAGGTGCGGa                 <  1:1497813/69‑1 (MQ=255)                                                gACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGGTGTTACGGGAGGTGCGGa                 <  1:251715/69‑1 (MQ=255)                                                      ttGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGgtg                                 <  1:154242/47‑1 (MQ=255)                                                        gATGTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGGTGTTACGGGAGGTGCGGAAACCGt           >  1:759228/1‑67 (MQ=255)                                                        gATCTATTTGTGCAGTGCTGAATTCGGATTTGACGACGGCGGGTGTTACGGGAGGTGCGGAAACCGt           >  1:2023997/1‑67 (MQ=255)                                                        gATCTATTTGTGCAGTGCTGAATTCGGATTTGACGACGGCGGGTGTTACGGGAGGTGCGGAAACCGt           >  1:2027833/1‑67 (MQ=255)                                                        gATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGGTGTTACGGGAgg                       <  1:3378385/55‑1 (MQ=255)                                                        gATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGGTGTTACGGGAgg                       <  1:2366743/55‑1 (MQ=255)                                                        gATCTAATTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGGTGTTACGGGAGGTGCGGAAACCGt           >  1:434905/1‑67 (MQ=255)                                                            tATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGGTGTTACGGGAGGTGCGGAAACCGTTGCa       >  1:2205049/1‑67 (MQ=255)                                                            tATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGGTGTTACGGGAGGTGCGGAAACCGTTGCAg      >  1:2772564/1‑68 (MQ=255)                                                            tATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGGTGTTACGGGAGGTGCGGAAACCGTTGCAg      >  1:2232161/1‑68 (MQ=255)                                                            tATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGGTGTTACGGGAGGTGCGGAAACCGTTGCAg      >  1:1530052/1‑68 (MQ=255)                                                            tATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGGTGTTACGGGAGGTGCGGAAACCGTTGCAg      >  1:3042950/1‑68 (MQ=255)                                                            tATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGGTGTTACGGGAGGTGCGGAAACCGTTGCAg      >  1:3038028/1‑68 (MQ=255)                                                            tATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGGTGTTACGGGAGGTGCGGAAACCGTTGCAg      >  1:3018932/1‑68 (MQ=255)                                                            tATTTGAGCAGTGCTGAATGCGGATTTGACGACGGCGGgtg                                 <  1:2907041/41‑1 (MQ=255)                                                                tgtgCAGTGCTGAATGCGGATTTGACGACGGCGGGTGTTACGGGAGGTGCGGAAACCGTTGCAGGaaa  >  1:3293121/1‑68 (MQ=255)                                                                tgtgCAGTGCTGAATGCGGATTTGACGACGGCGGGTGTTACGGGAGGTGCGGAAACCGTTGCAGGaaa  >  1:1423216/1‑68 (MQ=255)                                                                tgtgCAGTGCTGAATGCGGATTTGACGACGGCGGGTGTTACGGGAGGTGCGGAAACCGTTGCAGGaaa  >  1:105132/1‑68 (MQ=255)                                                                tgtgCAGTGCTGAATGCGGATTTGACGACGGCGGGTGTTACGGGAGGTGCGGAAACCGTTGCAGGaaa  >  1:558708/1‑68 (MQ=255)                                                                tgtgCAGTGCTGAATGCGGATTTGACGACGGCGGGTGTTACGGGAGGTGCGGAAACCGTTGCAGGaaa  >  1:1552640/1‑68 (MQ=255)                                                                tgtgCAGTGCTGAATGCGGATTTGACGACGGCGGGTGTTACGGGAGGTGCGGAAACCGTTGCAGGaaa  >  1:1607419/1‑68 (MQ=255)                                                                tgtgCAGTGCTGAATGCGGATTTGACGACGGCGGGTGTTACGGGAGGTGCGGAAACCGTTGCAGGaaa  >  1:2056085/1‑68 (MQ=255)                                                                tgtgCAGTGCTGAATGCGGATTTGACGACGGCGGGTGTTACGGGAGGTGCGGAAACCGTTGCAGGaaa  >  1:1917696/1‑68 (MQ=255)                                                                  |                                                                  ATCAGCACCTGCGAAAGTAGGGCGTCGGGATACAGCGCGACGGGCGCGACCCATTGATCTATTTGTGCAGTGCTGAATGCGGATTTGACGACGGCGGGTGTTACGGGAGGTGCGGAAACCGTTGCAGGAAA  >  minE/115708‑115838

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