Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,282,611	G→A	61.1%	A1464T (GCG→ACG)	yeeJ →	adhesin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,282,611	0	G	A	61.1%	43.2 / 58.6	90	A1464T (GCG→ACG)	yeeJ	adhesin
Reads supporting (aligned to +/- strand):  ref base G (31/3);  new base A (19/35);  total (51/39)
Fisher's exact test for biased strand distribution p-value = 1.46e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATTGACGGCAACCCTCGAAAATGGCGATTCGATGCAACAAACAGTGACCTATGTGCCGAACGTCGCGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGACAATAACGATCTCAC  >  minE/1282546‑1282675                                                                  |                                                                 aaTTGACGGCAACCCTCGAAAATGGCGATTCGATGCAACAAACAGTGACCTATGTGCCGAACGTcgcg                                                                <  1:2275308/68‑1 (MQ=255) aaTTGACGGCAACCCTCGAAAATGGCGATTCGATGCAACAAACAGTGACCTATGTGCCGAACGTcgcg                                                                <  1:1174994/68‑1 (MQ=255)    tGACGGCAACCCTCGAAAATGGCTATTCGATGCAACAAACAGTGAACTATTTGCCGAACGTCACGAATg                                                            <  1:3089547/69‑1 (MQ=255)    tGACGGCAACCCTCGAAAATGGCGATTCGATGCAACAAACAGTGACCTATGTGCCGAACGTCTCGAATg                                                            <  1:340765/69‑1 (MQ=255)    tGACGGCAACCCTCGAAAATGGCGATTCGATGCAACAAACAGTGAACTATGTGCCGAACGTCACGAATg                                                            <  1:2707996/69‑1 (MQ=255)    tGACGGCAACCCTCGAAAATGGCGATTCGATGCAACAAACAGTGAACTATGTGCCGAACGTCACGAATg                                                            <  1:914053/69‑1 (MQ=255)    tGACGGCAACCCTCGAAAATGGCGATTCGATGCAACAAACAGTGAACTATGTGCCGAACGTCACGAATg                                                            <  1:2152283/69‑1 (MQ=255)    tGACGGCAACCCTCGAAAATGGCGATTCGATGCAACAAACAGTGAACTATGTGCCGAACGTCACGAATg                                                            <  1:2208849/69‑1 (MQ=255)    tGACGGCAACCCTCGAAAATGGCGATTCGATGCAACAAACAGTGAACTATGTGCCGAACGTCACGAATg                                                            <  1:2463371/69‑1 (MQ=255)    tGACGGCAACCCTCGAAAATGGCGATTCGATGCAACAAACAGTGAACTATGTGCCGAACGTCACGAATg                                                            <  1:2615641/69‑1 (MQ=255)    tGACGGCAACCCTCGAAAATGGCGATTCGATGCAACAAACAGTGAACTATGTGCCGAACGTCACGAATg                                                            <  1:27040/69‑1 (MQ=255)    tGACGGCAACCCTCGAAAATGGCGATTCGATGCAACAAACAGTGAACTATGTGCCGAACGTCACGAATg                                                            <  1:383341/69‑1 (MQ=255)    tGACGGCAACCCTCGAAAATGGCGATTCGATGCAACAAACAGTGAACTATGTGCCGAACGTCACGAATg                                                            <  1:2897093/69‑1 (MQ=255)    tGACGGCAACCCTCGAAAATGGCGATTCGATGCAACAAACAGTGAACTATGTGCCGAACGTCACGAATg                                                            <  1:469775/69‑1 (MQ=255)    tGACGGCAACCCTCGAAAATGGCGATTCGATGCAACAAACAGTGAACTATGTGCCGAACGTCACGAATg                                                            <  1:2989170/69‑1 (MQ=255)    tGACGGCAACCCTCGAAAATGGCGATTCGATGCAACAAACAGTGAACTATGTGCCGAACGTCACGAATg                                                            <  1:310516/69‑1 (MQ=255)    tGACGGCAACCCTCGAAAATGGCGATTCGATGCAACAAACAGTGAACTATGTGCCGAACGTCACGAATg                                                            <  1:3128549/69‑1 (MQ=255)    tGACGGCAACCCTCGAAAATGGCGATTCGATGCAACAAACAGTGAACTATGTGCCGAACGTCACGAATg                                                            <  1:1026740/69‑1 (MQ=255)    tGACGGCAACCCTCGAAAATGGCGATTCGATGCAACAAACAGTGAACTATGTGCCGAACGTCACGAATg                                                            <  1:1573739/69‑1 (MQ=255)    tGACGGCAACCCTCGAAAATGGCGATTCGATGCAACAAACAGTGAACTATGTGCCGAACGTCACGAATg                                                            <  1:129426/69‑1 (MQ=255)    tGACGGCAACCCTCGAAAATGGCGATTCGATGCAACAAACAGTGAACTATGTGCCGAACGTCACGAATg                                                            <  1:1583164/69‑1 (MQ=255)     gACGGCAACCCTCGAAAATGGCGATTCGATGCAACAAACAGTGAACTATGTGCCGAACGTCACGAATg                                                            <  1:3223947/68‑1 (MQ=255)         gCAACCCTCGAAAATGGCGATTCGATGCAACAAACAGTGACCTATGTGCCGAACGTCGCGAATGCTGaa                                                       >  1:966652/1‑69 (MQ=255)              ccTCGAAAATGGCGATTCGATGCAACAAACAGTGACCTATGTGCCGAACGTCGCGAATGCTGAAATCAc                                                  <  1:3211464/69‑1 (MQ=255)              ccTCGAAAATGGCGATTCGATGCAACAAACAGTGACCTATGTGCCGAACGTCGCGAATGCTGAAATCAc                                                  <  1:2744188/69‑1 (MQ=255)              ccTCGAAAATGGCGATTCGATGCAACAAACAGTGACCTATGTGCCGAACGTCGCGAATGCTGAAATCAc                                                  <  1:1546648/69‑1 (MQ=255)                        ggCGATTCGATGCAACAAACAGTGACCTATGTGCCGAACGTCGCGAATGCTGAAATCACGCTGGCAGcc                                        >  1:3374559/1‑69 (MQ=255)                         gCGATTCGATGCAACAAACAGTGACCTATGTGCCGAACGTCGCGAATGCTGAAATCACGCTGGCAGcc                                        >  1:627322/1‑68 (MQ=255)                         gCGATTCGATGCAACAAACAGTGACCTATGTGCCGAACGTCGCGAATGCTGAAATCACGCTGGCAGcc                                        >  1:391035/1‑68 (MQ=255)                         gCGATTCGATGCAACAAACAGTGACCTATGTGCCGAACGTCGCGAATGCTGAAATCACGCTGGCAGcc                                        >  1:1119232/1‑68 (MQ=255)                         gCGATTCGATGCAACAAACAGTGACCTATGTGCCGAACGTCGCGAATGCTGAAATCACGCTGGCAGcc                                        >  1:2521530/1‑68 (MQ=255)                         gCGATTCGATGCAACAAACAGTGACCTATGTGCCGAACGTCGCGAATGCTGAAATCACGCTGGCAGcc                                        >  1:2517446/1‑68 (MQ=255)                         gCGATTCGATGCAACAAACAGTGACCTATGTGCCGAACGTCGCGAATGCTGAAATCACGCTGGCAGcc                                        >  1:2429303/1‑68 (MQ=255)                         gCGATTCGATGCAACAAACAGTGACCTATGTGCCGAACGTCGCGAATGCTGAAATCACGCTGGCAGcc                                        >  1:336143/1‑68 (MQ=255)                         gCGATTCGATGCAACAAACAGTGACCTATGTGCCGAACGTCGCGAATGCTGAAATCACGCTGGCAGcc                                        >  1:1847548/1‑68 (MQ=255)                         gCGATTCGATGCAACAAACAGTGACCTATGTGCCGAACGTCGCGAATGCTGAAATCACGCTGGCAGcc                                        >  1:1847612/1‑68 (MQ=255)                         gCGATTCGATGCAACAAACAGTGACCTATGTGCCGAACGTCGCGAATGCTGAAATCACGCTGGCAGcc                                        >  1:1998420/1‑68 (MQ=255)                                   gCAACAAACAGTGACCTATGTGCCGAACGTCGCGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATcc                             >  1:1078692/1‑69 (MQ=255)                                   gCAACAAACAGTGACCTATGTGCCGAACGTCGCGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATcc                             >  1:3390811/1‑69 (MQ=255)                                   gCAACAAACAGTGACCTATGTGCCGAACGTCGCGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATcc                             >  1:482771/1‑69 (MQ=255)                                    gaacaaacaGTGACCTATGTGCCGAACGT‑‑CGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGt                          >  1:3109378/2‑69 (MQ=255)                                               gaacTATGTGCCGAACGTCACGAATGCTTAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGa                 <  1:1039658/66‑1 (MQ=255)                                               gaacTATGTGCCGAACGTCACGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGa                 <  1:2359569/66‑1 (MQ=255)                                               gaacTATGTGCCGAACGTCACGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGa                 <  1:314637/66‑1 (MQ=255)                                               gaacTATGTGCCGAACGTCACGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGa                 <  1:1028994/66‑1 (MQ=255)                                               gaacTATGTGCCGAACGTCACGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGa                 <  1:281006/66‑1 (MQ=255)                                               gaacTATGTGCCGAACGTCACGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGa                 <  1:2741501/66‑1 (MQ=255)                                               gaacTATGTGCCGAACGTCACGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGa                 <  1:1230475/66‑1 (MQ=255)                                               gaacTATGTGCCGAACGTCACGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGa                 <  1:2093387/66‑1 (MQ=255)                                               gaacTATGTGCCGAACGTCACGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGa                 <  1:3250090/66‑1 (MQ=255)                                               gaacTATGTGCCGAACGTCACGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGa                 <  1:1531654/66‑1 (MQ=255)                                               gaacTATGTGCCGAACGTCACGAATGCAGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGa                 <  1:2927394/66‑1 (MQ=255)                                                 acTATGTGCCGAACGTCACGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGa                 <  1:435935/66‑1 (MQ=255)                                                 acTATGTGCCGAACGACACGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGa                 <  1:2064672/66‑1 (MQ=255)                                                   tATGTGCCGAACGTCAAGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGa                 <  1:763319/65‑1 (MQ=255)                                                       tgCCGAACGTCACGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGa                 <  1:2491489/61‑1 (MQ=255)                                                         ccGAACGTCGCGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGACAATAACGa        >  1:2609860/1‑68 (MQ=255)                                                         ccGAACGTCGCGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGACAATAACGAt       >  1:1716825/1‑69 (MQ=255)                                                         ccGAACGTCGCGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGACAATAACGAt       >  1:1177787/1‑69 (MQ=255)                                                         ccGAACGTCGCGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGACAATAACGAt       >  1:3176995/1‑69 (MQ=255)                                                         ccGAACGTCGCGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGACAATAACGAt       >  1:1748584/1‑69 (MQ=255)                                                         ccGAACGTCGCGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGACAATAACGAt       >  1:342260/1‑69 (MQ=255)                                                         ccGAACGTCGCGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGACAATAACGAt       >  1:2242179/1‑69 (MQ=255)                                                         ccGAACGTCGCGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGACAATAACGAt       >  1:2141751/1‑69 (MQ=255)                                                         ccGAACGTCACGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTg                         >  1:390670/1‑51 (MQ=255)                                                         ccGAACGTCACGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGACAATAACGa        >  1:2884670/1‑68 (MQ=255)                                                         ccGAACGTCACGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGACAATAACGa        >  1:1019403/1‑68 (MQ=255)                                                         ccGAACGTCACGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGACAATAACGAt       >  1:551442/1‑69 (MQ=255)                                                         ccGAACGTCACGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGACAATAACGAt       >  1:855848/1‑69 (MQ=255)                                                         ccGAACGTCACGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGACAATAACGAt       >  1:866092/1‑69 (MQ=255)                                                         ccGAACGTCACGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGACAATAACGAt       >  1:88903/1‑69 (MQ=255)                                                         ccGAACGTCACGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGACAATAACGAt       >  1:967918/1‑69 (MQ=255)                                                         ccGAACGTCACGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGACAATAACGAt       >  1:985400/1‑69 (MQ=255)                                                         ccGAACGTCACGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGACAATAACGAt       >  1:988895/1‑69 (MQ=255)                                                         ccGAACGTCACGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGACAATAACGAt       >  1:1026501/1‑69 (MQ=255)                                                         ccGAACGTCACGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGACAATAACGAt       >  1:1419522/1‑69 (MQ=255)                                                         ccGAACGTCACGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGACAATAACGAt       >  1:1492083/1‑69 (MQ=255)                                                         ccGAACGTCACGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGACAATAACGAt       >  1:2551802/1‑69 (MQ=255)                                                         ccGAACGTCACGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGACAATAACGAt       >  1:2425960/1‑69 (MQ=255)                                                         ccGAACGTCACGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGACAATAACGAt       >  1:155478/1‑69 (MQ=255)                                                         ccGAACGTCACGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGACAATAACGAt       >  1:2174987/1‑69 (MQ=255)                                                         ccGAACGTCACGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGACAATAACGAt       >  1:2123665/1‑69 (MQ=255)                                                         ccGAACGACACGAATGCTGAAATCACGCAGGCAGCCTCGa                                    >  1:2798527/1‑40 (MQ=37)                                                          cGAACGTCACGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGa                 <  1:460484/58‑1 (MQ=255)                                                              cGTCGCGAATGCTGAAATCACGCTGGCAGCCTCGaa                                   >  1:605392/1‑36 (MQ=255)                                                              cGTCGCGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGACAATAACGATCTCAc  >  1:2329788/1‑69 (MQ=255)                                                              cGTCGCGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGACAATAACGATCTCAc  >  1:1950489/1‑69 (MQ=255)                                                              cGTCGCGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGACAATAACGATCTCAc  >  1:90184/1‑69 (MQ=255)                                                              cGTCGCGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGACAATAACGATCTCAc  >  1:3054557/1‑69 (MQ=255)                                                              cGTCGCGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGACAATAACGATCTCAc  >  1:2566434/1‑69 (MQ=255)                                                              cGTCGCGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGACAATAACGATCTCAc  >  1:2962071/1‑69 (MQ=255)                                                              cGTCGCGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGACAATAACGATCTCAc  >  1:2973356/1‑69 (MQ=255)                                                                  |                                                                 AATTGACGGCAACCCTCGAAAATGGCGATTCGATGCAACAAACAGTGACCTATGTGCCGAACGTCGCGAATGCTGAAATCACGCTGGCAGCCTCGAAGGATCCGGTGATTGCCGACAATAACGATCTCAC  >  minE/1282546‑1282675

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