Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,238,049	0	A	G	44.7%	46.6 / 10.7	38	E139G (GAA→GGA)	yjbH	DUF940 family extracellular polysaccharide protein
Reads supporting (aligned to +/- strand):  ref base A (5/16);  new base G (16/1);  total (21/17)
Fisher's exact test for biased strand distribution p-value = 1.32e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAGCAGCGTCGAAGCGTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGTACGGGGCTGTTTGATGCGGAATA‑TCTTGTTGCCAGC‑AAAG‑CC‑TGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCT‑GGGGGTA‑TTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGG  >  NC_000913/4237902‑4238199                                                                                                                                                     |                                                                                                                                                             aCAGCAGCGTCGAAGCGTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGTACGGGGCTGTTTGATGCGGAat                                                                                                                                                            <  1:511427/150‑1 (MQ=255)      agcGTCGAAGCGTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGTACGGGGCTGTTTGATGCGGAATA‑TCttg                                                                                                                                                     <  2:294681/151‑1 (MQ=255)            gAAGCGTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGTACGGGGCTGTTTGATGCGGAATA‑TCTTGTTGCCa                                                                                                                                               <  2:425924/151‑1 (MQ=255)                             tCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGGTATCGGC‑GGGTCGGGGGTGGTTTATGCGGGATA‑ATTTGTTTCCCGC‑AAAA‑CC‑Cggggggc                                                                                                                            >  2:481815/1‑148 (MQ=255)                             tCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGGAACGGGGTGTTTTATGCGGGAAA‑TCTTTTTGCCAGC‑AAAA‑AC‑Cggggg                                                                                                                              >  1:372754/1‑148 (MQ=255)                                  cGTATAAAGGTAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGTACGGGGCTGTTTGATGCGGAATA‑TCTTGTTGCCAGC‑AAAG‑CC‑TGGGGGCCGTTCGa                                                                                                                      <  2:353684/151‑1 (MQ=255)                                    tataAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGGAAGGGGGTGGTTTATGCGGGAAA‑TCTTTTTTCCCGCAAAAC‑CC‑TGGGGGCGGttgatt                                                                                                                     >  1:404773/1‑147 (MQ=255)                                        aagataaagCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTTGGGCGCGGGATATCGGCGGGAAGGGGGCGTTTTTATGGGGAA‑A‑ACTTGGTTTGAGC‑AAAA‑CCGTGGGGGCCGTTGGtttt                                                                                                                   >  1:139563/1‑149 (MQ=255)                                         agataaagCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGTACGGGGGTGTTTTTTGCGGGAAA‑TTTTGTTGGCAGC‑AAAA‑GC‑TGGGGGCCGGTTGAtttt                                                                                                                  >  1:254376/1‑148 (MQ=255)                                               agCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGTACCGGGGTGTTTGGAGGGGGAAA‑TATTTTTTGCCGG‑AAAG‑GC‑TGGGGGCCGGTCGCtttttaactcgg                                                                                                          >  1:135492/1‑142 (MQ=255)                                                   ttCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGGAAGGGGGTGTTTTTTGCGGGAAA‑TATTGTTGTCAGC‑AAAA‑CC‑CGGGGGGCGGTTGATTTTTCGCTCGGcc                                                                                                        >  1:666582/1‑148 (MQ=255)                                                    tCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGTACGGGGCTGTTTGATGCGGAATA‑TCTTGTTGCCAGC‑AAAG‑CC‑TGGGGGCCGTTCGATTTTACGCTCGGCctg                                                                                                      <  1:586819/149‑1 (MQ=255)                                                    tCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGGACGGGGGTGTTTTGTGGGGGAAA‑TCTTTTTTGCCGC‑AAAA‑CC‑TGGGGGCCGGTTGGTTTTTCGCCCGGccctgg                                                                                                    >  1:208037/1‑147 (MQ=255)                                                    tCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGGACGGGGGCGTTTTATTCGGGATA‑TTTTTTTGTCCGC‑AAAG‑CC‑CGGGGGGCGGTTGGTTTTTCGCTCGGCctggg                                                                                                    >  2:317316/1‑151 (MQ=255)                                                     cGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGGACGGGGGTGTTTTGTGGGGGAAA‑TCTTTTTTGCCGC‑AAAA‑CC‑CGGGGGGCGGTTGATTTTTAGCCCGGCccggg                                                                                                    >  1:651437/1‑150 (MQ=255)                                                         ctcAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGGAAGGGGGTGTTTTATGCGGGAAA‑TCTTTTTGCCAGC‑CAAA‑AC‑TGGGGGGCGGTTGGTTTTTCGCCCCGCccggggcggg                                                                                               >  1:328761/1‑146 (MQ=255)                                                           cAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGTACGGGGCTGTTTGATGCGGAATA‑TCTTGTTGCCAGC‑AAAG‑CC‑TGGGGGCCGTTCGATTTTACGCTCGGCctgggct‑ggg                                                                                              <  1:370590/149‑1 (MQ=255)                                                                  cGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGGAAGGGGGGGTTTTGTGCGGGAAA‑TTTTGTTGGCCGC‑CAAA‑AC‑TGGGGGGCCGTCGATTTTTCGCCCCGGCCGGGGT‑gggggggttt                                                                                       >  1:631458/1‑144 (MQ=255)                                                                  cGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCAGC‑GGTACGGGGCTGTTTGATGCGGAATA‑TCTTGTTGCCAGC‑AAAG‑CC‑TGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCT‑GGGGGTA‑TTTgg                                                                                    <  1:265036/151‑1 (MQ=255)                                                                     ttGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTTGTGCGCGGGATATCGGC‑GGGACGGGGGCGTTTTATGCGGGATA‑ACTTTTTGCCCGC‑AAAA‑CC‑TGGGGGGCCGTCCCTTTTTCCCCCGGCCTCGGCC‑gggggggtt                                                                                        >  2:509321/1‑141 (MQ=255)                                                                       gtgGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGGACGGGGGCGTTTTATGGGGGAAA‑TCTTGTTGGCCGC‑AAAGACC‑TGGGGG‑CGTTCGGTTTTTCGCCCCGCCCGGGGT‑GGGGGGA‑TTTGGGGAcc                                                                               >  1:463272/1‑151 (MQ=255)                                                                                gagTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGTAGGGGGTCGTTTGATGGGGGATA‑TTTTCTTTCCCCC‑CAAA‑CC‑CCGGGGGCGCTCTATTTTTATCTCCGCCTCGGCTCGGGGGGA‑ttttgcaccaaccgt                                                                          >  2:538285/1‑136 (MQ=255)                                                                                    tACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGTACGGGGGTGTTTTATGGGGAAAA‑TCTTTTTGCCAGC‑CAAA‑CC‑TGGGGGGCGGTCGGTTTTAACCTCGGCCTGGGGT‑GGGGGGATTTTGGGCA‑CAGCGGGAAtgt                                                                    >  1:406359/1‑149 (MQ=255)                                                                                      cTGGCTGCCGCAAGTGGCGGTTGGTGGGCGGGGTATCGGC‑GGGTCGGGGGTGGTTGGTGGGGGAGATTTTTTTTGCAAGA‑AAAA‑CC‑GGGGGGCGGTTGGTTTTTAGCCTCGGCC‑GGGGT‑GGGGGTG‑TTTTGGGGCCAGCGGGGATGTGaa                                                                 >  2:247042/1‑150 (MQ=255)                                                                                      cTGGCTGCCGCAAGTGGCGGTTGGGGGGCGGGATATCGGC‑GGGAAGGGGGTGTTTGGTGGGGGAAA‑TTTTGTTGCCAGC‑AAAA‑GC‑TGGGGGGCGGTTGATTTTTAGCTCGGGCCGGGGT‑GGGGGGG‑TTTTGGGACCAGCGGGAAagtg                                                                   >  1:316450/1‑148 (MQ=255)                                                                                        ggCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGTACGGGGCTGTTTGATGCGGAATA‑TCTTGTTGCCAGC‑AAAG‑CC‑TGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCT‑GGGGGTA‑TTTGGGCACCAGCGGTAATGTGaaa                                                                <  2:676467/149‑1 (MQ=255)                                                                                            gCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGTACGGGGCTGTTTGATGCGGAATA‑TCTTGTTGCCAGC‑AAAG‑CC‑TGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCT‑GGGGGTA‑TTTGGGCACCAGCGGTAATGTGAAAAATcc                                                           <  1:563342/150‑1 (MQ=255)                                                                                                   gTGGCGGTTGGTGCGCGGGATATCGGC‑GGTACGGGGCTGTTTGATGCGGAATA‑TCTTGTTGCCAGC‑AAAG‑CC‑TGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCT‑GGGGGTA‑TTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGtt                                                   <  1:457/151‑1 (MQ=255)                                                                                                     ggCGGTTGGTGCGCGGGATATCGGC‑GGTACGGGGCTGTTTGATGCGGAATA‑TCTTGTTGCCAGC‑AAAG‑CC‑TGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCT‑GGGGGTA‑TTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTc                                                  <  1:555160/150‑1 (MQ=255)                                                                                                         gTTGGTGCGCGGGATATCGGC‑GGTACGGGGCTGTTTGATGCGGAATA‑TCTTGTTGCCAGC‑AAAG‑CC‑TGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCT‑GGGGGTA‑TTTGGGCACCAGCGGTAATGTGAAAAATCCGCTATGTTCAGCCa                                             <  2:65233/151‑1 (MQ=255)                                                                                                           tGGTGCGCGGGATATCGGC‑GGTACGGGGCTGTTTGATGCGGAATA‑TCTTGTTGCCAGC‑AAAG‑CC‑TGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCT‑GGGGGTA‑TTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAg                                            <  1:391558/150‑1 (MQ=255)                                                                                                                 gcgGGATATCGGC‑GGTACGGGGCTGTTTGATGCGGAATA‑TCTTGTTGCCAGC‑AAAG‑CC‑TGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCT‑GGGGGTA‑TTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATaaa                                     <  1:583826/151‑1 (MQ=255)                                                                                                                    ggATATCGGC‑GGTACGGGGCTGGTTGATGCGGAATA‑TCTTGTTGCCAGC‑AAAG‑CC‑TGGGGGCCGTTCGATTTTACGCCCGGCCTGGGCT‑GGGGGGA‑TTTGGGCACCAGCGGGAATGGGAAAAATCCGCTCTGTTCAGCCAGGGATaaaaa                                   >  2:342629/1‑148 (MQ=255)                                                                                                                          cggc‑ggTACGGGGCTGTTTGATGCGGAATA‑TCTTGTTGCCAGC‑AAAG‑CC‑TGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCT‑GGGGGTA‑TTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTa                             <  1:654764/150‑1 (MQ=255)                                                                                                                                gTACGGGGCTGTTTGATGCGGAATA‑TCTTGTTGCCAGC‑AAAG‑CC‑TGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCT‑GGGGGTA‑TTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATcgcg                        >  2:144299/1‑150 (MQ=255)                                                                                                                                gTACGGGGCTGTTTGATGCGGAATA‑TCTTGTTGCCAGC‑AAAG‑CC‑TGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCT‑GGGGGTA‑TTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATcgcg                        >  1:666303/1‑150 (MQ=255)                                                                                                                                          gTTTGATGCGGAATA‑TCTTGTTGCCAGC‑AAAG‑CC‑TGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCT‑GGGGGTA‑TTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCt               <  2:431408/149‑1 (MQ=255)                                                                                                                                            ttGATGCGGGATA‑TCTTGTTGCCAGC‑AAAG‑CC‑TGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCT‑GGGGGTA‑TTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTa              <  1:571977/148‑1 (MQ=255)                                                                                                                                                     aaTA‑TCTTGTTGCCAGC‑AAAG‑CC‑TGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCT‑GGGGGTA‑TTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCggg  >  2:100701/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                             ACAGCAGCGTCGAAGCGTTCTCTGGCGATCAAACGTATAAAGATAAAGCCTTCGATCTCAAACTGCGTTTGTGGGAAGAGAGTTACTGGCTGCCGCAAGTGGCGGTTGGTGCGCGGGATATCGGC‑GGTACGGGGCTGTTTGATGCGGAATA‑TCTTGTTGCCAGC‑AAAG‑CC‑TGGGGGCCGTTCGATTTTACGCTCGGCCTGGGCT‑GGGGGTA‑TTTGGGCACCAGCGGTAATGTGAAAAATCCGCTCTGTTCAGCCAGTGATAAATATTGCTATCGCGATAACAGCTACAAACAGGCGGG  >  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 < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑