Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	579,680	G→A	66.7%	intergenic (+371/‑377)	ybcW → / → nohB	hypothetical protein/DNA packaging protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	579,680	0	G	A	66.7%	23.3 / 23.3	30	intergenic (+371/‑377)	ybcW/nohB	hypothetical protein/DNA packaging protein
Reads supporting (aligned to +/- strand):  ref base G (6/4);  new base A (9/11);  total (15/15)
Fisher's exact test for biased strand distribution p-value = 7.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGTTTTTACGTT  >  W3110S.gb/579616‑579725                                                                 |                                              ggAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCGTATTc                                          <  1:775177/70‑1 (MQ=35) ggAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCGTATTc                                          <  1:3291452/70‑1 (MQ=35) ggAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCGTATTc                                          <  1:2718927/70‑1 (MQ=35) ggAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCGTATTc                                          <  1:229624/70‑1 (MQ=35)                cGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGGCCCATATTCAGTGTCGCTGATTTg                           <  1:757627/70‑1 (MQ=255)                cGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTg                           <  1:1143167/70‑1 (MQ=255)                cGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTg                           <  1:1551445/70‑1 (MQ=255)                cGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTg                           <  1:700335/70‑1 (MQ=255)                cGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTg                           <  1:2321522/70‑1 (MQ=255)                 gggTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTg                           <  1:1506107/69‑1 (MQ=255)                         gttCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTGTATTGTCTg                  >  1:2429250/1‑70 (MQ=35)                         gttCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTGTATTGTCTg                  >  1:843136/1‑70 (MQ=35)                         gttCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTGTATTGTCTg                  >  1:2822924/1‑70 (MQ=35)                         gttCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTGTATTGTCTg                  >  1:213447/1‑70 (MQ=35)                         gttCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTGTATTGTCTg                  >  1:1904440/1‑70 (MQ=35)                         gttCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTGTATTGTCTg                  >  1:2222487/1‑70 (MQ=35)                                    tctgttctTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTGTATTGTCTGAAg               >  1:2005528/1‑62 (MQ=25)                                    tctgttctTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTGTATTGTCTGAAg               >  1:900384/1‑62 (MQ=25)                                    tctgttctTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTGTATTGTCTGAAg               >  1:1607129/1‑62 (MQ=25)                                    tctgttctTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTGTATTGTCTGAAg               >  1:756198/1‑62 (MQ=25)                                      tgttctTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGtttt        <  1:3233929/67‑1 (MQ=25)                                      tgttctTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGtttt        <  1:2122715/67‑1 (MQ=25)                                      tgttctTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGtttt        <  1:2666132/67‑1 (MQ=25)                                      tgttctTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGtttt        <  1:1735845/67‑1 (MQ=25)                                      tgttctTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGtttt        <  1:2380185/67‑1 (MQ=25)                                          ctTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGTTTTTACGtt  >  1:3324252/1‑69 (MQ=25)                                          ctTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGTTTTTACGtt  >  1:429040/1‑69 (MQ=25)                                          ctTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGTTTTTACGtt  >  1:1162917/1‑69 (MQ=25)                                          ctTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGTTTTTACGt   >  1:3200746/1‑68 (MQ=25)                                          ctTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGTTTTTACGt   >  1:1127738/1‑68 (MQ=25)                                                                 |                                              GGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGTTTTTACGTT  >  W3110S.gb/579616‑579725

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