Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	907,580	0	T	G	41.7%	46.2 / 11.4	36	I235L (ATC→CTC)	ybjT	putative NAD‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (16/5);  new base G (0/15);  total (16/20)
Fisher's exact test for biased strand distribution p-value = 3.22e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGTTGTGGGATGAGTGCACGTAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGG‑A‑TCAACCAGCGGCGCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTT  >  NC_000913/907433‑907729                                                                                                                                                    |                                                                                                                                                         gCCGTTGTGGGATGAGTGCACGTAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGggggatggg                                                                                                                                                       <  1:215601/151‑1 (MQ=255)  ccGTTGTGGGATGAGTGCACGTAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGggggatgggg                                                                                                                                                      <  2:723059/151‑1 (MQ=255)   ctttGTGGGATGAGTGCACGTAGCGCGGTATCATCCGCCTGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGggggatgt                                                                                                                                                        <  2:308719/146‑2 (MQ=255)            aTGAGGGCACGTAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAAACACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGG‑A‑TCAACCAgc                                                                                                                                          >  2:345835/1‑151 (MQ=255)               aGTGCACGTAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGG‑A‑TCAACCAgcggc                                                                                                                                       <  1:138892/151‑1 (MQ=255)                   cACGTAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGG‑A‑TCAACCAGCGGCGCt                                                                                                                                    <  1:515485/150‑1 (MQ=255)                    aCGTAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGG‑A‑TCAACCAGCGGCGCtt                                                                                                                                   >  1:473199/1‑150 (MQ=255)                                           tcagATCGTGTTTCAGCCCCTGAATCACCGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGG‑A‑TCAACCAGCGGCGCTTA‑CCGCTACCCGCCATAAAATGTTc                                                                                                          >  2:417131/2‑151 (MQ=255)                                            cagATCGGGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGG‑A‑TCAACCAGCGGCGCTTA‑CCGCTCACCGCCATAAAATGTTc                                                                                                          >  2:96347/1‑150 (MQ=255)                                                    gtTTCACCCCCTAAATCACCCCCCGGGCGGTGGGGGGCGGTACGGAAGTAATCACATTGAAAAACCACCCCGAAACCCAGCGGGTGGGGAGGGGGAGGGGG‑A‑TCAACCGGCGGCGCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTg                                                                                                    <  2:43885/148‑1 (MQ=255)                                                     ttccAGCCCCGGAACCACCGCCCGGGGGGGGGTGGGGGGTACGGAATTAATCACATTGAAAACCCCCCCCGAATTCGAGCGGGGGGGGAGGGGGGGGGGG‑A‑TAACCCGGAGGCGCTCA‑CCGCTCACCGCCATAAAATGTTCAAACTg                                                                                                    <  2:457610/144‑1 (MQ=255)                                                        cAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGG‑A‑TCAACCAGCGGCGCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTg                                                                                              >  2:294010/1‑150 (MQ=255)                                                         aGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGG‑A‑TCAACCAGCGGCGCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGa                                                                                             >  2:144589/1‑150 (MQ=255)                                                          gCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGG‑A‑TCAACCAGCGGCGCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGAt                                                                                            >  1:369705/1‑150 (MQ=255)                                                           ccccTAAACCACCCCCCTGGGGGGGGGGGGCGGTAGGAAAGTAATCAAATGGAAAACCCCCACCGAAACCCGGGGGGGGGGGAGGGGGAGGGGA‑A‑TCAGCCGGCGGCGCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATaa                                                                                          <  2:218145/151‑1 (MQ=255)                                                           ccccAGAAACACCGCCCGGGGGGGGGGGGGGGGTAAGAAAGTAACCAAATGGAAACCCCAC‑CCTAAAACCCGGGGGGGGGGGGGGGGATGGGG‑ACCCAGCCAGCGGGGTTTC‑CCTTGCTCCGCCATAAAATGTTCAAACTGTTGCTGa                                                                                             <  2:517645/148‑1 (MQ=255)                                                               tgaataacgccccTGGCGGTGGGGGGCGGTACGGAATAAATCACATTAAAAAACCCCCCCAAAACCCGGCGGGTGGTGAGGGGGAGGGTG‑A‑TCAACCAGCGGCGTTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATa                                                                                           <  1:524716/136‑1 (MQ=255)                                                                  aTCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGG‑A‑TCAACCAGCGGCGCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTg                                                                                       >  2:25653/1‑147 (MQ=255)                                                                       cGCCCGGGGGGTGGGGGGGGGTACGGAAGTAATCAAATTAAAAACCCCCACCAAAACCCAGCGGGTGGGGGGGGGGAGGGGA‑A‑TCACCCGGCGGCGCTTC‑CCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACctct                                                                              <  2:180454/151‑1 (MQ=255)                                                                       cGCCCGGGGGGGGGGGGGGGGGAGGGAATTAAACACATAAAAAACCCACCCAGAAACCCGGGGGGGGGGAAGGGGGTGGGGG‑A‑TAAACCAGCGGCGCTTC‑CCCTCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACct                                                                                <  2:445314/149‑1 (MQ=255)                                                                            tGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGG‑A‑TCAACCAGCGGCGCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTcc                                                                         >  1:94221/1‑151 (MQ=255)                                                                            tGGCGGTGGGGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGGGGGGAGGGGGATGGGG‑A‑TCAACCAGCGGCGCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTcc                                                                         >  2:322895/1‑151 (MQ=255)                                                                             ggCGGTGGTGGGGGGTACGGAAGTAATCACATGGAAAACCCCCCCCAAAACCCGGGGGGGGGGGGGGGGGAGGGGA‑A‑TCAACCAGCGGCGCTTC‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcg                                                                        <  1:332662/151‑1 (MQ=255)                                                                             ggCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGG‑A‑TCAACCAGCGGCGCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcg                                                                        >  1:633495/1‑151 (MQ=255)                                                                                         cGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGG‑A‑TCAACCAGCGGCGCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGa                                                              >  2:572145/1‑149 (MQ=255)                                                                                              cGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGG‑A‑TCAACCAGCGGCGCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCg                                                       >  1:358244/1‑151 (MQ=255)                                                                                               gaaaGAAACCACATTGAAAAACCACACCGAAATCCAGGGGGGGGGGAGGGGGATGGGA‑A‑TCAACCAGCGGCCCTTA‑CCGCTCCCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGAt                                                          <  2:671149/145‑1 (MQ=255)                                                                                                     aaTCCCATGGAAAAACCACACCAAAATCCAGGGGGGGGGGGGGGGGAGGGGG‑A‑CCACCCGGCGGC‑CTTACCCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGtt                                                  <  1:244378/149‑1 (MQ=255)                                                                                                            ttGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGG‑A‑TCAACCAGCGGCGCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGcc                                          >  1:664423/1‑150 (MQ=255)                                                                                                                 aaaCCACCCCGAAACCCGGGGGGGGGGGGGGGGGAGGGGA‑A‑TCACCCGGCGGCGCTTC‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGgatg                                     <  1:74514/150‑1 (MQ=255)                                                                                                                       cacCGAAATCCAGCGGGTGGGGAGGGGGATGGGG‑A‑TCAACCAGCGGCGCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTaac                              >  1:600739/1‑151 (MQ=255)                                                                                                                         ccgaaccccgggggggggggggggggATGGGG‑A‑TCAACGGGCGGCTCTTC‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTaa                               <  2:417000/136‑1 (MQ=255)                                                                                                                            aaacccggggtggggggAGGGGGGGGGGG‑A‑TCACCCAGCGTCTCTTC‑CCGCTCTCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACg                          <  2:681559/137‑1 (MQ=255)                                                                                                                            aaaccagggggggggggggggggAGGGGA‑A‑TAACCCGGCGCCCCTTA‑CCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTaa                               <  1:226801/136‑1 (MQ=255)                                                                                                                            aaaccaaGCGGGGGGGGGGGGGGAGGGGG‑A‑TCAACCGGCGGCGCTTA‑CCGCAAACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTaacaa                            <  1:783931/142‑1 (MQ=255)                                                                                                                                    ggggTGGGGAGGGGGAGGGGG‑A‑CCAACAGGCGCCGTTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGAt                 <  1:357294/150‑1 (MQ=255)                                                                                                                                     ggggggggggggggAGGGGG‑A‑TCCACCAGCGGCGCTTC‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATa                <  1:43958/151‑1 (MQ=255)                                                                                                                                       ggggggggggggAGGGGGAA‑CCAACCGGCGGCGCTTA‑CCCCCCCCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAg               <  2:184947/149‑1 (MQ=255)                                                                                                                                          gggAGGGGGATGGGG‑A‑TAAACCAGCGGCGCTTA‑CCGCTCACCGCCATAAAATGTTCAAACGGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTgc            >  2:85511/1‑150 (MQ=255)                                                                                                                                              gggggaggggg‑a‑TCAACCAGCGGCGCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCaa       <  1:596528/151‑1 (MQ=255)                                                                                                                                                 ggatgggg‑a‑tCAACCAGCGGCGCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGtt    >  1:705837/1‑151 (MQ=255)                                                                                                                                                   atgggg‑a‑tCAACCAGCGGCGCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGtttt  >  1:81773/1‑151 (MQ=255)                                                                                                                                                    tgggt‑a‑tCAACCCGCGTCGCTCC‑CCGCTCGCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGtttt  <  2:542455/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                         GCCGTTGTGGGATGAGTGCACGTAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGG‑A‑TCAACCAGCGGCGCTTA‑CCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTT  >  NC_000913/907433‑907729

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