Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,001,164	0	A	C	46.0%	73.5 / 20.5	50	T274P (ACC→CCC)	xdhA	xanthine dehydrogenase, molybdenum binding subunit
Reads supporting (aligned to +/- strand):  ref base A (24/3);  new base C (0/23);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 2.43e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATT‑CCGGTGAAAGTTT‑CCCTTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTC  >  NC_000913/3001022‑3001313                                                                                                                                                 |                                                                                                                                                      aCGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATT‑CCGGTGAAAGTTT‑CCCTTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTAcc                                                                                                                                               <  1:355324/151‑1 (MQ=255)  cGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATT‑CCGGTGAAAGTTT‑CCCTTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTAccc                                                                                                                                              <  2:717451/151‑1 (MQ=255)    aGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATT‑CCGGTGAAAGTTT‑CCCTTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGc                                                                                                                                            >  2:890969/1‑151 (MQ=255)       catcaAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATT‑CCGGTGAAAGTTT‑CCCTTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCAc                                                                                                                                         >  2:302901/1‑151 (MQ=255)       catcaAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATT‑CCGGTGAAAGTTT‑CCCTTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCAc                                                                                                                                         >  2:153085/1‑151 (MQ=255)           aaccaaTTTGTCGGTGGCGGTTTTGGAAAAAAACAGGTTGTACGGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATT‑CCGGTAAAAGTTT‑CCCTTAGCCGTAAAGAGTTTTCCCCCCCAACCCGTCCCCGCCACGCtt                                                                                                                                     <  2:549740/146‑1 (MQ=255)           aaaCCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATT‑CCGGTGAAAGTTT‑CCCTTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCtt                                                                                                                                     >  1:822223/1‑151 (MQ=255)                                tttGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATT‑CCGGTGAAAGTTT‑CCCTTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATgg                                                                                                                >  1:518531/1‑151 (MQ=255)                                tttGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATT‑CCGGTGAAAGTTT‑CCCTTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATgg                                                                                                                >  1:554621/1‑151 (MQ=255)                                 ttGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATT‑CCGGTGAAAGTTT‑CCCTTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATggg                                                                                                               >  1:286347/1‑151 (MQ=255)                                  tGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATT‑CCGGTGAAAGTTT‑CCCTTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGc                                                                                                              >  2:513821/1‑151 (MQ=255)                                     taataaACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATT‑CCGGTGAAAGTTT‑CCCTTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTg                                                                                                           >  2:450399/1‑151 (MQ=255)                                      aaaaaaCAGGATGTACGGGAAGAGCCAAGGGCGGCTTCCCTGCCCAGCAAGCTGGGCGGCATT‑CCGG‑GAAACTTTCCCCTTAGCCGTGAGGGGTGTTCCCCCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGa                                                                                                          <  2:146934/148‑1 (MQ=255)                                       aaaaaCAGGATGTACTGGAAGAGCCAGGGGCGGCTTTCCTGACCAGCAAGCTGGGCGGCATT‑CCGGTAAAAGTTT‑CCCTTGGCCGTGAAGAGTTTTTCCCCGCACCCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGaa                                                                                                         <  1:143867/149‑1 (MQ=255)                                          aaCAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATT‑CCGGTGAAAGTTT‑CCCTTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACcg                                                                                                      >  2:633659/1‑151 (MQ=255)                                           aCAGGATGTACGGAAAGAGCAAATGGCGGCTTCCCTAACCAGAAAGCTGGGGGCCATT‑CCGGTAAGATTTC‑CCCTGGGCCGTGAAGGGTTTTTCCCCCCCCCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACcgc                                                                                                     <  1:427859/151‑1 (MQ=255)                                              ggATGTGCTGAAAGAGCAAATGGCGGCATTCCTGACCAGAAGGCTTGGCGGCATT‑CCGGTAAAATTTC‑CCCTGAGCCGTGAGGGGTGTTTCCTCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACcg                                                                                                      <  2:358026/147‑1 (MQ=255)                                                  gtagggaaaGAGCAAATGGCGCCATCCCGGACCAGAAAGCTGGGCGGCTTC‑CCGGTGAAAGTTC‑CCCTTAGCCGTGAAGGGTTTTCCCCCGCACCCCGCCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGa                                                                                               <  1:505691/143‑1 (MQ=255)                                                    acgggAAGGGCAAGTGGCGGCATTCCGGCCCAGCAGGCTGGGCGGCTTT‑CGGGAGAAAGTTT‑CCTTTCGCCGTAAAGGGTTTTCCCCCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAAc                                                                                             <  2:538763/147‑1 (MQ=255)                                                     gagaaaGACCAAATGGCGGCCTCCCTCACCCGCAAGGTTGGGGTCATC‑CCGGAAAACTTTC‑CCCTTACCGGTGAGGAGTGTTCCCACGCACCCCCCCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACa                                                                                            <  1:192117/146‑1 (MQ=255)                                                     cTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATT‑CCGGTGAAAGTTT‑CCCTTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACAt                                                                                           >  1:327895/1‑151 (MQ=255)                                                         aaGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATT‑CCGGTGAAAGTTT‑CCCTTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGa                                                                                        >  2:653848/1‑150 (MQ=255)                                                                aaTGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATT‑CCGGTGAAAGTTT‑CCCTTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTtat                                                                                >  2:460862/1‑151 (MQ=255)                                                                aaTGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATT‑CCGGTGAAAGTTT‑CCCTTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTtat                                                                                >  1:815831/1‑151 (MQ=255)                                                                   ggcgccATCCCTGACCAGAAATCTGGGCGCCATC‑CCGGTAAAATTTC‑CCCTGAGCCGAAAAGTGTGTTTCCCCCCACCCCGACCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAg                                                                              <  2:490915/150‑1 (MQ=255)                                                                       gcATTCCTGACCAGCAAGCTTGGCGGCATT‑CCGGTGAAAGTTT‑CCCTTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTgg                                                                         >  1:507852/1‑151 (MQ=255)                                                                            tccgccccggcaggttGGGGGCCTT‑CCGGGTAAATTTTCCCCTTAGCGGGGAAGAGTGTTCCCCCCCACCCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCt                                                                           <  2:803372/131‑1 (MQ=255)                                                                            ccggACCAGCAGGTTTGGCGGATTC‑CCGGTAAAATTTC‑CCCTTGGCGGAAAAGAGTGTTCCCCCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGtt                                                                    <  2:869611/148‑1 (MQ=255)                                                                            ccTGCCCAGCAAGCTTGGCGGCATT‑CCGGTGAAATTTC‑CCCTTAGCCGTGAAGATTTTTCCCCCCCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGtt                                                                    <  1:221701/151‑1 (MQ=255)                                                                            ccTGACCAGCAGGCTGGGCGCCTTC‑CGGGTAAAATTTC‑CCCTTGGCCGTAAAGAGTTTTTCCCCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATg                                                                      <  1:621003/149‑1 (MQ=255)                                                                              tGACCAGCAAGCTTGGCGGCATT‑CCGGTGAAAGTTT‑CCCTTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCt                                                                  >  1:851859/1‑151 (MQ=255)                                                                                  ccggaaagttGGCCGCCTTCCCGGTAAAATTTT‑CCCTTGGCGGTGAAGGGTTTTTCCCCACACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTg                                                                 <  2:555514/141‑1 (MQ=255)                                                                                      aaGCTTGGCGGCATC‑CCGGTGAAAGTTC‑CCCTTAGCCGTGAAGAGTTTTCCCCCGCACCCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTaa                                                            <  1:212438/149‑1 (MQ=255)                                                                                          ttGGCGGCATT‑CCGGTGAAAGTTT‑CCCTTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCgg                                                      >  2:134566/1‑151 (MQ=255)                                                                                                   atccgggTAAATGTTC‑CCCTTGGCGGTAAAGTGTTTTCCCTCGCACCCCGCCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTa                                                <  2:148396/144‑1 (MQ=255)                                                                                                    tcccGGTAAAAGTTC‑CCTTTAGCCGGGAAGGGTTTTTCCCCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGgcgc                                                   <  2:324153/143‑1 (MQ=255)                                                                                                            aaaTTTC‑CCCTTAGCCGGGAAGGGTGTTCCCCCCCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCAtctc                                        <  1:653932/148‑1 (MQ=255)                                                                                                                 tt‑CCCTTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCAc                                >  2:23585/1‑151 (MQ=255)                                                                                                                    ccctgggCCGGGAAGTGTTTTCCCCCCCACCCCTTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACt                               <  2:417348/144‑1 (MQ=255)                                                                                                                           ccGTGAAGGTTGTTCCCGCGAACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCt                       <  1:891090/151‑1 (MQ=255)                                                                                                                           ccGTGAAGATTGTTCCCCCGCACCCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATc                          <  1:1965/148‑1 (MQ=255)                                                                                                                            cGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCtt                      >  2:316521/1‑151 (MQ=255)                                                                                                                                aaGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGc                  >  1:259996/1‑151 (MQ=255)                                                                                                                                 agagTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCg                 >  2:332507/1‑151 (MQ=255)                                                                                                                                   agTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCggg               >  1:809835/1‑151 (MQ=255)                                                                                                                                     ttttcccccGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCg                 <  1:176085/139‑1 (MQ=255)                                                                                                                                      gtTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCgggggg            >  1:6861/1‑151 (MQ=255)                                                                                                                                          ccccGCACCCCTTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGa          <  1:467376/146‑1 (MQ=255)                                                                                                                                          ccTCCCACCCCGCCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCgggggg            <  2:199947/147‑1 (MQ=255)                                                                                                                                             cGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATaa      >  1:660742/1‑150 (MQ=255)                                                                                                                                                aaCCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTc  >  1:722550/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                                      ACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATT‑CCGGTGAAAGTTT‑CCCTTAGCCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTC  >  NC_000913/3001022‑3001313

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