Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,677,626	C→T	16.7%	A102T (GCA→ACA)	yhgA ←	predicted transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,677,626	0	C	T	16.7%	226.6 / 33.5	90	A102T (GCA→ACA)	yhgA	predicted transposase
Reads supporting (aligned to +/- strand):  ref base C (33/42);  new base T (9/6);  total (42/48)
Fisher's exact test for biased strand distribution p-value = 2.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

TGGTAAAACAACAATGGCACCACCATCGGCAACGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGCGAAATGCGATCAGTTTGTTTGAGGTGCTTTGATGTTCAATCAGA  >  minE/2677559‑2677693                                                                    |                                                                    tGGTAAAACAACAATGGCACCACCATCGGCAACGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGc                                                                  <  1:1181448/71‑1 (MQ=255)   gTAAAACAACAATGGCACCACCATCGGCAACGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGCa                                                                 <  1:200521/70‑1 (MQ=255)     aaaaCAACAATGGCACCACCATCGGCAACGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGTGGCaa                                                                >  1:1874990/1‑69 (MQ=255)     aaaaCAACAATGGCACCACCATCGGCAACGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGTGGCaa                                                                >  1:1513951/1‑69 (MQ=255)      aaaCAACAATGGCACCACCATCGGCAACGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGCAATa                                                              <  1:2110012/70‑1 (MQ=255)         caacaaTGGCACCACCATCGGCAACGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGt                                                          <  1:396277/71‑1 (MQ=255)         caacaaTGGCACCACCATCGGCAACGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGt                                                          <  1:1900760/71‑1 (MQ=255)           acaaTGGCACCACCATCGGCAACGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGTaa                                                        <  1:489076/71‑1 (MQ=255)           acaaTGGCACCACCATCGGCAACGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGTaa                                                        <  1:1212329/71‑1 (MQ=255)             aaTGGCACCACCATCGGCAACGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGTGGCAATAGCGTAACg                                                      <  1:421032/71‑1 (MQ=255)                 gCACCACCATCGGCAACGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGTGGCAATAGCGTAACGcat                                                   >  1:1461310/1‑70 (MQ=255)                 gCACCACCATCGGCAACGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGTAACGcatc                                                  <  1:1180034/71‑1 (MQ=255)                  caccacCATCGGCAACGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGTAACGcatca                                                 <  1:1166663/71‑1 (MQ=255)                   accaccATCGGCAACGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGTAACGcatca                                                 <  1:1918963/70‑1 (MQ=255)                    ccaccaTCGGCAACGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGc                                                                  <  1:1932583/52‑1 (MQ=255)                    ccaccaTCGGCAACGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGc                                                                  <  1:342087/52‑1 (MQ=255)                     caccaTCGGCAACGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGCaa                                                                >  1:1811435/1‑53 (MQ=255)                     caccaTCGGCAACGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGCaa                                                                >  1:1497922/1‑53 (MQ=255)                      accaTCGGCAACGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGc                                              >  1:2136593/1‑70 (MQ=255)                      accaTCGGCAACGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGc                                              >  1:969958/1‑70 (MQ=255)                      accaTCGGCAACGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGc                                              >  1:755800/1‑70 (MQ=255)                      accaTCGGCAACGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGc                                              >  1:1336797/1‑70 (MQ=255)                      accaTCGGCAACGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGc                                              >  1:838938/1‑70 (MQ=255)                      accaTCGGCAACGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGc                                              >  1:1282571/1‑70 (MQ=255)                      accaTCGGCAACGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGCg                                             <  1:307349/71‑1 (MQ=255)                       ccaTCGGCAACGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGCAATAGc                                                            >  1:1544443/1‑55 (MQ=255)                       ccaTCGGCAACGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGTAACGca                                                    <  1:664292/63‑1 (MQ=255)                       ccaTCGGCAACGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGCg                                             >  1:713212/1‑70 (MQ=255)                        caTCGGCAACGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGCa                                                                 <  1:1386116/49‑1 (MQ=255)                        caTCGGCAACGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGCGa                                            <  1:1415468/70‑1 (MQ=255)                        caTCGGCAACGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGCGa                                            <  1:1487288/70‑1 (MQ=255)                        caTCGGCAACGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGTAACGCATCATACGa                                            <  1:1246226/70‑1 (MQ=255)                         aTCGGCAACGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGCGaaa                                          >  1:2146830/1‑71 (MQ=255)                          tCGGCAACGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGTGGCAATAGCGTAACGCATCATGCGAAAt                                         <  1:1443918/71‑1 (MQ=255)                           cGGCAACGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGCGAAATg                                        >  1:1622004/1‑71 (MQ=255)                               aaCGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGCGa                                            <  1:1181803/63‑1 (MQ=255)                               aaCGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGCGAAATGCGAt                                    >  1:1548839/1‑71 (MQ=255)                               aaCGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGCGAAATGCGAt                                    >  1:274004/1‑71 (MQ=255)                                aCGTTTTGTATCCAGCATCTAGGTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGCGAAATGCGATc                                   <  1:825125/71‑1 (MQ=255)                                  gTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGCaa                                                                >  1:98847/1‑40 (MQ=255)                                  gTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGCAATAg                                                             >  1:616619/1‑43 (MQ=255)                                  gTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGCAATAg                                                             >  1:112007/1‑43 (MQ=255)                                  gTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGCGAAATGCGATCAg                                 >  1:1025719/1‑71 (MQ=255)                                  gTTTTGTATCAAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGCGAAATGCGATCAg                                 >  1:880851/1‑71 (MQ=255)                                   ttttGTATCCAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGTaa                                                        >  1:306106/1‑47 (MQ=255)                                   ttttGTATCCAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGCGAAATGCGAt                                    <  1:1135006/67‑1 (MQ=255)                                   ttttGTATCCAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGCGAAATGCGATCAgt                                <  1:33469/71‑1 (MQ=255)                                     ttGTATCCAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGCGAAATgag                                      >  1:756127/1‑61 (MQ=255)                                      tGTATCCAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGCGAAATGCGATCAgttt                              <  1:795960/70‑1 (MQ=255)                                       gTATCCAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGCGAAATGCGATCAgtttgt                            >  1:1468477/1‑71 (MQ=255)                                        tATCCAGCATCAAGGTGATTTTGCATTGTGGCAATAGCGTAACGCATCATGCGAAATGCGATCAgtttgtt                           <  1:1509160/71‑1 (MQ=255)                                        tATCCAGCATCAAGGTGATTTTGCATTGTGGCAATAGCGTAACGCATCATGCGAAATGCGATCAgtttgtt                           <  1:1305451/71‑1 (MQ=255)                                        tATCCAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGCGAAATGCGATCAgtttgt                            <  1:1578243/70‑1 (MQ=255)                                        tATCCAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGCGAAATGCGATCAgtttgt                            <  1:1615702/70‑1 (MQ=255)                                          tCCAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGCGAAATGCGATCAGtttgtttg                         <  1:830080/71‑1 (MQ=255)                                          tCCAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGCGAAATGCGATCAGtttgtttg                         <  1:617625/71‑1 (MQ=255)                                            cAGCATCAAGGTGATTTTGCATTGTGGCAATAGCGTAACGCATCATGCGAAATGCGATCAGTTTGTTTGAg                       >  1:103533/1‑71 (MQ=255)                                            cAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGCGAAATGCGATCAGTTTGTTTGAg                       >  1:1252048/1‑71 (MQ=255)                                                  cAAGGTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGCGAAATGCGATCAGTTTGTTTGAGGTGCtt                 >  1:1498209/1‑71 (MQ=255)                                                   aaGGTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGCGAAATGCGATCAGTTTGTTTGAgg                      >  1:1703968/1‑65 (MQ=255)                                                   aaGGTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGCGAAATGCGATCAGTTTGTTTGAgg                      >  1:1029681/1‑65 (MQ=255)                                                   aaGGTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGCGAAATGCGATCAGTTTGTTTGAGGTGCttt                <  1:313636/71‑1 (MQ=255)                                                   aaGGTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGCGAAATGCGATCAGTTTGTTTGAGGTGCttt                <  1:1600168/71‑1 (MQ=255)                                                   aaGGTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGCGAAATGCGATCAGTTTGTTTGAGGTGCttt                >  1:154184/1‑71 (MQ=255)                                                     ggTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGCGAAATGCGATCAgtttgtt                           >  1:2130713/1‑58 (MQ=255)                                                      gTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGCGAAATGCGATCAgtttg                             <  1:1630405/55‑1 (MQ=255)                                                      gTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGCGAAATGCGATCAgtttg                             <  1:1664547/55‑1 (MQ=255)                                                         aTTTTGCATTGTGGCAATAGCGTAACGCATCATGCGAAATGCGATCAGTTTGTTTGAg                       >  1:1010190/1‑58 (MQ=255)                                                          ttttGCATTGCGGCAATAGCGTAACGCATCATGCGAAATGCGATCAGTTTGTTTGAGGTGCttt                <  1:349354/64‑1 (MQ=255)                                                           tttGCATTGCGGCAATAGCGTAACGCATCATGCGAAATGCGATCAGTTTGTTTGAGGTGCTTTGAt             >  1:627018/1‑66 (MQ=255)                                                            ttGCATTGCGGCAATAGCGTAACGCATCATGCGAAATGCGATCAGTTTGTTTGAGGTGCTTTGATGTTCa        <  1:362821/70‑1 (MQ=255)                                                            ttGCATTGCGGCAATAGCGTAACGCATCATGCGAAATGCGATCAGTTTGTTTGAGGTGCTTTGATGTTCa        <  1:87741/70‑1 (MQ=255)                                                              gCATTGTGGCAATAGCGTAACGCATCATGCGAAATGCGATCAgtttgtt                           >  1:591387/1‑49 (MQ=255)                                                              gCATTGTGGCAATAGCGTAACGCATCATGCGAAATGCGATCAgtttgtt                           >  1:1091860/1‑49 (MQ=255)                                                              gCATTGTGGCAATAGCGTAACGCATCATGCGAAATGCGATCAgtttgtt                           >  1:2174817/1‑49 (MQ=255)                                                              gCATTGTGGCAATAGCGTAACGCATCATGCGAAATGCGATCAgtttgtt                           >  1:2170666/1‑49 (MQ=255)                                                               cATTGCGGCAATAGCGTAACGCATCATGCGAAATGCGATCAGTTTGTTTGAGGTGCTTTGATGTTCAATCa    <  1:979838/71‑1 (MQ=255)                                                               cATTGCGGCAATAGCGTAACGCATCATGCGAAATGCGATCAGTTTGTTTGAGGTGCTTTGATGTTCAATCa    <  1:766271/71‑1 (MQ=255)                                                                aTTGCGGCAATAGCGTAACGCATCATGCGAAATGCGATCAGTTTGTTTTAGGTGCttt                <  1:2030293/58‑1 (MQ=255)                                                                aTTGCGGCAATAGCGTAACGCATCATGCGAAATGCGATCAGTTTGTTTGAGGTGCttt                <  1:194278/58‑1 (MQ=255)                                                                aTTGCGGCAATAGCGTAACGCATCATGCGAAATGCGATCAGTTTGTTTGAGGTGCTTTGATGTTCAAt      >  1:2057236/1‑68 (MQ=255)                                                                aTTGCGGCAATAGCGTAACGCATCATGCGAAATGCGATCAGTTTGTTTGAGGTGCTTTGATGTTCAATCAg   <  1:316585/71‑1 (MQ=255)                                                                 ttGTGGCAATAGCGTAACGCATCATGCGAAATGCGATCAgtttg                             <  1:594120/44‑1 (MQ=255)                                                                 ttGTGGCAATAGCGTAACGCATCATGCGAAATGCGATCAgtttg                             <  1:1257944/44‑1 (MQ=255)                                                                 ttGCGGCAATAGCGTAACGCATCATGCGAAATGCGAt                                    <  1:566356/37‑1 (MQ=255)                                                                 ttGCGGCAATAGCGTAACGCATCATGCGAAATGCGAt                                    <  1:1743630/37‑1 (MQ=255)                                                                 ttGCGGCAATAGCGTAACGCATCATGCGAAATGCGATCAGTTTGTTTGAGGTGCttt                <  1:2173930/57‑1 (MQ=255)                                                                  tGCGGCAATAGCGTAACGCATCATGCGAAATGCGATCAGTTTGTTTGAGGTGCTTTGa              >  1:926009/1‑58 (MQ=255)                                                                  tGCGGCAATAGCGTAACGCATCATGCGAAATGCGATCAGTTTGTTTGAGGTGCTTTGATGTTCAATCAGa  >  1:692702/1‑70 (MQ=255)                                                                   gCGGCAATAGCGTAACGCATCATGCGAAATGCGATCa                                  >  1:976847/1‑37 (MQ=255)                                                                    |                                                                    TGGTAAAACAACAATGGCACCACCATCGGCAACGTTTTGTATCCAGCATCAAGGTGATTTTGCATTGCGGCAATAGCGTAACGCATCATGCGAAATGCGATCAGTTTGTTTGAGGTGCTTTGATGTTCAATCAGA  >  minE/2677559‑2677693

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