Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,865,246	0	A	G	41.9%	66.2 / 12.0	44	E49G (GAA→GGA)	ygbN	putative transporter
Reads supporting (aligned to +/- strand):  ref base A (1/24);  new base G (18/0);  total (20/24)
Fisher's exact test for biased strand distribution p-value = 3.12e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATGTCCACAATTACATTGTTATGCATTGCGTTAGCTGGCGTAATCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCG‑CCGGGATGGGCGGTGT‑GCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGCTGGGGCGTATGATCGAACACTCTGGTGGTGCAGAGTCACTGGCTAATTATTTCAGTCGCAAGTTAGGTG  >  NC_000913/2865100‑2865392                                                                                                                                                   |                                                                                                                                                     cATGTCCACAATTACATTGTTATGCATTGCGTTAGCTGGCGTAATCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGGAggg                                                                                                                                                  >  1:474301/1‑149 (MQ=255)  aTGTCCACAATTACATTGTTATGCATTGCGTTAGCTGGCGTAATCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGGAgggg                                                                                                                                                 >  1:435893/1‑151 (MQ=255)   tGTCCACAATTACATTGTTATGCATTGCGTTAGCTGGCGTAATCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGGAGgggg                                                                                                                                                >  2:326689/1‑151 (MQ=255)    gTCCACAATTACATTGTTATGCATTGCGTTAGCTGGCGTAATCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGGGGAGggggg                                                                                                                                               >  2:338006/1‑150 (MQ=255)    gTCCACAATTACATTGTTATGCATTGCGTTAGCTGGCGTAATCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGTggg                                                                                                                                                <  2:708762/150‑1 (MQ=255)               cATTGTTATGCATTGCGTTAGCTGGCGTAATCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGGAGggggg                                                                                                                                               >  2:394724/1‑139 (MQ=255)                      aTGCATTGCGTTAGCTGGCGTAATCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCG‑CCggg                                                                                                                            <  2:24038/151‑1 (MQ=255)                        gCATTGCGTTAGCTGGCGTAATCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCG‑CCGGGa                                                                                                                           <  2:126194/150‑1 (MQ=255)                          aTTGCGTTAGCTGGCGTAATCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAGCCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCG‑CCGGGa                                                                                                                           <  1:26367/148‑1 (MQ=255)                               gTTAGCTGGCGTAATCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGGAGGGGGGTAAGAGGTGATCG‑CCggggggggggg                                                                                                                    >  1:323073/1‑142 (MQ=255)                                 tAGCTGGCGTAATCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGGGGAGTGGGGGAAGTGGTGATCG‑CCgggggtggggggg                                                                                                                  >  2:42130/1‑140 (MQ=255)                                   gCTGGCGTAATCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGGAGGGGGGAAAAGGATGGTCG‑CCggggggggggggggg                                                                                                                >  1:150431/1‑138 (MQ=255)                                     tGGCGTAATCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGGAGGGGGGAAAAGGATGGTCG‑CCggggggggggggggggct                                                                                                             >  1:301368/1‑136 (MQ=255)                                      ggCGTAATCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGGAGGGGGGAAAGTGATTATCG‑CCGGGGGGGGGGGGGG‑GCtt                                                                                                            >  1:432298/1‑150 (MQ=255)                                      ggCGTAATCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCG‑CCGGGATGGGCGGTGT‑GCTTg                                                                                                           <  1:141592/151‑1 (MQ=255)                                         gTAATCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGGGGAGGGGGGGAAGGGATGAACG‑CCCGGGTggggggggggggttg                                                                                                           >  1:694373/1‑140 (MQ=255)                                              cATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCG‑CCGGGATGGGCGGTGT‑GCTTGGTTCCGt                                                                                                    <  2:406690/150‑1 (MQ=255)                                                 gCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGGGGAGTGGGGGAAGAGATGATCG‑CCGGGGTGGGGGGGGT‑GGTTTGGTTCGTCa                                                                                                  >  1:335031/1‑149 (MQ=255)                                                 gCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCG‑CCGGGATGGGCGGTGT‑GCTTGGTTCCGTCACt                                                                                                <  2:701681/151‑1 (MQ=255)                                                   tGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGGAGGGGGGGAAAAGGTGATTG‑CCGGGGGGGGGGGGGG‑GCTTgggt                                                                                                        >  1:698763/1‑139 (MQ=255)                                                    gTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCG‑CCGGGATGGGCGGTGT‑GCTTGGTTCCGTCACtatt                                                                                             <  1:317192/151‑1 (MQ=255)                                                          gctgGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCG‑CCGGGATGGGCGGTGT‑GCTTGGTTCCGTCACTattattatt                                                                                       <  2:123130/151‑1 (MQ=255)                                                                       gCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGGAGTGGGGAAAGAGATGATCG‑CCCGGGGGGGCGGGGTGGCTTGGTT‑CGTCACTatttttatt                                                                                       >  1:298795/1‑138 (MQ=255)                                                                            ggTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCG‑CCGGGATGGGCGGTGT‑GCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGCTg                                                                     <  2:277594/151‑1 (MQ=255)                                                                               aCAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCG‑CCGGGATGGGCGGTGT‑GCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGCTggg                                                                   <  1:649287/150‑1 (MQ=255)                                                                                 aaCCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCG‑CCGGGATGGGCGGTGT‑GCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGCTGGGGCg                                                                <  1:162360/151‑1 (MQ=255)                                                                                    cATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGGGGAGAGGGGTAAAAGAAGATCG‑GCGGGGTGGGGGGGGG‑GGTTTGGTTCGTACCTATTTTTTTTTGTTTTGGGGCCATGCTgggggg                                                                >  1:289/1‑146 (MQ=255)                                                                                       tCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCG‑CCGGGATGGGCGGTGT‑GCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGCTGGGGCGTATGa                                                           <  1:236873/150‑1 (MQ=255)                                                                                         gTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGGAGGGGGGAAAAGGATGATTG‑CCGGGGTGGGGGGGGG‑GGTTTGTTTCGGCCCTATTTTTATTTGGTTGGGGGCTATGCTgggggg                                                                >  1:489277/1‑141 (MQ=255)                                                                                           tgctttgctCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCG‑CCGGGATGGGCGGTGT‑GCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGCTGGGGCGTATGATCGa                                                       <  2:208304/150‑1 (MQ=255)                                                                                                       cGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGGAGGGGGGGAAAGGGTTATCGCCCGGGGGGGGGGGGGG‑GCTTGGTT‑CGTCTCCCTTTTTTTTTTTTTGGGGGGCATGCTGGGGGGGGTGATCGAAAACCCtggtgg                                           >  2:648461/1‑150 (MQ=255)                                                                                                        gTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGCAGGGGGGAAAGAGATGATCG‑CCCGGGGGGGCCGGGT‑GCTTTGTTCCCGCACTTTTTTTTTTTGGGTGGGGGCTATTCTGGGGGGGATGATCGAACACTCctgtgggg                                         >  1:86166/1‑149 (MQ=255)                                                                                                             ccTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCG‑CCGGGATGGGCGGTGT‑GCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGCTGGGGCGTATGATCGAACACTCTGGTGGTGCagag                                    <  1:584142/151‑1 (MQ=255)                                                                                                                gTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGGAGGGGGGGAAGGGATGAACG‑ACGGGGGGGGCGGGGG‑GGTTTGTTCCGTCCCCATTTTTTTTTGTCTGGGGGCTCTGCTGGGGGGTATggtgcgacactcc                                                >  2:646412/1‑124 (MQ=255)                                                                                                                   aGTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCG‑CCGGGATGGGCGGTGT‑GCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGCTGGGGCGTATGATCGAACACTCTGGTGGTGCAGAGTCACt                               <  2:489274/150‑1 (MQ=255)                                                                                                                    gTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCG‑CCGGGATGGGCGGTGT‑GCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGCTGGGGCGTATGATCGAACACTCTGGTGGTGCAGAGTCACTgg                             <  2:258867/151‑1 (MQ=255)                                                                                                                    gTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCG‑CCGGGATGGGCGGTGT‑GCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGCTGGGGCGTATGATCGAACACTCTGGTGGTGCAGAGTCACTgg                             <  1:169441/151‑1 (MQ=255)                                                                                                                      cGCACTTGCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCG‑CCGGGATGGGCGGTGT‑GCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGCTGGGGCGTATGATCGAACACTCTGGTGGTGCAGAGTCACTGGc                            <  2:340301/150‑1 (MQ=255)                                                                                                                                  aGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCG‑CCGGGATGGGCGGTGT‑GCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGCTGGGGCGTATGATCGAACACTCTGGTGGTGCAGAGTCACTGGCTAATTATTTCAGt               <  2:323070/151‑1 (MQ=255)                                                                                                                                   ggTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCG‑CCGGGATGGGCGGTGT‑GCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGCTGGGGCGTATGATCGAACACTCTGGTGGTGCAGAGTCACTGGCTAATTATTTCAg                <  2:52961/149‑1 (MQ=255)                                                                                                                                    gTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCG‑CCGGGATGGGCGGTGT‑GCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGCTGGGGCGTATGATCGAACACTCTGGTGGTGCAGAGTCACTGGCTAATTATTTCAGTCg             <  1:235868/151‑1 (MQ=255)                                                                                                                                          cggcggGCGAAGTGGGTAAAGTGATGATCG‑CCGGGATGGGCGGTGT‑GCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGCTGGGGCGTATGATCGAACACTCTGGTGGTGCAGAGTCACTGGCTAATTATTTCAGTCGCAAGtt       <  1:326692/151‑1 (MQ=255)                                                                                                                                              gggcgAAGTGGGTAAAGTGATGATCG‑CCGGGATGGGCGGTGT‑GCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGCTGGGGCGTATGATCGAACACTCTGGTGGTGCAGAGTCACTGGCTAATTATTTCAGTCGCAAGTTAGGt   <  2:708934/151‑1 (MQ=255)                                                                                                                                               ggcgAAGTGGGTAAAGTGATGATCG‑CCGGGATGGGCGGTGT‑GCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGCTGGGGCGTATGATCGAACACTCTGGTGGTGCAGAGTCACTGGCTAATTATTTCAGTCGCAAGTTAGGTg  >  2:718312/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                     CATGTCCACAATTACATTGTTATGCATTGCGTTAGCTGGCGTAATCATGCTGTTGCTGCTGGTCATCAAGGCAAAGGTACAACCATTCGTTGCTTTGCTCCTCGTCAGCCTGTTAGTCGCACTTGCGGCAGGTATACCGGCGGGCGAAGTGGGTAAAGTGATGATCG‑CCGGGATGGGCGGTGT‑GCTTGGTTCCGTCACTATTATTATTGGTCTGGGCGCTATGCTGGGGCGTATGATCGAACACTCTGGTGGTGCAGAGTCACTGGCTAATTATTTCAGTCGCAAGTTAGGTG  >  NC_000913/2865100‑2865392

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