Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	1,477,935	G→A	32.4%	E360K (GAA→AAA)	ydbD →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,477,935	0	G	A	32.4%	94.5 / 61.6	74	E360K (GAA→AAA)	ydbD	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (26/24);  new base A (16/8);  total (42/32)
Fisher's exact test for biased strand distribution p-value = 3.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCA  >  W3110S.gb/1477868‑1477992                                                                    |                                                          tGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTaaa                                                         >  1:1519150/1‑70 (MQ=255) tGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTaaa                                                         >  1:1773199/1‑70 (MQ=255) tGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTaaa                                                         >  1:3081563/1‑70 (MQ=255) tGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTaaa                                                         >  1:2127140/1‑70 (MQ=255) tGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTaaa                                                         >  1:2083927/1‑70 (MQ=255) tGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCATGCGCTACGTaaa                                                         >  1:2475876/1‑70 (MQ=255)   tCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTaaa                                                         <  1:37194/68‑1 (MQ=255)   tCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTaaa                                                         <  1:1600543/68‑1 (MQ=255)   tCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTaaa                                                         <  1:289083/68‑1 (MQ=255)   tCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTaaa                                                         <  1:3315967/68‑1 (MQ=255)   tCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTaaa                                                         <  1:1287460/68‑1 (MQ=255)   tCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTaaa                                                         <  1:1272643/68‑1 (MQ=255)   tCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTaaa                                                         <  1:1153344/68‑1 (MQ=255)   tCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTaaa                                                         <  1:1012131/68‑1 (MQ=255)        ataatGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTAAAAAtct                                                    >  1:396086/1‑68 (MQ=255)        ataatGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTAAAAATCTcc                                                  >  1:3314647/1‑70 (MQ=255)                    aTGAGTCGTCTAGCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGt                                     <  1:293446/71‑1 (MQ=255)                    aTGAGTCGTCTACCTGATGGTTTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGt                                     <  1:1600350/71‑1 (MQ=255)                    aTGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGt                                     <  1:2562582/71‑1 (MQ=255)                    aTGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGt                                     <  1:972110/71‑1 (MQ=255)                    aTGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGt                                     <  1:2151782/71‑1 (MQ=255)                    aTGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGt                                     <  1:1676182/71‑1 (MQ=255)                    aTGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGt                                     <  1:2220941/71‑1 (MQ=255)                           gtcTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGt                                     <  1:3047353/64‑1 (MQ=255)                                ccTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTgg                          >  1:1849673/1‑70 (MQ=255)                                ccTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTgg                          >  1:1292170/1‑70 (MQ=255)                                ccTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGc                         >  1:3035067/1‑71 (MQ=255)                                ccTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGc                         >  1:2002531/1‑71 (MQ=255)                                ccTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGc                         >  1:2996285/1‑71 (MQ=255)                                   gATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCg                        >  1:172145/1‑69 (MQ=255)                                   gATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGa                       >  1:1467982/1‑70 (MQ=255)                                   gATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTAATGACTCTGGCGa                       >  1:2482531/1‑70 (MQ=255)                                   gATGGTGTCGGTCCGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGa                       >  1:1192359/1‑70 (MQ=255)                                      ggTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGaa                      >  1:3261315/1‑68 (MQ=255)                                      ggTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGaa                      >  1:593934/1‑68 (MQ=255)                                      ggTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGaa                      >  1:1071790/1‑68 (MQ=255)                                         gtCGGTTCGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGaaaa                    >  1:2013824/1‑67 (MQ=255)                                         gtCGGTCCGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGaaaa                    >  1:2470785/1‑67 (MQ=255)                                         gtCGGTCCGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGaaaa                    >  1:917002/1‑67 (MQ=255)                                         gtCGGTCCGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGaaaa                    >  1:2840896/1‑67 (MQ=255)                                               ccGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGTCGAAAAGAGATTGGGa          >  1:2945015/1‑71 (MQ=255)                                               ccGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGa          >  1:3329602/1‑71 (MQ=255)                                               ccGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGa          >  1:1116603/1‑71 (MQ=255)                                                cGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGAc                               >  1:418455/1‑49 (MQ=255)                                                cGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGa          >  1:41139/1‑70 (MQ=255)                                                cGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGa          >  1:549366/1‑70 (MQ=255)                                                cGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGa          >  1:541071/1‑70 (MQ=255)                                                cGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGa          >  1:1567759/1‑70 (MQ=255)                                                cGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGa          >  1:441288/1‑70 (MQ=255)                                                    gAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATg      <  1:1548236/70‑1 (MQ=255)                                                    gAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATg      <  1:36094/70‑1 (MQ=255)                                                    gAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATg      <  1:2870550/70‑1 (MQ=255)                                                    gAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATg      <  1:2120793/70‑1 (MQ=255)                                                    gAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATg      <  1:669923/70‑1 (MQ=255)                                                    gAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATg      <  1:1380759/70‑1 (MQ=255)                                                    gAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATg      <  1:2305584/70‑1 (MQ=255)                                                    gAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATg      <  1:2164672/70‑1 (MQ=255)                                                       tCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGAt    >  1:3315688/1‑69 (MQ=255)                                                       tCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCa  >  1:3305922/1‑71 (MQ=255)                                                       tCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCa  >  1:3302003/1‑71 (MQ=255)                                                       tCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCa  >  1:439564/1‑71 (MQ=255)                                                       tCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCa  >  1:448029/1‑71 (MQ=255)                                                       tCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCa  >  1:1094663/1‑71 (MQ=255)                                                       tCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCa  >  1:738921/1‑71 (MQ=255)                                                       tCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCa  >  1:76340/1‑71 (MQ=255)                                                       tCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCa  >  1:1448435/1‑71 (MQ=255)                                                       tCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCa  >  1:2942720/1‑71 (MQ=255)                                                       tCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCa  >  1:2253502/1‑71 (MQ=255)                                                       tCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCa  >  1:2461040/1‑71 (MQ=255)                                                       tCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCa  >  1:284435/1‑71 (MQ=255)                                                       tCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCa  >  1:2629866/1‑71 (MQ=255)                                                       tCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTAGGATATGATCa  >  1:1481437/1‑71 (MQ=255)                                                                   tGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCa  <  1:2782701/59‑1 (MQ=255)                                                                   tGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCa  <  1:1035393/59‑1 (MQ=255)                                                                   tGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCa  <  1:1065295/59‑1 (MQ=255)                                                                   tGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCa  <  1:605271/59‑1 (MQ=255)                                                                   tGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCa  <  1:2136728/59‑1 (MQ=255)                                                                   tGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCa  <  1:542708/59‑1 (MQ=255)                                                                   tGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCa  <  1:529453/59‑1 (MQ=255)                                                                   tGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCa  <  1:2757826/59‑1 (MQ=255)                                                                   tGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCa  <  1:2172162/59‑1 (MQ=255)                                                                   tGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCa  <  1:2017243/59‑1 (MQ=255)                                                                   tGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCa  <  1:1860389/59‑1 (MQ=255)                                                                   tGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCa  <  1:1360195/59‑1 (MQ=255)                                                                   tGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCa  <  1:1695339/59‑1 (MQ=255)                                                                   tGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCa  <  1:1377535/59‑1 (MQ=255)                                                                   tGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCa  <  1:1385629/59‑1 (MQ=255)                                                                   tGAAAATCTCCCTCTGGATATCGTAATGACTCTGGCGAAAAGAGATTGGGATATGATCa  <  1:20931/59‑1 (MQ=255)                                                                    gAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCa  <  1:2797233/58‑1 (MQ=255)                                                                    gAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCa  <  1:2155009/58‑1 (MQ=255)                                                                    |                                                          TGTCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCA  >  W3110S.gb/1477868‑1477992

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