Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	3,059,788	C→T	36.6%	L95L (CTG→TTG)	yliK →	methylmalonyl‑CoA mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,059,788	0	C	T	36.6%	39.8 / 39.0	41	L95L (CTG→TTG)	yliK	methylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand):  ref base C (19/7);  new base T (5/10);  total (24/17)
Fisher's exact test for biased strand distribution p-value = 2.12e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAA  >  W3110S.gb/3059718‑3059858                                                                       |                                                                       aCCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAAcc                                                                        <  1:2728618/71‑1 (MQ=255) aCCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAAcc                                                                        <  1:4544559/71‑1 (MQ=255)   cGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACTTg                                                                      <  1:2643467/71‑1 (MQ=255)     tGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACTTGgc                                                                    <  1:4555130/71‑1 (MQ=255)      ggACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGgcc                                                                   >  1:666994/1‑71 (MQ=255)          cATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGccgccg                                                               >  1:2517842/1‑71 (MQ=255)              cGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGc                                                            >  1:1144583/1‑70 (MQ=255)                tCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACTTGgccgc                                                                 >  1:3184471/1‑63 (MQ=255)                tCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACTTGgccgc                                                                 >  1:3107066/1‑63 (MQ=255)                  aGTATGCTGGTTTTTCAACAGCAAAAGAGTCCATCGCTTTTTATCGCCGTAACTTGGCCGCCCGGCaaaa                                                        <  1:1316133/70‑1 (MQ=255)                  aGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACTTGGCCGCCGGGCaaaa                                                        <  1:4054377/70‑1 (MQ=255)                   gTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACTTGGCCGCCGGGCAAAAAg                                                      <  1:309663/71‑1 (MQ=255)                    tATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACTTGGCCGCCGGGCAAAAAgg                                                     <  1:4513312/71‑1 (MQ=255)                      tGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGTCGGGCAAAAAGGTc                                                   >  1:534411/1‑71 (MQ=255)                      tGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTc                                                   >  1:2389174/1‑71 (MQ=255)                           gTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACTTGGCCGCCGGGCAAAAAgg                                                     >  1:3692119/1‑64 (MQ=255)                           gTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACTTGGCCGCCGGGCAAAAAgg                                                     >  1:2991765/1‑64 (MQ=255)                             ttttCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGt                                            >  1:1009235/1‑71 (MQ=255)                             ttttCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGt                                            >  1:1546409/1‑71 (MQ=255)                             ttttCAACAGCAAAAGAGTCAAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAg                                                      >  1:2357981/1‑61 (MQ=255)                              tttCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAg                                                      >  1:3065728/1‑60 (MQ=255)                              tttCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGtt                                           >  1:1997515/1‑71 (MQ=255)                                tCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGc                                         >  1:1107317/1‑71 (MQ=255)                                     aGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACTTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTg                                    <  1:703037/71‑1 (MQ=255)                                        aaaaGAGTCCAACGCTTTTTATCGCCGTAACTTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTg                                    <  1:3883154/68‑1 (MQ=255)                                         aaaGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCt                                <  1:4327075/71‑1 (MQ=255)                                          aaGAGTCCAACGCTTTTTATCGCCGTAACTTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCtt                               <  1:368966/71‑1 (MQ=255)                                              gTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGcca                           >  1:4630743/1‑71 (MQ=255)                                                 cAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGccacc                         >  1:3001662/1‑70 (MQ=255)                                                   aCGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGccacccac                      >  1:1464316/1‑71 (MQ=255)                                                   aCGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGccaccca                       <  1:3712085/70‑1 (MQ=255)                                                     gCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCt                                                  >  1:2859570/1‑41 (MQ=255)                                                     gCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCt                                                  >  1:273329/1‑41 (MQ=255)                                                      cTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGt                   >  1:3127600/1‑71 (MQ=255)                                                       tttttATCGCCGTAACTTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTg                  >  1:1791216/1‑71 (MQ=255)                                                        ttttATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCg                                             <  1:814721/43‑1 (MQ=255)                                                        ttttATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTg                  >  1:4431494/1‑70 (MQ=255)                                                             tCGCCGTAACTTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTAc             <  1:5028458/70‑1 (MQ=255)                                                                     aCCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGa     >  1:3961696/1‑70 (MQ=255)                                                                     aCCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGAt    <  1:3973753/71‑1 (MQ=255)                                                                     aCCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGAt    <  1:3795566/71‑1 (MQ=255)                                                                     aCCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGAt    <  1:3033356/71‑1 (MQ=255)                                                                     aCCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGAt    <  1:2355307/71‑1 (MQ=255)                                                                       cTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATaa  <  1:3566280/71‑1 (MQ=255)                                                                       |                                                                       ACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCAAAAGAGTCCAACGCTTTTTATCGCCGTAACCTGGCCGCCGGGCAAAAAGGTCTTTCCGTTGCGTTTGACCTTGCCACCCACCGTGGCTACGACTCCGATAA  >  W3110S.gb/3059718‑3059858

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