Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	572,219	A→G	45.7%	I176I (ATT→ATC)	aqpZ ←	aquaporin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	572,219	0	A	G	45.7%	‑0.6 / 27.0	35	I176I (ATT→ATC)	aqpZ	aquaporin
Reads supporting (aligned to +/- strand):  ref base A (4/15);  new base G (13/3);  total (17/18)
Fisher's exact test for biased strand distribution p-value = 6.11e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGGAAGATAGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTG  >  minE/572155‑572285                                                                 |                                                                   ccTGGAAGATAGCAACCGCGGTGCTGCGCGCCTGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTa                                                                <  1:1079448/69‑1 (MQ=255) ccTGGAAGATAGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTa                                                                <  1:3280701/69‑1 (MQ=255) ccTGGAAGATAGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTa                                                                <  1:2343845/69‑1 (MQ=255)               aCCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTa                                                                <  1:523685/55‑1 (MQ=255)               aCCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGtt                                                  <  1:1134956/69‑1 (MQ=255)               aCCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGtt                                                  <  1:657155/69‑1 (MQ=255)               aCCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGtt                                                  <  1:2321927/69‑1 (MQ=255)               aCCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGtt                                                  <  1:2249836/69‑1 (MQ=255)               aCCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGtt                                                  <  1:1749768/69‑1 (MQ=255)               aCCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGtt                                                  <  1:1565933/69‑1 (MQ=255)               aCCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGtt                                                  <  1:155337/69‑1 (MQ=255)                ccGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTAACCGGAATACTAATTAAGTGAATCAGGGtt                                                  <  1:1932223/68‑1 (MQ=255)                                   gTTAACGGAAGTGTTAGTCACCGGAATACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAAt                              <  1:1157276/69‑1 (MQ=255)                                   gTTAACGGAAGTGTTAGTCACCGGAATACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAAt                              <  1:1428513/69‑1 (MQ=255)                                   gTTAACGGAAGTGTTAGTCACCGGAATACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAAt                              <  1:2581928/69‑1 (MQ=255)                                    ttAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGAt                              >  1:3208379/1‑68 (MQ=255)                                    ttAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATc                             >  1:333286/1‑69 (MQ=255)                                               gtTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCAGCGCAAAAcc                  >  1:1980053/1‑69 (MQ=255)                                                           aaTACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGaa      <  1:1948636/69‑1 (MQ=255)                                                           aaTACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGaa      <  1:1435480/69‑1 (MQ=255)                                                           aaTACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGaa      <  1:2623768/69‑1 (MQ=255)                                                            aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGaa      >  1:948065/1‑68 (MQ=255)                                                            aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGaa      >  1:3009856/1‑68 (MQ=255)                                                            aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGaa      >  1:2063640/1‑68 (MQ=255)                                                            aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGaa      >  1:1992461/1‑68 (MQ=255)                                                            aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGAAt     >  1:2804589/1‑69 (MQ=255)                                                            aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGAAt     >  1:2518255/1‑69 (MQ=255)                                                            aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGAAt     >  1:1209389/1‑69 (MQ=255)                                                            aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGAAt     >  1:1186120/1‑69 (MQ=255)                                                            aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGAAt     >  1:3360277/1‑69 (MQ=255)                                                            aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGAAt     >  1:444566/1‑69 (MQ=255)                                                            aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGAAt     >  1:1178516/1‑69 (MQ=255)                                                            aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGAAt     >  1:857902/1‑69 (MQ=255)                                                            aTACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAATTGGCGCAAAACCAGCAGGAGCGAAt     >  1:1480765/1‑69 (MQ=255)                                                               cTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTg  >  1:2918117/1‑69 (MQ=255)                                                                 |                                                                   CCTGGAAGATAGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTG  >  minE/572155‑572285

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