Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	3,998,703	A→G	31.4%	Y285H (TAC→CAC)	yhiP ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,998,703	0	A	G	31.4%	112.6 / 75.7	86	Y285H (TAC→CAC)	yhiP	predicted transporter
Reads supporting (aligned to +/- strand):  ref base A (30/29);  new base G (16/11);  total (46/40)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.26e-01

AAACCGAGAATTTCATGATGCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACAAACATTTTATTGCGCCCGGTTTT  >  W3110S.gb/3998636‑3998772                                                                    |                                                                      aaaCCGAGAATTTCATGATGCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTagag                                                                    <  1:3004505/71‑1 (MQ=255) aaaCCGAGAATTTCATGATGCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTagag                                                                    <  1:682202/71‑1 (MQ=255) aaaCCGAGAATTTCATGATGCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTagag                                                                    <  1:5122912/71‑1 (MQ=255)     cGAGAATTTCATGATGCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTGGAGAAt                                                                 >  1:2226457/1‑70 (MQ=255)      gagaATTTCATGATGCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTGGAGAATGt                                                               <  1:4148115/71‑1 (MQ=255)       agaATTTCATGATGCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGt                                                               >  1:281374/1‑70 (MQ=255)        gaATTTCATGATGCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTGGAGAATGTaa                                                             <  1:2062588/71‑1 (MQ=255)        gaATTTCATGATGCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTaa                                                             <  1:1229480/71‑1 (MQ=255)        gaATTTCATGATGCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTaa                                                             <  1:3688373/71‑1 (MQ=255)        gaATTTCATGATGCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTaa                                                             <  1:4116826/71‑1 (MQ=255)        gaATTTCATGATGCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTaa                                                             <  1:4852178/71‑1 (MQ=255)         aaTTTCATGATGCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTGGAGAATGTaaa                                                            <  1:2899569/71‑1 (MQ=255)         aaTTTCATGATGCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTGGAGAATGTaaa                                                            <  1:86680/71‑1 (MQ=255)             tCATGATGCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAacac                                                        >  1:1124322/1‑71 (MQ=255)             tCATGATGCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAacac                                                        >  1:2187631/1‑71 (MQ=255)              cATGATGCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTGGAGAATGTaaaa                                                           >  1:797802/1‑67 (MQ=255)              cATGATGCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTGGAGAATGTaaaa                                                           >  1:3465235/1‑67 (MQ=255)              cATGATGCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAacac                                                        <  1:3018332/70‑1 (MQ=255)                tgatgCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAAcactac                                                     >  1:3265434/1‑71 (MQ=255)                  atgCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACaccacc                                                    <  1:2845545/70‑1 (MQ=255)                   tgCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTaa                                                             <  1:3362737/60‑1 (MQ=255)                   tgCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTaa                                                             >  1:2679603/1‑60 (MQ=255)                    gCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTGGAGAATGTAAAACACCACCGCt                                                 >  1:4754944/1‑71 (MQ=255)                     cACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCtt                                                >  1:2744869/1‑71 (MQ=255)                     cACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCtt                                                >  1:4534864/1‑71 (MQ=255)                      aCGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTaaaa                                                           >  1:3438016/1‑59 (MQ=255)                      aCGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTc                                               >  1:1312923/1‑71 (MQ=255)                       cGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTGGAGAATGTAAAACACCACCGCTTc                                               >  1:421936/1‑70 (MQ=255)                       cGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTGGAGAATGTAAAACACCACCGCTTc                                               >  1:2034518/1‑70 (MQ=255)                       cGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTGGAGAATGTAAAACACCACCGCTTCg                                              <  1:1598897/71‑1 (MQ=255)                       cGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTGGAGAATGTAAAACACCACCGCTTCg                                              <  1:3224737/71‑1 (MQ=255)                        gttgttGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTGgaga                                                                   <  1:215788/49‑1 (MQ=255)                           gttgATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTGGAGAAt                                                                 >  1:2221614/1‑48 (MQ=255)                            ttgATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTGGAGAATGTAAAACACCACCGCTTCGAGCAt                                         >  1:4842461/1‑71 (MQ=255)                             tgATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCAt                                         >  1:2010927/1‑70 (MQ=255)                                tGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTGGAGAATGTaaaa                                                           >  1:588498/1‑49 (MQ=255)                                tGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTGGAGAATGTaaa                                                            >  1:4394793/1‑48 (MQ=255)                                tGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTGGAGAATGTAAAACACCACCGCTTCGAGCATCAg                                      <  1:2061400/70‑1 (MQ=255)                                tGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTGGAGAATGTAAAACACCACCGCTTCGAGCATCAg                                      <  1:432026/70‑1 (MQ=255)                                tGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAacac                                                        <  1:4476689/52‑1 (MQ=255)                                   cAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACa                                  <  1:4666483/71‑1 (MQ=255)                                   cAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAAAGTAAAACACCACCGCTTCGAGCATCAGGACa                                  <  1:1416119/71‑1 (MQ=255)                                    aaaGAAGTTCAGCGATGTTGGCATCTGGGCGTGGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACa                                  >  1:3616432/1‑70 (MQ=255)                                     aagaagTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACaaa                                <  1:1740254/71‑1 (MQ=255)                                     aagaagTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACaaa                                <  1:299103/71‑1 (MQ=255)                                     aagaagTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACaaa                                <  1:2803267/71‑1 (MQ=255)                                     aagaagTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACaaa                                <  1:3284490/71‑1 (MQ=255)                                      agaagTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAg                               >  1:573380/1‑71 (MQ=255)                                        aagTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGc                             <  1:2449118/71‑1 (MQ=255)                                        aagTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGc                             >  1:726337/1‑71 (MQ=255)                                            tCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAaca                                                         >  1:3244101/1‑39 (MQ=255)                                            tCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACGGCTTCGAGCATCAg                                      >  1:4915927/1‑58 (MQ=255)                                            tCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAg                                      >  1:2942688/1‑58 (MQ=255)                                             cAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACaccacc                                                    >  1:2178278/1‑43 (MQ=255)                                             cAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACaccacc                                                    >  1:4071081/1‑43 (MQ=255)                                             cAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACaacacc                                                    >  1:4759747/1‑43 (MQ=255)                                             cAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACaaa                                >  1:2614062/1‑63 (MQ=255)                                             cAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACaaa                                >  1:4814137/1‑63 (MQ=255)                                              aGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACaaa                       <  1:2247654/71‑1 (MQ=255)                                              aGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACaaa                       <  1:1511244/71‑1 (MQ=255)                                               gCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGGATCAGGACAAAGGCGACAAAc                      <  1:3181743/71‑1 (MQ=255)                                               gCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACAAAc                      <  1:1719095/71‑1 (MQ=255)                                                cGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACAAACa                     >  1:4105520/1‑71 (MQ=255)                                                cGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACAAACa                     >  1:4216217/1‑71 (MQ=255)                                                cGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACAAACa                     >  1:2571208/1‑71 (MQ=255)                                                  aTGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGa                                    <  1:1336704/54‑1 (MQ=255)                                                   tGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCt                                                 <  1:2298556/40‑1 (MQ=255)                                                     ttGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTc                                               >  1:200734/1‑40 (MQ=255)                                                     ttGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACAAACAtttt                 <  1:1925859/70‑1 (MQ=255)                                                     ttGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACAAACAtttt                 <  1:1715073/70‑1 (MQ=255)                                                       ggCATCTGGGCGTGGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACAAACAtttt                 <  1:3513416/68‑1 (MQ=255)                                                       ggCATCTGGGCGTGGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACAAACAtttt                 <  1:3461232/68‑1 (MQ=255)                                                       ggCATCTGGGCGTGGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACAAACATTTTa                >  1:362315/1‑69 (MQ=255)                                                         cATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACAAACATTTTATTgc            >  1:3275686/1‑71 (MQ=255)                                                           tCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCg                            <  1:857916/53‑1 (MQ=255)                                                             tGGGCGTGGAGAATGTAAAACACCACCGCTTCGAGc                                           >  1:2381044/1‑36 (MQ=255)                                                             tGGGCGTGGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACAAACa                     >  1:2860099/1‑58 (MQ=255)                                                             tGGGCGTGGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACAAACa                     >  1:1704676/1‑58 (MQ=255)                                                             tGGGCGTGGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACAAACa                     >  1:1872213/1‑58 (MQ=255)                                                             tGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACAAAc                      >  1:1323991/1‑57 (MQ=255)                                                             tGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACAAAc                      >  1:399526/1‑57 (MQ=255)                                                             tGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACAAACAttt                  <  1:4645315/61‑1 (MQ=255)                                                             tGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACAAACATTTTATTGCGccc        >  1:933964/1‑71 (MQ=255)                                                               ggCGTAGAGAATGTAAAACACCACCGCTTCGAGCTTCAGGACAAAGGCGACAAACATTTTATTGCGccc        <  1:1777724/69‑1 (MQ=255)                                                               ggCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACAAACATTTTATTGCGCCCg       <  1:3452445/70‑1 (MQ=255)                                                                gCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACAAACATTTTa                >  1:2073837/1‑60 (MQ=255)                                                                gCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACAAACATTTTa                >  1:3575202/1‑60 (MQ=255)                                                                   tAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACa                                  <  1:2566653/39‑1 (MQ=255)                                                                   tAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACAAACATTTTATTGCGCCCGGtttt  <  1:3060828/71‑1 (MQ=255)                                                                    agagAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGc                             <  1:4705251/43‑1 (MQ=255)                                                                    agagAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGc                             <  1:4569407/43‑1 (MQ=255)                                                                    |                                                                      AAACCGAGAATTTCATGATGCACGTTGTTGATGGCAAAGAAGTTCAGCGATGTTGGCATCTGGGCGTAGAGAATGTAAAACACCACCGCTTCGAGCATCAGGACAAAGGCGACAAACATTTTATTGCGCCCGGTTTT  >  W3110S.gb/3998636‑3998772

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