Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,001,297	1	.	G	32.1%	66.7 / 74.4	53	coding (953/2259 nt)	xdhA	xanthine dehydrogenase, molybdenum binding subunit
Reads supporting (aligned to +/- strand):  ref base . (19/17);  new base G (8/9);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 7.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCGTTGAGTC  >  NC_000913/3001162‑3001432                                                                                                                                         |                                                                                                                                        cAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑ggg                                                                                                                                      >  1:277344/1‑139 (MQ=255)  aaCCCTTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTGACGGGCACTCCATCGCTTCTGC‑gggg                                                                                                                                     <  2:265429/139‑1 (MQ=255)            cGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑GGGGGGGAATAAAg                                                                                                                           >  1:1148670/1‑139 (MQ=255)              ccACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑g                                                                                                                                        >  1:712140/1‑124 (MQ=255)              ccACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑g                                                                                                                                        <  2:712140/124‑1 (MQ=255)               cACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑GGGGGGGAATAAAGTCg                                                                                                                        >  2:603944/1‑139 (MQ=255)                              gACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCt                                                                                                          >  2:1296593/1‑139 (MQ=255)                                 gggCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGt                                                                                                       <  1:1296593/139‑1 (MQ=255)                                      attagggacgAACCGCGACCGAACACTGCCAGCTTATAGTCTGGATGTTCTGTCTAACACCGGCGGTTATGCATCTCACGGGCACTCCAGCGCTTCTGC‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGca                                                                                                 <  2:1430962/130‑2 (MQ=255)                                       aTGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGcc                                                                                                 >  2:861377/1‑139 (MQ=255)                                       aTGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGcc                                                                                                 <  1:467725/139‑1 (MQ=255)                                              tGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTa                                                                                         <  1:359208/139‑1 (MQ=255)                                                             cATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCt                                                                                                                      <  2:613365/96‑1 (MQ=255)                                                             cATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCt                                                                                                                      >  1:613365/1‑96 (MQ=255)                                                                tGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑GGTGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCta                                                                       <  1:798625/139‑1 (MQ=255)                                                                        tataGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAAc                                                                <  2:277344/139‑1 (MQ=255)                                                                         ataGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAAcc                                                               >  2:994760/1‑139 (MQ=255)                                                                         ataGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑GGGGGGGAATAAAGTCGCt                                                                                                                      <  2:856216/83‑1 (MQ=255)                                                                         ataGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑GGGGGGGAATAAAGTCGCt                                                                                                                      >  1:856216/1‑83 (MQ=255)                                                                              cTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAg                                                                               >  2:759259/1‑117 (MQ=255)                                                                              cTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAg                                                                               <  1:759259/117‑1 (MQ=255)                                                                                ggATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCg                                                       >  1:66614/1‑139 (MQ=255)                                                                                 gATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCATTCCATCGCTTCTGC‑g                                                                                                                                        >  1:1259333/1‑57 (MQ=255)                                                                                 gATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCATTCCATCGCTTCTGC‑g                                                                                                                                        <  2:1259333/57‑1 (MQ=255)                                                                                 gATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCa                                                                                  >  2:1129054/1‑111 (MQ=255)                                                                                 gATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCa                                                                                  <  1:1129054/111‑1 (MQ=255)                                                                                 gATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCt                                                                        >  2:1244863/1‑121 (MQ=255)                                                                                 gATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCt                                                                        <  1:1244863/121‑1 (MQ=255)                                                                                 gATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCgg                                                      >  2:1013757/1‑139 (MQ=255)                                                                                 gATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCgg                                                      >  1:1114215/1‑139 (MQ=255)                                                                                       cTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTg                                                <  1:1013757/139‑1 (MQ=255)                                                                                         gTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGa                                                                              >  1:1275702/1‑107 (MQ=255)                                                                                         gTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGa                                                                              <  2:1275702/107‑1 (MQ=255)                                                                                               cacCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑GGGGGGGAATAAAGTCGCt                                                                                                                      >  1:1306185/1‑61 (MQ=255)                                                                                               cacCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑GGGGGGGAATAAAGTCGCt                                                                                                                      <  2:1306185/61‑1 (MQ=255)                                                                                                    gcgcTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTGCCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATg                                           <  1:1191339/132‑1 (MQ=255)                                                                                                    gcgcTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATg                                           >  2:1191339/1‑132 (MQ=255)                                                                                                    gcgcTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGtt                                    >  2:468917/1‑139 (MQ=255)                                                                                                    gcgcTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGtt                                    >  2:236041/1‑139 (MQ=255)                                                                                                    gcgcTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑gggg                                                                                                                                     >  1:558568/1‑41 (MQ=255)                                                                                                    gcgcTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑gggg                                                                                                                                     <  2:558568/41‑1 (MQ=255)                                                                                                     cgcTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTa                                   <  1:994760/139‑1 (MQ=255)                                                                                                         tATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGgc                               <  1:1184006/139‑1 (MQ=255)                                                                                                               tctcACGGGCACTCCATCGCTTCTGC‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTc                                                             >  1:435038/1‑102 (MQ=255)                                                                                                               tctcACGGGCACTCCATCGCTTCTGC‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTc                                                             <  2:435038/102‑1 (MQ=255)                                                                                                                  cACGGGCACTCCATCGCTTCTGC‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAAc                                                                >  1:578509/1‑96 (MQ=255)                                                                                                                  cACGGGCACTCCATCGCTTCTGC‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAAc                                                                <  2:578509/96‑1 (MQ=255)                                                                                                                    cGGGCACTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGt                    >  1:24826/1‑139 (MQ=255)                                                                                                                          cTCCATCGCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTAtt              <  2:975682/139‑1 (MQ=255)                                                                                                                          cTCCATCGCTTCTGC‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGc                                                     >  1:188132/1‑99 (MQ=255)                                                                                                                          cTCCATCGCTTCTGC‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGc                                                     <  2:188132/99‑1 (MQ=255)                                                                                                                           tCCATCGCTTCTGC‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTg            >  1:35929/1‑139 (MQ=255)                                                                                                                           tCCATCGCTTCTGC‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTg            >  1:37438/1‑139 (MQ=255)                                                                                                                                cGCTTCTGC‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCGtt       <  2:35929/139‑1 (MQ=255)                                                                                                                                 gCTTCTGCGGGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCGtt       <  1:236041/139‑1 (MQ=255)                                                                                                                                    tCTGC‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCGTTGAGt   <  2:37438/139‑1 (MQ=255)                                                                                                                                     cTGC‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCg                                              >  2:1358657/1‑95 (MQ=255)                                                                                                                                     cTGC‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCg                                              <  1:1358657/95‑1 (MQ=255)                                                                                                                                     cTGC‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCGTTGAGTc  >  1:1285507/1‑139 (MQ=255)                                                                                                                                     cTGC‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCGTTGAGTc  <  2:990376/139‑1 (MQ=255)                                                                                                                                     cTGC‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCGTTGAGTc  >  1:990376/1‑139 (MQ=255)                                                                                                                                         |                                                                                                                                        CAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGC‑GGGGGGGAATAAAGTCGCTTACCTTTATCCTCGTTGTGCCTACGCTTACAGTTCAAAGACCTGCTATACCAACCTCCCCTCGGCTGGTGCGATGCGTGGTTATGGCGCGCCACAAGTCGTATTTGCCGTTGAGTC  >  NC_000913/3001162‑3001432

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