Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	1,477,935	G→A	73.6%	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	73.6%	82.9 / 36.5	53	E360K (GAA→AAA)	ydbD	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (4/10);  new base A (17/22);  total (21/32)
Fisher's exact test for biased strand distribution p-value = 3.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.93e-01

TCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGAGC  >  W3110S.gb/1477870‑1478001                                                                  |                                                                   tCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGaa                                                                  <  1:2559304/68‑1 (MQ=255) tCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGaa                                                                  <  1:377031/68‑1 (MQ=255) tCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGaa                                                                  <  1:2010445/68‑1 (MQ=255)    aaaTAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGaa                                                                  <  1:3203301/65‑1 (MQ=255)    aaaTAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAAtctc                                                            <  1:2389533/71‑1 (MQ=255)                gAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATc                                                <  1:1081998/71‑1 (MQ=255)                   tGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGt                                              >  1:3136707/1‑70 (MQ=255)                   tGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGt                                              >  1:176707/1‑70 (MQ=255)                        cgtcTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGAc                                        >  1:1849501/1‑71 (MQ=255)                         gtcTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGAc                                        <  1:2890289/70‑1 (MQ=255)                         gtcTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGAc                                        <  1:3227500/70‑1 (MQ=255)                         gtcTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGAc                                        <  1:2617896/70‑1 (MQ=255)                         gtcTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGAc                                        <  1:3310133/70‑1 (MQ=255)                         gtcTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGAc                                        <  1:3616862/70‑1 (MQ=255)                         gtcTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGAc                                        <  1:1118969/70‑1 (MQ=255)                              ccTGATTGTGTCGGTCCGTTGAATCTTGCGCTACGTGaaa                                                                 <  1:2089711/40‑1 (MQ=255)                              ccTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGaaa                                                                 <  1:2003739/40‑1 (MQ=255)                              ccTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGaaa                                                                 <  1:2991794/40‑1 (MQ=255)                              ccTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGaaa                                                                 <  1:946845/40‑1 (MQ=255)                               cTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGaaa                                                                 <  1:685315/39‑1 (MQ=255)                                tGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAAt                                                               >  1:730551/1‑40 (MQ=255)                                tGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAAt                                                               >  1:1345392/1‑40 (MQ=255)                                tGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAAt                                                               >  1:963070/1‑40 (MQ=255)                                tGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAAt                                                               >  1:780211/1‑40 (MQ=255)                                      tgtCGGTCCGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAagag                          <  1:2877533/71‑1 (MQ=255)                                      tgtCGGTCCGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAagag                          <  1:481649/71‑1 (MQ=255)                                       gtCGGTCCGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGaaaa                             >  1:1816121/1‑67 (MQ=255)                                       gtCGGTCCGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGaaaa                             >  1:3285521/1‑67 (MQ=255)                                       gtCGGTCCGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGaaaa                             >  1:2079174/1‑67 (MQ=255)                                              cGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAagag                          <  1:1338358/63‑1 (MQ=255)                                              cGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGa                   >  1:1032315/1‑70 (MQ=255)                                              cGTTGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGa                   >  1:2863455/1‑70 (MQ=255)                                                ttGAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGatat                <  1:738104/71‑1 (MQ=255)                                                  gAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGa              >  1:3023237/1‑71 (MQ=255)                                                  gAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGa              >  1:1060395/1‑71 (MQ=255)                                                  gAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGa              >  1:1078143/1‑71 (MQ=255)                                                  gAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGa              >  1:2789844/1‑71 (MQ=255)                                                  gAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGa              >  1:262157/1‑71 (MQ=255)                                                  gAATCTTGCGCTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGa              >  1:1752089/1‑71 (MQ=255)                                                         gcgcTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAgaga                         >  1:670713/1‑53 (MQ=255)                                                          cgcTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAg      <  1:2171067/71‑1 (MQ=255)                                                          cgcTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAg      <  1:1424725/71‑1 (MQ=255)                                                          cgcTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAg      <  1:1428829/71‑1 (MQ=255)                                                          cgcTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAg      <  1:2336031/71‑1 (MQ=255)                                                          cgcTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAg      <  1:2376153/71‑1 (MQ=255)                                                          cgcTACGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAg      <  1:2796107/71‑1 (MQ=255)                                                            cTACGTAAAAATCTCCATCTGGATATCGTCATGACTCTGGCGaa                               <  1:1025534/44‑1 (MQ=38)                                                              aCGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGAGc  >  1:649120/1‑71 (MQ=255)                                                              aCGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGAGc  >  1:3444942/1‑71 (MQ=255)                                                               cGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGAt             <  1:4098/59‑1 (MQ=255)                                                               cGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGAt             <  1:2475149/59‑1 (MQ=255)                                                               cGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGAt             <  1:3238668/59‑1 (MQ=255)                                                               cGTAAAAATCTCCCTCTGGATATCGTCATGACTCTGGCG‑AAATAGATTGGGATATGAt             <  1:3078555/58‑1 (MQ=38)                                                                  |                                                                   TCGAAATAATGACGCGAATGAGTCGTCTACCTGATGGTGTCGGTCCGTTGAATCTTGCGCTACGTGAAAATCTCCCTCTGGATATCGTCATGACTCTGGCGAAAAGAGATTGGGATATGATCATTCAGGAGC  >  W3110S.gb/1477870‑1478001

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