Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,717,014	0	T	G	46.9%	76.4 / 10.1	49	G156G (GGT→GGG)	yiaD	multicopy suppressor of bamB; outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (7/19);  new base G (23/0);  total (30/19)
Fisher's exact test for biased strand distribution p-value = 3.54e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAA  >  NC_000913/3716866‑3717163                                                                                                                                                     |                                                                                                                                                      tCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGGTAACGTGGTTGGGTATACCGACAGCACgggggg                                                                                                                                                     >  1:462327/1‑151 (MQ=255)       atatGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGGTTGGGTATACCGACAGCACGGGGGGTCAc                                                                                                                                                 >  1:345458/1‑149 (MQ=255)            ccGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGa                                                                                                                                          <  2:543377/151‑1 (MQ=255)             cGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGGTAACGTGGTTGGGTATACCGACAGCACGGGGGGGCGCGGCCTg                                                                                                                                           >  1:92566/1‑149 (MQ=255)              gAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAAc                                                                                                                                        <  2:822870/151‑1 (MQ=255)               aaCAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGGTTGGTTATAACGACAACACGGGGGGTCACGACCTGaa                                                                                                                                         >  2:842518/1‑149 (MQ=255)               aaCAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAAAGTGATTGGGTATAACGACAGCACGGGGGGGCACGACCTGaa                                                                                                                                         >  2:464900/1‑149 (MQ=255)                aCAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACAt                                                                                                                                      <  2:763264/151‑1 (MQ=255)                 cAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTTAACTGATTTGGTTTACCGGCAACCCGGGGGGGCACGCCCTGaaaa                                                                                                                                       >  2:401079/1‑147 (MQ=255)                 cAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATg                                                                                                                                     <  2:530222/151‑1 (MQ=255)                  aaTGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTTACGTGGTTGGTTTTACCGACAACACCGGGGGGCACGACCTCAACa                                                                                                                                       >  1:525317/1‑148 (MQ=255)                  aaTGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGGTAAAGTGATTTGGTATTCCGACCGCCCGGGGGGTCACGAccctaa                                                                                                                                         >  1:918831/1‑142 (MQ=255)                        aCCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAAAGTGATTGGTTATACCGACAGCACGGGGGGGCACGACCTGAACAAGCGtctc                                                                                                                               >  2:718142/1‑150 (MQ=255)                         ccTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGGTAAAGGGATTGGTTATACCCACAGCAAGGGGGGGCACGGCCCGAACATGc                                                                                                                                    >  1:777921/1‑144 (MQ=255)                          cTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGGTTAAGTGATTGGTTATACCGACAACACGGGTGGGTGCGGCCCGAAAATGCGtctctc                                                                                                                             >  1:151714/1‑150 (MQ=255)                            tCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGGTTAAGCGATTGGTTATACCGGCAGAACGGGGGGGCACGAACTGAACATGCGGCTCTccc                                                                                                                           >  1:783758/1‑150 (MQ=255)                             cGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAg                                                                                                                         <  2:356175/151‑1 (MQ=255)                             cGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGGGATTTGGTTTACCGACAACACGGGGGGGCACGACCTGAACATGCGTCTCTccc                                                                                                                           >  2:674505/1‑149 (MQ=255)                              gACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGATTGGTTATAACGACAGCACGGGGGGGCACGACCCGAAAAAGAGGCTCTccc                                                                                                                           >  1:911695/1‑148 (MQ=255)                                  gcagcagcGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGGGATTGGGTATACCGACAGCACGGGGGGTCACGACCTGAAAATGCGTCTCTCCCCGCAAc                                                                                                                     >  1:789565/1‑150 (MQ=255)                                  gcagcagcGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGGTAACGTGATTGGTTTTACCCACAGCACGGGGGGTCTCGAGCCGAAAATGCGTCTCCCCCacaaaa                                                                                                                     >  2:512419/1‑145 (MQ=255)                                   cagcagcGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGGTAACGTGGTTGGTTATACCGACAACACGGGGGGTCACGACCTGGACAAGCGGCTCCCCCCGCaa                                                                                                                      >  1:200209/1‑148 (MQ=255)                                       agcGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCg                                                                                                                <  1:685788/150‑1 (MQ=255)                                               ccTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGtt                                                                                                       <  1:466581/151‑1 (MQ=255)                                                  gAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTg                                                                                                      <  2:382258/149‑1 (MQ=255)                                                  gAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTg                                                                                                      <  2:878846/149‑1 (MQ=255)                                                        ggcgggcgCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGGTAACGTGATTGGGTATACCCACAGCACGGGGGGGCACGACCTGAACATGCGTCTCTCCCCGCAACGTGCGGGTTCCGTTGCCAGcgcg                                                                                              >  1:134685/1‑151 (MQ=255)                                                           gggcgCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGATTGGTTATACCGACAACACGGGGGGGCACGACCTGAACATGCGTCTCTCCCAGCAACGGGCGGATTCCGTTGCCAGcgcggt                                                                                            >  2:596045/1‑148 (MQ=255)                                                                 tAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGATTGGGTATACCGACAGCACGGGGGGTCACCACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCAcc                                                                                     >  2:648689/1‑151 (MQ=255)                                                                    cacCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGATTGGGTATACCGACAGCACGGGGGGTCACGACCTGAACATGCGTCTCTCCCAGCCAAGTGCGGATTCCGTTGCCAGCGCGTTGATCAccc                                                                                    >  1:891995/1‑149 (MQ=255)                                                                       ccTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGc                                                                               <  2:415656/151‑1 (MQ=255)                                                                              ggCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGGTTGGGTATACCGACAGCACGGGGGGTCACGACCTGAACATGCGTCTCTCCCCGCAACGGGCGGGTTCCGTTGCCAGCGCGTTGATCACCCCGGGCgggg                                                                           >  1:113878/1‑148 (MQ=255)                                                                                    gCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGcc                                                                  <  1:132463/151‑1 (MQ=255)                                                                                          gTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATcc                                                            <  1:488868/151‑1 (MQ=255)                                                                                               gAAAGAGTATCCGAAAACGGCGGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTa                                                         <  2:659815/149‑1 (MQ=255)                                                                                                  agagTATCCGAAAACGGCGGTTAACGTGATTGGTTATACCGACAGCACGGGGGGTCACGACCTGAACATGCGTCTCTCCCCGCAACGTGCGGATTCCGTTGCCAGCGCGGTGATTACCCAGGGCGTGGGCGCCAGCCGCATCCGTAccc                                                      >  1:430096/1‑147 (MQ=255)                                                                                                       aTCCGAAAACGGCGGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGcc                                                <  1:780208/150‑1 (MQ=255)                                                                                                         ccGAAAACGGCGGTTAACGTTATTGGGTATACCGACAGCACGGGGGGTCACGGCCTGAACATGGGTCTTTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGGTGATGACCCAGGGCGGGGACCCCAGCCGCATCCTTACTCAgg                                                   >  1:723068/1‑145 (MQ=255)                                                                                                         ccGAAAACGGCGGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCtt                                              <  2:728071/150‑1 (MQ=255)                                                                                                            aaaaCGGCGGCTAACGTGATTGGTTATACCCACAGCACCGGTGGGCACGACCTGAACAAGCGTCTCTCCCAGCAACGTGCGGATTCCCCTGCCAGCGCGGTGATCACCCAGGGCGTGGACGCCAGCCCCATCCGCACTCCGGGCCTTGGc                                           >  1:442295/1‑150 (MQ=255)                                                                                                               aCGGCGGTTAAAGTGATTTGTTATACCGACAGCACGGGGGGGCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTTCGTTGCCAGCGCGTTGATCACCCAGGGCGGGGGCGCCAGCCGCATCCCTACTCAGGGCCCTTGccc                                         >  1:821533/1‑149 (MQ=255)                                                                                                                cggcggTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGc                                      <  1:200055/151‑1 (MQ=255)                                                                                                                  gcggTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTa                                    <  1:446649/151‑1 (MQ=255)                                                                                                                       tAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAAccc                                >  1:784678/1‑150 (MQ=255)                                                                                                                            tGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGc                          >  1:62042/1‑151 (MQ=255)                                                                                                                                       aCCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCTCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGc                >  1:655354/1‑150 (MQ=255)                                                                                                                                          gACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCg            <  1:569749/151‑1 (MQ=255)                                                                                                                                              gCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGc           >  2:534146/1‑148 (MQ=255)                                                                                                                                               cACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCGGa        >  1:218237/1‑150 (MQ=255)                                                                                                                                                aCGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAg      <  1:527067/151‑1 (MQ=255)                                                                                                                                                    gtggtCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTaa  <  2:220869/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      TCCTCAATATGCCGAACAATGTGACCTTCGACAGCAGCAGCGCGACCCTGAAACCGGCGGGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAA  >  NC_000913/3716866‑3717163

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