Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,677,416	G→T	14.9%	R172S (CGT→AGT) ‡	yhgA ←	predicted transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,677,416	0	G	T	14.9%	240.5 / 19.8	87	R172S (CGT→AGT) ‡	yhgA	predicted transposase
Reads supporting (aligned to +/- strand):  ref base G (41/33);  new base T (2/11);  total (43/44)
Fisher's exact test for biased strand distribution p-value = 1.37e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

TGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAG  >  minE/2677349‑2677485                                                                    |                                                                      tGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCg                                                                    >  1:5485432/1‑71 (MQ=255) tGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCg                                                                    >  1:1750307/1‑71 (MQ=255) tGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCg                                                                    >  1:1580988/1‑71 (MQ=255) tGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCg                                                                    >  1:4488858/1‑71 (MQ=255) tGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGAGCCATACGTCg                                                                    >  1:4893237/1‑71 (MQ=255)   ctctACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCgg                                                                   <  1:422764/70‑1 (MQ=255)   ctctACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCgg                                                                   <  1:2703192/70‑1 (MQ=255)   ctctACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCgg                                                                   <  1:1545122/70‑1 (MQ=255)    tctACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtg                                                                 >  1:140882/1‑71 (MQ=255)        ccAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtgctg                                                              >  1:2235050/1‑70 (MQ=255)        ccAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtgctg                                                              >  1:2848944/1‑70 (MQ=255)        ccAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtgctg                                                              >  1:15778/1‑70 (MQ=255)          aGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtg                                                                 >  1:4615940/1‑65 (MQ=255)          aGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtg                                                                 >  1:2862998/1‑65 (MQ=255)           atcccAACAGATCGCGCTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCAt                                                           <  1:2602926/68‑1 (MQ=255)            tcccAACAGATCGCGCTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCAtaa                                                         >  1:5350984/2‑69 (MQ=255)              ccATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGAtt                                                       <  1:1087263/71‑1 (MQ=255)              ccATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGAtt                                                       <  1:4987377/71‑1 (MQ=255)                 tCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCa                                                            >  1:4551637/1‑63 (MQ=255)                 tCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTtca                                                    >  1:2979993/1‑71 (MQ=255)                  cAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTtca                                                    >  1:4320673/1‑70 (MQ=255)                   aGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTtcat                                                   >  1:3072811/1‑70 (MQ=255)                   aGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTtcat                                                   >  1:3514677/1‑70 (MQ=255)                   aGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTtcat                                                   >  1:1993096/1‑70 (MQ=255)                     gTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGt                                                                      <  1:4713278/49‑1 (MQ=255)                     gTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTtcatcat                                                <  1:2660088/71‑1 (MQ=255)                        gcgcTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCtcgtcatcc                                              <  1:1036944/70‑2 (MQ=255)                        gcgcTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCtcgtcatcc                                              <  1:4776253/70‑2 (MQ=255)                        gcgcTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCtcgtcatcc                                              <  1:4977092/70‑2 (MQ=255)                        gcgTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCgg                                                                   >  1:5460968/1‑49 (MQ=255)                         cgctGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCgg                                            <  1:5353350/68‑1 (MQ=255)                           ttGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTc                                                                     <  1:886975/44‑1 (MQ=255)                           ttGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCa                                          <  1:5202700/71‑1 (MQ=255)                            tGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCAt                                         >  1:3925169/1‑71 (MQ=255)                             gACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCAt                                         >  1:1599291/1‑70 (MQ=255)                              aCGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCAt                                         <  1:921756/69‑1 (MQ=255)                              aCGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGa                                       <  1:5650889/71‑1 (MQ=255)                                gAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACg                                     >  1:2702693/1‑71 (MQ=255)                                gAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCAta                                                          <  1:3023090/50‑2 (MQ=255)                                gAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCAta                                                          <  1:5350784/50‑2 (MQ=255)                                gAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCAta                                                          <  1:3992153/50‑2 (MQ=255)                                 aaTATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtg                                                                 <  1:2917339/42‑1 (MQ=255)                                  atatGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGt                                   >  1:5192344/1‑71 (MQ=255)                                  atatGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGt                                   <  1:2240814/71‑1 (MQ=255)                                  atatGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGt                                   >  1:1906093/1‑71 (MQ=255)                                  atatGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGt                                   >  1:1891867/1‑71 (MQ=255)                                  atatGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGt                                   >  1:665995/1‑71 (MQ=255)                                   tatGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGa                                                         >  1:1334354/1‑48 (MQ=255)                                   tatGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGt                                   <  1:889489/70‑1 (MQ=255)                                     tGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGAc                                <  1:389629/71‑1 (MQ=255)                                     tGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGAc                                <  1:1823691/71‑1 (MQ=255)                                     tGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGAATGATTTCATCATCAGGCATGACGGTGAc                                <  1:2451865/71‑1 (MQ=255)                                      gTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACa                               >  1:1055009/1‑71 (MQ=255)                                       tttttGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAg                                             <  1:5547609/56‑1 (MQ=255)                                       tttttGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAg                                             <  1:3642943/56‑1 (MQ=255)                                        ttttGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACat                              >  1:2630057/1‑70 (MQ=255)                                        ttttGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACat                              >  1:2246251/1‑70 (MQ=255)                                         tttGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACatc                             <  1:3539959/70‑1 (MQ=255)                                         tctGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCtcgtcatc                                               >  1:336148/3‑52 (MQ=255)                                         tctGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCCGCACCACGGTAATatc                             <  1:5120915/68‑1 (MQ=255)                                          ctGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACGGTAATatcc                            <  1:1557292/69‑2 (MQ=255)                                           tGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACatcaat                          >  1:2548112/1‑71 (MQ=255)                                            gAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCGGCACCACGGTAATatc                             <  1:2017244/67‑1 (MQ=255)                                             aaTTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACatcaatca                        <  1:5275029/71‑1 (MQ=255)                                                tAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTtcatcat                                                >  1:3982781/1‑44 (MQ=255)                                                 aaCTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCg                     >  1:4759618/1‑70 (MQ=255)                                                  aCTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGa                   <  1:2503480/71‑1 (MQ=255)                                                    tCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGaa                  >  1:3936979/1‑70 (MQ=255)                                                     ccAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAAt                >  1:3397653/1‑71 (MQ=255)                                                      cagcagCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGaaa                 >  1:4642543/1‑69 (MQ=255)                                                      cagcagCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATg               >  1:3098297/1‑71 (MQ=255)                                                      cagcagCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATg               >  1:391025/1‑71 (MQ=255)                                                       agcagcGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATg               >  1:4353244/1‑70 (MQ=255)                                                         cagcGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCgg            <  1:4105215/71‑1 (MQ=255)                                                         cagcGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCg             >  1:4473905/1‑70 (MQ=255)                                                          agcGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGa           <  1:494097/71‑1 (MQ=255)                                                          agcGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGa           <  1:4718488/71‑1 (MQ=255)                                                          agcGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGa           <  1:3833914/71‑1 (MQ=255)                                                          agcGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGa           <  1:3949063/71‑1 (MQ=255)                                                            cgcCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAgg         <  1:4749342/71‑1 (MQ=255)                                                             gcCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGc        <  1:4308832/71‑1 (MQ=255)                                                                aTACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATa     <  1:5114754/71‑1 (MQ=255)                                                                aTACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATa     <  1:1096705/71‑1 (MQ=255)                                                                 tACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATaa    <  1:1019130/71‑1 (MQ=255)                                                                 tACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATaa    <  1:5646981/71‑1 (MQ=255)                                                                   cGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAg  >  1:1007875/1‑71 (MQ=255)                                                                   cGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAg  >  1:4786541/1‑71 (MQ=255)                                                                    |                                                                      TGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAG  >  minE/2677349‑2677485

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