Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	132,239	C→A	26.4%	L178M (CTG→ATG) ‡	yadD →	predicted transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	132,239	0	C	A	26.4%	170.7 / 42.6	83	L178M (CTG→ATG) ‡	yadD	predicted transposase
Reads supporting (aligned to +/- strand):  ref base C (25/36);  new base A (14/8);  total (39/44)
Fisher's exact test for biased strand distribution p-value = 8.42e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGG  >  minE/132173‑132309                                                                   |                                                                       gtggATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCa                                                                    <  1:5028866/71‑1 (MQ=255)   ggATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCaaa                                                                  <  1:879856/71‑1 (MQ=255)    gATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCaaaa                                                                 <  1:3186170/71‑1 (MQ=255)    gATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCaaaa                                                                 <  1:2104096/71‑1 (MQ=255)    gATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCaaa                                                                  <  1:5527922/70‑1 (MQ=255)     atatCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCaaaa                                                                 <  1:2659004/70‑1 (MQ=255)        tCACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCa                                                              >  1:4331714/1‑70 (MQ=255)         cACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCata                                                            >  1:4456591/1‑71 (MQ=255)         cACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCat                                                             >  1:5139991/1‑70 (MQ=255)         cACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCat                                                             >  1:3151006/1‑70 (MQ=255)         cACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACata                                                            >  1:137935/1‑71 (MQ=255)          aCCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCata                                                            <  1:1682564/70‑1 (MQ=255)          aCCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACatat                                                           <  1:5411175/71‑1 (MQ=255)             aTCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCg                                                        <  1:4109234/71‑1 (MQ=255)              tCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGc                                                       >  1:1724271/1‑71 (MQ=255)              tCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGc                                                       <  1:3733406/71‑1 (MQ=255)              tCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGc                                                       <  1:2811198/71‑1 (MQ=255)              tCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGc                                                       <  1:1990572/71‑1 (MQ=255)                acacCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGcc                                                      <  1:5225852/70‑1 (MQ=255)                acacCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGcc                                                      <  1:1456262/70‑1 (MQ=255)                  acCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCaaaa                                                                 <  1:4922131/57‑1 (MQ=255)                  acCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCaaaa                                                                 <  1:1591358/57‑1 (MQ=255)                    cGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAgcgc                                                 >  1:5603944/1‑71 (MQ=255)                    cGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAgcgc                                                 <  1:1982600/71‑1 (MQ=255)                    cGGATGACGAAATCATGCAACATAGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAgcg                                                  >  1:4314735/1‑70 (MQ=255)                     ggATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACa                                                              <  1:4411448/57‑1 (MQ=255)                     ggATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACa                                                              <  1:5499548/57‑1 (MQ=255)                       aTGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCaaat                                                                 <  1:1765041/52‑2 (MQ=255)                        tGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGc                                                       <  1:5480886/61‑1 (MQ=255)                        tGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACt                                             <  1:5105414/71‑1 (MQ=255)                          aCGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACtt                                            >  1:4311932/1‑70 (MQ=255)                             agaTCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAAtc                                        >  1:5264174/3‑70 (MQ=255)                              aaTCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGt                                       <  1:4346016/71‑1 (MQ=255)                               aTCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACcgcg                                                <  1:1981471/61‑1 (MQ=255)                               aTCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAAtcgg                                      <  1:3144933/71‑4 (MQ=255)                                tCATGCAACATCGGCGTATAGCTCTGCTTGAACTGATTCAAAAGCATATTCGCGACCGCGat                                              <  1:437409/62‑2 (MQ=255)                                tCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGat                                              <  1:3389478/62‑2 (MQ=255)                                tCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGat                                              <  1:5381492/62‑2 (MQ=255)                                tCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGAttt                                            <  1:4614380/64‑1 (MQ=255)                                tCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTAAtcgg                                      <  1:4086207/70‑4 (MQ=255)                                tCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCaaaaa                                                                <  1:2386059/44‑1 (MQ=255)                                tCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCaaaaa                                                                <  1:3964826/44‑1 (MQ=255)                                tCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCaaaaa                                                                <  1:5304759/44‑1 (MQ=255)                                 cATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCata                                                            >  1:1096838/1‑47 (MQ=255)                                 cATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCat                                                             >  1:3643736/1‑46 (MQ=255)                                 cATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCat                                                             >  1:3007904/1‑46 (MQ=255)                                 cATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCat                                                             >  1:3949494/1‑46 (MQ=255)                                 cATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGat                                              >  1:2371657/1‑60 (MQ=255)                                 cATGCAACATAGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGat                                              >  1:3248330/1‑60 (MQ=255)                                  aTGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACa                                                              >  1:176267/1‑44 (MQ=255)                                   tGCAACATCGGCGGATTGCGATTCTCGAACTACTGCaaa                                                                  <  1:202621/39‑1 (MQ=255)                                    gCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGcga                                                     >  1:2947347/1‑49 (MQ=255)                                     cAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAgcg                                                  >  1:1678405/1‑53 (MQ=255)                                     cAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCt                                <  1:4014833/71‑1 (MQ=255)                                     cAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCt                                >  1:1755295/1‑71 (MQ=255)                                     cAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCt                                >  1:121453/1‑71 (MQ=255)                                      aaCATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTa                                           >  1:2550432/1‑59 (MQ=255)                                      aaCATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGATTTa                                           >  1:2398642/1‑59 (MQ=255)                                       aCATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAgcgc                                                 <  1:825848/52‑1 (MQ=255)                                         aTCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTaa                                          >  1:5389881/1‑57 (MQ=255)                                         aTCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTaa                                          >  1:4756399/1‑57 (MQ=255)                                         aTCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAg                            >  1:3152257/1‑71 (MQ=255)                                           cggcggATTGCGATTCTCGAACTACTGCAAAAACatat                                                           >  1:2487730/1‑38 (MQ=255)                                           cggcggATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAAt                                         >  1:4727016/1‑56 (MQ=255)                                           cggcggATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAAt                                         >  1:4733468/1‑56 (MQ=255)                                           cggcggATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAAt                                         >  1:5388259/1‑56 (MQ=255)                                            ggcggATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCaa                         <  1:1932043/71‑1 (MQ=255)                                            ggcggATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCaa                         <  1:1726524/71‑1 (MQ=255)                                            ggcggATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCaa                         >  1:5025523/1‑71 (MQ=255)                                            ggcggATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCaa                         >  1:549698/1‑71 (MQ=255)                                               ggATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACt                       <  1:1854585/70‑1 (MQ=255)                                                      gATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCtt                               <  1:1875035/55‑1 (MQ=255)                                                      gATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCtt                               <  1:1957148/55‑1 (MQ=255)                                                        ttCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCt              >  1:5450095/1‑70 (MQ=255)                                                        ttCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTg             <  1:2192961/71‑1 (MQ=255)                                                         tctcGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGa            >  1:1205416/1‑71 (MQ=255)                                                         tctcGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGa            >  1:3072992/1‑71 (MQ=255)                                                          ctcGAACTACTGCAAAAACATATTCGCCAGCGCGACTTaa                                          >  1:1582883/1‑40 (MQ=255)                                                          ctcGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCaa                         >  1:4439828/1‑57 (MQ=255)                                                          ctcGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCaa                         >  1:3476834/1‑57 (MQ=255)                                                             gAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTAt                                    >  1:1217928/1‑43 (MQ=255)                                                             gAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGAt           <  1:498423/68‑1 (MQ=255)                                                             gAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATc          <  1:2636217/69‑1 (MQ=255)                                                               actactGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATcgcc       >  1:4945734/1‑68 (MQ=255)                                                               actactGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATcgac       <  1:27477/70‑1 (MQ=255)                                                               actactGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATcgac       >  1:4301254/1‑70 (MQ=255)                                                                  actGCAAAAACATATTCGCCAGCGCGACTTTATGt                                       <  1:2813211/35‑1 (MQ=255)                                                                  actGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAg   >  1:4110214/1‑71 (MQ=255)                                                                   ctGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAAc                        >  1:2780334/1‑49 (MQ=255)                                                                   ctGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAAc                        >  1:820289/1‑49 (MQ=255)                                                                   ctGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATc          <  1:1448313/63‑1 (MQ=255)                                                                   ctGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGaa    >  1:723550/1‑69 (MQ=255)                                                                   ctGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAgg  >  1:3767974/1‑71 (MQ=255)                                                                   |                                                                       GTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAGG  >  minE/132173‑132309

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