Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,097,480	0	A	G	65.4%	24.9 / 18.4	26	intergenic (+651/+85)	ycdU/serX	putative inner membrane protein/tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (5/4);  new base G (7/10);  total (12/14)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

GTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTC‑GGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCG  >  NC_000913/1097331‑1097623                                                                                                                                                      |                                                                                                                                                  gtgCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAAc                                                                                                                                                  <  2:458996/151‑1 (MQ=31)  tgCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGaa                                                                                                                                                   <  2:530260/149‑1 (MQ=11)   gCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTGACTCGCTTCCCTCCCCCCGCGGGGAGCCCCCCCACTTCGTTCACCGTCTGTCCCACTGGCCGTCGCCCGTTGTCCAAccc                                                                                                                                                >  2:814364/1‑151 (MQ=11)   gCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAAcc                                                                                                                                                 <  1:375585/150‑1 (MQ=31)     aaTATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTGACTCGCTTCGCTCGCCCCGCGGGCAGCCCCCCCACTGCGTTCCCGGGCTGTCCAACTGGCTGTCGCCCGTTGTCGAAaccc                                                                                                                                               >  2:1662824/1‑150 (MQ=11)        ataCGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCCCTCACTGCGTTCACGGTCTGTCCCACTGGCTGTCGCCAGTTGTCGAAcccccg                                                                                                                                             >  2:552118/1‑147 (MQ=9)        ataCGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTGGAACCCCGGTc                                                                                                                                           <  1:198026/151‑1 (MQ=14)        ataCGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGt                                                                                                                                            <  2:1001911/150‑1 (MQ=31)                                                                                                          cACTGCGTTCACGGTCTGGCCAACTGGGTGGCGCCCGTTGTTGGGGCCCGGTCGGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTc                                          <  2:683042/150‑1 (MQ=11)                                                                                                               cGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTC‑GGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGtttacactcc                                   >  1:370080/1‑143 (MQ=0)                                                                                                               cGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTC‑GGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAgggggggggttactcgc                                    >  1:557447/1‑140 (MQ=2)                                                                                                                   cACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAGCCCCGGTCGGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAAcccccc                                 >  2:898722/1‑149 (MQ=255)                                                                                                                    aCGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAGCCCCGGTCGGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGa                                <  1:1662716/150‑1 (MQ=255)                                                                                                                      ggTCTGTCCAACTGGCTGTCGCCAGTTGTCGAGCCCCGGTCGGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGGGGGGGGGGGTTGAAaccccccaa                              >  2:1158272/1‑146 (MQ=1)                                                                                                                           gTCCAACTGGCTGTCGCCAGTTGTCGAGCCCCGGTCGGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAACCCCCCATACgct                          >  2:214536/1‑147 (MQ=11)                                                                                                                           gTCCAACTGGCTGTCGCCAGTTGTCGAGCCCCGGTCGGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTg                         <  1:566261/150‑1 (MQ=255)                                                                                                                                       gTCGCCAGTTGTCGAGCCCCGGTCGGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATacaca              <  2:638150/149‑1 (MQ=255)                                                                                                                                        tCGCCAGTTGTCGAGCCCCGGTCGGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑Ct           <  2:1194821/150‑1 (MQ=255)                                                                                                                                         cGCCAGTTGTCGAGCCCCGGTCGGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑Cttt         <  2:687674/151‑1 (MQ=255)                                                                                                                                         cGCCAGTTGTCGAGCCCCGGTCGGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGGGGGGGGGGGTCTAACCCCCCGCACCTTGCCGTATACACA‑Cttt         >  1:183758/1‑151 (MQ=14)                                                                                                                                           ccAGTTGTCGAGCCCCGGTCGGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGGTCGAAACCCCCATACGTTGCCGTATACACA‑Cctt         >  1:1119297/1‑149 (MQ=14)                                                                                                                                           ccAGTTGTCGAGCCCCGGTCGGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTcc       <  2:96335/151‑1 (MQ=255)                                                                                                                                           ccAGTTGTCGAGCCCCGGTCGGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTc        <  1:202410/150‑1 (MQ=255)                                                                                                                                            cAGTTGTCGAGCCCCGGTCGGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCa      <  1:1630340/151‑1 (MQ=255)                                                                                                                                              gTTGTCGAGCCCCGGTCGGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTTGTACCCCCCATACCCTTCCCTATACCCACCTTTccgg     >  1:547417/1‑149 (MQ=11)                                                                                                                                               ttGTCGAGCCCCGGTCGGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGTTCGGAACCCCCCTACGGTGTCGTATAAACA‑Cttttccgg    >  1:1079625/1‑145 (MQ=14)                                                                                                                                                 gTCGAGCCCCGGTCGGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCg  <  2:1294028/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                  GTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTC‑GGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCG  >  NC_000913/1097331‑1097623

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