Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,677,423	C→T	17.1%	Q169Q (CAG→CAA)	yhgA ←	predicted transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,677,423	0	C	T	17.1%	171.2 / 23.6	76	Q169Q (CAG→CAA)	yhgA	predicted transposase
Reads supporting (aligned to +/- strand):  ref base C (34/29);  new base T (2/11);  total (36/40)
Fisher's exact test for biased strand distribution p-value = 1.43e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTG  >  minE/2677356‑2677493                                                                    |                                                                       ccAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtgctg                                                                      >  1:2235050/1‑70 (MQ=255) ccAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtgctg                                                                      >  1:2848944/1‑70 (MQ=255) ccAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtgctg                                                                      >  1:15778/1‑70 (MQ=255)    atcccAACAGATCGCGCTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCAt                                                                   <  1:2602926/68‑1 (MQ=255)     tcccAACAGATCGCGCTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCAtaa                                                                 >  1:5350984/2‑69 (MQ=255)       ccATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGAtt                                                               <  1:4987377/71‑1 (MQ=255)       ccATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGAtt                                                               <  1:1087263/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:1993096/1‑70 (MQ=255)            aGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTtcat                                                           >  1:3514677/1‑70 (MQ=255)              gTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTtcatcat                                                        <  1:2660088/71‑1 (MQ=255)                 gcgcTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCtcgtcatcc                                                      <  1:4776253/70‑2 (MQ=255)                 gcgcTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCtcgtcatcc                                                      <  1:1036944/70‑2 (MQ=255)                 gcgcTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCtcgtcatcc                                                      <  1:4977092/70‑2 (MQ=255)                  cgctGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAATCTCGTCATCCgg                                                    <  1:5353350/68‑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:3992153/50‑2 (MQ=255)                         gAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCAta                                                                  <  1:3023090/50‑2 (MQ=255)                         gAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCAta                                                                  <  1:5350784/50‑2 (MQ=255)                           atatGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGt                                           <  1:2240814/71‑1 (MQ=255)                           atatGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGt                                           >  1:1891867/1‑71 (MQ=255)                           atatGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGt                                           >  1:5192344/1‑71 (MQ=255)                           atatGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGt                                           >  1:1906093/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:1823691/71‑1 (MQ=255)                              tGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGAc                                        <  1:389629/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:2246251/1‑70 (MQ=255)                                 ttttGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACat                                      >  1:2630057/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:4718488/71‑1 (MQ=255)                                                   agcGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGa                   <  1:494097/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:1096705/71‑1 (MQ=255)                                                         aTACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATa             <  1:5114754/71‑1 (MQ=255)                                                          tACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATaa            <  1:1019130/71‑1 (MQ=255)                                                          tACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATaa            <  1:5646981/71‑1 (MQ=255)                                                            cGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAg          >  1:4786541/1‑71 (MQ=255)                                                            cGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAg          >  1:1007875/1‑71 (MQ=255)                                                               cGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGaaa                         >  1:4045948/1‑53 (MQ=255)                                                                ggTGCTGCATGCTTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGc                      >  1:2603163/1‑55 (MQ=255)                                                                ggTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGc                      >  1:2956120/1‑55 (MQ=255)                                                                ggTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGc                      >  1:2020387/1‑55 (MQ=255)                                                                ggTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGc                      >  1:726427/1‑55 (MQ=255)                                                                ggTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGc      <  1:3930425/71‑1 (MQ=255)                                                                 gTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGcc     <  1:5686244/71‑1 (MQ=255)                                                                   gctgcATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCtt   >  1:1512572/1‑71 (MQ=255)                                                                   gctgcATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCtt   >  1:115905/1‑71 (MQ=255)                                                                    ctgcATGATTTCATCATCAGGCATGACGGTGACatca                                    <  1:1581167/37‑1 (MQ=255)                                                                    ctgcATGATTTCATCATCAGGCATGACGGTGACatca                                    <  1:1296667/37‑1 (MQ=255)                                                                    ctgcATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATaa            >  1:949821/1‑61 (MQ=255)                                                                    ctgcATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATaa            >  1:106189/1‑61 (MQ=255)                                                                    ctgcATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTg  >  1:4132067/1‑71 (MQ=255)                                                                    ctgcATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTg  >  1:4157304/1‑71 (MQ=255)                                                                    ctgcATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATACGGAGGCAt              >  1:4245261/1‑59 (MQ=255)                                                                    |                                                                       CCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTG  >  minE/2677356‑2677493

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