Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,677,369	G→A	7.7%	D187D (GAC→GAT)	yhgA ←	predicted transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,677,369	0	G	A	7.7%	278.2 / 9.2	91	D187D (GAC→GAT)	yhgA	predicted transposase
Reads supporting (aligned to +/- strand):  ref base G (45/39);  new base A (4/3);  total (49/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

TTTGATTTGTCTGTCATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCAT  >  minE/2677300‑2677436                                                                      |                                                                    tttGATTTGTCTGTCATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGt                                                                    >  1:2225561/1‑71 (MQ=255)   tGATTTGTCTGTCATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTcg                                                                  >  1:2652157/1‑71 (MQ=255)   tGATTTGTCTGTCATTAGCGTACCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTcg                                                                  >  1:3460570/1‑71 (MQ=255)     aTTTGTCTGTCATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTcgcg                                                                <  1:1070768/71‑1 (MQ=255)       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:2824569/71‑1 (MQ=255)                      gTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGttttt                                               <  1:3890914/71‑1 (MQ=255)                      gTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGttttt                                               <  1:3886274/71‑1 (MQ=255)                      gTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGttttt                                               <  1:1832287/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:3210283/71‑1 (MQ=255)                                                tttGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGt                     <  1:2162599/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:3896256/1‑71 (MQ=255)                                                       tACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGtgct              >  1:99037/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:418647/1‑71 (MQ=255)                                                           aGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCAt          >  1:269866/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:2803257/71‑1 (MQ=255)                                                            gCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATg         <  1:2960844/71‑1 (MQ=255)                                                            gCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATg         <  1:3542177/71‑1 (MQ=255)                                                             tcccAACAGATCGCGCTGACGAATATGTTTCTGAATTAACTCCAACAGCGCCATTTTGCGGTGTTGCAtaa        <  1:3508644/70‑3 (MQ=255)                                                             ccccATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGa        >  1:3398659/1‑71 (MQ=255)                                                               ccAACAGATCGCGCTGACGAATATGTTTCTGAATTaa                                        >  1:1024970/1‑37 (MQ=21)                                                               ccAACAGATCGCGCTGACGAATATGTTTCTGAATTaa                                        >  1:1909450/1‑37 (MQ=21)                                                               ccAACAGATCGCGCTGACGAATATGTTTCTGAATTaa                                        >  1:2166678/1‑37 (MQ=21)                                                               ccAACAGATCGCGCTGACGAATATGTTTCTGAATTaa                                        >  1:22136/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)                                                                      |                                                                    TTTGATTTGTCTGTCATTAGCGTATCCACTACTTAATAAGCAGGCCATTTGCTCTACCAGCCCCATCAGGTCGCGTTGACGAATATGTTTTTGAATTAACTCCAGCAGCGCCATACGTCGGTGCTGCATGATTTCAT  >  minE/2677300‑2677436

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