Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	810,589	T→G	M275R (ATG→AGG)	bioB →	biotin synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	810,589	0	T	G	79.1%	120.8 / 13.0	43	M275R (ATG→AGG)	bioB	biotin synthase
Reads supporting (aligned to +/- strand):  ref base T (0/9);  new base G (26/8);  total (26/17)
Fisher's exact test for biased strand distribution p-value = 4.31e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.35e-01

ATGCCAACCTCTTACGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGG  >  W3110S.gb/810525‑810654                                                                 |                                                                  atgCCAACCTCTTACGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCTATGTGCtt                                                             <  1:1704722/71‑1 (MQ=255) atgCCAACCTCTTACGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCtt                                                             <  1:4017700/71‑1 (MQ=255) atgCCAACCTCTTACGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCtt                                                             <  1:384905/71‑1 (MQ=255) atgCCAACCTCTTACGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGAGGTGCtt                                                             <  1:2243889/71‑1 (MQ=255)                               ggACGCGAGCAGATGAACGAACAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGATTTTc                                    >  1:1431911/1‑66 (MQ=255)                                gACGCGAGCAGATGAACGAACAGACTCAGGCGAggtg                                                                >  1:2039426/1‑37 (MQ=255)                                gACGCGAGCAGATGAACGAACAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGAtt                                       >  1:293030/1‑62 (MQ=255)                                gACGCGAGCAGATGAACGAACAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGAtt                                       >  1:648828/1‑62 (MQ=255)                                gACGCGAGCAGATGAACGAACAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGAtt                                       >  1:563519/1‑62 (MQ=255)                                gACGCGAGCAGATGAACGAACAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGAtt                                       >  1:3366004/1‑62 (MQ=255)                                gACGCGAGCAGATGAACGAACAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGAtt                                       >  1:1077317/1‑62 (MQ=255)                                gACGCGAGCAGATGAACGAACAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGAtt                                       >  1:230322/1‑62 (MQ=255)                                gACGCGAGCAGATGAACGAACAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGAtt                                       >  1:2093817/1‑62 (MQ=255)                                gACGCGAGCAGATGAACGAACAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGAtt                                       >  1:204096/1‑62 (MQ=255)                                    cgAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCa                          <  1:3512744/71‑1 (MQ=255)                                    cgAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCa                          <  1:3913483/71‑1 (MQ=255)                                    cgAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCa                          <  1:3209491/71‑1 (MQ=255)                                    cgAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCa                          <  1:454175/71‑1 (MQ=255)                                    cgAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCa                          <  1:2678802/71‑1 (MQ=255)                                    cgAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCa                          <  1:281021/71‑1 (MQ=255)                                       gCAGATGAACGAACAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCaaa                        >  1:2607947/1‑70 (MQ=255)                                       gCAGATGAACGAACAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCaa                         >  1:3147513/1‑69 (MQ=255)                                       gCAGATGAACGAACAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGATTTTATACGGTTGCaaa                        >  1:3555484/1‑70 (MQ=255)                                             gaacgaacAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGa                 >  1:3202159/1‑71 (MQ=255)                                             gaacgaacAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGa                 >  1:3946608/1‑71 (MQ=255)                                             gaaagaacAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAActg                     >  1:2387007/1‑67 (MQ=255)                                                    cAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGcc          >  1:1532011/1‑71 (MQ=255)                                                    cAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGcc          >  1:1993294/1‑71 (MQ=255)                                                    cAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGcc          >  1:181437/1‑71 (MQ=255)                                                    cAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGc           >  1:397705/1‑70 (MQ=255)                                                    cAGACTCAGGCGAGGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGc           >  1:1121430/1‑70 (MQ=255)                                                            ggCGAGGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCgg  <  1:604154/71‑1 (MQ=255)                                                            ggCGAGGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCgg  <  1:2260896/71‑1 (MQ=255)                                                            ggCGAGGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCgg  <  1:3384466/71‑1 (MQ=255)                                                            ggCGAGGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCgg  <  1:3284101/71‑1 (MQ=255)                                                            ggCGAGGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCgg  <  1:2673578/71‑1 (MQ=255)                                                            ggCGAGGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCgg  <  1:2615583/71‑1 (MQ=255)                                                            ggCGAGGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCgg  <  1:2525426/71‑1 (MQ=255)                                                             gCGAGGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCgg  >  1:3554053/1‑70 (MQ=255)                                                             gCGAGGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCgg  >  1:2160021/1‑70 (MQ=255)                                                             gCGAGGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAAACCgg  >  1:4030288/1‑70 (MQ=255)                                                              cGAGGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCgg  >  1:3447105/1‑69 (MQ=255)                                                              cGAGGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCgg  >  1:308406/1‑69 (MQ=255)                                                              cGAGGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCgg  >  1:109200/1‑69 (MQ=255)                                                                 |                                                                  ATGCCAACCTCTTACGTGCGCCTTTCTGCCGGACGCGAGCAGATGAACGAACAGACTCAGGCGATGTGCTTTATGGCAGGCGCAAACTCGATTTTCTACGGTTGCAAACTGCTGACCACGCCGAATCCGG  >  W3110S.gb/810525‑810654

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