Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	393,168	T→C	82.4%	G310G (GGT→GGC)	yaiU →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	393,168	0	T	C	82.4%	122.0 / 18.5	51	G310G (GGT→GGC)	yaiU	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (3/6);  new base C (20/22);  total (23/28)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

TCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCA  >  W3110S.gb/393101‑393235                                                                    |                                                                    tCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGCACt                                                                  <  1:3480569/71‑1 (MQ=255) tCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGCACt                                                                  <  1:829665/71‑1 (MQ=255) tCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGCACt                                                                  <  1:602382/71‑1 (MQ=255) tCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGCACt                                                                  <  1:523062/71‑1 (MQ=255)  cTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGCACt                                                                  <  1:3266982/70‑1 (MQ=255)   tttGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGCACtt                                                                 >  1:2333474/1‑70 (MQ=255)   tttGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGCACTTa                                                                >  1:3572402/1‑71 (MQ=255)   tttGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGCACTTa                                                                >  1:457077/1‑71 (MQ=255)     tgttgATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTAc                                                               >  1:3479949/1‑70 (MQ=255)     tgttgATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACg                                                              >  1:1943162/1‑71 (MQ=255)          aTGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGAc                                                         <  1:1013977/71‑1 (MQ=255)          aTGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGAc                                                         <  1:361216/71‑1 (MQ=255)             gTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACgg                                                       <  1:531003/70‑1 (MQ=255)             gTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACgg                                                       <  1:3016556/70‑1 (MQ=255)             gTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGCACTTACGTTGACgg                                                       <  1:2467106/70‑1 (MQ=255)             gTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGCACTTACGTTGACgg                                                       <  1:392605/70‑1 (MQ=255)             gTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGCACTTACGTTGACgg                                                       <  1:2865082/70‑1 (MQ=255)             gTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGCACTTACGTTGACgg                                                       >  1:2942404/1‑70 (MQ=255)                gCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGCACTTACGTTGACGGTAg                                                    <  1:447348/70‑1 (MQ=255)                gCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGCACTTACGTTGACGGTAg                                                    <  1:3361492/70‑1 (MQ=255)                    gAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGCACTTACGTTGACGGTAGCAcc                                                >  1:2251719/1‑70 (MQ=255)                    gAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGCACTTACGTTGACGGTAGCAcc                                                <  1:3617792/70‑1 (MQ=255)                    gAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGCACTTACGTTGACGGTAGCAcc                                                <  1:162687/70‑1 (MQ=255)                             tcACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACg                                      >  1:1215885/1‑71 (MQ=255)                             tcACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGCACTTACGTTGACGGTAGCACCAACTCCGACg                                      >  1:2026746/1‑71 (MQ=255)                             tcACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGCACTTACGTTGACGGTAGCACCAACTCCGACg                                      >  1:400465/1‑71 (MQ=255)                                  tCAACAACGACCTGTCTGCAACCATGAGCAACGGCACTTACGTTGACGGTAGCACCAACTCCGACGCTTgg                                 >  1:1230095/1‑71 (MQ=255)                                   caacaaCGACCTGTCTGCAACCATGAGCAACGGCACTTACGTTGTCGGTAGCACCAACTCCGACGCTTggg                                >  1:1978706/1‑71 (MQ=255)                                   caacaaCGACCTGTCTGCAACCATGAGCAACGGCACTTACGTTGACGGTAGCACCAACTCCGACGCTTggg                                >  1:2940990/1‑71 (MQ=255)                                      caacGACCTGTCTGCAACCATGAGCAACGGCACTTAc                                                               <  1:519164/37‑1 (MQ=255)                                      caacGACCTGTCTGCAACCATGAGCAACGGCACTTAc                                                               <  1:1018823/37‑1 (MQ=255)                                      caacGACCTGTCTGCAACCATGAGCAACGGCACTTAc                                                               <  1:1524590/37‑1 (MQ=255)                                      caacGACCTGTCTGCAACCATGAGCAACGGCACTTAc                                                               <  1:1800329/37‑1 (MQ=255)                                        acGACCTGTCTGCAACCATGAGCAACGGCACTTACGt                                                             >  1:857747/1‑37 (MQ=255)                                          gACCTGTCTGCAACCATGAGCAACGGCACTTACGTTGAc                                                         <  1:918512/39‑1 (MQ=255)                                          gACCTGTCTGCAACCATGAGCAACGGCACTTACGTTGAc                                                         <  1:2922343/39‑1 (MQ=255)                                                    cAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGcc               <  1:1655732/71‑1 (MQ=255)                                                    cAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGcc               <  1:722400/71‑1 (MQ=255)                                                     aaCCATGAGCAACGGCACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCg              <  1:2304941/71‑1 (MQ=255)                                                     aaCCATGAGCAACGGCACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCg              <  1:3505193/71‑1 (MQ=255)                                                        cATGAGCAACGGCACTTACGTTGACGGTAGCACCAACTCCGGCGCtt                                   >  1:3635441/1‑47 (MQ=39)                                                        cATGAGCAACGGCACTTACGTTGACGGTAGCACCAACTCCGACGCtt                                   >  1:2695524/1‑47 (MQ=255)                                                        cATGAGCAACGGCACTTACGTTGACGGTAGCACCAACTCCGACGCtt                                   >  1:1814429/1‑47 (MQ=255)                                                        cATGAGCAACGGCACTTACGTTGACGGTAGCACCAACTCCGACGCtt                                   >  1:285491/1‑47 (MQ=255)                                                        cATGAGCAACGGCACTTACGTTGACGGTAGCACCAACTCCGACGCtt                                   >  1:1588665/1‑47 (MQ=255)                                                        cATGAGCAACGGCACTTACGTTGACGGTAGCACCAACTCCGACGCtt                                   >  1:599698/1‑47 (MQ=255)                                                            aGCAACGGCACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGa       >  1:1909521/1‑71 (MQ=255)                                                            aGCAACGGCACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGa       >  1:574524/1‑71 (MQ=255)                                                            aGCAACGGCACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGa       >  1:358294/1‑71 (MQ=255)                                                                 cGGCACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCa  <  1:157436/71‑1 (MQ=255)                                                                 cGGCACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGCCTTCGGTTTGAAAGCCGGTTACGACTTCa  <  1:2850335/71‑1 (MQ=255)                                                                    |                                                                    TCTTTGTTGATGGTAGCTTGAGCTACTCTCACTTCAACAACGACCTGTCTGCAACCATGAGCAACGGTACTTACGTTGACGGTAGCACCAACTCCGACGCTTGGGGCTTCGGTTTGAAAGCCGGTTACGACTTCA  >  W3110S.gb/393101‑393235

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