Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	45,813	C→T	56.5%	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	56.5%	14.3 / 56.1	46	P3S (CCG→TCG)	yaaU	predicted transporter
Reads supporting (aligned to +/- strand):  ref base C (12/8);  new base T (10/16);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 2.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATT  >  W3110S.gb/45746‑45863                                                                    |                                                   gCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAATCgg                                                 <  1:1390630/71‑2 (MQ=255) gCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAATCGt                                                 <  1:516754/71‑1 (MQ=255) gCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAATCGt                                                 <  1:104878/71‑1 (MQ=255) gCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAATCGt                                                 <  1:2498788/71‑1 (MQ=255) gCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAATCGt                                                 <  1:1863399/71‑1 (MQ=255) gCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAATCGt                                                 <  1:3004543/71‑1 (MQ=255) gCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAATCGt                                                 <  1:501929/71‑1 (MQ=255) gCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAATCGt                                                 <  1:722728/71‑1 (MQ=255)      tGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGa                                            <  1:1183270/71‑1 (MQ=255)                     tAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAATCGTCCAGAAACTTTGACGAtct                              >  1:3086807/1‑70 (MQ=255)                       aCGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCa                            <  1:1260437/70‑1 (MQ=255)                       aCGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCa                            <  1:1185169/70‑1 (MQ=255)                       aCGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCa                            <  1:1821572/70‑1 (MQ=255)                       aCGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCa                            <  1:2628968/70‑1 (MQ=255)                       aCGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCa                            <  1:251538/70‑1 (MQ=255)                       aCGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCa                            <  1:1132550/70‑1 (MQ=255)                         gccgcACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCaaa                          >  1:682184/1‑70 (MQ=255)                         gccgcACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCaaa                          >  1:2386637/1‑70 (MQ=255)                         gccgcACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCaaa                          >  1:744434/1‑70 (MQ=255)                         gccgcACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCaaa                          >  1:1840956/1‑70 (MQ=255)                         gccgcACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTAGACGATCTCaaa                          >  1:2621278/1‑70 (MQ=255)                               cTGACTCTCATTGCAAAAAACAGGAATAACCATGCAATCGt                                                 <  1:1294582/41‑1 (MQ=255)                                 gACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTctcctc                  >  1:1513648/1‑70 (MQ=255)                                  aCTCTCATTGCAAAAAACAGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATTctcctc                  <  1:3052554/69‑1 (MQ=255)                                  aCTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTa                <  1:2343628/71‑1 (MQ=255)                                   ctctCATTGCAAAAAACAGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATtctc                     <  1:843626/65‑1 (MQ=255)                                   ctctCATTGCAAAAAACAGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATtctc                     <  1:1139384/65‑1 (MQ=255)                                   ctctCATTGCAAAAAACAGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATtctc                     <  1:521853/65‑1 (MQ=255)                                   ctctCATTGCAAAAAACAGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATtctc                     <  1:1946088/65‑1 (MQ=255)                                   ctctCATTGCAAAAAACAGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATtctc                     <  1:298622/65‑1 (MQ=255)                                   ctctCATTGCAAAAAACAGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATtctc                     <  1:1520626/65‑1 (MQ=255)                                           gCAAAAAACAGGAATAACCATGCAATCGTCCAGAAACTTTGACGAt                                >  1:331867/1‑46 (MQ=255)                                           gCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCAccgcc       >  1:2097577/1‑71 (MQ=255)                                           gCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCAccgcc       >  1:1466764/1‑71 (MQ=255)                                           gCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCAccgcc       >  1:2696533/1‑71 (MQ=255)                                           gCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCAccgcc       >  1:250106/1‑71 (MQ=255)                                           gCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCAccgc        >  1:2412781/1‑70 (MQ=255)                                           gCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCAccgc        >  1:2098128/1‑70 (MQ=255)                                                aaaCAGGAATAACCATGCAATCTTCCAGAAACTTTGACg                                  >  1:1420969/1‑39 (MQ=38)                                                aaaCAGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCAtt  >  1:1784660/1‑71 (MQ=255)                                                aaaCAGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCAtt  >  1:324595/1‑71 (MQ=255)                                                aaaCAGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCAtt  >  1:1830860/1‑71 (MQ=255)                                                aaaCAGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCAtt  >  1:668155/1‑71 (MQ=255)                                                aaaCAGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCAtt  >  1:1984489/1‑71 (MQ=255)                                                aaaCAGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCAtt  >  1:75034/1‑71 (MQ=255)                                                aaaCAGGAATAACCATGCAATCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCAtt  >  1:2004935/1‑71 (MQ=255)                                                                    |                                                   GCTGATGTTGGTTTGATACGTAACGCCGCACTGACTCTCATTGCAAAAAACAGGAATAACCATGCAACCGTCCAGAAACTTTGACGATCTCAAATTCTCCTCTATTCACCGCCGCATT  >  W3110S.gb/45746‑45863

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