Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	579,680	G→A	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	83.3%	59.1 / 8.1	30	intergenic (+371/‑377)	ybcW/nohB	hypothetical protein/DNA packaging protein
Reads supporting (aligned to +/- strand):  ref base G (5/0);  new base A (15/10);  total (20/10)
Fisher's exact test for biased strand distribution p-value = 1.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

GTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTGTATTGTCTGAAGT  >  W3110S.gb/579612‑579713                                                                     |                                  gTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCATa                                 <  1:606086/71‑1 (MQ=21) gTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCATa                                 <  1:3614985/71‑1 (MQ=21)    aGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCATAtt                               >  1:1584534/1‑70 (MQ=25)          gCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCATATTCAGTGTc                        >  1:970100/1‑71 (MQ=21)          gCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCATATTCAGTGTc                        >  1:461079/1‑71 (MQ=21)          gCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCATATTCAGTGTc                        >  1:455818/1‑71 (MQ=21)          gCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCATATTCAGTGTc                        >  1:1018204/1‑71 (MQ=21)          gCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCATATTCAGTGTc                        >  1:3918744/1‑71 (MQ=21)          gCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCATATTCAGTGTc                        >  1:2735079/1‑71 (MQ=21)          gCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCATATTCAGTGTc                        >  1:2785245/1‑71 (MQ=21)                   gCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGAttt                <  1:4125906/70‑1 (MQ=255)                   gCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGAttt                <  1:282233/70‑1 (MQ=255)                   gCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTg               >  1:4123375/1‑71 (MQ=255)                   gCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTg               >  1:3898238/1‑71 (MQ=255)                    cGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTg               <  1:3707724/70‑1 (MQ=255)                    cGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGACCCATATTCAGTGTCGCTGATTTGt              <  1:1897476/71‑1 (MQ=255)                    cGGGTTGTTGTTCTGCGGGTTCTGTTCTCAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTGt              <  1:1778393/71‑1 (MQ=11)                      ggTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCATa                                 <  1:2599455/50‑1 (MQ=255)                             gttCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGAttt                >  1:2870563/1‑60 (MQ=35)                             gttCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTGTATTGTCTg      >  1:3091661/1‑70 (MQ=35)                             gttCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTGTATTGTCTg      >  1:2170767/1‑70 (MQ=35)                             gttCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTGTATTGTCTg      >  1:871815/1‑70 (MQ=35)                             gttCTGCGGGATCTGTTCTTAGTTGACATGAGGTTGCCCAGTATTCAGTGTCGCTGATTTGTa             >  1:1259725/1‑63 (MQ=11)                                cTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTg               >  1:159042/1‑58 (MQ=255)                                cTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTg               >  1:1111733/1‑58 (MQ=255)                                cTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTGTATTGTCTGAAgt  >  1:3841424/1‑71 (MQ=21)                                 tGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTGTATTg          <  1:3644045/62‑1 (MQ=255)                                 tGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTGTATTg          <  1:2207121/62‑1 (MQ=255)                                   cGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTGTATTGTc        >  1:2430686/1‑62 (MQ=255)                                   cGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCATATTCAGTGTCGCTGATTTGTATTGTc        >  1:2127605/1‑62 (MQ=255)                                                                     |                                  GTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTAGTTGACATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTGTATTGTCTGAAGT  >  W3110S.gb/579612‑579713

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