Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,003,633	0	C	.	66.7%	47.3 / 62.7	39	coding (195/516 nt)	ycbV	putative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base C (6/7);  new base . (14/12);  total (20/19)
Fisher's exact test for biased strand distribution p-value = 7.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGGGGCTGTGAGCTATTTGCCCATGATGGCACGGTCAACATTAGCGGATCGTTTCGCCGTAATACATGCGTGCTGGCACAGGATAGCAAGCAAATCAATGTGCAGTTAGGCGATGTCTCGCTGACACGTTTTTCTCATGGGAATTATGGCCCGGAGAAATCTTTCATCATTAATTTGCAGGATTGCGGAACAGATGTCAGCACGGTTGACGTCACCTTTTCAGGAACCCCCGACGGAGTGCAGAGCGAAATGTTGTCCATTGAAAGTGGTACGGATGCAGCCAGCGGGCTGGCAATTGCG  >  NC_000913/1003483‑1003783                                                                                                                                                       |                                                                                                                                                       gTGGGGCTGTGAGCTATTTGCCCATGATGGCACGGTCAACATTAGCGGATCGTTTCGCCGTAATACATGCGTGCTGGCACAGGATAGCAAGCAAATCAATGTGCAGTTAGGCGATGTCTCGCTGACACGTTTTTCTCATGGGAATTATGGc                                                                                                                                                        <  2:37651/151‑1 (MQ=255)     ggCTGTGAGCTATTTGCCCATGATGGCACGGTCAACATTAGCGGATCGTTTCGCCGTAATACATGCGTGCTGGCACAGGATAGCAAGCAAATCAATGTGCAGTTAGGCGATGTCTCGCTGACACGTTTTTCTCATGGGAATTATGGCCCgg                                                                                                                                                    <  2:441653/151‑1 (MQ=255)         gtgAGCTATTTGCCCATGATGGCACGGTCAACATTAGCGGATCGTTTCGCCGTAATACATGCGTGCTGGCACAGGATAGCAAGCAAATCAATGTGCAGTTAGGCGATGTCTCGCTGACACGTTTTTCTCATGGGAATTATGG‑CCGGAGaaa                                                                                                                                               <  2:256106/151‑1 (MQ=255)            aGCTATTTGCCCATGATGGCACGGTCAACATTAGCGGATCGTTTCGCCGTAATACATGCGTGCTGGCACAGGATAGCAAGCAAATCAATGTGCAGTTAGGCGATGTCTCGCTGACACGTTTTTCTCATGGGAATTATGG‑CCGGAGAAATCt                                                                                                                                            >  2:541596/1‑151 (MQ=255)                    gCCCATGATGGCACGGTCAACATTAGCGGATCGTTTCGCCGTAATACATGCGTGCTGGCACAGGATAGCAAGCAAATCAATGTGCAGTTAGGCGATGTCTCGCTGACACGTTTTTCTCATGGGAATTATGGCCCGGAGAAATCTTtcatca                                                                                                                                     <  1:336043/151‑1 (MQ=255)                      ccATGATGGCACGGTCAACATTAGCGGATCGTTTCGCCGTAATACATGCGTGCTGGCACAGGATAGCAAGCAAATCAATGTGCAGTTAGGCGATGTCTCGCTGACACGTTTTTCTCATGGGAATTATGG‑CCGGAGAAATCTTTCATCAtt                                                                                                                                   >  1:23921/1‑150 (MQ=255)                       cATGATGGCACGGTCAACATTAGCGGATCGTTTCGCCGTAATACATGCGTGCTGGCACAGGATAGCAAGCAAATCAATGTGCAGTTAGGCGATGTCTCGCTGACACGTTTTTCTTATGGGAATTATGG‑CCGGAGAAATCTTTCATCATTaa                                                                                                                                 <  2:122137/151‑1 (MQ=255)                         tgatgGCACGGTCAACATTAGCGGATCGTTTCGCCGTAATACATGCGTGCTGGCACAGGATAGCAAGCAAATCAATGTGCAGTTAGGCGATGTCTCGCTGACACGTTTTTCTCATGGGAATTATGG‑CCGGAGAAATCTTTCATCATTAAtt                                                                                                                               >  2:309327/1‑151 (MQ=255)                          gatgGCACGGTCAACATTAGCGGATCGTTTCGCCGTAATACATGCGTGCTGGCACAGGATAGCAAGCAAATCAATGTGCAGTTAGGCGATGTCTCGCTGACACGTTTTTCTCATGGGAATTATGG‑CCGGAGAAATCTTTCATCATTAAttt                                                                                                                              <  2:411220/151‑1 (MQ=255)                               cACGGTCAACATTAGCGGATCGTTTCGCCGTAATACATGCGTGCTGGCACAGGATAGCAAGCAAATCAATGTGCAGTTAGGCGATGTCTCGCTGACACGTTTTTCTCATGGGAATTATGG‑CCGGAGAAATCTTTCATCATTAATTTGCAgg                                                                                                                         >  2:268681/1‑151 (MQ=255)                                aCGGTCAACATTAGCGGATCGTTTCGCCGTAATACATGCGTGCTGGCACAGGATAGCAAGCAAATCAATGTGCAGTTAGGCGATGTCTCGCTGACACGTTTTTCTCATGGGAATTATGG‑CCGGAGAAATCTTTCATCATTAATTTGCAGGa                                                                                                                        <  1:335786/151‑1 (MQ=255)                                   gTCAACATTAGCGGATCGTTTCGCCGTAATAAATGCGTGCTGGCACAGGATAGCAAGCAAATCAATGTGCAGTTAGGCGATGTCTCGCTGACACGTTTTTCTCATGGGAATTATGGCCCGGAGAAATCTTTCATCATTAATTTGCAGGAtt                                                                                                                      >  2:465582/1‑151 (MQ=255)                                      aaCATTAGCGGATCGTTTCGCCGTAATACATGCGTGCTGGCACAGGATAGCAAGCAAATCAATGTGCAGTTAGGCGATGTCTCGCTGACACGTTTTTCTCATGGGAATTATGGCCCGGAGAAATCTTTCATCATTAATTTGCAGGATTGCg                                                                                                                   <  1:295278/151‑1 (MQ=255)                                      aaCATTAGCGGATCGTTTCGCCGTAATACATGCGTGCTGGCACAGGATAGCAAGCAAATCAATGTGCAGTTAGGCGATGTCTCGCTGACACGTTTTTCTCATGGGAATTATGG‑CCGGAGAAATCTTTCATCATTAATTTGCAGGATTGCgg                                                                                                                  >  1:210167/1‑151 (MQ=255)                                           tAGCGGATCGTTTCGCCGTAATACATGCGTGCTGGCACAGGATAGCAAGCAAATCAATGTGCAGTTAGGCGATGTCTCGCTGACACGTTTTTCTCATGGGAATTATGG‑CCGGAGAAATCTTTCATCATTAATTTGCAGGATTGCGGAAc                                                                                                               >  2:536210/1‑149 (MQ=255)                                            aGCGGATCGTTTCGCCGTAATACATGCGTGCTGGCACAGGATAGCAAGCAAATCAATGTGCAGTTAGGCGATGTCTCGCTGACACGTTTTTCTCATGGGAATTATGGCCCGGAGAAATCTTTCATCATTAATTTGCAGGATTGCGGAACa                                                                                                              <  1:604322/150‑1 (MQ=255)                                             gCGGATCGTTTCGCCGTAATACATGCGTGCTGGCACAGGATAGCAAGCAAATCAATGTGCAGTTAGGCGATGTCTCGCTGACACGTTTTTCTCATGGGAATTATGG‑CCGGAGAAATCTTTCATCATTAATTTGCAGGATTGCGGAACAGAt                                                                                                           <  1:391729/151‑1 (MQ=255)                                                  tCGTTTCGCCGTAATACATGCGTGCTGGCACAGGATAGCAAGCAAATCAATGTGCAGTTAGGCGATGTCTCGCTGACACGTTTTTCTCATGGGAATTATGG‑CCGGAGAAATCTTTCATCATTAATTTGCAGGATTGCGGAACAGATGTc                                                                                                        >  1:403759/1‑149 (MQ=255)                                                   cGTTTCGCCGTAATACATGCGTGCTGGCACAGGATAGCAAGCAAATCAATGTGCAGTTAGGCGATGTCTCGCTGACACGTTTTTCTCATGGGAATTATGG‑CCGGAGAAATCTTTCATCATTAATTTGCAGGATTGCGGAACAGATGTCAGc                                                                                                     >  2:156060/1‑151 (MQ=255)                                                               aTACATGCGTGCTGGCACAGGATAGCAAGCAAATCAATGTGCAGTTAGGCGATGTCTCGCTGACACGTTTTTCTCATGGGAATTATGGCCCGGAGAAATCTTTCATCATTAATTTGCAGGATTGCGGAACAGATGTCAGCACGGTTGAc                                                                                            >  1:208177/1‑149 (MQ=255)                                                               aTACATGCGTGCTGGCACAGGATAGCAAGCAAATCAATGTGCAGTTAGGCGATGTCTCGCTGACACGTTTTTCTCATGGGAATTATGG‑CCGGAGAAATCTTTCATCATTAATTTGCAGGATTGCGGAACAGATGTCAGCACGGTTGACGTc                                                                                         <  2:571574/151‑1 (MQ=255)                                                                           tGGCACAGGATAGCAAGCAAATCAATGTGCAGTTAGGCGATGTCTCGCTGACACGTTTTTCTCATGGGAATTATGGCCCGGAGAAATCTTTCATCATTAATTTGCAGGATTGCGGAACAGATGTCAGCACGGTTGACGTCACCTTTTCAgg                                                                              >  2:525205/1‑151 (MQ=255)                                                                            ggCACAGGATAGCAAGCAAATCAATGTGCAGTTAGGCGATGTCTCGCTGACACGTTTTTCTCATGGGAATTATGGCCCGGAGAAATCTTTCATCATTAATTTGCAGGATTGCGGAACAGATGTCAGCACGGTTGACGTCACCTTTTCAg                                                                               <  1:524108/149‑1 (MQ=255)                                                                               acaGGATAGCAAGCAAATCAATGTGCAGTTAGGCGATGTCTCGCTGACACGTTTTTCTCATGGGAATTATGG‑CCGGAGAAATCTTTCATCATTAATTTGCAGGATTGCGGAACAGATGTCAGCACGGTTGACGTCACCTTTTCAGGAAccc                                                                         >  2:617780/1‑151 (MQ=255)                                                                                      agcaagcaAATCAATGTGCAGTTAGGCGATGTCTCGCTGACACGTTTTTCTCATGGGAATTATGGCCCGGAGAAATCTTTCATCATTAATTTGCAGGATTGCGGAACAGATGTCAGCACGGTTGACGTCACCTTTTCAGGAACCCCCGAc                                                                    <  2:420965/150‑1 (MQ=255)                                                                                       gcaagcaaATCAATGTGCAGTTAGGCGATGTCTCGCTGACACGTTTTTCTCATGGGAATTATGG‑CCGGAGAAATCTTTCATCATTAATTTGCAGGATTGCGGAACAGATGTCAGCACGGTTGACGTCACCTTTTCAGGAACCCCCGACgg                                                                  >  2:250409/1‑150 (MQ=255)                                                                                       gcaagcaaATCAATGTGCAGTTAGGCGATGTCTCGCTGACACGTTTTTCTCATGGGAATTATGG‑CCGGAGAAATCTTTCATCATTAATTTGCAGGATTGCGGAACAGATGTCAGCACGGTTGACGTCACCTTTTCAGGAACCCCCGACGGa                                                                 <  1:425407/151‑1 (MQ=255)                                                                                              aaTCAATGTGCAGTTAGGCGATGTCTCGCTGACACGTTTTTCTCATGGGAATTATGG‑CCGGAGAAATCTTTCATCATTAATTTGCAGGATTGCGGAACAGATGTCAGCACGGTTGACGTCACCTTTTCAGGAACCCCCGACGGAGTGCaga                                                          <  2:124941/151‑1 (MQ=255)                                                                                                   aTGTGCAGTTAGGCGATGTCTCGCTGACACGTTTTTCTCATGGGAATTATGG‑CCGGAGAAATCTTTCATCATTAATTTGCAGGATTGCGGAACAGATGTCAGCACGGTTGACGTCACCTTTTCAGGAACCCCCGACGGAGTGCAGAGCGa                                                      >  1:571032/1‑150 (MQ=255)                                                                                                     gtgCAGTTAGGCGATGTCTCGCTGACACGTTTTTCTCATGGGAATTATGG‑CCGGAGAAATCTTTCATCATTAATTTGCAGGATTGCGGAACAGATGTCAGCACGGTTGACGTCACCTTTTCAGGAACCCCCGACGGAGTGCAGAGCGAAAt                                                   <  2:210168/151‑1 (MQ=255)                                                                                                        cAGTTAGGCGATGTCTCGCTGACACGTTTTTCTCATGGGAATTATGG‑CCGGAGAAATCTTTCATCATTAATTTGCAGGATTGCGGAACAGATGTCAGCACGGTTGACGTCACCTTTTCAGGAACCCCCGACGGAGTGCAGAGCGAAAtgtt                                                <  1:23244/151‑1 (MQ=255)                                                                                                                            gACACGTTTTTCTCATGGGAATTATGG‑CCGGAGAAATCTTTCATCATTAATTTGCAGGATTGCGGAACAGATGTCAGCACGGTTGACGTCACCTTTTCAGGAACCCCCGACGGAGTGCAGAGCGAAATGTTGTCCATTGAAAGTGGTACgg                            >  2:171043/1‑151 (MQ=255)                                                                                                                                cGTTTTTCTCATGGGAATTATGGCCCGGAGAAATCTTTCATCATTAATTTGCAGGATTGCGGAACAGATGTCAGCACGGTTGACGTCACCTTTTCAGGAACCCCCGACGGAGTGCAGAGCGAAATGTTGTCCATTGAAAGTGGTACGGAt                          >  1:209324/1‑150 (MQ=255)                                                                                                                                cGTTTTTCTCATGGGAATTATGGCCCGGAGAAAGCTTTCATCATTAATTTGCAGGATTGCGGAACAGATGTCAGCACGATTGACGTGACCTTTTCAGGAACCCCCGACGGAGTGCAGAGCGAAATGTTGTCCATTGAAAGTGGTACGGAt                          <  1:270712/150‑1 (MQ=255)                                                                                                                                  tttttCTCATGGGAATTATGGCCCGGAGAAATCTTTCATCATTAATTTGCAGGATTGCGGAACAGATGTCAGCACGGTTGACGTCACCTTTTCAGGAACCCCCGACGGAGTGCAGAGCGAAATGTTGTCCATTGAAAGTGGTACGGATGc                        >  1:584890/1‑150 (MQ=255)                                                                                                                                   ttttCTCATGGGAATTATGG‑CCGGAGAAATCTTTCATCATTAATTTGCAGGATTGCGGAACAGATGTCAGCACGGTTGACGTCACCTTTTCAGGAACCCCCGACGGAGTGCAGAGCGAAATGTTGTCCATTGAAAGTGGTACGGATGcag                      <  2:571694/150‑1 (MQ=255)                                                                                                                                       ctcATGGGAATTATGG‑CCGGAGAAATCTTTCATCATTAATTTGCAGGATTGCGGAACAGATGTCAGCACGGTTGACGTCACCTTTTCAGGAACCCCCGACGGAGTGCAGAGCGAAATGTTGTCCATTGAAAGTGGTACGGATGcagccagc                 >  2:639712/1‑151 (MQ=255)                                                                                                                                        tcATGGGAATTATGG‑CCGGAGAAATCTTTCATCATTAATTTGCAGGATTGCGGAACAGATGTCAGCACGGTTGACGTCACCTTTTCAGGAACCCCCGACGGAGTGCAGAGCGAAATGTTGTCCATTGAAAGTGGTACGGATGCAGCCAGCg                >  2:116578/1‑151 (MQ=255)                                                                                                                                         cATGGGAATTATGG‑CCGGAGAAATCTTTCATCATTAATTTGCAGGATTGCGGAACAGATGTCAGCACGGTTGACGTCACCTTTTCAGGAACCCCCGACGGAGTGCAGAGCGAAATGTTGTCCATTGAAAGTGGTACGGATGcagccagc                 <  2:53615/149‑1 (MQ=255)                                                                                                                                                      gCCCGGAGAAATCTTTCATCATTAATTTGCAGGATTGCGGAACAGATGTCAGCACGGTTGACGTCACCTTTTCAGGAACCCCCGACGGAGTGCAGAGCGAAATGTTGTCCATTGAAAGTGGTACGGATGCAGCCAGCGGGCTGGCAATTGc   >  2:583404/1‑151 (MQ=255)                                                                                                                                                       cccGGAGAAATCTTTCATCATTAATTTGCAGGATTGCGGAACAGATGTCAGCACGGTTGACGTCACCTTTTCAGGAACCCCCGACGGAGTGCAGAGCGAAATGTTGTCCATTGAAAGTGGTACGGATGCAGCCAGCGGGCTGGCAATTGCg  <  2:169143/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                       GTGGGGCTGTGAGCTATTTGCCCATGATGGCACGGTCAACATTAGCGGATCGTTTCGCCGTAATACATGCGTGCTGGCACAGGATAGCAAGCAAATCAATGTGCAGTTAGGCGATGTCTCGCTGACACGTTTTTCTCATGGGAATTATGGCCCGGAGAAATCTTTCATCATTAATTTGCAGGATTGCGGAACAGATGTCAGCACGGTTGACGTCACCTTTTCAGGAACCCCCGACGGAGTGCAGAGCGAAATGTTGTCCATTGAAAGTGGTACGGATGCAGCCAGCGGGCTGGCAATTGCG  >  NC_000913/1003483‑1003783

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