Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,469,153	A→G	31.2%	S857P (TCC→CCC)	yfaL ←	adhesin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,469,153	0	A	G	31.2%	71.6 / 40.5	83	S857P (TCC→CCC)	yfaL	adhesin
Reads supporting (aligned to +/- strand):  ref base A (31/26);  new base G (9/17);  total (40/43)
Fisher's exact test for biased strand distribution p-value = 1.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTGATCGCTTACCGAGTCA  >  minE/1469089‑1469220                                                                 |                                                                    ggATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGttt                                                                 <  1:2998298/69‑2 (MQ=255) ggATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTa                                                                 <  1:2628656/69‑1 (MQ=255) ggATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTa                                                                 <  1:2191267/69‑1 (MQ=255) ggATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTa                                                                 <  1:489211/69‑1 (MQ=255) ggATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTa                                                                 <  1:2867406/69‑1 (MQ=255) ggATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTa                                                                 <  1:416512/69‑1 (MQ=255) ggATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTa                                                                 <  1:2202329/69‑1 (MQ=255) ggATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTa                                                                 <  1:1698733/69‑1 (MQ=255) ggATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTa                                                                 <  1:2573127/69‑1 (MQ=255)   aTCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAac                                                               <  1:1041341/69‑1 (MQ=255)   aTCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAac                                                               <  1:636658/69‑1 (MQ=255)   aTCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAac                                                               <  1:2866891/69‑1 (MQ=255)          gCGAAATCAACTACTTTAATACCTGTAGATGTCGGCTCACCAATACCCGTAATGGGGTTAATtaccacg                                                        >  1:1283206/1‑68 (MQ=255)          gCGAAATCAACTACTTTAATACCTGTAGATGTCGGCTCACCAATACCCGTAATGGGGTTAATtaccacg                                                        >  1:2868881/1‑68 (MQ=255)          gCGAAATCAACTACTTTAATACCTGTAGATGTCGGCTCACCAATACCCGTAATGGGGTTAATtaccacg                                                        >  1:301844/1‑68 (MQ=255)          gCGAAATCAACTACTTTAATACCTGTAGATGTCGGCTCACCAATACCCGTAATGGGGTTAATtaccacg                                                        >  1:2515264/1‑68 (MQ=255)          gCGAAATCAACTACTTTAATACCTGTAGATGTCGGCTCACCAATACCCGTAATGGGGTTAATtaccacg                                                        >  1:1309246/1‑68 (MQ=255)          gCGAAATCAACTACTTTAATACCTGTAGATGTCGGCTCACCAATACCCGTAATGGGGTTAATtaccacg                                                        >  1:1211409/1‑68 (MQ=255)          gCGAAATCAACTACTTTAATACCTGTAGATGTCGGCTCACCAATACCCGTAATGGGGTTAATtaccacg                                                        >  1:444036/1‑68 (MQ=255)          gCGAAATCAACTACTTTAATACCTGTAGATGTCGGCTCACCAATACCCGTAATGGGGTTAATtaccacg                                                        >  1:618791/1‑68 (MQ=255)          gCGAAATCAACTACTTTAATACCTGTAGATGTCGGCTCACCAATACCCGTAATGGGGTTAATtaccacg                                                        >  1:705875/1‑68 (MQ=255)             aaaTCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGt                                                      <  1:1833292/68‑1 (MQ=255)             aaaTCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGt                                                      <  1:1905776/68‑1 (MQ=255)             aaaTCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGt                                                      <  1:2675991/68‑1 (MQ=255)             aaaTCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGt                                                      <  1:304852/68‑1 (MQ=255)             aaaTCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGt                                                      <  1:490794/68‑1 (MQ=255)             aaaTCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAAACACCACAGt                                                      <  1:1322917/68‑1 (MQ=255)                     cACTTTAATGCCTGTCGATGTCGGCTCA‑CAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGcc                                             >  1:1264997/1‑68 (MQ=255)                      aCTTTAATGCCTGTCGATGTCGGCTCAGCAATCCCTGTAATGGAGTTa                                                                 <  1:1420643/48‑1 (MQ=255)                            aTGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCcagcag                                        <  1:3428359/67‑1 (MQ=255)                            aTGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCcagcag                                        <  1:974050/67‑1 (MQ=255)                            aTGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGt                                       <  1:3279426/68‑1 (MQ=255)                            aTGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGt                                       <  1:1014206/68‑1 (MQ=255)                            aTGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGt                                       <  1:1443720/68‑1 (MQ=255)                                              tCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCaa                     >  1:406401/1‑68 (MQ=255)                                              tCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCaa                     >  1:393406/1‑68 (MQ=255)                                              tCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCaa                     >  1:2511025/1‑68 (MQ=255)                                              tCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCaa                     >  1:1192780/1‑68 (MQ=255)                                              tCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCaa                     >  1:1856608/1‑68 (MQ=255)                                              tCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCaa                     >  1:3348112/1‑68 (MQ=255)                                              tCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCaa                     >  1:2099172/1‑68 (MQ=255)                                              tCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTCCCAGCAGTATTACCGTTCATCACCaa                     >  1:2105174/1‑68 (MQ=255)                                                 ccAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCcagcag                                        <  1:3432029/46‑1 (MQ=255)                                                    aTCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCcagcat                                        <  1:704828/43‑2 (MQ=255)                                                    aTCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTGATCg              >  1:1392995/1‑69 (MQ=255)                                                    aTCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTGATCg              >  1:391902/1‑69 (MQ=255)                                                    aTCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTGATCg              >  1:3283608/1‑69 (MQ=255)                                                    aTCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTGATCg              >  1:2865773/1‑69 (MQ=255)                                                    aTCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTGATCg              >  1:1807154/1‑69 (MQ=255)                                                    aTCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTGATCg              >  1:2639565/1‑69 (MQ=255)                                                    aTCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTGATCg              >  1:589577/1‑69 (MQ=255)                                                    aTCCCTGAAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTGATCg              >  1:1081747/1‑69 (MQ=255)                                                           tattGGGGTTAATTACCACGGCCGTGTTGCCAGCAGTATTACCGTTCAATACCAGTTGGTc               <  1:245139/58‑1 (MQ=255)                                                           tAATTGGGTTAATTACCACGGCCGTGTTGCCAGCAGTATTACCGTTCAATACCAGTTGGTc               <  1:409562/61‑1 (MQ=255)                                                           tAATGGGGTTAATTACCACGGCCGTGTTGCCATCAGTATTACCGTTCAATACCAGTTGGTc               <  1:3432322/61‑1 (MQ=255)                                                           tAATGGGGTTAATTACCACGGCCGTGTTGCCAGTAGTATTACCGTTCAATACCAGTTGGTc               <  1:1644129/61‑1 (MQ=255)                                                           tAATGGGGTTAATTACCACGGCCGTGTTGCCAGCAGTATTACCGTTCAATACCAGTTGGTc               <  1:1004468/61‑1 (MQ=255)                                                           tAATGGGGTTAATTACCACGGCCGTGTTGCCAGCAGTATTACCGTTCAATACCAGTTGGTc               <  1:1147310/61‑1 (MQ=255)                                                           tAATGGGGTTAATTACCACGGCCGTGTTGCCAGCAGTATTACCGTTCAATACCAGTTGGTc               <  1:1138577/61‑1 (MQ=255)                                                           tAATGGGGTTAATTACCACGGCCGTGTTGCCAGCAGTATTACCGTTCAATACCAGTTGGTc               <  1:558650/61‑1 (MQ=255)                                                           tAATGGGGTTAATTACCACGGCCGTGTTGCCAGCAGTATTACCGTTCAATACCAGTTGGTc               <  1:19297/61‑1 (MQ=255)                                                           tAATGGGGTTAATTACCACGGCCGTGTTGCCAGCAGTATTACCGTTCAATACCAGTTGGTc               <  1:3179100/61‑1 (MQ=255)                                                           tAATGGGGTTAATTACCACGGCCGTGTTGCCAGCAGTATTACCGTTCAATACCAGTTGGTc               <  1:1560909/61‑1 (MQ=255)                                                           tAATGGGGTTAATTACCACGGCCGTGTTGCCAGCAGTATTACCGTTCAATACCAGTTGGTc               <  1:1621669/61‑1 (MQ=255)                                                           tAATGGGGTTAATTACCACGGCCGTGTTGCCAGCAGTATTACCGTTCAATACCAGTTGGTc               <  1:2816263/61‑1 (MQ=255)                                                           tAATGGGGTTAATTACCACGGCCGTGTTGCCAGCAGTATTACCGTTCAATACCAGTTGGTc               <  1:1838768/61‑1 (MQ=255)                                                           tAATGGGGTTAATTACCACGGCCGTGTTGCCAGCAGTATTACCGTTCAATACCAGTTGGTc               <  1:2264990/61‑1 (MQ=255)                                                           tAATGGGGTAAATTACCACGGCCGTGTTGCCAGCAGTATTACCGTTCAATACCAGTTGGTc               <  1:1090492/61‑1 (MQ=255)                                                             aTGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTg                  >  1:2099020/1‑56 (MQ=255)                                                             aTGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTg                  >  1:2122644/1‑56 (MQ=255)                                                             aTGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTg                  >  1:2173930/1‑56 (MQ=255)                                                             aTGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTg                  >  1:1796662/1‑56 (MQ=255)                                                             aTGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTg                  >  1:163178/1‑56 (MQ=255)                                                             aTGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTg                  >  1:1442680/1‑56 (MQ=255)                                                              tGGGGTTAATTACCACGGCCGTGTTGCCAGCAGTATTACCGTTCAATACCAGTTGGTc               <  1:1612916/58‑1 (MQ=255)                                                                gAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTGATCGCTTACCGAGTc   >  1:1112648/1‑68 (MQ=255)                                                                gAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTGATCGCTTACCGAGTCa  >  1:1953842/1‑69 (MQ=255)                                                                gAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTGATCGCTTACCGAGTCa  >  1:833991/1‑69 (MQ=255)                                                                gAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTGATCGCTTACCGAGTCa  >  1:1931694/1‑69 (MQ=255)                                                                gAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTGATCGCTTACCGAGTCa  >  1:203304/1‑69 (MQ=255)                                                                gAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTGATCGCTTACCGAGTCa  >  1:2429039/1‑69 (MQ=255)                                                                gAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTGATCGCTTACCGAGTCa  >  1:297439/1‑69 (MQ=255)                                                                gAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTGATCGCTTACCGAGTCa  >  1:1425488/1‑69 (MQ=255)                                                                 |                                                                    GGATCAGCTGCGAAATCAACCACTTTAATGCCTGTCGATGTCGGCTCACCAATCCCTGTAATGGAGTTAACCACCACAGTTGTGTTGCCAGCAGTATTACCGTTCATCACCAATTGATCGCTTACCGAGTCA  >  minE/1469089‑1469220

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