Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,137,615	0	T	G	46.3%	43.8 / 12.2	41	T15P (ACC→CCC)	yijF	DUF1287 family protein
Reads supporting (aligned to +/- strand):  ref base T (5/17);  new base G (19/0);  total (24/17)
Fisher's exact test for biased strand distribution p-value = 3.15e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGATAAGTAAGCTGCACATACGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGT‑GAAGCTTTCATCTGGTTCG‑TCCGTG‑AAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGAAGCGGATAAAAAATTCCTGAACAACGTGGG  >  NC_000913/4137468‑4137764                                                                                                                                                    |                                                                                                                                                         cAGGATAAGTAAGCTGCACATCCGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGCCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCTCGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGTaag                                                                                                                                                       <  2:422152/151‑2 (MQ=255)  aGGATAAGTAAGCTGCACATACGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGGGAATTTAAGGGGTGGGTGGGaaa                                                                                                                                                       >  2:168291/1‑150 (MQ=255)   ggATAAGTAAGCTGCACATACGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGCATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTATGGTGTGGCTGGTgaa                                                                                                                                                       >  2:241528/1‑149 (MQ=255)         gtaagCTGCACATACGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGTAAATGCCGTTa                                                                                                                                               <  2:495206/151‑1 (MQ=255)           aagCTGCACATACGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGGGATTTTAAGGGGGGGCTGGGaaaatccggttg                                                                                                                                              >  1:625158/1‑141 (MQ=255)           aagCTGCACATACGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGTAAATGCCGTTAg                                                                                                                                              <  2:516205/150‑1 (MQ=255)               tGCACATACGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGAGTTTTAGGAGGGGCGGGGaaaagggcgt                                                                                                                                                >  2:560177/1‑137 (MQ=255)                   caTACGCCGGATCGTAAAATAGCGTGCTGCCAATTTGGTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGTCTGGTAAATGCCGTTAGCAGGCTGag                                                                                                                                     <  2:353695/151‑1 (MQ=255)                   caTACGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTTAGGAGTGGCTGGGAAATGACGGTAGCAGGGTg                                                                                                                                       >  2:352249/1‑149 (MQ=255)                      aCGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGagca                                                                                                                                   <  1:192743/150‑1 (MQ=255)                         ccGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAgcgc                                                                                                                               <  2:746370/151‑1 (MQ=255)                            gATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGGAAATGCCGTTTGCAGGGTTGGGAGCGcc                                                                                                                              >  1:420362/1‑149 (MQ=255)                             aTCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGGGATTTTTAGGAGGGGGCGGTAAATAGCGGTAGGAGGGTGAGGAgggc                                                                                                                               >  2:322434/1‑147 (MQ=255)                             aTCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGGTTTTTAGGAGGGGGTGGGAAATGGCGTTTGCAGGCCGAGGAGCGcc                                                                                                                              >  1:770947/1‑148 (MQ=255)                             aTCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGGTGGGAAAAGCCGTTAGCAGGCTGAGGAGCGcc                                                                                                                              >  1:388321/1‑148 (MQ=255)                                         cGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTTAGGAGTGGGTGGGGAAAGCCGGTTGCAAGGGGGGCAGCGCCCGG‑GAAGGTTTCAt                                                                                                               >  2:370572/1‑150 (MQ=255)                                                  aaTTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGGTGGGTAATGCCGGTAGCAGGGTGAGCACCGCCCGT‑GAAGCTTTTCTCTTGTTCG‑t                                                                                                     >  2:135303/1‑150 (MQ=255)                                                       gCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGGTGGGAAAAGCCGTTAGCAGGGTGAGCAGCGCCAGT‑GAAGATTTTATCTGGGTTG‑TCCGTg                                                                                                >  1:423255/1‑150 (MQ=255)                                                               gCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGT‑GAAGCTTTCATCTGGTTCG‑TCCGTG‑AAAATAaga                                                                                      <  2:379546/151‑1 (MQ=255)                                                                cACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTTATGAGTGGGTGGGGAATTGCGTTAACAGGGTGAGGAGGGCCGGT‑GAAGCGTTTTTTTTGTTTG‑TTCGTGTAAAATa                                                                                         >  1:242512/1‑148 (MQ=255)                                                                cACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGGTGGGAAAATGCGGTAGCAGGCTGGGCAGGGCCAGT‑GAAAGTTTCCACTTGTTTGTTCCTTG‑aaaaaa                                                                                         >  1:632993/1‑146 (MQ=255)                                                                  ccGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGT‑GAAGCTTTCATCTGGTTCG‑TCCGTG‑AAAATAAGAGt                                                                                    <  2:433146/150‑1 (MQ=255)                                                                    gTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTTAGGGGTGGGTGGGTAAAGCCGTTAACAGGGGGAGGCGGGCTCGTGGAGGCTTTC‑TCTGGGGTG‑TCCGGG‑AAAATTagagcacc                                                                                 >  2:143416/1‑147 (MQ=255)                                                                                  ggTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGGTGGGAAATTCCGGTAGCAGGGTGGGGAGCGGCCGT‑GAAGCTTTCATTTGGTTCG‑TCCGGG‑GAAATAAGAGGCACAGATTTGACCtt                                                                     >  2:469654/1‑149 (MQ=255)                                                                                     tGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGGGATTTTAAGGAGGGGCTGGGAAAATGCGGTAGCAGGGCGAGGAGCGCCCGG‑GGAGGTTTTATCTTGTTTG‑TCCGGG‑GAAAAAAAAGGCACAGAATCGACCCTcccg                                                                 >  1:38128/1‑148 (MQ=255)                                                                                         aTTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGT‑GAAGCTTTCATCTGGTTCG‑TCCGTG‑AAAATAAGAGTCACAGATTCGACCTTCCGGGCTaa                                                            <  1:562692/151‑1 (MQ=255)                                                                                                 ggATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTTAGGAGTGGCTGGGAAATGCCGTTTGCAGGCTGAGCAGCGCCAGT‑GAAGCTTTTATCTGGTTCG‑TCCGTG‑GAAATAAGAGGCACAGATTCGACCTTCCGGGCTAACGATTGtt                                                    >  1:563763/1‑151 (MQ=255)                                                                                                    tttGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGGGGCTGGGAAATGGCGTTAGCAGGCTGAGCAGCGCCCGG‑GAAGCTTTCATCTGGGTCG‑TCCGTG‑GAAAAAAAAGGCAAAGATTCGACCTTCCGGGCTAACGATTGtt                                                    >  1:255546/1‑148 (MQ=255)                                                                                                      tGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGT‑GAAGCTTTCATCTGGTTCG‑TCCGTG‑AAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTg                                                 <  2:219234/149‑1 (MQ=255)                                                                                                          ccgcGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGT‑GAAGCTTTCATCTGGTTCG‑TCCGTG‑AAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTagag                                            <  2:725020/147‑1 (MQ=255)                                                                                                           cACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGGAAATGCCGTTAGCAGGCTGAGCAGCGCCCGT‑GAAGCTTTTATCTGGTTTG‑TCCGTG‑GAAATAAGAGTCACAGATTTGACCTTCCGGGCTAACGATTGTTCTGTagag                                            >  1:387317/1‑149 (MQ=255)                                                                                                            aCGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGTCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGT‑GAAGCTTTCATCTGGTTCG‑TCCGTG‑AAAATAAGAGACACAGATTCGACCTTCCTGGCTAACGTTTGTTCTGTaga                                             <  2:722483/147‑1 (MQ=255)                                                                                                              ggTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGT‑GAAGCTTTCATCTGGTTCG‑TCCGTG‑AAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAAt                                       <  1:565773/151‑1 (MQ=255)                                                                                                                   gtggAACGGCGGGAGATTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGT‑GAAGCTTTCATCTGGTTCG‑TCCGTG‑AAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCa                                  <  2:709281/151‑1 (MQ=255)                                                                                                                          ggcggGAGATTTTAAGGAGTGGCTGGTAAAAGCCCGTAGCAGGCTGAGCAGCGCCAGT‑GAAGCTTTCATCTGGTTCG‑TCCGTG‑AAAATAAGAGTCACAGATTCGACCTTTCGGGCTAACGATTGTTTTGTAGAGAAAATGTCCAGAAGCgg                           >  2:535427/1‑151 (MQ=255)                                                                                                                           gcggGAGATTTTAAGGAGTGGCTGGGAAATGCCGTTAGCAGGGTGAGCAGCGCCAGT‑GAAGCTTTCCTCTGGTTCG‑TCCGTG‑AAAATAAGAGGCCCAGAATCGACCTTCCGGGGTAACGATTGTTCTGTAGAGAAAATGGCCAGAAGCgg                           >  2:180370/1‑150 (MQ=255)                                                                                                                             gggAGATTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGT‑GAAGCTTTCATCTGGTTCG‑TCCGTG‑AAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGAAGCGGATa                        <  1:240987/151‑1 (MQ=255)                                                                                                                                 gaTTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGT‑GAAGCTTTCATCTGGTTCG‑TCCGTG‑AAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGAAGCGGATaaaaa                    <  2:38123/151‑1 (MQ=255)                                                                                                                                    tttAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGT‑GAAGCTTTCATCTGGTTCG‑TCCGTG‑AAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGAAGCGGATAAAAAAt                  <  1:570875/150‑1 (MQ=255)                                                                                                                                                 tGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGT‑GAAGCTTTCATCTGGTTCG‑TCCGTG‑AAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGAAGCGGATAAAAAATTCCTGAACAACGTg    >  1:572251/1‑151 (MQ=255)                                                                                                                                                 tGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGT‑GAAGCTTTCATCTGGTTCG‑TCCGTG‑AAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGAAGCGGATAAAAAATTCCTGAACAACGTg    >  1:436006/1‑151 (MQ=255)                                                                                                                                                    tAAATGCCGTTAGCAGGCTGAGCAGCGCCAGT‑GAAGCTTTCATCTGGTTCG‑TCCGTG‑AAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGCTTGTTCTGTAGAGAAAATGTCCAGAAGCGGATAAAAAATTCCTGAACAACGTggg  >  2:193246/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                         CAGGATAAGTAAGCTGCACATACGCCGGATCGTAAAATAGCGTGCTGCCAATTTGCTGTCTGGCACCGTCTGCGATGGCAAGGTTGGTATTGGCCTGGATTTGTACCACGGTTGGTGGAACGGCGGGAGATTTTAAGGAGTGGCTGGTAAATGCCGTTAGCAGGCTGAGCAGCGCCAGT‑GAAGCTTTCATCTGGTTCG‑TCCGTG‑AAAATAAGAGTCACAGATTCGACCTTCCGGGCTAACGATTGTTCTGTAGAGAAAATGTCCAGAAGCGGATAAAAAATTCCTGAACAACGTGGG  >  NC_000913/4137468‑4137764

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