Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	572,195	A→G	33.3%	S184S (TCT→TCC)	aqpZ ←	aquaporin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	572,195	0	A	G	33.3%	26.1 / 17.9	39	S184S (TCT→TCC)	aqpZ	aquaporin
Reads supporting (aligned to +/- strand):  ref base A (8/18);  new base G (10/3);  total (18/21)
Fisher's exact test for biased strand distribution p-value = 1.54e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATC  >  minE/572133‑572258                                                               |                                                                aGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGCAACCGCAGTACTGCGCGCCGGGTTAACGGAAgt                                                            >  1:671215/1‑68 (MQ=255) aGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGCAACCGCAGTACTGCGCGCCGGGTTAACGGAAgt                                                            >  1:2853452/1‑68 (MQ=255) aGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGCAACCGCAGTACTGCGCGCCGGGTTAACGGAAgt                                                            >  1:290551/1‑68 (MQ=255) aGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGCAACCGCAGTACTGCGCGCCGGGTTAACGGAAgt                                                            >  1:2296941/1‑68 (MQ=255) aGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGCAACCGCAGTACTGCGCGCCGGGTTAACGGAAgt                                                            >  1:3322624/1‑68 (MQ=255) aGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGCAACCGCAGTACTGCGCGCCGGGTTAACGGAAgt                                                            >  1:1610671/1‑68 (MQ=255) aGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGCAACCGCAGTACTGCGCGCCGGGTTAACGGAAgt                                                            >  1:507714/1‑68 (MQ=255) aGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGCAACCGCAGTACTGCGCGCCGGGTTAACGGAAgt                                                            >  1:515047/1‑68 (MQ=255) aGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGCAACCGCAGTACTGCGCGCCGGGTTAACGGAAgt                                                            >  1:1403411/1‑68 (MQ=255) aGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGCAACCGCAGTACTGCGCGCCGGGTTAACGGAAgt                                                            >  1:1213585/1‑68 (MQ=255)           aTGCCCAGCCGCCCTGGAAGATAGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCAcc                                                  >  1:2783685/1‑68 (MQ=255)           aTGCCCAGCCGCCCTGGAAGATAGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCAcc                                                  >  1:1637576/1‑68 (MQ=255)           aTGCCCAGCCGCCCTGGAAGATAGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCAcc                                                  >  1:2474718/1‑68 (MQ=255)           aTGCCCAGCCGCCCTGGAAGATAGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCAcc                                                  >  1:1753789/1‑68 (MQ=255)           aTGCCCAGCCGCCCTGGAAGATAGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCAcc                                                  >  1:1186698/1‑68 (MQ=255)             gCCCAGCCGCCCTGGAAGATAGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCgg                                                >  1:2408818/1‑68 (MQ=255)              cccAGCCGCCCTGGAAGATAGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGaa                                              <  1:2035157/69‑1 (MQ=255)              cccAGCCGCCCTGGAAGATAGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGaa                                              <  1:2054175/69‑1 (MQ=255)              cccAGCCGCCCTGGAAGATAGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGaa                                              <  1:1847001/69‑1 (MQ=255)              cccAGCCGCCCTGGAAGATAGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGaa                                              <  1:3358777/69‑1 (MQ=255)              cccAGCCGCCCTGGAAGATAGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGaa                                              <  1:1504695/69‑1 (MQ=255)               ccAGCCGCCCTGGAAGATAGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGaa                                              <  1:1163275/68‑1 (MQ=255)                       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:657155/69‑1 (MQ=255)                                     aCCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGtt                       <  1:2249836/69‑1 (MQ=255)                                     aCCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGtt                       <  1:2321927/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)                                     aCCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGtt                       <  1:1134956/69‑1 (MQ=255)                                      ccGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTAACCGGAATACTAATTAAGTGAATCAGGGtt                       <  1:1932223/68‑1 (MQ=255)                                                         gTTAACGGAAGTGTTAGTCACCGGAATACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAAt   <  1:2581928/69‑1 (MQ=255)                                                         gTTAACGGAAGTGTTAGTCACCGGAATACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAAt   <  1:1428513/69‑1 (MQ=255)                                                         gTTAACGGAAGTGTTAGTCACCGGAATACTGATTAAGTGAATCAGGGTTAAGGCCAGACCAATGGCAAt   <  1:1157276/69‑1 (MQ=255)                                                          ttAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGAt   >  1:3208379/1‑68 (MQ=255)                                                          ttAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATc  >  1:333286/1‑69 (MQ=255)                                                               |                                                                AGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATC  >  minE/572133‑572258

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