Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	132,238	A→G	26.5%	L177L (CTA→CTG)	yadD →	predicted transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	132,238	0	A	G	26.5%	125.9 / 51.7	83	L177L (CTA→CTG)	yadD	predicted transposase
Reads supporting (aligned to +/- strand):  ref base A (25/36);  new base G (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

CTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAG  >  minE/132170‑132308                                                                     |                                                                       cTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAacta                                                                        >  1:1291540/1‑69 (MQ=255) cTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAacta                                                                        >  1:5320433/1‑69 (MQ=255)    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)           taaCCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAacta                                                                        >  1:3366443/3‑59 (MQ=255)           tCACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCa                                                             >  1:4331714/1‑70 (MQ=255)           tCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAacta                                                                        >  1:1715133/1‑59 (MQ=255)           tCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAacta                                                                        >  1:5194850/1‑59 (MQ=255)           tCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAacta                                                                        >  1:3654363/1‑59 (MQ=255)            cACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCata                                                           >  1:4456591/1‑71 (MQ=255)            cACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCat                                                            >  1:3151006/1‑70 (MQ=255)            cACCATCGTACCTGACGATGAGATCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCat                                                            >  1:5139991/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)                 tCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAacta                                                                        <  1:1364574/53‑1 (MQ=255)                 tCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGc                                                      <  1:3733406/71‑1 (MQ=255)                 tCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGc                                                      <  1:2811198/71‑1 (MQ=255)                 tCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGc                                                      <  1:1990572/71‑1 (MQ=255)                 tCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGc                                                      >  1:1724271/1‑71 (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:1982600/71‑1 (MQ=255)                       cGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAgcgc                                                >  1:5603944/1‑71 (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:5381492/62‑2 (MQ=255)                                   tCATGCAACATCGGCGTATAGCTCTGCTGGAACTGATTCAAAAGCATATTCGCGACCGCGat                                             <  1:3389478/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:1755295/1‑71 (MQ=255)                                        cAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCt                               <  1:4014833/71‑1 (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:549698/1‑71 (MQ=255)                                               ggcggATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCaa                        >  1:5025523/1‑71 (MQ=255)                                               ggcggATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCaa                        <  1:1932043/71‑1 (MQ=255)                                               ggcggATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCaa                        <  1:1726524/71‑1 (MQ=255)                                                  ggATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACt                      <  1:1854585/70‑1 (MQ=255)                                                         gATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCtt                              <  1:1957148/55‑1 (MQ=255)                                                         gATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCtt                              <  1:1875035/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:3476834/1‑57 (MQ=255)                                                             ctcGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCaa                        >  1:4439828/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:4301254/1‑70 (MQ=255)                                                                  actactGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATcgac      <  1:27477/70‑1 (MQ=255)                                                                     actGCAAAAACATATTCGCCAGCGCGACTTTATGt                                      <  1:2813211/35‑1 (MQ=255)                                                                     actGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAg  >  1:4110214/1‑71 (MQ=255)                                                                     |                                                                       CTGGTGGATATCACCATCACACCGGATGACGAAATCATGCAACATCGGCGGATTGCGATTCTCGAACTACTGCAAAAACATATTCGCCAGCGCGACTTAATGTTATTGCTTGAGCAACTGGTCACGCTGATCGACGAAG  >  minE/132170‑132308

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