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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,677,375	0	T	C	7.8%	250.7 / 9.8	90	Q185Q (CAA→CAG)	yhgA	predicted transposase
Reads supporting (aligned to +/- strand):  ref base T (43/40);  new base C (4/3);  total (47/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

TTGTCTGTCATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGG  >  minE/2677306‑2677443                                                                      |                                                                     ttGTCTGTCATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGt                                                                      >  1:2785377/1‑70 (MQ=255)  tgtctgtcATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGtt                                                                     >  1:2247712/1‑70 (MQ=255)   gtctgtcATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTg                                                                    <  1:3826493/70‑1 (MQ=255)   gtctgtcATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTg                                                                    <  1:3830991/70‑1 (MQ=255)      tgtcATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACg                                                                 >  1:2071921/1‑70 (MQ=255)      tgtcATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACg                                                                 >  1:1207380/1‑70 (MQ=255)        tcATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGa                                                                <  1:2272721/69‑1 (MQ=255)        tcATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAat                                                              <  1:452830/71‑1 (MQ=255)           ttAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATg                                                           <  1:3720596/71‑1 (MQ=255)             aGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGt                                                          <  1:289529/70‑1 (MQ=255)              gCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGtt                                                         <  1:1957323/70‑1 (MQ=255)                gTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGttttt                                                      <  1:1832287/71‑1 (MQ=255)                gTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGttttt                                                      <  1:3886274/71‑1 (MQ=255)                gTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGttttt                                                      <  1:2824569/71‑1 (MQ=255)                gTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGttttt                                                      <  1:3890914/71‑1 (MQ=255)                  aTCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGttt                                                        >  1:2200781/1‑67 (MQ=255)                  aTCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGttt                                                        <  1:3770208/67‑1 (MQ=255)                   tCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGaa                                                   >  1:1232981/1‑71 (MQ=255)                   tCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGa                                                    <  1:511991/70‑1 (MQ=255)                       ctactTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGaa                                                   >  1:1417220/1‑67 (MQ=255)                       ctactTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTaa                                               >  1:2723191/1‑71 (MQ=255)                        tactTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGAc                                                                  >  1:930502/1‑51 (MQ=255)                        tactTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAAc                                              >  1:447278/1‑71 (MQ=255)                        tactTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAAc                                              >  1:632769/1‑71 (MQ=255)                        tactTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTAGACGAATATGTTTTTGAATTa                                                >  1:3738356/1‑69 (MQ=255)                             aataaGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAAtt                                                 <  1:1032001/63‑1 (MQ=255)                             aataaGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTcc                                           <  1:1573181/69‑1 (MQ=255)                             aataaGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTcc                                           <  1:294496/69‑1 (MQ=255)                              ataaGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTcc                                           <  1:620420/68‑1 (MQ=255)                               taaGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGttt                                                        >  1:1955693/1‑54 (MQ=255)                               taaGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGttt                                                        >  1:2528655/1‑54 (MQ=255)                               taaGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCcag                                         >  1:1006126/1‑69 (MQ=255)                               taaGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCcag                                         >  1:1400013/1‑69 (MQ=255)                                 aGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCagcagc                                     >  1:659153/1‑71 (MQ=255)                                    aGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAgcgc                                   >  1:2557141/1‑70 (MQ=255)                                     ggCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCa                                 <  1:1225575/71‑1 (MQ=255)                                     ggCCATTTGCTATACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCa                                 >  1:2141876/1‑71 (MQ=255)                                      gCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCcagca                                       <  1:2668072/64‑1 (MQ=255)                                      gCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCa                                 >  1:2080422/1‑70 (MQ=255)                                      gCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCa                                 >  1:2771993/1‑70 (MQ=255)                                       ccATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGaa                                                               <  1:1324491/39‑1 (MQ=255)                                        cATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAAtt                                                 <  1:570166/52‑1 (MQ=255)                                        cATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATAc                              >  1:3814505/1‑71 (MQ=255)                                         aTTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAatat                                                            <  1:1439093/40‑1 (MQ=255)                                          tttGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAAt                                                  <  1:3856267/49‑1 (MQ=255)                                          tttGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAAt                                                  <  1:1020589/49‑1 (MQ=255)                                          tttGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGt                            <  1:2162599/71‑1 (MQ=255)                                          tttGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGt                            <  1:3210283/71‑1 (MQ=255)                                           ttGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTc                           >  1:2944899/1‑71 (MQ=255)                                              ctctACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGt                        <  1:700351/71‑1 (MQ=255)                                               tctACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTg                                                     <  1:1340377/41‑1 (MQ=255)                                               tctACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTcc                                           <  1:666806/51‑1 (MQ=255)                                                ctACCAGCCCCATCAGGTCGCGTTGACGAATATGttttt                                                      >  1:2213496/1‑39 (MQ=255)                                                ctACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTaa                                               >  1:2489470/1‑46 (MQ=255)                                                 tACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTc                           >  1:694205/1‑65 (MQ=255)                                                 tACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtgct                     >  1:99037/1‑71 (MQ=255)                                                 tACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtgct                     >  1:3896256/1‑71 (MQ=255)                                                 tACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtgct                     >  1:1822442/1‑71 (MQ=255)                                                 tACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtgct                     >  1:2217154/1‑71 (MQ=255)                                                 tACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtgc                      >  1:1251752/1‑70 (MQ=255)                                                   ccAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtgctg                    >  1:2229720/1‑70 (MQ=255)                                                   ccAGCCCCATCAGGTCGCGTTGACGAATATGATTTTGAATTAACTCCAGCAGCGCCATACGTCGGtgctg                    >  1:3393685/1‑70 (MQ=255)                                                    cAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCa                  <  1:3872632/71‑1 (MQ=255)                                                     aGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCg                          >  1:221171/1‑62 (MQ=255)                                                     aGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCAt                 >  1:269866/1‑71 (MQ=255)                                                     aGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCAt                 >  1:418647/1‑71 (MQ=255)                                                      gCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAAtt                                                 <  1:1573021/38‑1 (MQ=255)                                                      gCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtg                       <  1:229142/64‑1 (MQ=255)                                                      gCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCAt                 >  1:1292942/1‑70 (MQ=255)                                                      gCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATg                <  1:3542177/71‑1 (MQ=255)                                                      gCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATg                <  1:2960844/71‑1 (MQ=255)                                                      gCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATg                <  1:2803257/71‑1 (MQ=255)                                                       tcccAACAGATCGCGCTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCAtaa               <  1:3508644/70‑3 (MQ=255)                                                       ccccATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGa               >  1:3398659/1‑71 (MQ=255)                                                         ccAACAGATCGCGCTGACGAATATGTTTCTGAATTaa                                               >  1:22136/1‑37 (MQ=21)                                                         ccAACAGATCGCGCTGACGAATATGTTTCTGAATTaa                                               >  1:1909450/1‑37 (MQ=21)                                                         ccAACAGATCGCGCTGACGAATATGTTTCTGAATTaa                                               >  1:2166678/1‑37 (MQ=21)                                                         ccAACAGATCGCGCTGACGAATATGTTTCTGAATTaa                                               >  1:1024970/1‑37 (MQ=21)                                                         ccAACAGATCGCGCTGACGAATATGTTTCTGAATTAACTcc                                           <  1:605972/41‑1 (MQ=25)                                                         ccAACAGATCGCGCTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAAtc             <  1:1252999/71‑2 (MQ=255)                                                          cATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGAttt            >  1:3907307/1‑71 (MQ=255)                                                          cATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGAtt             >  1:1785837/1‑70 (MQ=255)                                                           aTCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTtc           <  1:2651186/71‑1 (MQ=255)                                                            tCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCcagcag                                      >  1:1485487/1‑43 (MQ=255)                                                            tCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGt                            <  1:1651933/53‑1 (MQ=255)                                                            tCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTtca          >  1:430257/1‑71 (MQ=255)                                                              aGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTtcat         <  1:3591437/70‑1 (MQ=255)                                                                 tCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtgct                     >  1:1834015/1‑55 (MQ=255)                                                                   gcgcTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCATAAtctc           >  1:1056936/1‑63 (MQ=255)                                                                   gcgTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAg   >  1:564050/1‑71 (MQ=255)                                                                    cgTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAgg  <  1:2194101/71‑1 (MQ=255)                                                                     gTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTtcat         >  1:878175/1‑63 (MQ=255)                                                                     gTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTtca          <  1:1240339/62‑1 (MQ=255)                                                                      ttGAGGAATATGTTTTTGAATTAACTCCAGCAgcgc                                   <  1:3678936/36‑1 (MQ=255)                                                                      |                                                                     TTGTCTGTCATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCATCATCAGG  >  minE/2677306‑2677443

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