Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	45,813	C→T	66.7%	P3S (CCG→TCG)	yaaU →	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	45,813	0	C	T	66.7%	27.1 / 23.0	27	P3S (CCG→TCG)	yaaU	predicted transporter
Reads supporting (aligned to +/- strand):  ref base C (7/2);  new base T (3/15);  total (10/17)
Fisher's exact test for biased strand distribution p-value = 3.65e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGG  >  W3110S.gb/45757‑45880                                                         |                                                                    tttGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACttt                                                       <  1:1334821/71‑1 (MQ=255) tttGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACttt                                                       <  1:139842/71‑1 (MQ=255)            aCGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCa                                             <  1:248836/70‑1 (MQ=255)             cgccgcACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCaaa                                           <  1:36434/71‑1 (MQ=255)             cgccgcACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCaaa                                           <  1:1190734/71‑1 (MQ=255)             cgccgcACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCaaa                                           <  1:1192996/71‑1 (MQ=255)             cgccgcACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCaaa                                           <  1:1689969/71‑1 (MQ=255)             cgccgcACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCaaa                                           <  1:1454876/71‑1 (MQ=255)                      gACTCTCATTGTAAAAAACAGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCaaa                                           >  1:1477768/1‑62 (MQ=255)                       aCTCTCATTGCAAAAAACAGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATTctcctc                                   <  1:908777/69‑1 (MQ=255)                       aCTCTCATTGCAAAAAACAGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATTctcctc                                   <  1:414749/69‑1 (MQ=255)                            cATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATtc                                        >  1:1376526/1‑59 (MQ=255)                              ttGCAAAAAACAGGAATAACCATGCAATCGTCCAGAAACTTTGAc                                                    <  1:1295331/45‑1 (MQ=255)                              ttGCAAAAAACAGGAATAACCATGCAATCGTCCAGAAACTTTGAc                                                    <  1:1096693/45‑1 (MQ=255)                                gCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCAccgcc                        >  1:1214243/1‑71 (MQ=255)                                gCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCAccgcc                        >  1:226818/1‑71 (MQ=255)                                gCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCAccgcc                        >  1:1121668/1‑71 (MQ=255)                                        cAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTg                >  1:1567000/1‑71 (MQ=255)                                        cAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTg                >  1:1141155/1‑71 (MQ=255)                                        cAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTg                >  1:878193/1‑71 (MQ=255)                                         aGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGc               <  1:1479303/71‑1 (MQ=255)                                         aGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGc               <  1:1695361/71‑1 (MQ=255)                                         aGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGc               <  1:438756/71‑1 (MQ=255)                                         aGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGc               <  1:815925/71‑1 (MQ=255)                                         aGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGc               <  1:1015935/71‑1 (MQ=255)                                                      cAATCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCgg  >  1:1373015/1‑71 (MQ=255)                                                      cAATCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCg   >  1:328078/1‑70 (MQ=255)                                                         |                                                                    TTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATTTTGCTGTGGGGAAGCGG  >  W3110S.gb/45757‑45880

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