Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	984,541	C→T	60.5%	S203N (AGT→AAT)	ynfC ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	984,541	0	C	T	60.5%	36.4 / 59.0	86	S203N (AGT→AAT)	ynfC	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (19/15);  new base T (24/28);  total (43/43)
Fisher's exact test for biased strand distribution p-value = 5.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

CGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTACATTACCAACCCGTTGATTATTCAGTAAAGTA  >  minE/984479‑984597                                                               |                                                         cGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCat                                                    <  1:249999/69‑1 (MQ=255) cGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCat                                                    <  1:2916391/69‑1 (MQ=255)       cggGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTcac                                              <  1:1066128/69‑1 (MQ=255)       cggGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTcac                                              <  1:2214899/69‑1 (MQ=255)       cggGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTcac                                              <  1:1149105/69‑1 (MQ=255)       cggGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTcac                                              <  1:3364587/69‑1 (MQ=255)        gggTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTcac                                              <  1:1397194/68‑1 (MQ=255)           ttttACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGct                                          >  1:2867468/1‑69 (MQ=255)           ttttACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGc                                           >  1:1931957/1‑68 (MQ=255)               aCTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTg                                       <  1:277640/68‑1 (MQ=255)                     tcatcaACAATGATTATCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTa                                  <  1:37699/67‑1 (MQ=255)                     tcatcaACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTa                                  <  1:3378546/67‑1 (MQ=255)                     tcatcaACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTa                                  <  1:37730/67‑1 (MQ=255)                     tcatcaACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTa                                  <  1:405760/67‑1 (MQ=255)                     tcatcaACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTa                                  <  1:1798517/67‑1 (MQ=255)                       atcaACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTa                                  <  1:2539783/65‑1 (MQ=255)                         caacaaTGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTACATTa                             <  1:3242995/68‑1 (MQ=255)                             aaTGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTACATTAcc                           >  1:979248/1‑66 (MQ=255)                             aaTGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTACATTAcc                           >  1:2823010/1‑66 (MQ=255)                             aaTGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTACATTAcc                           >  1:2716676/1‑66 (MQ=255)                             aaTGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTACATTAcc                           >  1:1822611/1‑66 (MQ=255)                             aaTGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTACATTAcc                           >  1:1425573/1‑66 (MQ=255)                             aaTGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTACATTAcc                           >  1:1372636/1‑66 (MQ=255)                             aaTGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTACATTAcc                           >  1:13207/1‑66 (MQ=255)                                 aTTAGCTGACAGTCCACCGGATTGGCGTGATTGTCATATTCACAGCTCTGtttt                                   >  1:1488789/1‑54 (MQ=255)                                 aTTAGCTGACAGTCCACCGGATTGGCGTGATTGTCATATTCACAGCTCTGTTTTa                                  >  1:2944604/1‑55 (MQ=255)                                 aTTAGCTGACAGTCCACCGGATTGGCGTGATTGTCATATTCACAGCTCTGTTTTa                                  >  1:567020/1‑55 (MQ=255)                                 aTTAGCTGACAGTCCACCGGATTGGCGTGATTGTCATATTCACAGCTCTGTTTTa                                  >  1:2708958/1‑55 (MQ=255)                                 aTTAGCTGACAGTCCACCGGATTGGCGTGATTGTCATATTCACAGCTCTGTTTTa                                  >  1:1448393/1‑55 (MQ=255)                                 aTTAGCTGACAGTCCACCGGATTGGCGTGATTGTCATATTCACAGCTCTGTTTTa                                  >  1:2798012/1‑55 (MQ=255)                                 aTTAGCTGACAGTCCACCGGATTGGCGTGATTGTCATATTCACAGCTCTGTTTTa                                  >  1:1574649/1‑55 (MQ=255)                                 aTTAGCTGACAGTCCACCGGATTGGCGTGATTGTCATATTCACAGCTCTGTTTTa                                  >  1:1186564/1‑55 (MQ=255)                                 aTTAGCTGACAGTCCACCGGATTGGCGTGATTGTCATATTCACAGCTCTGTTTTa                                  >  1:3095754/1‑55 (MQ=255)                                 aTTAGCTGACAGTCCACCGGATTGGCGTGATTGTCATATTCACAGCTCTGTTTTa                                  >  1:2092890/1‑55 (MQ=255)                                 aTTAGCTGACAGTCCACCGGATTGGCGTGATTGTCATATTCACAGCTCTGTTTTa                                  >  1:2029581/1‑55 (MQ=255)                                 aTTAGCTGACAGTCCACCGGATTGGCGTGATTGTCATATTCACAGCTCTGTTTTa                                  >  1:1774971/1‑55 (MQ=255)                                 aTTAGCTGACAGTCCACCGGATTGGCGTGATTGTCATATTCACAGCTCTGTTTTa                                  >  1:3312930/1‑55 (MQ=255)                                    aGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTACATTACCAACCCGTTg                  >  1:2230866/1‑68 (MQ=255)                                    aGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTACATTACCAACCCGTTg                  >  1:2320955/1‑68 (MQ=255)                                    aGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTACATTACCAACCCGTTg                  >  1:2366039/1‑68 (MQ=255)                                    aGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTACATTACCAACCCGTTg                  >  1:3177095/1‑68 (MQ=255)                                    aGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTACATTACCAACCCGTTg                  >  1:255475/1‑68 (MQ=255)                                                   ggATTAGCGTGACTGTCATATTCACAGCTCTGTTTTACATTACCAACCCGTTGATTATTCagtaaagta  >  1:1712197/1‑69 (MQ=255)                                                   ggATTAGCGTGACTGTCATATTCACAGCTCTGTTTTACATTACCAACCCGTTGATTATTCagtaaagta  >  1:3264719/1‑69 (MQ=255)                                                   ggATTAGCGTGACTGTCATATTCACAGCTCTGTTTTACATTACCAACCCGTTGATTATTCagtaaagta  >  1:2616298/1‑69 (MQ=255)                                                   ggATTAGCGTGACTGTCATATTCACAGCTCTGTTTTACATTACCAACCCGTTGATTATTCagtaaagta  >  1:811973/1‑69 (MQ=255)                                                   ggATTAGCGTGACTGTCATATTCACAGCTCTGTTTTACATTACCAACCCGTTGATTATTCagtaaagta  >  1:882477/1‑69 (MQ=255)                                                    gATTGGCGTGATTGTCATATTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtaaagta  >  1:1211449/1‑68 (MQ=255)                                                    gATTGGCGTGATTGTCATATTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtaaagta  >  1:787155/1‑68 (MQ=255)                                                    gATTGGCGTGATTGTCATATTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtaaagta  >  1:2040011/1‑68 (MQ=255)                                                    gATTGGCGTGATTGTCATATTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtaaagta  >  1:1279261/1‑68 (MQ=255)                                                    gATTGGCGTGATTGTCATATTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtaaagta  >  1:3133574/1‑68 (MQ=255)                                                    gATTGGCGTGATTGTCATATTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtaaagta  >  1:214146/1‑68 (MQ=255)                                                    gATTGGCGTGATTGTCATATTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtaaagta  >  1:35404/1‑68 (MQ=255)                                                    gATTGGCGTGATTGTCATATTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtaaagta  >  1:1171162/1‑68 (MQ=255)                                                    gATTGGCGTGATTGTCATATTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtaaagta  >  1:2338312/1‑68 (MQ=255)                                                    gATTGGCGTGATTGTCATATTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtaaagta  >  1:863770/1‑68 (MQ=255)                                                    gATTGGCGTGATTGTCATATTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtaaagt   >  1:2284392/1‑67 (MQ=255)                                                       tggCGTGATTGTCATTTTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtat      <  1:971934/59‑2 (MQ=255)                                                       tggCGTGATTGTCATATTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtaa      <  1:1143679/59‑1 (MQ=255)                                                       tggCGTGATTGTCATATTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtaa      <  1:390160/59‑1 (MQ=255)                                                       tggCGTGATTGTCATATTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtaa      <  1:109940/59‑1 (MQ=255)                                                       tggCGTGATTGTCATATTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtaa      <  1:914786/59‑1 (MQ=255)                                                       tggCGTGATTGTCATATTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtaa      <  1:426243/59‑1 (MQ=255)                                                       tggCGTGATTGTCATATTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtaa      <  1:875034/59‑1 (MQ=255)                                                       tggCGTGATTGTCATATTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtaa      <  1:1007696/59‑1 (MQ=255)                                                       tggCGTGATTGTCATATTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtaa      <  1:337709/59‑1 (MQ=255)                                                       tggCGTGATTGTCATATTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtaa      <  1:470669/59‑1 (MQ=255)                                                       tggCGTGATTGTCATATTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtaa      <  1:2681641/59‑1 (MQ=255)                                                       tggCGTGATTGTCATATTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtaa      <  1:1981069/59‑1 (MQ=255)                                                       tggCGTGATTGTCATATTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtaa      <  1:2083869/59‑1 (MQ=255)                                                       tggCGTGATTGTCATATTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtaa      <  1:211343/59‑1 (MQ=255)                                                       tggCGTGATTGTCATATTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtaa      <  1:2304998/59‑1 (MQ=255)                                                       tggCGTGATTGTCATATTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtaa      <  1:2339272/59‑1 (MQ=255)                                                       tggCGTGATTGTCATATTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtaa      <  1:1563347/59‑1 (MQ=255)                                                       tggCGTGATTGTCATATTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtaa      <  1:2539130/59‑1 (MQ=255)                                                       tggCGTGATTGTCATATTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtaa      <  1:2678003/59‑1 (MQ=255)                                                       tggCGTGATTGTCATATTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtaa      <  1:331545/59‑1 (MQ=255)                                                       tggCGTGATTGTCATATTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtaa      <  1:2738879/59‑1 (MQ=255)                                                       tggCGTGATTGTCATATTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtaa      <  1:2810282/59‑1 (MQ=255)                                                       tggCGTGATTGTCATATTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtaa      <  1:2841180/59‑1 (MQ=255)                                                       tggCGTGATTGTCATATTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtaa      <  1:2844847/59‑1 (MQ=255)                                                       tggCGTGATTGTCATATTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtaa      <  1:28633/59‑1 (MQ=255)                                                       tggCGTGATTGTCATATTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtaa      <  1:2864787/59‑1 (MQ=255)                                                       tggCGTGATTGTCATATTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtaa      <  1:1348402/59‑1 (MQ=255)                                                       tggCGTGATTGTCATATTCACAGCTCTGTTTTACATTTCCAACCCGTTGATTATTCagtaa      <  1:3085286/59‑1 (MQ=255)                                                               |                                                         CGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTACATTACCAACCCGTTGATTATTCAGTAAAGTA  >  minE/984479‑984597

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