Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	4,096,796	G→A	6.4%	S152S (TCC→TCT)	yhgA ←	predicted transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,096,796	0	G	A	6.4%	286.2 / 5.9	94	S152S (TCC→TCT)	yhgA	predicted transposase
Reads supporting (aligned to +/- strand):  ref base G (34/53);  new base A (3/3);  total (38/56)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.53e-01

GCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGAAACAATCCAGCCAACACAGCGAATAGGGAT  >  W3110S.gb/4096731‑4096863                                                                  |                                                                    gcCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAgg                                                                 >  1:1964311/1‑70 (MQ=255)    aTACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCAt                                                              >  1:2389438/1‑70 (MQ=255)     tACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATaa                                                            <  1:580108/71‑1 (MQ=255)      aCGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATa                                                             >  1:1798204/1‑69 (MQ=255)       cGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGa                                                                   >  1:2735718/1‑62 (MQ=255)       cGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCa                                                               >  1:2763772/1‑66 (MQ=255)       cGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCa                                                               >  1:3128310/1‑66 (MQ=255)       cGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAg                                                          >  1:95725/1‑71 (MQ=255)           ggTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTg                                                       <  1:2912639/70‑1 (MQ=255)           ggTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGc                                                      <  1:3396722/71‑1 (MQ=255)                tgcATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATaaa                                                           <  1:180619/61‑1 (MQ=255)                  cATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCTGAGGCATAAAGCTGCCTTGc                                                 >  1:3366495/1‑69 (MQ=255)                  cATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTGc                                                 <  1:2328146/69‑1 (MQ=255)                  cATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTGc                                                 <  1:1468613/69‑1 (MQ=255)                  cATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTGc                                                 <  1:3477288/69‑1 (MQ=255)                    tGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATa                                                             <  1:2246213/55‑1 (MQ=255)                    tGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATa                                                             <  1:3837275/55‑1 (MQ=255)                      aTTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGttt                                           <  1:3426420/71‑1 (MQ=255)                       tttCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCg                                                                     <  1:2320029/44‑1 (MQ=255)                       tttCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCg                                                                     <  1:3144362/44‑1 (MQ=255)                       tttCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCg                                                                     <  1:3089971/44‑1 (MQ=255)                       tttCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCg                                                                     <  1:453477/44‑1 (MQ=255)                       tttCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCa                                                               <  1:1032895/50‑1 (MQ=255)                       tttCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCa                                                               <  1:1238496/50‑1 (MQ=255)                       tttCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCa                                                               <  1:1242636/50‑1 (MQ=255)                       tttCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTg                                          >  1:1269697/1‑71 (MQ=255)                       tttCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTg                                          >  1:3306439/1‑71 (MQ=255)                         tcatcatcaGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATa                                                             >  1:3473147/1‑50 (MQ=255)                         tcatcatcaGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGttt                                           >  1:3662852/1‑68 (MQ=255)                         tcatcatcaGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGttt                                           <  1:827871/68‑1 (MQ=255)                         tcatcatcaGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTgg                                         >  1:3052541/1‑70 (MQ=255)                         tcatcatcaGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTgg                                         <  1:2570192/70‑1 (MQ=255)                          catcatcaGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTggg                                        <  1:1193868/70‑1 (MQ=255)                           atcatcaGGCATGACGGTGACATCAATCAGCGGAAATGCg                                                                     <  1:553883/40‑1 (MQ=255)                              atcaGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCg                                    <  1:1443363/70‑1 (MQ=255)                              atcaGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCg                                    <  1:1073346/70‑1 (MQ=255)                                 aGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGa                                <  1:590938/71‑1 (MQ=255)                                 aGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGa                                <  1:1769061/71‑1 (MQ=255)                                   gCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGt                                             <  1:1755955/56‑1 (MQ=255)                                     aTGACTGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGAAAc                             >  1:971665/1‑70 (MQ=255)                                      tGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTggg                                        <  1:403558/58‑1 (MQ=255)                                      tGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTggg                                        <  1:3337455/58‑1 (MQ=255)                                      tGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTggg                                        <  1:483909/58‑1 (MQ=255)                                      tGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGa                                <  1:741370/66‑1 (MQ=255)                                      tGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGa                                <  1:1095146/66‑1 (MQ=255)                                        aCGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGcc                                                     <  1:177733/43‑1 (MQ=255)                                        aCGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGaaa                              >  1:2537563/1‑66 (MQ=255)                                        aCGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGAAACAATc                         <  1:520567/71‑1 (MQ=255)                                        aCGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGAAACAATc                         >  1:991999/1‑71 (MQ=255)                                        aCGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGAAACAATc                         <  1:3047554/71‑1 (MQ=255)                                        aCGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGAAACAATc                         >  1:2491579/1‑71 (MQ=255)                                        aCGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTCCCTTGCCAg                                              >  1:2408021/1‑50 (MQ=255)                                          ggTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGAAACAATc                         >  1:1433533/1‑69 (MQ=255)                                           gTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGAAACAATCCAg                      <  1:2039260/71‑1 (MQ=255)                                           gTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGAAACAATCCAg                      <  1:703227/71‑1 (MQ=255)                                            tGACATCAATCAGCGGAAATGCAGAGGCATAAAGCTGCCTTGc                                                 <  1:2832599/43‑1 (MQ=255)                                            tGACATCAATCAGCGGAAATGCAGAGGCATAAAGCTGCCTTGc                                                 <  1:2149069/43‑1 (MQ=255)                                            tGACATCAATCAGCGGAAATGCAGAGGCATAAAGCTGCCTTGc                                                 <  1:1642034/43‑1 (MQ=255)                                             gACATCAATCAGCGGAAATGCAGAGGCATGAAGCTGCCTTGCCAGGTTGGGATCGGCGAAACAATCCa                       >  1:912779/1‑68 (MQ=255)                                             gACATCAATCAGCGGAAATGCAGAGGCATAAAGCTGCCTTGCCAGGTTGGGATCGGCGAAACAATCCa                       >  1:522884/1‑68 (MQ=255)                                             gACATCAATCAGCGGAAATGCAGAGGCATAAAGCTGCCTTGCCAGGTTGGGATCGGCGAAACAATCCa                       >  1:3803931/1‑68 (MQ=255)                                              aCATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTGc                                                 <  1:2231714/41‑1 (MQ=255)                                              aCATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTgg                                         >  1:3013177/1‑49 (MQ=255)                                              aCATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTgg                                         >  1:214332/1‑49 (MQ=255)                                               cATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGAAACAATCCAGCCaa                  <  1:461153/71‑1 (MQ=255)                                                 tcaatcaGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGtt                                            <  1:3734313/43‑1 (MQ=255)                                                 tcaatcaGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGtt                                            <  1:3738962/43‑1 (MQ=255)                                                  caatcaGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGAAACAAt                          >  1:3713377/1‑60 (MQ=255)                                                  caatcaGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGAAACAAt                          >  1:1268152/1‑60 (MQ=255)                                                   aatcaGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTg                                          >  1:188553/1‑43 (MQ=255)                                                   aatcaGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGa                                <  1:761008/53‑1 (MQ=255)                                                   aatcaGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGa                                <  1:1645600/53‑1 (MQ=255)                                                   aatcaGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGAAACAATCCAGCCAaca                <  1:2392937/69‑1 (MQ=255)                                                   aatcaGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGAAACAATCCAGCCAaca                <  1:24060/69‑1 (MQ=255)                                                   aatcaGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGAAACAATCCAGCCAaca                >  1:3546157/1‑69 (MQ=255)                                                   aatcaGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGAAACAATCCAGCCAaca                >  1:2629512/1‑69 (MQ=255)                                                   aatcaGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGAAACAATCCAGCCAAcaca              <  1:2882765/71‑1 (MQ=255)                                                       aGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGAAACAATCCAGCCAACACAGCg           >  1:894199/1‑70 (MQ=255)                                                       aGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGAAACAATCCAGCCAACACAGCg           >  1:1128733/1‑70 (MQ=255)                                                       aGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGAAACAATCCAGCCAACACAGCg           >  1:1946339/1‑70 (MQ=255)                                                       aGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGAAACAATCCAGCCAACACAGCg           >  1:2890539/1‑70 (MQ=255)                                                       aGCGG‑AATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGAAACAATCCAGCCAAAACAGCGaa         >  1:2271989/1‑71 (MQ=255)                                                         cGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGAAACAATCCAGCCAACACAGCGaa         <  1:3416936/70‑1 (MQ=255)                                                         cGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGAAACAATCCAGCCAACACAGCGAat        <  1:1974876/71‑1 (MQ=255)                                                         cGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGAAACAATCCAGCCAACACAGCGAat        <  1:734675/71‑1 (MQ=255)                                                          ggAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATc                                     <  1:103069/41‑1 (MQ=255)                                                          ggAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATc                                     <  1:2885044/41‑1 (MQ=255)                                                          ggAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATc                                     <  1:2973656/41‑1 (MQ=255)                                                          ggAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATc                                     <  1:566245/41‑1 (MQ=255)                                                          ggAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGAAACAATc                         >  1:312200/1‑53 (MQ=255)                                                           gAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGAAACAATCCAGCCAACACAGCGAat        >  1:442464/1‑69 (MQ=255)                                                            aaaTGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGAAACAATCCAGCCAACACAGCGAatagg     <  1:916454/71‑1 (MQ=255)                                                            aaaTGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGAAACAATCCAGCCAACACAGCGAatagg     <  1:1984400/71‑1 (MQ=255)                                                             aaTGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGAAACAATCCAGCCAACACAGCGAata       >  1:350882/1‑68 (MQ=255)                                                             aaTGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGAAACAATCCAGCCAACACAGCGAata       >  1:24752/1‑68 (MQ=255)                                                              aTGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGAAACAATCCAGCCAACACAGCGAataggga   <  1:2347946/71‑1 (MQ=255)                                                               tGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGAAACAATCCAGCCAACACAGCg           <  1:2470520/62‑1 (MQ=255)                                                               tGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGAAACAATCCAGCCAACACAGCGAatagggat  <  1:1849802/71‑1 (MQ=255)                                                                 cGGAGTCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGAAACAATCCAGCCAACACAGCGAat        <  1:2057749/63‑1 (MQ=255)                                                                  |                                                                    GCCATACGTCGGTGCTGCATGATTTCATCATCAGGCATGACGGTGACATCAATCAGCGGAAATGCGGAGGCATAAAGCTGCCTTGCCAGTTTGGGATCGGCGAAACAATCCAGCCAACACAGCGAATAGGGAT  >  W3110S.gb/4096731‑4096863

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