Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,097,480	0	A	G	68.4%	12.3 / 10.5	19	intergenic (+651/+85)	ycdU/serX	putative inner membrane protein/tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (2/4);  new base G (4/9);  total (6/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.90e-01

GTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGC  >  NC_000913/1097331‑1097626                                                                                                                                                      |                                                                                                                                                   gtgCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCCTCGCCCGCCCTGCGGGCAGCCCGCTCACTGCGTTCCCGGTCTGTCCAACTGGGTGTCGCCCGTTGTCGAAc                                                                                                                                                   >  2:593789/1‑151 (MQ=255)    cAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAAcccc                                                                                                                                                <  1:263275/151‑1 (MQ=31)    cAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAAccc                                                                                                                                                 <  1:185743/150‑1 (MQ=31)        ataCGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGt                                                                                                                                             <  2:3809/150‑1 (MQ=31)         taCGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTc                                                                                                                                            <  1:630034/150‑1 (MQ=11)                                                                                                              gCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAGCCCCGGTCTGTGCTTTTCATCCCTCCGGTGTGTGAAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAAccc                                      <  1:523968/151‑1 (MQ=11)                                                                                                                   cACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGa                                  <  2:92747/150‑1 (MQ=38)                                                                                                                           gTCCAACTGGCTGTCGCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGcc                         <  2:379664/151‑1 (MQ=255)                                                                                                                           gTCCAACTGGCTGTCGCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGGGGGGGGGGATCGAAACCCCCATACCtt                            >  1:84909/1‑148 (MQ=11)                                                                                                                           gTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTAATccccccgccggccc                           >  1:35407/1‑140 (MQ=1)                                                                                                                                   ggCTGTCGCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATac                   <  2:84910/149‑1 (MQ=255)                                                                                                                                    gCTGTCGCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAAACCCCCATACGTTGCCGTATaca                  >  2:125400/1‑149 (MQ=255)                                                                                                                                       gTCGCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATacacac               <  2:70287/149‑1 (MQ=255)                                                                                                                                       gTCGCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACt              <  1:204676/150‑1 (MQ=255)                                                                                                                                          gCCAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTcc          <  2:274589/151‑1 (MQ=255)                                                                                                                                           ccAGTTGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCa         <  2:74807/151‑1 (MQ=255)                                                                                                                                                tGTCGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGGGGGGGGGGATTGAACCCCCGATACGTTGCCGTATAAACACTTTTCCGGCg     >  2:443654/1‑150 (MQ=17)                                                                                                                                                   cGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGc  <  2:35406/150‑1 (MQ=255)                                                                                                                                                   cGAGCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGGGGGGGGGATTTGAAACCCCCCTACCTTTCCGTCTACCCACTTTcc          >  1:254874/1‑142 (MQ=14)                                                                                                                                                      |                                                                                                                                                   GTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGC  >  NC_000913/1097331‑1097626

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑