Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	3,998,196	A→G	35.1%	Y454H (TAT→CAT)	yhiP ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,998,196	0	A	G	35.1%	73.3 / 73.3	74	Y454H (TAT→CAT)	yhiP	predicted transporter
Reads supporting (aligned to +/- strand):  ref base A (23/25);  new base G (23/3);  total (46/28)
Fisher's exact test for biased strand distribution p-value = 8.72e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

ACCATCAACAGCATCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGC  >  W3110S.gb/3998129‑3998261                                                                    |                                                                  aCCATCAACAGCATCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGAc                                                                <  1:1981696/71‑1 (MQ=255)  ccATCAACAGCATCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGtata                                                                   <  1:2840259/67‑1 (MQ=255)  ccATCAACAGCATCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATGGa                                                                 >  1:2889281/1‑69 (MQ=255)  ccATCAACAGCATCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATGGa                                                                 >  1:227131/1‑69 (MQ=255)  caaTCAACAGCATCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATGGa                                                                 >  1:967920/3‑69 (MQ=255)    aTCAACAGCATCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACggg                                                             >  1:2413364/1‑71 (MQ=255)    aTCAACAGCATCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACggg                                                             >  1:3665765/1‑71 (MQ=255)     tCAACAGCATCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGc                                                            <  1:1996529/71‑1 (MQ=255)     tCAACAGCATCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGc                                                            <  1:1165486/71‑1 (MQ=255)     tCAACAGCATCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGc                                                            <  1:3206581/71‑1 (MQ=255)     tCAACAGCATCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGc                                                            <  1:1486286/71‑1 (MQ=255)      cAACAGCATCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACATGTTGGTATAGACGGGc                                                            <  1:2740796/70‑1 (MQ=255)       aaCAGCATCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATGGACGGGCaa                                                          >  1:948424/1‑71 (MQ=255)       aaCAGCATCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACAAAACACGTTGGTATGGACGGGCaa                                                          >  1:810286/1‑71 (MQ=255)             aTCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGtata                                                                   >  1:2131422/1‑56 (MQ=255)             aTCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGtata                                                                   >  1:3474942/1‑56 (MQ=255)             aTCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGtata                                                                   >  1:164189/1‑56 (MQ=255)             aTCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGtata                                                                   >  1:3760855/1‑56 (MQ=255)             aTCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGtctc                                                    <  1:1045468/71‑1 (MQ=255)             aTCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGtctc                                                    <  1:1069248/71‑1 (MQ=255)             aTCACTAATGCAACGCCCAGCGTGACCAGCCCAATCTTACCAAACACGTTGGtata                                                                   >  1:2516754/1‑56 (MQ=255)              tCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATGGACGGGCAACGtctc                                                    >  1:1366627/1‑70 (MQ=255)              tCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATGGACGGGCAACGtctc                                                    >  1:2347767/1‑70 (MQ=255)              tCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATGGACGGGCAACGtctc                                                    >  1:1307738/1‑70 (MQ=255)              tCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATGGACGGGCAACGtctc                                                    >  1:3657632/1‑70 (MQ=255)              tCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCa                                                   <  1:1894251/71‑1 (MQ=255)                actactGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATGGACGGGCAACGTCTCAAg                                                 >  1:1840906/1‑71 (MQ=255)                actactGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAg                                                 >  1:2429573/1‑71 (MQ=255)                actactGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAg                                                 >  1:3192188/1‑71 (MQ=255)                actactGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAg                                                 >  1:42923/1‑71 (MQ=255)                    ctGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACTGGCAACGTCTCAAGCgg                                              <  1:1800599/70‑1 (MQ=255)                    ctGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGt                                                       >  1:1188598/1‑61 (MQ=255)                       cAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATGGACGGGCAACGTCTCAAGCGGATc                                           <  1:194004/70‑1 (MQ=255)                        aaCGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATGGACGGGCAACGTCTCAAGCGGATc                                           <  1:1103093/69‑1 (MQ=255)                         aCGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGt                                        <  1:3643393/71‑1 (MQ=255)                         aCGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGt                                        <  1:1834625/71‑1 (MQ=255)                           gcacAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTaa                                      <  1:1733515/68‑1 (MQ=255)                           gCCCAGCGTGACCAGACCTATCTTACCAAACACGTTGGTATGGACGGGCAACGTCTCAAGCGGATCGGTaa                                      >  1:1383211/1‑71 (MQ=255)                           gCCCAGCGTGACCAGACCAATCTTACCAAACACTTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTaa                                      <  1:2846625/71‑1 (MQ=255)                           gCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATGGACGGGCAACGTCTCAAGCGGATCGGTaa                                      >  1:3997997/1‑71 (MQ=255)                           gCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATGGACGGGCAACGTCTCAAGCGGATCGGTaa                                      >  1:358935/1‑71 (MQ=255)                               aGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAAtgtt                                  >  1:1938883/1‑71 (MQ=255)                                  gTGACCAGACCAATCTTACCAAACACGTTGGTATAGa                                                                 >  1:1913742/1‑37 (MQ=255)                                  gTGACCAGACCAATCTTACCAAACACGTTGGTATAGa                                                                 >  1:2461097/1‑37 (MQ=255)                                  gTGACCAGACCAATCTTACCAAACACGTTGGTATAGa                                                                 >  1:3097037/1‑37 (MQ=255)                                  gTGACCAGACCAATCTTACCAAACACGTTGGTATAGa                                                                 >  1:2146169/1‑37 (MQ=255)                                  gTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTc                               >  1:1077529/1‑71 (MQ=255)                                   tGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTcc                              <  1:3336883/71‑1 (MQ=255)                                     accagaccaATCTTACCAAACACGTTGGTATGGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCgg                            >  1:1068363/1‑71 (MQ=255)                                     accagaccaATCTTACCAAACACGTTGGTATGGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCgg                            >  1:2377206/1‑71 (MQ=255)                                     accagaccaATCTTACCAAACACGTTGGTATGGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCgg                            >  1:2405196/1‑71 (MQ=255)                                     accagaccaATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAAtgttg                                 <  1:3643301/66‑1 (MQ=255)                                     accagaccaATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCgg                            <  1:2576250/71‑1 (MQ=255)                                      ccagaccaATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGc                           <  1:1297987/71‑1 (MQ=255)                                           ccaATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCAc                         >  1:2965256/1‑68 (MQ=255)                                           ccaATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACcgc                      >  1:2869234/1‑71 (MQ=255)                                           ccaATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACcgc                      >  1:66858/1‑71 (MQ=255)                                              aTCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAAtg                                    >  1:307376/1‑54 (MQ=255)                                              aTCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCgg                    >  1:686741/1‑70 (MQ=255)                                              aTCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCgg                    >  1:478073/1‑70 (MQ=255)                                              aTCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCgg                    >  1:444934/1‑70 (MQ=255)                                                 ttACCAAACACGTTGGTATGGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTaaa                >  1:1491647/1‑71 (MQ=255)                                                 ttACCAAACACGTTGGTATGGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTaaa                >  1:1321450/1‑71 (MQ=255)                                                 ttACCAAACACGTTGGTATAGACGGGCAACGTCTCAAg                                                 >  1:2187536/1‑38 (MQ=255)                                                 ttACCAAACACGTTGGTATAGACGGGCAACGTCTCAAg                                                 >  1:2139188/1‑38 (MQ=255)                                                   aCCAAACACGTTGGTATGGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATg              >  1:859649/1‑71 (MQ=255)                                                   aCCAAACACGTTGGTATGGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATg              >  1:2075744/1‑71 (MQ=255)                                                   aCCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTaaa                <  1:400473/69‑1 (MQ=255)                                                   aCCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTa                  <  1:130774/67‑1 (MQ=255)                                                   aCCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAAt               <  1:2581956/70‑1 (MQ=255)                                                    ccAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGt             <  1:3945610/71‑1 (MQ=255)                                                     cAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGtt            <  1:1944796/71‑1 (MQ=255)                                                         cacGTTGGTATGGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAAt               >  1:3009679/1‑64 (MQ=255)                                                         cacGTTGGTATGGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAAt               >  1:2931770/1‑64 (MQ=255)                                                         cacGTTGGTATGGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCca        <  1:2177785/71‑1 (MQ=255)                                                         cacGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCgg                            >  1:3241997/1‑51 (MQ=255)                                                         cacGTTGGTATAGACGGGCAACGTCTCAAGAGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCca        >  1:4074426/1‑71 (MQ=255)                                                          acGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCcac       <  1:1155521/71‑1 (MQ=255)                                                             ttGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTaaa                <  1:1244511/59‑1 (MQ=255)                                                               ggTATGGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAgc  >  1:1272586/1‑71 (MQ=255)                                                                    |                                                                  ACCATCAACAGCATCACTACTGCAACGCCCAGCGTGACCAGACCAATCTTACCAAACACGTTGGTATAGACGGGCAACGTCTCAAGCGGATCGGTAATGTTGTCCGGCACCGCGGTAAATGTTGCCACATAGC  >  W3110S.gb/3998129‑3998261

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