Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	563,959	A→G	22.8%	A314A (GCT→GCC)	ybjL ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	563,959	0	A	G	22.8%	106.7 / 25.4	83	A314A (GCT→GCC)	ybjL	predicted transporter
Reads supporting (aligned to +/- strand):  ref base A (18/46);  new base G (6/13);  total (24/59)
Fisher's exact test for biased strand distribution p-value = 7.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGACGCGGTTAAGGAAGCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACGCATGTCGAGAAGGTCACGATCGAA  >  minE/563896‑564024                                                                |                                                                  aaTGACGCGGTTAAGGAAGCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGt                                                               <  1:1961869/68‑1 (MQ=255) aaTGACGCGGTTAAGGAAGCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGt                                                               <  1:2120073/68‑1 (MQ=255) aaTGACGCGGTTAAGGAAGCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGt                                                               <  1:1509623/68‑1 (MQ=255) aaTGACGCGGTTAAGGAAGCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGt                                                               <  1:1252723/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:714912/1‑68 (MQ=255)            tAAGGAAGCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTa                                                    >  1:1422601/1‑68 (MQ=255)            tAAGGAAGCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTa                                                    >  1:562943/1‑68 (MQ=255)             aaGGAAGCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAAc                                                  <  1:1234861/69‑1 (MQ=255)             aaGGAAGCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAAc                                                  <  1:1467601/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:1670857/69‑1 (MQ=255)                   gCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCAc                                            <  1:1561458/69‑1 (MQ=255)                   gCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCAc                                            >  1:1416453/1‑69 (MQ=255)                      aCCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTGCCCACGGCGTTATGGttttt                                                     <  1:463194/57‑1 (MQ=255)                      aCCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTGCCCACGGCGTTATGGttttt                                                     <  1:2834615/57‑1 (MQ=255)                      aCCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTGCCCACGGCGTTATGGttttt                                                     <  1:2787635/57‑1 (MQ=255)                      aCCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTGCCCACGGCGTTATGGttttt                                                     <  1:1275692/57‑1 (MQ=255)                      aCCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTGCCCACGGCGTTATGGttttt                                                     <  1:2643931/57‑1 (MQ=255)                      aCCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTGCCCACGGCGTTATGGttttt                                                     <  1:1313381/57‑1 (MQ=255)                      aCCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTGCCCACGGCGTTATGGttttt                                                     <  1:431755/57‑1 (MQ=255)                      aCCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTGCCCACGGCGTTATGGttttt                                                     <  1:2481670/57‑1 (MQ=255)                      aCCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTGCCCACGGCGTTATGGttttt                                                     <  1:622355/57‑1 (MQ=255)                      aCCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTGCCCACGGCGTTATGGttttt                                                     <  1:924322/57‑1 (MQ=255)                      aCCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTGCCCACGGCGTTATGGttttt                                                     <  1:120833/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:1407351/1‑69 (MQ=255)                                   aaCTTCAGTTGTGCGAGACGTTTGCCCACGGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATa                            >  1:2297244/1‑69 (MQ=255)                                     cTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGATGATACg                          <  1:437493/69‑1 (MQ=255)                                     cTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACg                          <  1:413117/69‑1 (MQ=255)                                     cTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACg                          <  1:56471/69‑1 (MQ=255)                                     cTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACg                          <  1:3313751/69‑1 (MQ=255)                                     cTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACg                          <  1:3171179/69‑1 (MQ=255)                                     cTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACg                          <  1:1981943/69‑1 (MQ=255)                                     cTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACg                          <  1:634016/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)                                     cTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACg                          <  1:2898586/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:2978238/68‑1 (MQ=255)                                               gCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACGCATGTCgag                 <  1:2851913/68‑1 (MQ=255)                                               gCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACGCATGTCgag                 <  1:1392232/68‑1 (MQ=255)                                               gCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACGCATGTCgag                 <  1:2013289/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)                                                             tACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACGCATGTCGAGAAGGTCACGATCGaa  >  1:408761/1‑69 (MQ=255)                                                                |                                                                  AATGACGCGGTTAAGGAAGCAACCGTGATCGGTCAACTTCAGTTGTGCGAGACGTTTACCTACAGCGTTATGGTTTTTAACGACCACTTCTTCAGTGACGATACGCATGTCGAGAAGGTCACGATCGAA  >  minE/563896‑564024

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