Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	4,501,102	T→C	41.7%	noncoding (16/85 nt)	leuX →	tRNA‑Leu

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,501,102	0	T	C	41.7%	14.2 / 42.4	36	noncoding (16/85 nt)	leuX	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base T (6/15);  new base C (11/4);  total (17/19)
Fisher's exact test for biased strand distribution p-value = 1.66e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

GTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTT  >  W3110S.gb/4501037‑4501163                                                                  |                                                              gtggCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCgg                                                            <  1:2628155/69‑1 (MQ=255) gtggCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCgg                                                            <  1:2583318/69‑1 (MQ=255) gtggCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCgg                                                            <  1:2695294/69‑1 (MQ=255) gtggCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCgg                                                            <  1:1774291/69‑1 (MQ=255)   ggCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCgg                                                            <  1:2587609/67‑1 (MQ=255)   ggCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAACCGGTa                                                          <  1:1445121/69‑1 (MQ=255)    gCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAACCGGTa                                                          <  1:257387/68‑1 (MQ=255)              ggTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCgg                                                            <  1:3412163/56‑1 (MQ=255)              ggTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCgg                                                            <  1:2154815/56‑1 (MQ=255)              ggTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCgg                                                            <  1:1874629/56‑1 (MQ=255)                    aTCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCaaa                                         >  1:2917018/1‑69 (MQ=255)                      ccACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAACCGGTAGACGCAGTTGATTCaaaa                                        >  1:3252659/1‑68 (MQ=255)                      ccACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAACCGGTAGACGCAGTTGATTCaaaa                                        >  1:505050/1‑68 (MQ=255)                      ccACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAACCGGTAGACGCAGTTGATTCaaaa                                        >  1:1229123/1‑68 (MQ=255)                      ccACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAACCGGTAGACGAAGTTGATTCaaaa                                        >  1:105864/1‑68 (MQ=255)                              ttttCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGt                                >  1:2278563/1‑68 (MQ=255)                              ttttCCGCATACCTCTTCAGTGCCGAAGTGGCGAAACCGGTAGACGCAGTTGATTCAAAATCAACCGt                                >  1:3548193/1‑68 (MQ=255)                                   cGCATACCTCTTCAGTGCCGAAGTGGCGAAACCGGTAGACGCAGTTGATTCAAAATCAACCGTAGaaa                           <  1:3428598/68‑1 (MQ=255)                                     cATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATa                         <  1:1614412/68‑1 (MQ=255)                                         ccTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGa                             <  1:2411663/60‑1 (MQ=255)                                         ccTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGa                             <  1:1461893/60‑1 (MQ=255)                                         ccTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGa                             >  1:557762/1‑60 (MQ=255)                                         ccTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGc                    >  1:2432667/1‑69 (MQ=255)                                            ctTCAGTGCCGAAGTGGCGAAACCGGTAGACGCAGTTGATTCAAAATCa                                     <  1:1457949/49‑1 (MQ=255)                                               cAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCaa                                    <  1:1948998/47‑1 (MQ=255)                                               cAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCaa                                    <  1:1762990/47‑1 (MQ=255)                                               cAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCaa                                    <  1:1236650/47‑1 (MQ=255)                                               cAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCaa                                    <  1:710069/47‑1 (MQ=255)                                               cAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGt                >  1:3530321/1‑67 (MQ=255)                                               cAGTGCCGAAGTGGCGAAACCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGt                >  1:228218/1‑67 (MQ=255)                                                         gTGGCGAAACCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCg       >  1:1548754/1‑66 (MQ=255)                                                         gTGGCGAAACCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGcc    >  1:2281315/1‑69 (MQ=255)                                                         gTGGCGAAACCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGcc    >  1:559190/1‑69 (MQ=255)                                                         gTGGCGAAACCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGc     >  1:1975662/1‑68 (MQ=255)                                                         gTGGCGAAACCGGTAGACGCAGTTGATTCAAAATCAACAGTAGAAATACGTGCCGGTTCGAGTCCGGcc    >  1:274034/1‑69 (MQ=255)                                                            gCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCtt  >  1:875729/1‑68 (MQ=255)                                                                  |                                                              GTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTT  >  W3110S.gb/4501037‑4501163

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