Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	563,953	A→G	22.2%	G316G (GGT→GGC)	ybjL ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	563,953	0	A	G	22.2%	112.8 / 25.5	85	G316G (GGT→GGC)	ybjL	predicted transporter
Reads supporting (aligned to +/- strand):  ref base A (20/46);  new base G (6/13);  total (26/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.74e-01

CTGGCTACGAATGACGCGGTTAAGGAAGCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACGCATGTCGAGAAGGTCACGA  >  minE/563887‑564019                                                                   |                                                                   cTGGCTACGAATGACGCGGTTAAGGAAGCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTAcc                                                                  >  1:1608171/1‑69 (MQ=255) cTGGCTACGAATGACGCGGTTAAGGAAGCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTAcc                                                                  >  1:2865057/1‑69 (MQ=255) cTGGCTACGAATGACGCGGTTAAGGAAGCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTAc                                                                   >  1:701158/1‑68 (MQ=255)          aaTGACGCGGTTAAGGAAGCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGt                                                          <  1:1961869/68‑1 (MQ=255)          aaTGACGCGGTTAAGGAAGCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGt                                                          <  1:2120073/68‑1 (MQ=255)          aaTGACGCGGTTAAGGAAGCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGt                                                          <  1:1252723/68‑1 (MQ=255)          aaTGACGCGGTTAAGGAAGCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGt                                                          <  1:1509623/68‑1 (MQ=255)          aaTGACGCGGTTAAGGAAGCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGt                                                          <  1:751678/68‑1 (MQ=255)                   gTTAAGGAAGCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGttttt                                                <  1:200300/69‑1 (MQ=255)                   gTTAAGGAAGCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGttttt                                                <  1:831181/69‑1 (MQ=255)                    ttAAGGAAGCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGttttt                                                <  1:2610642/68‑1 (MQ=255)                     tAAGGAAGCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTa                                               >  1:1284520/1‑68 (MQ=255)                     tAAGGAAGCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTa                                               >  1:1422601/1‑68 (MQ=255)                     tAAGGAAGCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTa                                               >  1:714912/1‑68 (MQ=255)                     tAAGGAAGCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTa                                               >  1:562943/1‑68 (MQ=255)                      aaGGAAGCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAAc                                             <  1:1467601/69‑1 (MQ=255)                      aaGGAAGCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAAc                                             <  1:1234861/69‑1 (MQ=255)                      aaGGAAGCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAAc                                             <  1:675339/69‑1 (MQ=255)                      aaGGAAGCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAAc                                             <  1:1208182/69‑1 (MQ=255)                      aaGGAAGCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAAc                                             <  1:2037208/69‑1 (MQ=255)                      aaGGAAGCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAAc                                             <  1:2006353/69‑1 (MQ=255)                      aaGGAAGCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAAc                                             <  1:3333843/69‑1 (MQ=255)                         gAAGCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGAc                                          <  1:1613261/69‑1 (MQ=255)                            gCAACCGTGATCGGTCAACTTCAGTTTTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGTCCAc                                       <  1:227716/69‑1 (MQ=255)                            gCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCAc                                       <  1:740491/69‑1 (MQ=255)                            gCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCAc                                       <  1:231999/69‑1 (MQ=255)                            gCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCAc                                       <  1:1192223/69‑1 (MQ=255)                            gCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCAc                                       >  1:1416453/1‑69 (MQ=255)                            gCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCAc                                       <  1:1561458/69‑1 (MQ=255)                            gCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCAc                                       <  1:1670857/69‑1 (MQ=255)                               aCCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTGCCCACGGCGTTATGGttttt                                                <  1:463194/57‑1 (MQ=255)                               aCCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTGCCCACGGCGTTATGGttttt                                                <  1:1313381/57‑1 (MQ=255)                               aCCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTGCCCACGGCGTTATGGttttt                                                <  1:1275692/57‑1 (MQ=255)                               aCCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTGCCCACGGCGTTATGGttttt                                                <  1:2834615/57‑1 (MQ=255)                               aCCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTGCCCACGGCGTTATGGttttt                                                <  1:2787635/57‑1 (MQ=255)                               aCCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTGCCCACGGCGTTATGGttttt                                                <  1:2643931/57‑1 (MQ=255)                               aCCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTGCCCACGGCGTTATGGttttt                                                <  1:120833/57‑1 (MQ=255)                               aCCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTGCCCACGGCGTTATGGttttt                                                <  1:2481670/57‑1 (MQ=255)                               aCCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTGCCCACGGCGTTATGGttttt                                                <  1:622355/57‑1 (MQ=255)                               aCCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTGCCCACGGCGTTATGGttttt                                                <  1:431755/57‑1 (MQ=255)                               aCCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTGCCCACGGCGTTATGGttttt                                                <  1:924322/57‑1 (MQ=255)                                ccGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTGCCCACGGCGTTATGGttttt                                                <  1:2085886/56‑1 (MQ=255)                                  gTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCACttcttc                                 <  1:2566978/69‑1 (MQ=255)                                   tGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGt                                                          <  1:3382980/43‑1 (MQ=255)                                            aaCTTCAGTTGTGCGAGACGTTTGCCCACGGCGTTa                                                        >  1:2238669/1‑36 (MQ=37)                                            aaCTTCAGTTGTGCGAGACGTTTGCCCACGGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATa                       >  1:857078/1‑69 (MQ=255)                                            aaCTTCAGTTGTGCGAGACGTTTGCCCACGGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATa                       >  1:1188261/1‑69 (MQ=255)                                            aaCTTCAGTTGTGCGAGACGTTTGCCCACGGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATa                       >  1:688560/1‑69 (MQ=255)                                            aaCTTCAGTTGTGCGAGACGTTTGCCCACGGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATa                       >  1:2297244/1‑69 (MQ=255)                                            aaCTTCAGTTGTGCGAGACGTTTGCCCACGGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATa                       >  1:1407351/1‑69 (MQ=255)                                              cTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGATGATACg                     <  1:437493/69‑1 (MQ=255)                                              cTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACg                     <  1:3171179/69‑1 (MQ=255)                                              cTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACg                     <  1:1981943/69‑1 (MQ=255)                                              cTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACg                     <  1:56471/69‑1 (MQ=255)                                              cTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACg                     <  1:634016/69‑1 (MQ=255)                                              cTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACg                     <  1:413117/69‑1 (MQ=255)                                              cTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACg                     <  1:3313751/69‑1 (MQ=255)                                              cTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACg                     <  1:2898586/69‑1 (MQ=255)                                              cTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACg                     <  1:309408/69‑1 (MQ=255)                                              cTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACg                     <  1:3062235/69‑1 (MQ=255)                                              cTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACg                     <  1:2995880/69‑1 (MQ=255)                                              cTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACg                     <  1:2930703/69‑1 (MQ=255)                                               ttCAGTTGTGCGAGACGTTTGCCCACGGCGTTATGGttttt                                                <  1:1828643/41‑1 (MQ=255)                                                tCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACGc                    >  1:1269622/1‑68 (MQ=255)                                                       tgtgagaCGTTTACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACGCATGTCgag            <  1:2323882/66‑1 (MQ=255)                                                        gCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACGCATGTCgag            <  1:657411/68‑1 (MQ=255)                                                        gCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACGCATGTCgag            <  1:1392232/68‑1 (MQ=255)                                                        gCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACGCATGTCgag            <  1:2978238/68‑1 (MQ=255)                                                        gCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACGCATGTCgag            <  1:2013289/68‑1 (MQ=255)                                                        gCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACGCATGTCgag            <  1:2851913/68‑1 (MQ=255)                                                           agaCGTTTACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACGCATGTCGAGAAg         <  1:92560/68‑1 (MQ=255)                                                           agaCGTTTACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACGCATGTCGAGAAg         <  1:3024983/68‑1 (MQ=255)                                                           agaCGTTTACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACGCATGTCGAGAAg         <  1:289986/68‑1 (MQ=255)                                                           agaCGTTTACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACGCATGTCGAGAAg         <  1:2345251/68‑1 (MQ=255)                                                                  tacatacaGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACGCATGTCGAGAAGGTCACGa  >  1:554871/1‑68 (MQ=255)                                                                  tACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACGCATGTCGAGAAGGTCACGa  >  1:562081/1‑68 (MQ=255)                                                                  tACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACGCATGTCGAGAAGGTCACGa  >  1:54629/1‑68 (MQ=255)                                                                  tACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACGCATGTCGAGAAGGTCACGa  >  1:540997/1‑68 (MQ=255)                                                                  tACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACGCATGTCGAGAAGGTCACGa  >  1:3156521/1‑68 (MQ=255)                                                                  tACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACGCATGTCGAGAAGGTCACGa  >  1:306408/1‑68 (MQ=255)                                                                  tACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACGCATGTCGAGAAGGTCACGa  >  1:2745268/1‑68 (MQ=255)                                                                  tACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACGCATGTCGAGAAGGTCACGa  >  1:2640529/1‑68 (MQ=255)                                                                  tACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACGCATGTCGAGAAGGTCACGa  >  1:1563647/1‑68 (MQ=255)                                                                  tACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACGCATGTCGAGAAGGTCACGa  >  1:1207484/1‑68 (MQ=255)                                                                  tACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGAGACGATACGCATGTCga             >  1:1953331/1‑57 (MQ=255)                                                                  tACCTACAGAGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACGCATGTCGAGAAGGTCACGa  >  1:2957314/1‑68 (MQ=255)                                                                   |                                                                   CTGGCTACGAATGACGCGGTTAAGGAAGCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACGCATGTCGAGAAGGTCACGA  >  minE/563887‑564019

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