Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,677,183	T→C	5.8%	R249R (AGA→AGG) ‡	yhgA ←	predicted transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,677,183	0	T	C	5.8%	249.7 / 4.2	86	R249R (AGA→AGG) ‡	yhgA	predicted transposase
Reads supporting (aligned to +/- strand):  ref base T (49/32);  new base C (2/3);  total (51/35)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

AGGAACTCCGGATTCAAGCATTATTTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCGT  >  Exported/2677114‑2677253                                                                      |                                                                       agGAACTCCGGATTCAAGCATTATTTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATc                                                                       >  1:4063120/1‑71 (MQ=255) agGAACTCCGGATTCAAGCATTATTTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATc                                                                       >  1:3406544/1‑71 (MQ=255)        ccGGATTCAAGCATTATTTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCg                                                                 >  1:2227312/1‑70 (MQ=255)        ccGGATTCAAGCATTATTTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCg                                                                 >  1:2233742/1‑70 (MQ=255)        ccGGATTCAAGCATTATTTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCg                                                                 >  1:4432514/1‑70 (MQ=255)         cGGATTCAAGCATTATTTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATc                                                                       >  1:469898/1‑63 (MQ=255)         cGGATTCAAGCATTATTTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATc                                                                       >  1:3592952/1‑63 (MQ=255)         cGGATTCAAGCATTATTTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGc                                                                >  1:3978426/1‑70 (MQ=255)          ggATTCAAGCATTATTTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCaa                                                              <  1:2413551/71‑1 (MQ=255)              tCAAGCATTATTTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTc                                                                   >  1:477769/1‑62 (MQ=255)              tCAAGCATTATTTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGt                                                          >  1:4851246/1‑71 (MQ=255)                  gCATTATTTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCAt                                                       >  1:1679226/1‑70 (MQ=255)                   cATTATTTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTa                                                     <  1:3722993/71‑1 (MQ=255)                    attattTTGGCTATATGCAGGGC‑TTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTaaa                                                   >  1:1993270/1‑71 (MQ=255)                     ttattTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTaaa                                                   <  1:3674607/71‑1 (MQ=255)                     ttattTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTaaa                                                   <  1:2999804/71‑1 (MQ=255)                       attTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATc                                                                       >  1:497589/1‑49 (MQ=255)                       attTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAAct                                                 >  1:1442851/1‑71 (MQ=255)                       attTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAAc                                                  >  1:1614431/1‑70 (MQ=255)                          ttGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAActct                                               <  1:4030237/70‑1 (MQ=255)                             gCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCtt                                           <  1:2238824/71‑1 (MQ=255)                                atatGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCaa                                                              >  1:2286552/1‑49 (MQ=255)                                atatGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTtgtg                                        >  1:3018182/1‑71 (MQ=255)                                 tatGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCtt                                                                     <  1:865956/41‑1 (MQ=255)                                 tatGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTgtgt                                       >  1:2640550/1‑71 (MQ=255)                                   tGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAAt                                                                        <  1:985862/36‑1 (MQ=255)                                   tGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTa                                                     >  1:663721/1‑55 (MQ=255)                                   tGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAActc                                                <  1:2322354/60‑1 (MQ=255)                                   tGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGttt                                     >  1:5674497/1‑71 (MQ=255)                                   tGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGttt                                     >  1:458701/1‑71 (MQ=255)                                   tGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGttt                                     >  1:2039009/1‑71 (MQ=255)                                    gCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCttt                                                                    <  1:4568646/39‑1 (MQ=255)                                    gCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTc                                    >  1:2855571/1‑71 (MQ=255)                                     cAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCg                                   >  1:2114348/1‑71 (MQ=255)                                      aGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTcc                                                                  >  1:5476994/1‑39 (MQ=255)                                      aGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCgg                                  <  1:3123443/71‑1 (MQ=255)                                       gggCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTc                                                         <  1:1161130/47‑1 (MQ=255)                                       gggCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGttt                                     <  1:4605188/67‑1 (MQ=255)                                        ggCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGt                                 >  1:2284974/1‑70 (MQ=255)                                          cTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGTAATAGTCATTAAACTCTCCTTGTGTTTCGGTGaa                              >  1:665471/1‑71 (MQ=255)                                          cTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTaaa                                                   >  1:3730373/1‑50 (MQ=255)                                          cTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTg                                <  1:2654664/69‑1 (MQ=255)                                          cTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGaa                              >  1:5566757/1‑71 (MQ=255)                                           tttGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTa                                                     <  1:1186435/47‑1 (MQ=255)                                            ttGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGaa                              >  1:1695325/1‑69 (MQ=255)                                             tGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGt                           >  1:2339278/1‑71 (MQ=255)                                             tGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGt                           >  1:3608937/1‑71 (MQ=255)                                             tGGATTGTTCCACCTCCTGCCGCAATATTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGt                           >  1:3721111/1‑71 (MQ=255)                                               gATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAg                                                           >  1:1901807/1‑37 (MQ=255)                                               gATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGtt                          <  1:2158191/70‑1 (MQ=255)                                               gATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTc                         >  1:1466136/1‑71 (MQ=255)                                                 ttGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCa                                                        <  1:1809776/38‑1 (MQ=255)                                                 ttGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCgg                       <  1:3575475/71‑1 (MQ=255)                                                 ttGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCgg                       <  1:2933328/71‑1 (MQ=255)                                                  tGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTaaa                                                   >  1:1461287/1‑42 (MQ=255)                                                   gTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCAtt                                                      >  1:4177365/1‑38 (MQ=255)                                                   gTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTa                                                     <  1:4398553/39‑1 (MQ=255)                                                   gTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGc                      <  1:2483895/70‑1 (MQ=255)                                                   gTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGc                      <  1:1473989/70‑1 (MQ=255)                                                    ttCCCACTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCa                     <  1:1604583/70‑1 (MQ=255)                                                     tCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTaaa                                                   >  1:2550433/1‑39 (MQ=255)                                                     tCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCg                        >  1:4082537/1‑66 (MQ=255)                                                     tCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCg                        >  1:783954/1‑66 (MQ=255)                                                     tCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCaa                    >  1:5549790/1‑70 (MQ=255)                                                     tCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAAc                   >  1:3991758/1‑71 (MQ=255)                                                     tCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAAc                   >  1:2568676/1‑71 (MQ=255)                                                       ccccTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACg                  <  1:1969183/70‑1 (MQ=255)                                                        cccTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTACTTGTGTTTCGGTg                                <  1:4651551/55‑1 (MQ=255)                                                          ttcctGCCGTAGCCGTTCAGCAATAGTCATTAATCTCTCCTtttg                                        <  1:4111483/44‑1 (MQ=12)                                                          ttcctGCCGTAGCCGTTCAGCAATAGTCATTAATCTCTCCTTttgta                                      >  1:3008479/2‑46 (MQ=12)                                                          ctcctGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCaa                    <  1:2492474/65‑1 (MQ=255)                                                          ctcctGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCaa                    <  1:4730658/65‑1 (MQ=255)                                                          ctcctGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCaa                    <  1:1389454/65‑1 (MQ=255)                                                          ctcctGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCaa                    <  1:3778954/65‑1 (MQ=255)                                                          ctcctGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGt              >  1:1520338/1‑71 (MQ=255)                                                           tcctGCCGTAGCCGTTCAGCAATAGTCATTAATCTCTCCTTttgt                                       <  1:5020024/45‑1 (MQ=12)                                                           tcctGCCGTAGCCGTTCAGCAATAGTCATTAATCTCTCCTTttgt                                       <  1:1771849/45‑1 (MQ=12)                                                           tcctGCCGCAGCCGTTCAGCAATAGTCATTAATCTCTCCt                                            >  1:2142348/1‑40 (MQ=16)                                                               gCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGAt          >  1:5014097/1‑70 (MQ=255)                                                               gCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGAt          >  1:3202811/1‑70 (MQ=255)                                                                ccGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTgtgt                                       <  1:1110434/40‑1 (MQ=255)                                                                ccGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATaa        <  1:5446177/71‑1 (MQ=255)                                                                  gCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATaaaa      >  1:1918759/1‑71 (MQ=255)                                                                    aaTCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATc    <  1:2978835/71‑1 (MQ=255)                                                                    aaTCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATc    <  1:1136212/71‑1 (MQ=255)                                                                     aTCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCg   >  1:4838314/1‑71 (MQ=255)                                                                     aTCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCg   >  1:1451415/1‑71 (MQ=255)                                                                      tCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTg                                <  1:1744585/41‑1 (MQ=255)                                                                      tCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATaaaa      <  1:4637351/67‑1 (MQ=255)                                                                      tCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCGt  >  1:799715/1‑71 (MQ=255)                                                                      |                                                                       AGGAACTCCGGATTCAAGCATTATTTTGGCTATATGCAGGGCTTTGGATTGTTCCCCCTCCTGCCGCAATCTTTCCGCAATAGTCATTAAACTCTCCTTGTGTTTCGGTGAACGTTCGGCAACGCCGTCGATAAAATCGT  >  Exported/2677114‑2677253

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑