Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	158,626	G→A	S167S (AGC→AGT)	pcnB ←	poly(A) polymerase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	158,626	0	G	A	100.0%	187.4 / NA	64	S167S (AGC→AGT)	pcnB	poly(A) polymerase I
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (43/21);  total (43/21)

CGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCA  >  W3110S.gb/158559‑158693                                                                    |                                                                    cGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGACTg                                                                  <  1:1476053/71‑1 (MQ=255)   tCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGACTGtt                                                                <  1:125007/71‑1 (MQ=255)     cTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGACTGTTg                                                               >  1:3209431/1‑70 (MQ=255)     cTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGACTGTTGa                                                              >  1:3224719/1‑71 (MQ=255)     cTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGACTGTTGa                                                              <  1:4238805/71‑1 (MQ=255)     cTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGACTGTTGa                                                              >  1:4394336/1‑71 (MQ=255)       tCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGACTGTTGa                                                              >  1:2211640/1‑69 (MQ=255)        cAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGACTGTTGAt                                                             >  1:155036/1‑69 (MQ=255)        cAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGACTGTTGATAg                                                           >  1:3855554/1‑71 (MQ=255)             ccTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGACTGTTGATAGTGaaa                                                      >  1:1587694/1‑71 (MQ=255)             ccTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGACTGTTGATAGTGaaa                                                      >  1:3991735/1‑71 (MQ=255)             ccTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGACTGTTGATAGTGaaa                                                      >  1:1308828/1‑71 (MQ=255)             ccTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGACTGTTGATAGTGaaa                                                      >  1:2775821/1‑71 (MQ=255)             ccTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGACTGTTGATAGTGaa                                                       >  1:658476/1‑70 (MQ=255)             ccTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGACTGTTGATAGTGaa                                                       >  1:1428728/1‑70 (MQ=255)             ccTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGACTGTTGATAGTGaa                                                       >  1:2401992/1‑70 (MQ=255)             ccTTCAAGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGACTGTTGATAGTGaaa                                                      >  1:4971/1‑71 (MQ=255)                    gccgccAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGACTGTTGATAGTGAAATCgcgg                                                >  1:3759949/1‑70 (MQ=255)                     ccgccAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGACTGTTGATAGTGAAATCgcggc                                               <  1:3596569/70‑1 (MQ=255)                     ccgccAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGACTGTTGATAGTGAAATCgcggc                                               <  1:4060303/70‑1 (MQ=255)                     ccgccAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGACTGTTGATAGTGAAATCgcggc                                               <  1:27061/70‑1 (MQ=255)                             gTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGACTGTTGATAGTGAAATcgcg                                                 >  1:2398786/1‑60 (MQ=255)                             gTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGACTGTTGATAGTGAAATCgcgg                                                >  1:3752797/1‑61 (MQ=255)                             gTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGACTGTTGATAGTGAAATCgcgg                                                >  1:3866484/1‑61 (MQ=255)                             gTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGACTGTTGATAGTGAAATCGCGGCGCTGGGCg                                       >  1:3841260/1‑70 (MQ=255)                             gTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGt                                      >  1:1643475/1‑71 (MQ=255)                             gTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGACTGGTGATAGTGAAATCGCGGCGCTGGGCg                                       >  1:4297883/1‑70 (MQ=255)                              tAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGtc                                     <  1:2226207/71‑1 (MQ=255)                                aTCACGGACGGTAAAATCCGCTACGCTGTAATACAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGtctt                                   >  1:810486/1‑71 (MQ=255)                                aTCACGGACGGTAAAATCCGCTACGCTGTAATACAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGtct                                    <  1:4130624/70‑1 (MQ=255)                                aTCACGGACGGTAAAATCCGCTACGCTGTAATACAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGtct                                    <  1:3029812/70‑1 (MQ=255)                                   acggacggTAAAATCCGCTACGCTGTAATACAGACTGTTGATAGTGaaa                                                      >  1:698103/1‑49 (MQ=255)                                   acggacggTAAAATCCGCTACGCTGTAATACAGACTGTTGATAGTGAAATcgc                                                  >  1:2481381/1‑53 (MQ=255)                                   acggacggTAAAATCCGCTACGCTGTAATACAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGtcttct                                 >  1:3638192/1‑70 (MQ=255)                                   acggacggTAAAATCCGCTACGCTGTAATACAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTcttctt                                >  1:3290511/1‑71 (MQ=255)                                    cggacggTAAAATCCGCTACGCTGTAATACAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCttcttc                               <  1:396479/71‑1 (MQ=255)                                    cggacggTAAAATCCGCTACGCTGTAATACAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCttcttc                               <  1:3446804/71‑1 (MQ=255)                                     ggacggTAAAATCCGCTACGCTGTAATACAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCttcttc                               <  1:4522898/70‑1 (MQ=255)                                      gacggTAAAATCCGCTACGCTGTAATACAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGa                             <  1:2132598/71‑1 (MQ=255)                                            aaaaTCCGCTACGCTGTAATACAGACTGTTTATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGc                       >  1:11996/1‑71 (MQ=255)                                            aaaaTCCGCTACGCTGTAATACAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGc                       >  1:4310908/1‑71 (MQ=255)                                            aaaaTCCGCTACGCTGTAATACAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGc                       >  1:3168991/1‑71 (MQ=255)                                                tCCGCTACGCTGTAATACAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGa                    >  1:3221867/1‑70 (MQ=255)                                                   gCTACGCTGTAATACAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAAt                <  1:4591294/71‑1 (MQ=255)                                                       cGCTGTAATACAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAAtgtt             <  1:2889780/70‑1 (MQ=255)                                                       cGCTGTAATACAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAAtgtt             <  1:3194905/70‑1 (MQ=255)                                                             aaTACAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTcgcg       >  1:3434800/1‑70 (MQ=255)                                                              aTACAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTcgcg       >  1:4121858/1‑69 (MQ=255)                                                              aTACAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTcgcg       >  1:3556533/1‑69 (MQ=255)                                                              aTACAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTcgcg       >  1:1408918/1‑69 (MQ=255)                                                              aTACAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTcgcg       >  1:4276360/1‑69 (MQ=255)                                                              aTACAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGACGAAAATGTTGTcgcg       >  1:2635757/1‑69 (MQ=255)                                                              aTACAGACTGTTGATAGTCAAATCGCGGCGCTGGGCGTCTTCTTAGATGGAGCCGAAAATGTTGTcgcg       >  1:4126867/1‑69 (MQ=255)                                                                 cAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGa                             >  1:3089584/1‑44 (MQ=255)                                                                 cAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCgcagca  <  1:3466822/71‑1 (MQ=255)                                                                 cAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCgcagc   <  1:3478032/70‑1 (MQ=255)                                                                 cAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCgcagc   >  1:2040927/1‑70 (MQ=255)                                                                 cAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCgcagc   <  1:465989/70‑1 (MQ=255)                                                                 cAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCgcagc   >  1:3760937/1‑70 (MQ=255)                                                                 cAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCgcagc   >  1:1048495/1‑70 (MQ=255)                                                                 cAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCgcagc   >  1:4257376/1‑70 (MQ=255)                                                                 cAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCgcagc   <  1:340226/70‑1 (MQ=255)                                                                 cAGACTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTAGCgcagc   >  1:4293311/1‑70 (MQ=255)                                                                 cAGACTGTTGATAGTGAAATCGCGGCGCTGGCCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCgcagc   <  1:1430345/70‑1 (MQ=255)                                                                    |                                                                    CGTCCTTCAGATCCTTCATGCCGCCAACGTAATCACGGACGGTAAAATCCGCTACGCTGTAATACAGGCTGTTGATAGTGAAATCGCGGCGCTGGGCGTCTTCTTCGATGGAGCCGAAAATGTTGTCGCGCAGCA  >  W3110S.gb/158559‑158693

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