Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,001,164	0	A	C	40.8%	90.9 / 11.4	49	T274P (ACC→CCC)	xdhA	xanthine dehydrogenase, molybdenum binding subunit
Reads supporting (aligned to +/- strand):  ref base A (28/1);  new base C (0/20);  total (28/21)
Fisher's exact test for biased strand distribution p-value = 7.43e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATGCGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACA‑GGATGTACTGGAAGAGCCAATGGCGGCATTCCTG‑ACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TT‑AGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGG  >  NC_000913/3001015‑3001303                                                                                                                                                           |                                                                                                                                            cATGCGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACA‑GGATGTACTGGAAGAGCCAATGGCGGCATTCCTG‑ACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TT‑AGCCGTGAA‑GAGTGTTTCCTCGCaa                                                                                                                                             <  2:1109108/150‑1 (MQ=255)   tGCGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTTGGAATAAACA‑GGAGGTACTGGAAAAGCCAATGGGGGCATTCCTG‑ACCAGCAAGCTTGGGGGCATTCCGGTGAAAGTTTCCC‑TT‑AGCCGGGAA‑GAGTGTTTCCTCGCAAccc                                                                                                                                          >  2:732141/1‑151 (MQ=255)    gCGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACA‑GGATGTACTGGAAGAGCCAATGGCGGCATTCCTG‑ACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TT‑AGCCGTGAA‑GAGTGTTTCCTCGCAACCCg                                                                                                                                         >  1:987665/1‑151 (MQ=255)          gAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACA‑GGATGTACTGGAAGAGCCACTGGCGGCATTCCTG‑ACCAGCAAGCTTGGCGGCATTCCGGGGAAAGTTTCCC‑TT‑AGCCGTGCA‑GGGTGTTTCCTCGCAACCCGTAccc                                                                                                                                    >  1:544953/1‑150 (MQ=255)                 caAACCATTTGTCGGTGGCGGTTTTGGTAATAAACA‑GGATGTACTGGAAGAGCCAATGGCGGCATTCCTG‑ACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TT‑AGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGc                                                                                                                             <  1:33901/150‑1 (MQ=255)                        tttGTCGGTGGCGGTTTTGGTAATAAACA‑GTATGTACTGGAAGAGCCAATGGCGGCATTCCTG‑ACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TT‑AGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTAcc                                                                                                                      >  2:812519/1‑150 (MQ=255)                            gcgggtggggTTTTGGGAAAAAAAAGGGTTGAAGTGAAAGGGCCAAGGGCGTCCTCCCGG‑CCCAAAAGGCTTGGCGCCTTCCCGGTAAACGTTCCCC‑TT‑GGCCGGAAA‑GGGTTTTTCCACGCACCCCGTCCCCGCCACGCTTTTACCATTg                                                                                                                  <  2:508270/143‑1 (MQ=255)                             cGGTGGCGGTTTTGGTAATAAACA‑GGATGTTCTGGAAGAGCCAATGGCGGCATTCCTG‑ACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TT‑AGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGAc                                                                                                                >  1:11888/1‑151 (MQ=255)                                  ggggTTTTGGAAAAAAAAA‑GGGTTTAGGGGGAGGGGCCATGGGGGCCATCCGGAACCAGAAGGCTGGGC‑GCCTTCCGGGTAAAGTTTCCC‑CT‑TGGCTTAAAGGAGTTTTCCCCCGCACCCCGTCCCCGCCACGCTTTTACCATTGACggg                                                                                                             <  2:51575/148‑1 (MQ=255)                                        ttGGTAATAAACA‑GGATGTACTGGAAGAGCCAATGGCGGCATTCCTG‑ACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TT‑AGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATggg                                                                                                     >  2:884635/1‑151 (MQ=255)                                        ttGGTAAAAAACA‑GGATGTACTGGAAGAGCCAATGGCGGCATTCCTG‑ACCAGCAAGTTGGGCGGCATTCGGGTAAAAGTTTCCC‑TT‑AGCCGTGAA‑GAGTGTTCCCCCGCACCCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATgg                                                                                                      <  2:362135/150‑1 (MQ=255)                                             aataaACA‑GGATGTACTGGAAGAGCCAATGGCGGCATTCCTG‑ACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TT‑AGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTg                                                                                                 >  2:858582/1‑150 (MQ=255)                                               aaaaCG‑GGATGTACGGAAAGGGCAATGGGCGGTATCCCTG‑CCCGGCAAGTTGGGCGTTTTTCCGGTGAATGTTCCCC‑TT‑GGCGAGAAA‑GATTGTTCCCCCCCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAAc                                                                                              <  1:1095188/150‑1 (MQ=255)                                                aaaCA‑GGATGTACTGGAAGAGCCAATGGCGGCATTCCTG‑ACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TT‑AGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAAcc                                                                                             >  1:270837/1‑151 (MQ=255)                                                           tgaggaaaGAGCCATTGCCGTCTTTCCTG‑ACAAGCAAGCTCGGCCGCATCCCGGGAAAATTTTCCC‑TG‑GGCCGGAAA‑GGGTTTTTCCTCCCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACgg                                                                                      <  2:757512/142‑1 (MQ=255)                                                               ggaggggcaaagggcggcATTCCTG‑ACAAGCAAGCTGGGCGCCTTTCGGGAAAATTTTCCCC‑TT‑AGCGGGGAA‑GAGTGTTCCCCCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACAtt                                                                                <  2:1039484/138‑1 (MQ=255)                                                               ggagggcccAGTGGGGGCTTCCCTG‑CCCAGCAGGCTGGGCGGCATCCCGGTGAAATTTTCCC‑TT‑AGCGGAGAA‑GGGTTTTCCCTCGCACCCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACAtt                                                                                <  2:412792/143‑1 (MQ=255)                                                               ggAAGAGCCAATGGCGGCATTCCTG‑ACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TT‑AGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTg                                                                               >  1:602866/1‑151 (MQ=255)                                                                         aTGGCGGCATTCCTG‑ACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TT‑AGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTta                                                                       >  1:526561/1‑149 (MQ=255)                                                                           ggcggcTTTCCTG‑ACCAGCAAGCTGGGCGTCATCCGGGTAAATTTTTCCC‑TG‑AGCGGGAAA‑GTTTTTTCCCGCGCCCCCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAg                                                                    <  2:42116/150‑1 (MQ=255)                                                                           ggcggcATTCCTG‑ACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TT‑AGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGt                                                                   >  1:533197/1‑151 (MQ=255)                                                                           ggcggcATTCCTG‑ACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TT‑AGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGt                                                                   >  1:552860/1‑151 (MQ=255)                                                                           ggcggcATTCCTG‑ACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TT‑AGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGt                                                                   >  1:164549/1‑151 (MQ=255)                                                                              gtcttcccTG‑ACCAGCAAGGTTGGCGCCATCCCGGTAAATTTTCCCC‑TT‑AGCCGTAAA‑GGGTTTTCCCCCCCACCCCGCCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCt                                                                 <  2:836077/144‑1 (MQ=255)                                                                              ggcATTCCTG‑ACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TT‑AGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTg                                                                >  2:483083/1‑151 (MQ=255)                                                                                cATTCCTG‑ACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TT‑AGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTgg                                                               >  2:837905/1‑150 (MQ=255)                                                                                    ccTG‑CCCAGAAAGCTGGGCGGCATCCGGGTAAAATTTTCCC‑TG‑AGCGGTGAA‑GTTTTTTTCCCCCCACCCCTTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATg                                                            <  2:977876/149‑1 (MQ=255)                                                                                    ccTG‑ACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TT‑AGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGtt                                                          >  1:306437/1‑151 (MQ=255)                                                                                    ccTG‑ACCAGCAAGCTTGGCGGCATCCGGGTAAAATTTCCCC‑TT‑AGCCGGGAA‑GAGTTTTTCCCCACACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGAt                                                             <  1:36371/148‑1 (MQ=255)                                                                                    ccTG‑ACCAGAAGGCTGGGCGGCATTCGGGTAAAATTTCCCC‑TT‑AGCGGTGAA‑GGGTTTTTCCCCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATg                                                            <  1:626739/149‑1 (MQ=255)                                                                                      tG‑ACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TT‑AGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCt                                                        >  2:445713/1‑151 (MQ=255)                                                                                      cgacccAGCAGGCTTGGC‑GCCTTCCGGGTAAATTTTCCCTTT‑GGCCGTAAA‑GTTTTTTCCCCCCCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTc                                                         <  1:269473/147‑1 (MQ=255)                                                                                           cAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TT‑AGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCt                                                    >  2:797142/1‑151 (MQ=255)                                                                                               aaGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TT‑AGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAac                                                 >  2:101344/1‑150 (MQ=255)                                                                                                     ggcggcATTCCGGTGAAAGTTTCCC‑TT‑AGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGgcg                                          >  2:306797/1‑151 (MQ=255)                                                                                                     ggcggcA‑TCCGGGTAAATTTTCCC‑TTGAGCGGGGAA‑GGGTTTTCCCCCCCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGgc                                           <  2:953359/150‑1 (MQ=255)                                                                                                        ggcATTCCGGTGAAAGTTTCCC‑TT‑AGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCtt                                       >  2:171459/1‑151 (MQ=255)                                                                                                             tCCGGTGAAAGTTTCCC‑TT‑AGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGc                                   >  2:165873/1‑150 (MQ=255)                                                                                                                 gttaaaGTTTCCCTTT‑AGCGGTGAG‑GTGTTTTTCCTCACACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCa                                  <  2:815904/145‑1 (MQ=255)                                                                                                                  tGAAAGTTTCCC‑TT‑AGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCAtctc                              >  2:314278/1‑150 (MQ=255)                                                                                                                     aaGTTTCCC‑TT‑AGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACgg                          >  2:76885/1‑151 (MQ=255)                                                                                                                       aagttcccatgggCCGTGAA‑GGGTGTTCCCGCCCACCCCCTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCAtctc                              <  2:1075462/135‑1 (MQ=255)                                                                                                                                  aGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCg               >  2:842975/1‑151 (MQ=255)                                                                                                                                    ccGTGAA‑GGGTTTTCCCTCCCACCCCGCACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCAt                 <  2:860922/147‑1 (MQ=255)                                                                                                                                      gagAA‑GAGTTTTCCCCCGCACCCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCtt            <  1:475765/148‑1 (MQ=255)                                                                                                                                      gTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTc           >  1:258756/1‑151 (MQ=255)                                                                                                                                      gTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTc           >  1:808222/1‑151 (MQ=255)                                                                                                                                         aa‑GAGTGTTCCCTCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGc        <  1:466975/151‑1 (MQ=255)                                                                                                                                            gagTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCgg      >  2:358301/1‑151 (MQ=255)                                                                                                                                                 ttcccTCGCACCCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCgggggg  <  1:820065/147‑1 (MQ=255)                                                                                                                                                           |                                                                                                                                            CATGCGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACA‑GGATGTACTGGAAGAGCCAATGGCGGCATTCCTG‑ACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTTCCC‑TT‑AGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGG  >  NC_000913/3001015‑3001303

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