Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	907,602	0	A	C	53.1%	40.5 / 69.0	49	G227G (GGT→GGG)	ybjT	putative NAD‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (16/7);  new base C (0/26);  total (16/33)
Fisher's exact test for biased strand distribution p-value = 7.32e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.19e-07

AGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAG‑CGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCGTGGTGC  >  NC_000913/907456‑907750                                                                                                                                                    |                                                                                                                                                     aGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAG‑CGGCGCTTACCg                                                                                                                                                   <  2:764079/150‑1 (MQ=255) aGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAG‑CGGCGCTTACCg                                                                                                                                                   <  1:704532/150‑1 (MQ=255)  gcgcGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAG‑CGGCGCTTACCGCt                                                                                                                                                 <  1:349499/151‑1 (MQ=255)         atcatcCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAG‑CGGCGCTTACCGCTCACCGcc                                                                                                                                          >  2:47333/1‑151 (MQ=255)           catcCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAG‑CGGCGCTTACCGCTCACCGCCAt                                                                                                                                        >  2:842406/1‑151 (MQ=255)              ccGCGAGCAGATCGTGTTCCACCCCCTGAATCACCGCCCTGGCGGTGGTGGGCGGTACGAAAGTAATCACATGGAAAAACCACACCGAAACCCAGCGGGTGGGGAGGGGGAGGGGGTTCACCAAG‑CGGCGCTTCCCGCTCACCGCCATaa                                                                                                                                      <  1:403951/150‑1 (MQ=255)                          ggtgtTTCACCCCCAAAATCAACCCCCGGGGGGGGGGGGGCGGAACGGAAATAATCAAATTGAAAACCCCCACCGAAACCCAGGGGGTGGGGGGGGGGAGGGGGATCACCCAG‑CGGCGCTTCCCCCTCACCGCCATAAAATGTTCAAACTg                                                                                                                         <  1:107452/150‑1 (MQ=255)                               ttCACCCCCTGAATCAACGCCCGGGCGGGGGGGGGGGGTGGGGAAGTAATCAAATTGAAAACCCACCCCAAAATCCGGGGGGGGGGGGGGGGGAGGGGGATCAACCGG‑CGGCGTTCCCCCCTCACCGCCATAAAATGTTCAAACTGTTGc                                                                                                                     <  2:463389/150‑1 (MQ=255)                                tcccccccTGAATCACCGCCCGGGGGGTGGGGGGCGGTACGGAAGTAATCACATTGAAAACCCCCCCCAAAACCCGGGGGGTGGGGGGGGGGATGGGAATCACCCAG‑CGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCt                                                                                                                    <  1:411617/146‑1 (MQ=255)                                 cAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAG‑CGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTg                                                                                                                   >  2:504158/1‑150 (MQ=255)                                    ccctTGAACCACCGCCCGGGGGGTGGGGGGGGGTAGGAAAGTAATCACATTGAAAACCCACCCCAAAACCCGGGGGGGGGGGGGGGGGATGGGGATCACCCGG‑GGGCTTTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATa                                                                                                                <  2:303445/150‑1 (MQ=255)                                      ccTGAACCACCGCCCGGGGGGGGGGGGGGGGTAGGAAAGTAATCACATGGAAAACCCCCCCCAAAACCCAGGGGGGGGGGGGGGGGAGGGGGATCACCCAG‑CGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACt                                                                                                             <  2:276806/151‑1 (MQ=255)                                          aaccaccGCCCTGGGGGTGGTGGGGGGTACGGAATTAATCACATGGAAAACCCACCCCGAAACCCGGGGGGGGGGGGGGGGGAGGGGGATCAACCGG‑CGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGc                                                                                                         <  2:713211/145‑1 (MQ=255)                                          aaccaccGCCCGGGGGGTGGGGGGGGGTACGGAAGAAATCACTTGGAAAACCCCCCCCAAAACCCGGGGGGTGGGGGGGGGGAGGGGGATCAACCAG‑CGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGa                                                                                                           <  2:347432/143‑1 (MQ=255)                                          aaacAACCCCCGGGGGGGGGTGGGGGGTAGGAAAGTAATCACATTGAAAACCCACACCAAAACCCGGGGGGTGGGGGGGGGGAGGGGGATCACCCGG‑CGGCGCTTCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTg                                                                                                            <  1:22351/145‑1 (MQ=255)                                                cccccTGGGGGGGGTGGGGGGTACGGAAGTAATCACATTGAGAACCCCCACCAAAATCCGGGGGGTGGGGGGGGGGGTGGGGATCACCCGG‑GGGCGTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACctct                                                                                                   <  2:636936/149‑1 (MQ=255)                                                cGCCCTGGGGGGGGGGGGGGGTAGGGAAGTAATCACATTGAAAACCCACCCCAAAACCCGGGGGGGGGGGGGGGGGAGGGGGATCAACCGG‑CGGCGTTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACctct                                                                                                   <  2:734718/151‑1 (MQ=255)                                                cGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGGGGGGAGGGGGATGGGGATCAACCAG‑CGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACctc                                                                                                    >  1:194116/1‑150 (MQ=255)                                                  ccggggggggggggggCGGTGAGGAAGTAATCACTGGAAAAAACCACCCCGAAACCCAGGGGGGGGGGAGGGGGATGGGGATCAACCGG‑GGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTg                                                                                                  <  1:758820/137‑1 (MQ=255)                                                      gggggtgggGGGGGGTACGGAAATAATCACATGGAAAAACCACCCCAAAACCCGGGGGGGGGGGGGGGGGATGGGGATCACCCAG‑CGGCGCTTCCCCCTCCCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcg                                                                                             <  1:986613/148‑1 (MQ=255)                                                      ggggggggtggGCGGTGAGGAAGAAATCACATGGAAAAACCCCCCCAAAACCCGGGGGGTGGGGAGGGGGTGGGGGATCACCCAG‑CGGCGCTCCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcg                                                                                             <  1:47333/145‑1 (MQ=255)                                                      gggggggggggggggTAGGGAAGAAAACACATTGAAAACCCACACCAAAACCCGGGGGGGGGGGAGGGGGATGGGGACAACCCAG‑CGGCCCTTCCCGCTCTCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcg                                                                                             <  1:127854/142‑1 (MQ=255)                                                        gggtgggGGGCGGTGCGGAAGTAATCACATAGAAAACCCACCCCGAAACCCGGGGGGTGGGGAGGGGGAGGGGGATCAACCAG‑CGGCGCTTCCCTCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcggc                                                                                           <  2:735387/150‑1 (MQ=255)                                                                          aaGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAG‑CGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGgtg                                                                         >  1:675408/1‑151 (MQ=255)                                                                           aGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAG‑CGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGgtg                                                                         >  1:828563/1‑150 (MQ=255)                                                                              aaTCAAATTGAAAACCCACCCCAAAACCCGGGGGGGGGGGGGGGGGAGGGGAATCAACCGG‑CGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTc                                                                      <  2:34877/150‑1 (MQ=255)                                                                                 cacaTTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAG‑CGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTg                                                                  >  1:338343/1‑151 (MQ=255)                                                                                     ttGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAG‑CGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCg                                                              >  1:615116/1‑151 (MQ=255)                                                                                      tGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAG‑CGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCgg                                                             >  2:307226/1‑151 (MQ=255)                                                                                           aaCCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAG‑CGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGgatgat                                                        >  2:482231/1‑151 (MQ=255)                                                                                                   cGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAG‑CGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTaacaac                                                >  1:984539/1‑151 (MQ=255)                                                                                                    gAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAG‑CGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACg                                               >  2:78750/1‑151 (MQ=255)                                                                                                     aaacccggggggTGGGGGGGGGGAGGGGGATCAACCAG‑CGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACg                                               <  2:746289/141‑1 (MQ=255)                                                                                                         ccAGCGGGTGGGGAGGGGGATGGGGATCAACCAG‑CGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAAcc                                          >  1:863606/1‑151 (MQ=255)                                                                                                             gggggggggggggggggggGAATCACCCAG‑CGGCCCTCCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGAt                                      <  1:432737/150‑1 (MQ=255)                                                                                                              ggggggggggggggATGGGAATCACCCGG‑CGGCCCTCCCCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCa                                         <  1:348408/147‑1 (MQ=255)                                                                                                              ggggggggggggggAGGGGGATCAACCAG‑CGGCGTTTCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATa                                     <  2:707095/151‑1 (MQ=255)                                                                                                               gggggggggggggAGGGGAATCACCCAG‑CGGCGTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAg                                    <  2:232043/148‑1 (MQ=255)                                                                                                                  ggggggggggGTGGGAATCACCCAG‑CGGTCCTACCCCCCCCCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTg                                  <  1:783997/150‑1 (MQ=255)                                                                                                                     gggggggagggggaTCACCCGG‑GGGCGTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTgcagc                              <  1:685854/149‑1 (MQ=255)                                                                                                                                gtgacaaacgacgggccCTTCCCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGtttt                       <  2:170125/135‑1 (MQ=255)                                                                                                                                ggacCACCCGG‑CGCCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCGCCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAg                   <  2:799918/151‑1 (MQ=255)                                                                                                                                    gcaCCAGCCGGCGCTACCCCCCCCCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAgc                  <  1:402188/147‑1 (MQ=255)                                                                                                                                    cAACCAG‑CGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAgcgc                >  1:161015/1‑150 (MQ=255)                                                                                                                                       ccgg‑cggcgCTTACCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAgcg                 <  1:674537/143‑1 (MQ=255)                                                                                                                                            cggccCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATggg          <  2:585591/149‑1 (MQ=255)                                                                                                                                             ggcgCTCACCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGc         <  2:820964/149‑1 (MQ=255)                                                                                                                                                  ttACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCGTGGTGc  >  1:611037/1‑151 (MQ=255)                                                                                                                                                  ttACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCGTGGTGc  >  2:772853/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                     AGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGGATGGGGATCAACCAG‑CGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCGTGGTGC  >  NC_000913/907456‑907750

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