Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,535,789	0	C	A	42.5%	53.3 / 17.3	47	V24F (GTT→TTT)	yddE	PhzC‑PhzF family protein
Reads supporting (aligned to +/- strand):  ref base C (7/20);  new base A (20/0);  total (27/20)
Fisher's exact test for biased strand distribution p-value = 9.79e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑AAAAC‑AACACCGG‑CAG‑AATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAAT  >  NC_000913/1535640‑1535937                                                                                                                                                       |                                                                                                                                                         aCCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑AAaac‑a                                                                                                                                                        <  2:620654/151‑1 (MQ=255)  ccGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑AAaac‑aa                                                                                                                                                       <  1:155323/151‑1 (MQ=255)    gCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑AAAac‑aaca                                                                                                                                                     <  1:454230/151‑1 (MQ=255)      aaaTCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑Gaaaa‑aaaaac                                                                                                                                                   >  1:809760/1‑147 (MQ=255)          cGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑AAAAC‑AACACCGG‑c                                                                                                                                               <  2:895467/150‑1 (MQ=255)              aCTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑AAAAC‑AACACCGG‑CAG‑AAt                                                                                                                                         <  2:566605/151‑1 (MQ=255)              aCTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑AAAAC‑AACACCGG‑CAG‑AAt                                                                                                                                         <  1:596203/151‑1 (MQ=255)                  cAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑GAAAA‑AAAACCCG‑CAG‑aaat                                                                                                                                        >  1:11504/1‑146 (MQ=255)                       gTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑AAAAC‑AAAACCCG‑GAG‑aaat                                                                                                                                        >  1:594196/1‑141 (MQ=255)                          ggCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑AAAAC‑AACACCGG‑CAG‑AATTGCCGC‑GAAAc                                                                                                                             <  1:709498/150‑1 (MQ=255)                           gCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑GAAAA‑AAAACCCG‑CGG‑AATTGCCGCGGAAACg                                                                                                                            >  2:906390/1‑151 (MQ=255)                               aaaGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑GAAAA‑AAAACCGGCCGG‑AATTCCCCG‑GAAACGGt                                                                                                                          >  2:607522/1‑149 (MQ=255)                                aaGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑AAAAC‑AACACCGG‑CAG‑AATTGCCGC‑GAAACGGTtgtg                                                                                                                      <  2:437321/151‑1 (MQ=255)                                  gTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑AAAAC‑AACACCGG‑CAG‑AATTGCCGC‑GAAACGGTTGTGAg                                                                                                                    <  1:525896/151‑1 (MQ=255)                                       ggATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑GAAAA‑AAAACCCG‑CAG‑AATTGCCGC‑GGAACAGGTGTGGGGGaa                                                                                                                >  2:355377/1‑150 (MQ=255)                                         aTGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑GAAAA‑AACACCCG‑CAG‑AAATGCCGC‑GGAAAGGGTGTGGGGGaaaa                                                                                                              >  1:622226/1‑150 (MQ=255)                                          tGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑AAAAA‑AAAACCCG‑CAG‑AAATTCCGC‑GGAACAGGTGTggggggaaaa                                                                                                             >  2:637138/1‑144 (MQ=255)                                              cacaTCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑GAAAA‑AAAACCCC‑GAG‑AATTGTCGC‑GGAAAgggggtggggggaaa                                                                                                              >  2:552469/1‑132 (MQ=255)                                                 aTCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑AAAAC‑AACACCGG‑CAG‑AATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATCGACg                                                                                                     <  1:490379/151‑1 (MQ=255)                                                  tCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑AAAAC‑AACACCGG‑CAG‑AAATTCCGC‑GCAACGGTTGTGGGGGaaaaaaatc                                                                                                         >  1:705420/1‑141 (MQ=255)                                                         cgtcgCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑GAGAC‑AACACCCG‑CAG‑AAATGCCGC‑CCAACAGTAGTTAGGGGAAAACACCCACCGGGTaca                                                                                              >  2:65462/1‑150 (MQ=255)                                                         cgtcgCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑AAAAC‑AACACCGG‑CAG‑AATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATCGACGTGGTaca                                                                                              <  1:374036/150‑1 (MQ=255)                                                               tgtgCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑GAAAA‑AAAACCCG‑CAG‑AAATGCCGC‑GGAACGGGTGTGGGGGAAAAAAATCGACGGGGGACAcccccgg                                                                                       >  1:463475/1‑146 (MQ=255)                                                               tgtgCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑AAAAC‑AACACCAG‑CAG‑AAATGCCGC‑GAAAAGGGTGTGGGGGAAAAACATCGACGGGGTACACCcgcg                                                                                        >  1:561207/1‑150 (MQ=255)                                                                      cagAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑AAAAC‑AACACCGG‑CAG‑AATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCAtt                                                                                <  2:675070/151‑1 (MQ=255)                                                                                gTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑AAAAC‑AACACCGG‑CAG‑AATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTcc                                                                      <  1:210209/151‑1 (MQ=255)                                                                                  ttCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑GAAAA‑AAAACCCG‑CAGCAATTTCC‑C‑GCAACGGGTGGGGGGGGAAAGCAACAACCCGGGACACCCCCGGGGTCATTTAAGTtttc                                                                       >  2:533459/1‑148 (MQ=255)                                                                                      gAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑GAAAA‑AACACCCG‑CAG‑GAATGCCGC‑GAAAAGGGTGTGGGGGAAAAAAAACGACGGGGGACAACTGCGGGTTTAtttaaagtctcc                                                                      >  1:368302/1‑135 (MQ=255)                                                                                      gAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑AAAAA‑AAAACAGG‑CCG‑AAATGCCGC‑GGAACGGGTGTGGGGGGAAAAAATACACGTGGGGGACCCGCGGGTTTATTaaaggttttccctttc                                                                >  2:882712/1‑137 (MQ=255)                                                                                       aaTGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑GAAAA‑AAAACCCG‑CAG‑AATTGCCGC‑GAAACGGGTGTGGGGGAAAAAAAACCACGGGGGACACCTGCGGGTTCATTAAAGTTCTCCtgtg                                                                  >  2:831542/1‑148 (MQ=255)                                                                                         tGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑GAAAA‑AAAACCCG‑CAG‑AATTGCCGC‑GAAAAGGGTGTGAGGGAAAAACAACCACGGGGGACACCCGCGGGTTCATTAAAATTCTcc                                                                      >  2:801708/1‑142 (MQ=255)                                                                                          ggCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑GAAAA‑AACACCCG‑CAG‑AAATACCCC‑GGAAAGGGGGTGGGGGAAAAAAAACCACGGGGGACACCTCCGGGTTTATTAACATTTTCCTTTGAc                                                                >  2:408125/1‑147 (MQ=255)                                                                                               aaCTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑AAAAC‑AACACCGG‑CAG‑AATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAAc                                                        <  2:711359/150‑1 (MQ=255)                                                                                                   cgcgGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑AAAAC‑AACACCGG‑CAG‑AATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCa                                                   <  2:306649/151‑1 (MQ=255)                                                                                                      gggCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑AAAAC‑AACACCGG‑CAG‑AATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAg                                                <  1:567906/151‑1 (MQ=255)                                                                                                        gCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑GAAAA‑AAAAACCG‑CAG‑GATTGGCGC‑CGAAAGGGTGTGGGGGAAAAGCACAGACGCGGGACACCTCCGGGGTTATTAAAGGTCTCCCGTGGCTTGAGAAATATCAAAGcc                                              >  2:541361/1‑151 (MQ=255)                                                                                                          gATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑AAAAC‑AACACCGG‑CAG‑AATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTg                                            <  2:707620/151‑1 (MQ=255)                                                                                                            tAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑GAAAA‑AACACCGG‑CAG‑AATTGCCGC‑GGAACGGGTGTGGGGGAAAAGCATCGACGTGGGACACCTGCGGGTTCATTAACGTTCTCCTGTGACTGGAGAAATATCATAGcc                                              >  1:514984/1‑147 (MQ=255)                                                                                                            tAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑GAAAA‑AAAACCGG‑CAGAAATTGCC‑C‑GAAAAGGGTGGGGGGTGAAAAAAACGAAGGGGGACAACCGCGGGTTTATTTAAGGTTTCCTCTGACTGGGGAACTaacctaacc                                              >  2:888737/1‑139 (MQ=255)                                                                                                            tAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑GAAAA‑AAAACCCG‑CAG‑AATTGGCGC‑GCAACGGGTGTGGGGGGAAAACAACGACGTGGGACACCTGCGGGTTTATTAACATTCTCCTTTGGCTGGGGAAATATCATaaac                                              >  2:272041/1‑144 (MQ=255)                                                                                                             aaGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑AAAAC‑AACACCGG‑CAG‑AATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCaa                                         <  2:421266/151‑1 (MQ=255)                                                                                                               gCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑AAAAC‑AACACCGG‑CAG‑AATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGt                                       <  1:305822/151‑1 (MQ=255)                                                                                                                  gCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑AAAAC‑AACACCGG‑CAG‑AATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTg                                      <  2:465310/149‑1 (MQ=255)                                                                                                                  gCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑AAAAC‑AACACCCG‑CCG‑AATTGCCCC‑GGAACGGGTGTGGGGGAAAAAACACGCCGTGGTACACCCGCGGGTTTATTTACGGTCTCCCCTGTCTGGTGGAACATCATAGCCCGCAAATgg                                     >  2:298622/1‑150 (MQ=255)                                                                                                                    aTTTGCGCTTCGCTAAGATTATCGGCGGGG‑AAAAC‑AACACCGG‑CAG‑AATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCg                                  <  1:597113/151‑1 (MQ=255)                                                                                                                      ttGCGCTTCGCTAAGATTATCGGCGGGG‑GAAACAAAAACCGG‑C‑G‑AAATGCCGC‑GAAAAGGGTGTGAGGGAAAAAAAACGACGTGGGACACCTCCGGGTTTATTAAACTTCTCCTGTGGCTGGTGAAATATCATAGACTGCAAGAgg                                     >  2:647427/1‑146 (MQ=255)                                                                                                                            ttCGCTAAGATTATCGGCGGGGAAAAAA‑AAAACCGG‑AGG‑AATTTCGCG‑GAAACGGTGGGGGGGAAAAAAAACGGAGGGGGAAAACC‑GCGGGTTCATTAAAGTTTTCCTCTGACTGGGGAAAAATAATATCCCGCAAGAGGcccgggggg                            >  2:673535/1‑142 (MQ=255)                                                                                                                            ttCGCTAAGATTATCGGCGGGG‑GAAAA‑AACACCCG‑CAG‑AAATGCCGC‑GGAACGGGTGTGGGGGAAAAACAACGACGTGGGACACCCGCGGGTTTATTTAAGTTTTCCCTTGACTGGGGAACTATCATAGCCTGCAAGTGGCCGGGGAGCg                           >  1:330812/1‑150 (MQ=255)                                                                                                                                            gcggGG‑AAAAC‑AACACCGG‑CAG‑AATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAg          <  1:847598/151‑1 (MQ=255)                                                                                                                                                gg‑AAAAC‑AACACCGG‑CAG‑AATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATCGACGTGGTGCACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTg      >  2:488921/1‑151 (MQ=255)                                                                                                                                                   aaaaC‑AACACCGG‑CAG‑AATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGaa    <  1:462016/151‑1 (MQ=255)                                                                                                                                                     aac‑aacACCGG‑CAG‑AATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAAt  >  1:817154/1‑151 (MQ=255)                                                                                                                                                     aac‑aacACCGG‑CAG‑AATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAAt  >  2:406345/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                         ACCGCAAATCGGCACTTCAACCGTTGGCGTAAAGTAGCGGATGCGCACATCGCTGTCGTCGCTGTGCAGCAGAAAAGCAGTTTCCGAATGGCCTAACTCGCGGGCGATAAGCTGCATTTGCGCTTCGCTAAGATTATCGGCGGGG‑AAAAC‑AACACCGG‑CAG‑AATTGCCGC‑GAAACGGTTGTGAGGTAAAAGCATCGACGTGGTACACCTGCGGTTTCATTAACGTTCTCCTGTGACTGGAGAACTATCATAGCCTGCAAGTGGCCGGAGAGCGAAGGGCTATCCGGCCAGGGTGAAAT  >  NC_000913/1535640‑1535937

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: ‑