Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,238,049	0	A	G	46.9%	62.6 / 16.4	49	E139G (GAA→GGA)	yjbH	DUF940 family extracellular polysaccharide protein
Reads supporting (aligned to +/- strand):  ref base A (3/23);  new base G (23/0);  total (26/23)
Fisher's exact test for biased strand distribution p-value = 5.48e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAGCAGCGTCGAAGCGTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTGTTTGATGCGGAATATCTTGTT‑G‑CCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGT‑GAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGG  >  NC_000913/4237902‑4238199                                                                                                                                                     |                                                                                                                                                          aCAGCAGCGTCGAAGCGTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTGTTTGATGCGGAat                                                                                                                                                         <  1:235187/150‑1 (MQ=255)   agcagcGTCGAAGCGTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGGGCGCGGGATATCGGCGGGA‑CGGGGCGGTTTTGTGCGGGatat                                                                                                                                                       >  1:105789/1‑150 (MQ=255)           cGAAGCGTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTGTTTGATGCGGAATATCTTGTT‑G‑cc                                                                                                                                            <  2:595465/151‑1 (MQ=255)                cGTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGGTATCGGCGGGA‑CGGGGGCGTTTTTTGGGGGATATtttttttgc                                                                                                                                              >  1:152937/1‑143 (MQ=255)                 gTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTGTTTGATGCGGAATATCTTGTT‑G‑CCAGCaaa                                                                                                                                      <  1:150150/151‑1 (MQ=255)                  ttCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGA‑CGGGGGTGGTTTGTGCGGGATATTTTGTT‑G‑CCCGcca                                                                                                                                       >  2:693935/1‑147 (MQ=255)                   tctcTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTGTTTGATGCGGAATATCTTGTT‑G‑CCAGCAAAg                                                                                                                                     <  1:127446/150‑1 (MQ=255)                       tGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGA‑CGGGGGTGTTTTATGGGGGATATCTTGTT‑G‑CCCGCAAAACCTgg                                                                                                                                >  2:858113/1‑151 (MQ=255)                         gCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGGTGTTTGGTGGGGGAAATCTTGGT‑G‑CCCGCAAAGCCTgggg                                                                                                                              >  2:625023/1‑151 (MQ=255)                          cGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGA‑CGGGGGTGGTTTGTGCGGGATATTTTGTT‑G‑GCCGCCAAACCCggggg                                                                                                                             >  1:855513/1‑151 (MQ=255)                               aaaCGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGA‑CGGGGGTGGTTTATGGGGGAAATCTTGTT‑G‑GCCGCAAAACCCGGGGGGCggt                                                                                                                        >  2:35469/1‑149 (MQ=255)                               aaaCGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGA‑CGGGGCTGTTTTATGCGGGATATCTTGTT‑T‑CCCGCAAAACCCGGGGGcc                                                                                                                           >  2:218295/1‑148 (MQ=255)                                      taaagataaagCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGGGGTA‑CGGGGGTGTTTTATGGGGGATATCTTGTT‑G‑GCTGCAAAACCTGGGGGGCGGTTGGtttt                                                                                                                 >  1:29030/1‑151 (MQ=255)                                         agataaagCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTGTTTGATGCGGAATATCTTGTT‑G‑CCAGCAAAGCCTGGGGGCCGTTCGATTTTACg                                                                                                              <  1:539564/151‑1 (MQ=255)                                         agataaagCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGA‑CGGGGCTGTTTGTTGCGGAATATCTTGTT‑G‑GCAGCAAAGCCTTGGGGGCGGTCGGtttt                                                                                                                 >  1:644834/1‑148 (MQ=255)                                           ataaagCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTGTTTGATGCGGAATATCTTGTT‑G‑CCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCt                                                                                                            <  1:480613/151‑1 (MQ=255)                                           ataaagCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGA‑CGGGGGTGTTTTGTGCGGGATATCTTGTT‑T‑CCCGCAAAGCCTGGGGGCCGGTCGATTTTTCg                                                                                                              >  2:182934/1‑149 (MQ=255)                                            taaagCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGGGCGCGGGGTATCGGCGGGA‑AGGGGGTGTTTGATGCGGGAAATTTTGGT‑GCCCCGCAAACCCGGGGGGcgctcggtttt                                                                                                                  >  1:920841/1‑135 (MQ=255)                                             aaagCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTGTTTGATGCGGAATATCTTGTT‑G‑CCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTc                                                                                                           <  1:419647/150‑1 (MQ=255)                                                 ccTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGGTGTTTTATGCGGGATATCTTTTT‑T‑CCAGCCAAACCTGGGGGCCGGTCGGTTTTACGCCCGGcc                                                                                                       >  2:152629/1‑150 (MQ=255)                                                 ccTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTGTTTGATGCGGAATATCTTGTT‑G‑CCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCct                                                                                                      <  1:754255/151‑1 (MQ=255)                                                  cTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGGTATCGGCGGGA‑CCGGGGCGTTTTATGCGGGATATTTTGTT‑T‑TCCGCAAAGCCTGGGGGGCGTTTGAttttt                                                                                                                >  2:661793/1‑139 (MQ=255)                                                     cGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTGTTTGATGCGGAATATCTTGTT‑G‑CCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCctgg                                                                                                    <  1:371262/149‑1 (MQ=255)                                                        tctcAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGGTGTTTGGTGCGGGATATCTTTTT‑G‑GCAGCAAAGCCTGGGGGGCGGTCGATTTTTCCCTCGGccctgggtgg                                                                                               >  2:218039/1‑143 (MQ=255)                                                           cAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTGTTTGATGCGGAATATCTTGTT‑G‑CCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCctgggct                                                                                                 <  1:800147/146‑1 (MQ=255)                                                            aaaCTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTGTTTGATGCGGAATATCTTGTT‑G‑CCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGt                                                                                           <  1:402856/151‑1 (MQ=255)                                                                tGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGCTGGTGGGCGGGACATCGGCGGGG‑CGGGGCGGGTTGGTGGGGGATTATTTTTT‑G‑CCAGCCAAACCCGGGGGGCGGTCGCTTTTTACCGCGGGCCGGGGTggggggagtt                                                                                       >  1:53214/1‑146 (MQ=255)                                                                 gCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTGTTTGATGCGGAATATCTTGTT‑G‑CCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTAtt                                                                                        <  1:582579/149‑1 (MQ=255)                                                                 gCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGA‑CGGGGGTGTTTTTTGCGGGATATCTTGTT‑G‑CCCGCAAAACCTGGGGGGCGGTCGATTTTTAGCTCGGCCCGGGGTGGGGGGAttt                                                                                       >  2:803945/1‑150 (MQ=255)                                                                  cGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTGTTTGATGCGGAATATCTTGTT‑G‑CCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTgg                                                                                     <  2:152429/151‑1 (MQ=255)                                                                              gagagTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGA‑CGGGGGTGTTTTATGCGGGATATCTTTTT‑G‑CCCGCAAAACCTGGGGGGCGGTCGGTTTTTCGCTCGGCCTTGGCTGGGGGGATTTTGGGACCAGCgg                                                                           >  2:451190/1‑149 (MQ=255)                                                                                 agTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTGTTTGATGCGGAATATCTTGTT‑G‑CCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAAtg                                                                      <  1:262746/151‑1 (MQ=255)                                                                                 agTTACTGGCTGCCGCAAGTGGCGGTTGGGGCGCGGGATATCGGCGGGA‑CGGGGGCGTTTTATGCGGGATATTTTTTT‑G‑GCCGCCAAACCCGGGGGGCGGTCGATTTTTAGCTCGGCCCGGGGTGGGGGGGTTTTGGGACCCGCGGGAAt                                                                       >  1:731715/1‑150 (MQ=255)                                                                                     aCTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTGTTTGATGCGGAATATCTTGTT‑G‑CCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGT‑Ga                                                                  <  2:807939/150‑1 (MQ=255)                                                                                      cTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGA‑CGGGGCTGTTTGATGCGGGATATTTTGTT‑G‑CCAGCAAAACCCGGGGGCCGGTCGATTTTTAGCCCGGCCTGGGCTGGGGGGATTTTGGGCCCAGCGGGGATGT‑Gaa                                                                 >  2:707604/1‑150 (MQ=255)                                                                                           tGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGA‑AGGGGGTGTTTGGTGCGGGATATCTTTTT‑T‑CCCGCCAAAGCTGGGGGGCGTTTGGTTTTTCGCCCGGGCTTGGGTGGGGGGATTTTGGGCCCAGCGGGAATGT‑GAAAAAAcc                                                           >  2:440221/1‑151 (MQ=255)                                                                                            gCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGA‑AGGGGGTGTTTTATGCGGGATATCTTTTT‑G‑GCCGCAAAACCCGGGGGCCGTTCGATTTTTAGCTCGGGCCGGGCTGGGGGGATTTGGGGACCCACGGGAAAGT‑GAAAAAAcc                                                           >  2:391449/1‑150 (MQ=255)                                                                                               gCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTGTTTGATGCGGAATATCTTGTT‑G‑CCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGT‑GAAAAATCCGct                                                        <  1:576666/150‑1 (MQ=255)                                                                                                cAAGTGGCGGTTGGTGCGCGGGATATCGGCGGGACCGGGGCT‑TTTGGTGCGGGAAATCTTGTTGG‑CCGGCAAA‑CCTTGGGGGCGGTCGATTTTTCCCCCCGCCCGGGGTGGGGGGATTTGGGGACCCGCGGGGAAGTGGAAAAATccccc                                                        >  1:101260/1‑147 (MQ=255)                                                                                                 aaGTGGCGGTTGGTGCGCGGGATATCGGCGGGA‑CGGGGGTGTTTGATGCGGGATATCTTGTT‑G‑CCAGCAAAACCTGGGGGGCGGTCGATTTTTAGCCCGGGCCGGGGTGGGGGGATTTTGGGACCAGCGGGAATGG‑GGAAAAttc                                                           >  1:568194/1‑143 (MQ=255)                                                                                                      gCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTGTTTGATGCGGAATATCTTGTT‑G‑CCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGT‑GAAAAATCCGCTCTGTTCAg                                                <  2:849179/151‑1 (MQ=255)                                                                                                        ggTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTGTTTGATGCGGAATATCTTGTT‑G‑CCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGT‑GAAAAATCCGCTCTGTTCAGcc                                              <  1:566464/151‑1 (MQ=255)                                                                                                           tGGTGCGCGGGATATCGGCGGTA‑CGGGGCTGTTTGATGCGGAATATCTTGTT‑G‑CCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGT‑GAAAAATCCGCTCTGTTCAGCCAg                                            >  2:817734/1‑150 (MQ=255)                                                                                                               gcgcGGGATATCGGCGGTA‑CGGGGCTGTTTGATGCGGAATATCTTGTT‑G‑CCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGT‑GAAAAATCCGCTCTGTTCAGCCAGTGAt                                        <  1:734363/150‑1 (MQ=255)                                                                                                                         tCGGCGGTA‑CGGGGCTGTTTGATGCGGAATATCTTGTT‑G‑CCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGT‑GAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCt                              >  1:555110/1‑150 (MQ=255)                                                                                                                                          gTTTGATGCGGAATATCTTGTT‑G‑CCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGT‑GAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTa              <  1:500960/150‑1 (MQ=255)                                                                                                                                             tGATGCGGAATATCTTGTT‑G‑CCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGT‑GAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACaa           <  2:918648/150‑1 (MQ=255)                                                                                                                                                  cGGAATATCTTGTT‑G‑CCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGT‑GAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGc     <  1:415257/151‑1 (MQ=255)                                                                                                                                                    gAATATCTTGTT‑G‑CCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGT‑GAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCgg   <  1:182917/151‑1 (MQ=255)                                                                                                                                                     aaTATCTTGTT‑G‑CCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCAGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGT‑GAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCggg  >  2:142007/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                          ACAGCAGCGTCGAAGCGTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGCGGTA‑CGGGGCTGTTTGATGCGGAATATCTTGTT‑G‑CCAGCAAAGCCTGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCTGGGGGTATTTGGGCACCAGCGGTAATGT‑GAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGG  >  NC_000913/4237902‑4238199

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