Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	2,251,587	A→G	13.7%	F94L (TTT→CTT)	lysP ←	lysine transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,251,587	0	A	G	13.7%	207.6 / 24.2	80	F94L (TTT→CTT)	lysP	lysine transporter
Reads supporting (aligned to +/- strand):  ref base A (27/42);  new base G (6/5);  total (33/47)
Fisher's exact test for biased strand distribution p-value = 3.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

CAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCA  >  W3110S.gb/2251521‑2251657                                                                   |                                                                       cAACCAGGTCAACGGCGATAGTCACCGCCGAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCct                                                                    <  1:2139096/71‑1 (MQ=255) cAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGc                                                                      >  1:978979/1‑69 (MQ=255)    ccAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCaaaa                                                                        <  1:4519824/64‑1 (MQ=255)      aGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAg                                                                       >  1:1477856/1‑63 (MQ=255)      aGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGc                                                                      >  1:3158888/1‑64 (MQ=255)      aGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCctt                                                                   >  1:1575641/1‑67 (MQ=255)        gTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCaa                                                             <  1:1797575/71‑1 (MQ=255)        gTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCaa                                                             <  1:32878/71‑1 (MQ=255)        gTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCaa                                                             <  1:2313114/71‑1 (MQ=255)          cAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAGGCcttctt                                                                >  1:4791053/1‑66 (MQ=255)          cAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAAc                                                            >  1:3053007/1‑70 (MQ=255)          cAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACa                                                           >  1:4752720/1‑71 (MQ=255)           aaCGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACAt                                                          <  1:4908430/71‑1 (MQ=255)               gCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAGGCCTTCTTCAACATAGtt                                                      <  1:5011853/71‑1 (MQ=255)               gCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAGGCCTTCTTCAACATAGtt                                                      <  1:4921105/71‑1 (MQ=255)               gCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAGGCCTTCTTCAACATAGtt                                                      <  1:957125/71‑1 (MQ=255)                cGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTc                                                     <  1:777050/71‑1 (MQ=255)                cGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTc                                                     <  1:4271589/71‑1 (MQ=255)                 gATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCt                                                    <  1:135498/71‑1 (MQ=255)                 gATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCt                                                    <  1:3864450/71‑1 (MQ=255)                  aTAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTg                                                   <  1:1719899/71‑1 (MQ=255)                   tAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCaaat                                                                       >  1:2838685/1‑49 (MQ=255)                   tAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAg                                                                       >  1:4526804/1‑50 (MQ=255)                     gTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGAc                                                 <  1:1467292/70‑1 (MQ=255)                       cACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCg                                               >  1:2046154/1‑70 (MQ=255)                        aCCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCttcttc                                                               <  1:1122837/53‑1 (MQ=255)                        aCCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCttcttc                                                               <  1:5069461/53‑1 (MQ=255)                        aCCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGt                                              <  1:3421504/70‑1 (MQ=255)                        aCCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGt                                              <  1:3366307/70‑1 (MQ=255)                        aCCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGt                                              <  1:1618773/70‑1 (MQ=255)                         ccGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCt                                                    <  1:3044744/63‑1 (MQ=255)                          cGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAg                                           <  1:4823073/71‑1 (MQ=255)                              cAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCAAAAGCct                                                                    >  1:2867693/1‑42 (MQ=255)                               aGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCt                                                    <  1:2618066/57‑1 (MQ=255)                               aGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCt                                                    <  1:3956211/57‑1 (MQ=255)                               aGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCt                                                    <  1:3711902/57‑1 (MQ=255)                               aGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCa                                      <  1:2801014/71‑1 (MQ=255)                               aGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCa                                      <  1:4870706/71‑1 (MQ=255)                                 ttGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCttcttc                                                               <  1:4288261/44‑1 (MQ=255)                                 ttGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGa                                                  >  1:4164029/1‑57 (MQ=255)                                  tGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGAc                                                 >  1:4454723/1‑57 (MQ=255)                                   gTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAAc                                   >  1:397338/1‑70 (MQ=255)                                    tACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACg                                  <  1:2163741/70‑1 (MQ=255)                                    tACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACg                                  <  1:3624120/70‑1 (MQ=255)                                    tACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACg                                  <  1:4633088/70‑1 (MQ=255)                                       cAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCcttct                                                                 <  1:2946985/36‑1 (MQ=255)                                        agtagtTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAg                                           >  1:164612/1‑57 (MQ=255)                                        agtagtTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAg                                           >  1:2267097/1‑57 (MQ=255)                                        agtagtTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAg                                           >  1:2585507/1‑57 (MQ=255)                                           agtTCCAGCCCAGCGCGAAGCCAAGGCCTTCTTCAACATAGTTCTGa                                                  >  1:4738282/1‑47 (MQ=255)                                           agtTCCAGCCCAGCGCGAAGCCAAGGCCTTCTTCAACATAGTTCTGa                                                  >  1:2646733/1‑47 (MQ=255)                                           agtTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGa                           >  1:2407876/1‑70 (MQ=255)                                            gtTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGaaa                         >  1:5086816/1‑71 (MQ=255)                                            gtTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGaaa                         >  1:2546401/1‑71 (MQ=255)                                             ttCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAAc                        <  1:123148/71‑1 (MQ=255)                                               ccagcccagcGCGAAGCCAAGGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCg                      >  1:1039129/1‑71 (MQ=255)                                               ccagcccagcGCGAAGCCAAAGCCTTCTTCAACATAGtt                                                      >  1:357588/1‑39 (MQ=255)                                               ccagcccagcGCGAAGCCAAAGCCTTCTTCAACATAGtt                                                      >  1:2073609/1‑39 (MQ=255)                                               ccagcccagcGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCg                      >  1:4988401/1‑71 (MQ=255)                                               ccagcccagcGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCg                      >  1:1570992/1‑71 (MQ=255)                                               ccagcccagcGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCg                      >  1:2938551/1‑71 (MQ=255)                                                  gcccagcGCGAAGCCAAGGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCa                   >  1:3252839/1‑71 (MQ=255)                                                   cccagcGCGAAGCCAAGGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCa                   >  1:895987/1‑70 (MQ=255)                                                   cccagcGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAAc                        <  1:1726682/65‑1 (MQ=255)                                                   cccagcGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAAc                        <  1:2927326/65‑1 (MQ=255)                                                    cccgcgcGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCata                 <  1:1947510/68‑1 (MQ=255)                                                    ccagcGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCata                 <  1:153485/71‑1 (MQ=255)                                                    ccagcGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCata                 <  1:1924784/71‑1 (MQ=255)                                                    ccagcGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCat                  <  1:603229/70‑1 (MQ=255)                                                    ccagcGCGAACCCAAGGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCata                 <  1:2605795/71‑1 (MQ=255)                                                     cagcGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCatat                <  1:1075872/71‑1 (MQ=255)                                                     cagcGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCatat                <  1:706924/71‑1 (MQ=255)                                                     cagcGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCatat                <  1:614357/71‑1 (MQ=255)                                                     cagcGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCata                 >  1:3161216/1‑70 (MQ=255)                                                     cagcGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCata                 >  1:3372106/1‑70 (MQ=255)                                                      agcGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTgg                                        <  1:2247685/46‑1 (MQ=255)                                                       gcgcGAAGCCAAGGCCTTCTTCAACATAGTTCTGAcc                                                <  1:4178158/37‑1 (MQ=255)                                                       gcgcGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGa                                 >  1:2861194/1‑52 (MQ=255)                                                       gcgcGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGa                                 >  1:2636062/1‑52 (MQ=255)                                                          cgAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCa           <  1:1497677/71‑1 (MQ=255)                                                           gAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAg          <  1:2271711/71‑1 (MQ=255)                                                               ccAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCatat                <  1:3925105/61‑1 (MQ=255)                                                                 aaaGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCatatat              >  1:4826876/1‑61 (MQ=255)                                                                 aaaGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGt     >  1:2503106/1‑70 (MQ=255)                                                                   aGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCa  <  1:2348047/71‑1 (MQ=255)                                                                   aGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCa  <  1:2386429/71‑1 (MQ=255)                                                                   aGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCa  <  1:478798/71‑1 (MQ=255)                                                                   |                                                                       CAACCAGGTCAACGGCGATAGTCACCGCCCAGTTGTACCAGTAGTTCCAGCCCAGCGCGAAGCCAAAGCCTTCTTCAACATAGTTCTGACCGTAAGTGGCAAACGAACCGGAAACCGGCATATATGCAGCCAGTTCA  >  W3110S.gb/2251521‑2251657

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