Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	1,006,848	A→T	24.3%	K159* (AAA→TAA)	ycbW →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,006,848	0	A	T	24.3%	67.2 / 21.1	37	K159* (AAA→TAA)	ycbW	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (9/19);  new base T (3/6);  total (12/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

CACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTCTAA  >  W3110S.gb/1006780‑1006916                                                                     |                                                                     cACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTa                                                                      <  1:947426/69‑1 (MQ=255)   cTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTaaa                                                                    <  1:1815250/69‑1 (MQ=255)   cTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTaaa                                                                    <  1:2122406/69‑1 (MQ=255)     aTGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTa                                                                      >  1:1772468/1‑65 (MQ=255)     aTGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTa                                                                      >  1:200891/1‑65 (MQ=255)     aTGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAAt                                                                  >  1:3656518/1‑69 (MQ=255)     aTGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAAt                                                                  >  1:1847286/1‑69 (MQ=255)     aTGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAAt                                                                  >  1:2875332/1‑69 (MQ=255)     aTGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAAt                                                                  >  1:2659134/1‑69 (MQ=255)      tGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTtaa                                                                    >  1:597668/1‑66 (MQ=255)      tGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTtaa                                                                    >  1:1789667/1‑66 (MQ=255)      tGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTtaa                                                                    >  1:3377926/1‑66 (MQ=255)      tGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAATc                                                                 >  1:444616/1‑69 (MQ=255)      tGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAATc                                                                 >  1:2284530/1‑69 (MQ=255)              ggCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCTACGCCATTAAAATCATGAACGa                                                         <  1:1907115/69‑1 (MQ=255)                   aaCCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCATGAACGa                                                         <  1:1034446/64‑1 (MQ=255)                    aCCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCATGAACGACAGGc                                                    >  1:1939754/1‑68 (MQ=255)                    aCCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCATGAACGACAGGc                                                    >  1:921657/1‑68 (MQ=255)                    aCCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCATGAACGACAGGc                                                    >  1:2744387/1‑68 (MQ=255)                     cacTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCATGAACGa                                                         <  1:3122037/60‑1 (MQ=255)                      ccTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCATGAACGACAGGCTGa                                                 <  1:34827/69‑1 (MQ=255)                      ccTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCATGAACGACAGGCTGa                                                 <  1:3081334/69‑1 (MQ=255)                      ccTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCATGAACGACAGGATGa                                                 <  1:2795360/69‑1 (MQ=255)                      ccTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCATGAACGACAGGATGa                                                 <  1:2799621/69‑1 (MQ=255)                                 aTTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGtt                                      <  1:3284205/69‑1 (MQ=255)                                 aTTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGtt                                      <  1:2371536/69‑1 (MQ=255)                                 aTTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGtt                                      <  1:3057344/69‑1 (MQ=255)                                   tGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTTAAATCATGAACGACAGGCTGAAACCGCAGGTTa                                     <  1:3389814/68‑1 (MQ=255)                                    gCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTaa                                    <  1:2923304/68‑1 (MQ=255)                                       ggTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGtt                                <  1:102599/69‑1 (MQ=255)                                       ggTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGtt                                <  1:743031/69‑1 (MQ=255)                                       ggTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGtt                                <  1:841074/69‑1 (MQ=255)                                       ggTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGtt                                <  1:2035188/69‑1 (MQ=255)                                               gATGCAGTTGGGCGACGCCATTTAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCtt                        <  1:2841369/69‑1 (MQ=255)                                                        gggCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCg                >  1:420076/1‑68 (MQ=255)                                                             aCGCCATTTAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAgg          <  1:2357176/69‑1 (MQ=255)                                                             aCGCCATTTAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAgg          <  1:1352230/69‑1 (MQ=255)                                                             aCGCCATTTAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAgg          <  1:562477/69‑1 (MQ=255)                                                             aCGCCATTTAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAgg          <  1:416670/69‑1 (MQ=255)                                                               gCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCa        <  1:2009769/69‑1 (MQ=255)                                                                ccATTAAAATCATCAACGACAGGCTGAAACCGCAGGTTAATGtt                                <  1:924584/44‑1 (MQ=255)                                                                   ttAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTCt    <  1:180977/69‑1 (MQ=255)                                                                    tAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTCt    >  1:2176727/1‑68 (MQ=255)                                                                    tAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTCt    >  1:2025100/1‑68 (MQ=255)                                                                    tAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTCt    <  1:2239837/68‑1 (MQ=255)                                                                     aaaaTCATGAACGACAGGCTGAAACCGCAGGTTAATGtt                                <  1:3409868/39‑1 (MQ=255)                                                                     aaaaTCATGAACGACAGGCTGAAACCGCAGGTTAATGtt                                >  1:3177668/1‑39 (MQ=255)                                                                     aaaaTCATGAACGACAGGCTGAAACCGCAGGTTAATGtt                                >  1:2977230/1‑39 (MQ=255)                                                                     aaaaTCATGAACGACAGGCTGAAACCGCAGGTTAATGtt                                <  1:2120996/39‑1 (MQ=255)                                                                     aaaaTCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTCTaa  >  1:1813695/1‑69 (MQ=255)                                                                     |                                                                     CACTATGCCTGCTGGCACAACCCTGCGTTGTTATTGCGGGTCGCGCGATGCAGTTGGGCGACGCCATTAAAATCATGAACGACAGGCTGAAACCGCAGGTTAATGTTGACAGCTTCAGCCTCGAACAGGCAGTCTAA  >  W3110S.gb/1006780‑1006916

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑