Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	132,202	A→G	10.5%	E165E (GAA→GAG)	yadD →	predicted transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	132,202	0	A	G	10.5%	252.1 / 16.9	95	E165E (GAA→GAG)	yadD	predicted transposase
Reads supporting (aligned to +/- strand):  ref base A (53/32);  new base G (7/3);  total (60/35)
Fisher's exact test for biased strand distribution p-value = 7.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGT  >  minE/132135‑132272                                                                    |                                                                       aGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATc                                                                     <  1:1731507/71‑1 (MQ=255)  gCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCa                                                                    >  1:2408109/1‑71 (MQ=255)   cTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCAt                                                                   <  1:2266054/71‑1 (MQ=255)     ggCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATg                                                                  >  1:643475/1‑70 (MQ=255)     ggCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATg                                                                  >  1:750181/1‑70 (MQ=255)      gcgcGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATg                                                                  <  1:3744934/69‑1 (MQ=255)      gcgcGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGc                                                                 >  1:922081/1‑70 (MQ=255)      gcgcGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGc                                                                 <  1:1973193/70‑1 (MQ=255)      gcgcGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGc                                                                 <  1:1923041/70‑1 (MQ=255)       cgcgACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATc                                                                     >  1:3571361/1‑65 (MQ=255)        gcgACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCaa                                                               >  1:1598612/1‑70 (MQ=255)        gcgACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCaa                                                               >  1:2757408/1‑70 (MQ=255)        gcgACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAAc                                                              <  1:746966/71‑1 (MQ=255)        gcgACGCGTCTATAACAGTCCTTTCCCGATGGTGGATATCACCATCACACCGGATGACGAAATCATGCaa                                                               >  1:157050/1‑70 (MQ=255)         cgACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACa                                                             <  1:2993070/71‑1 (MQ=255)         cgACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACa                                                             >  1:729607/1‑71 (MQ=255)         cgACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGAAGAAATCATGCAACa                                                             >  1:1576723/1‑71 (MQ=255)           aCGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATc                                                           >  1:3718495/1‑71 (MQ=255)              cgTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATctgc                                                        <  1:1326904/71‑1 (MQ=255)                tCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCAt                                                                   <  1:2143457/58‑1 (MQ=255)                 cTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGa                                                     >  1:2375227/1‑71 (MQ=255)                    taACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGc                                                                 >  1:3311360/1‑56 (MQ=255)                    taACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGc                                                                 >  1:1131343/1‑56 (MQ=255)                     aaCAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGc                                                 >  1:1834634/1‑71 (MQ=255)                      aCAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGc                                                 <  1:3546386/70‑1 (MQ=255)                       cAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGa                                               >  1:2340633/1‑71 (MQ=255)                        aGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGAt                                              >  1:2146896/1‑71 (MQ=255)                          tCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATtc                                            <  1:1865246/71‑1 (MQ=255)                          tCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATtc                                            <  1:2289256/71‑1 (MQ=255)                            cTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTg                                                  >  1:2128209/1‑63 (MQ=255)                            cTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATtctc                                          <  1:1654283/71‑1 (MQ=255)                            cTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATtctc                                          >  1:1154686/1‑71 (MQ=255)                            cTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATtct                                           >  1:2360845/1‑70 (MQ=255)                             tttCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATcggc                                                        <  1:2305727/56‑1 (MQ=255)                               tCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCg                                         >  1:3672365/1‑69 (MQ=255)                                cccGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCg                                         <  1:1911822/68‑1 (MQ=255)                                  cGCTGGTGGATATCACCATCACACCGGATGACGAAATCATg                                                                  <  1:1888484/41‑1 (MQ=255)                                  cGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAacta                                    >  1:1559944/1‑71 (MQ=255)                                   gCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAactac                                   >  1:3282126/1‑71 (MQ=255)                                    cTGGTGGATATCACCATCACACCGGATGACGAAATCATGCa                                                                >  1:1649488/1‑41 (MQ=255)                                    atggtggATATCACCATCACACCGGATGACGAAATCATGCa                                                                >  1:3124053/2‑41 (MQ=255)                                     tggtggATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGaa                                       >  1:2689060/1‑65 (MQ=255)                                     tggtggATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTg                                 >  1:3699296/1‑71 (MQ=255)                                     tggtggATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTg                                 >  1:1378704/1‑71 (MQ=255)                                          gATATCACCATCGTACCTGACGATGAGATCATGCAACATcggcg                                                       >  1:3257111/1‑44 (MQ=255)                                          gATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCaaa                             <  1:2961241/70‑1 (MQ=255)                                           atatCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGa                                               >  1:788052/1‑51 (MQ=255)                                           atatCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGa                                               >  1:597024/1‑51 (MQ=255)                                            tatCACCATCACACCGGATGACGAAATCATGCAAc                                                              >  1:113110/1‑35 (MQ=255)                                            tatCACCATCACACCGGATGACGAAATCATGCAACATcggc                                                        >  1:2884706/1‑41 (MQ=255)                                             atCACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCa                         >  1:524861/1‑71 (MQ=255)                                             atCACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCa                         >  1:312576/1‑71 (MQ=255)                                             atCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCaaa                             <  1:232113/67‑1 (MQ=255)                                              tCACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCa                         >  1:3279009/1‑70 (MQ=255)                                              tCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACa                         >  1:1589625/1‑70 (MQ=255)                                               cACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCg                                         <  1:820672/53‑1 (MQ=255)                                                aCCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCatat                      <  1:3378901/71‑1 (MQ=255)                                                aCCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTAGAact                                     >  1:2996407/1‑56 (MQ=255)                                                aCCATCACACCGGATGACGAAATCATGCAACATcggcg                                                       >  1:3106062/1‑38 (MQ=255)                                                aCCATCACACCGGATGACGAAATCATGCAACATcggcg                                                       >  1:2493090/1‑38 (MQ=255)                                                aCCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACatat                      <  1:1866093/71‑1 (MQ=255)                                                aCCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACatat                      <  1:3916233/71‑1 (MQ=255)                                                 ccATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCaaaa                            >  1:501805/1‑64 (MQ=255)                                                 ccATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCatat                      >  1:369444/1‑70 (MQ=255)                                                 ccATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCatat                      <  1:895428/70‑1 (MQ=255)                                                 ccATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATAtt                     >  1:3018776/1‑71 (MQ=255)                                                  cATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTc                    >  1:3459369/1‑71 (MQ=255)                                                  cATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTc                    >  1:1094608/1‑71 (MQ=255)                                                  aaTCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTc                    >  1:2295944/2‑71 (MQ=255)                                                   aTCACACCGGATGACGAAATCATGCAACATCGGCGGATTg                                                  >  1:3632017/1‑40 (MQ=255)                                                   aTCACACCGGATGACGAAATCATGCAACATCGGCGGATTg                                                  >  1:3328833/1‑40 (MQ=255)                                                   aTCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTc                    >  1:1704122/1‑70 (MQ=255)                                                    tCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATtct                                           >  1:2133987/1‑46 (MQ=255)                                                    tCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCg                                         <  1:3474398/48‑1 (MQ=255)                                                    tCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACata                       >  1:413424/1‑66 (MQ=255)                                                    tCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACa                         <  1:3129806/64‑1 (MQ=255)                                                    tCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCg                   >  1:2795204/1‑70 (MQ=255)                                                    tCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGc                  <  1:822055/71‑1 (MQ=255)                                                     cgtacCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGc                  <  1:353302/67‑1 (MQ=255)                                                     cacaCCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACata                       >  1:1445872/1‑65 (MQ=255)                                                      acacCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGaa                                       >  1:1816522/1‑48 (MQ=255)                                                      acacCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGc                  >  1:938210/1‑69 (MQ=255)                                                      acacCGGATGACGAAATCATGCAACATCGGCGGATAGCGATTCTCg                                         >  1:538351/1‑46 (MQ=255)                                                       cacCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGcc                 <  1:1085283/69‑1 (MQ=255)                                                        acCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAactact                                  <  1:2242153/51‑1 (MQ=255)                                                        acCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAgc              >  1:1215807/1‑71 (MQ=255)                                                        acCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAgc              >  1:3892751/1‑71 (MQ=255)                                                          cGGATGACGAAATCATGCAACATCGGCGGATTGCg                                                >  1:192327/1‑35 (MQ=255)                                                           ggATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGcgcg           <  1:1739920/71‑1 (MQ=255)                                                              tGACGAAATCATGCAACATCGGCGGATTGCGATtct                                           <  1:1927416/36‑1 (MQ=255)                                                              tGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAgc              >  1:3322513/1‑65 (MQ=255)                                                               gACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCaaa                             <  1:2582517/49‑1 (MQ=255)                                                               gACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCaaa                             <  1:368264/49‑1 (MQ=255)                                                                 cGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGc                  <  1:2036055/58‑1 (MQ=255)                                                                 cGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGc                  <  1:3826494/58‑1 (MQ=255)                                                                    aaTCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAAc                          <  1:3338305/47‑1 (MQ=255)                                                                    aaTCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGt  >  1:74025/1‑71 (MQ=255)                                                                    |                                                                       AGCTGGCGCGACGCGTCTATAACAGTCCTTTCCCGCTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGT  >  minE/132135‑132272

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