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

TGCGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACC‑AGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAG‑CCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCG  >  NC_000913/3001017‑3001314                                                                                                                                                        |                                                                                                                                                       tGCGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACC‑AGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAG‑CCGTGAAGAGTGTTTCCTCGCAAccc                                                                                                                                                     >  2:108876/1‑151 (MQ=255)  gCGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACC‑AGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAG‑CCGTGAAGAGTGTTTCCTCGCAAccc                                                                                                                                                     <  2:185633/150‑1 (MQ=255)      aCGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACC‑AGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAG‑CCGTGAAGAGTGTTTCCTCGCAACCCGTAcc                                                                                                                                                <  2:305129/151‑1 (MQ=255)        gAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACC‑AGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAG‑CCGTGAAGAGTGTTTCCTCGCAACCCGTAccc                                                                                                                                               >  2:531559/1‑150 (MQ=255)         aGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACC‑AGCAAGCTTTGCGGCATTTCGGTGAAAGGTT‑CCC‑TTAG‑CCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGc                                                                                                                                             >  1:965946/1‑151 (MQ=255)               caAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACC‑AGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAG‑CCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGc                                                                                                                                        >  2:63206/1‑150 (MQ=255)                      tttgggggTGGGGTTTTTGAAAAAAAAGGGGATGTGGGGAAAGAGCAAAGGGCGGCATCCCGGACC‑AGAATGCTGGGGGCTATTCGGGT‑AAAGTTT‑CCCTTTAG‑CGGTGAAGAGTGTTTCCCCCCACCCCGTACCCGCCACGCTTTTAcc                                                                                                                                 <  1:751229/144‑1 (MQ=255)                            ggTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACC‑AGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAG‑CCGTGAAGAGTGTTTCCTCGCAAACCGTACCCGCCACGCTTTTACCATTGa                                                                                                                            >  2:493500/1‑149 (MQ=255)                                   gTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACC‑AGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAG‑CCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAAt                                                                                                                   >  2:514330/1‑151 (MQ=255)                                      ttGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACC‑AGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAG‑CCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATggg                                                                                                                >  2:684985/1‑151 (MQ=255)                                        ggTAAAAAAAGGGATGTAGTGGAAGGGCAATTGGCGGCCTTCTTGCCC‑AGAAAGTTCGGCGGTTTTCGGGTAAATTTTT‑CCC‑TTGG‑CCGGGAAGAGTTTTCCCACCCCCCCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGc                                                                                                               <  1:842980/150‑1 (MQ=255)                                           aaaaaaCGGGATGTAGGGAAAGACCCAAGGGC‑GCTTTCCCGCCCAACCAAGCTTGGCGG‑ATTCCGGGGAAAATTTCCCC‑TTAC‑CCGGGAAGGTTTTTCCCCCGCACCCCGTCGCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGa                                                                                                           <  1:204064/148‑1 (MQ=255)                                            ataaACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACC‑AGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAG‑CCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGa                                                                                                           >  2:156373/1‑150 (MQ=255)                                              aaaCAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACC‑AGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAG‑CCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAAcc                                                                                                        >  2:238781/1‑151 (MQ=255)                                                  aGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACC‑AGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAG‑CCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACcgcg                                                                                                     >  1:364623/1‑150 (MQ=255)                                                      tgtaggggAAGAGCAAGGGGCGTCTTCCCTGACC‑AGAAAGCTGGGCGGCTTCCCGGGGAAATTTT‑CCC‑TTAG‑CCGTAAAGAGTTTTCCCCCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACg                                                                                                  <  2:172335/143‑1 (MQ=255)                                                       gTACTGGAAGAGCCAATGGCGGCATTCCTGACC‑AGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAG‑CCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGa                                                                                                >  2:891688/1‑150 (MQ=255)                                                            ggaggaacctggggcggcTTTCCTTACC‑TGCAGGTGGGGCGGCATTCC‑GTCAAAGGTT‑CCCATTAG‑GGGTGCGTGTTGTTTCCCCGCAACCTGGCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGa                                                                                                <  1:748263/133‑1 (MQ=255)                                                            ggAAGAGCCAATGGCGGCATTCCTGACC‑AGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAG‑CCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTg                                                                                          >  2:99232/1‑151 (MQ=255)                                                                  gCCAATGGCGGCATTCCTGACC‑AGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAG‑CCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGt                                                                                    >  2:743817/1‑151 (MQ=255)                                                                      aTGGCGGCATTCCTGACC‑AGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAG‑CCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTtat                                                                                 >  1:149612/1‑150 (MQ=255)                                                                        ggcggcTTCCCTGCCC‑AGAAAGCTGGGCGGCATCCGGGAAAATGTTC‑CCC‑TTGG‑CGGGGAAGGGTTTTCCCCCCCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTtat                                                                                 <  1:415355/148‑1 (MQ=255)                                                                        ggcggcTTCCCTGACC‑AGCAAGCTTGGCGCTTTTCGGTTAAAATTTT‑CCC‑TTGG‑CGGTGAAGAGTTTTCCCCCACACCCCATCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTtata                                                                                <  1:841801/149‑1 (MQ=255)                                                                        ggcgccATTCCTGCCC‑AGCAAGCTGGGGGGCTTTCCGGTAAAAGTTT‑CCC‑TTAG‑CGGTAAGGATTGTTCCCCCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTtata                                                                                <  2:364829/149‑1 (MQ=255)                                                                        ggcgccATTCCTGACC‑AGAAAGTTGGGCGGCTTCCCGGTGAAAGTTT‑CCC‑TTGC‑CGGAGAAGGGTGTTTCCCCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGt                                                                                    <  1:341732/145‑1 (MQ=255)                                                                               ttCCTGCCA‑AGAAAGCTGGGCGGC‑TTTCGGGAAAATTTT‑CCC‑TTTGCCCGGAAAGGGTTTTCCCCCCCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGAt                                                                        <  2:512853/150‑1 (MQ=255)                                                                               tccgtGACC‑AGCAAGGTGGGGGCGTTTCGGGTAAAATTTC‑CCC‑TTAG‑CTGTGAAGTTTTTTCCCCCGCACCCCTTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGa                                                                         <  1:908291/145‑1 (MQ=255)                                                                                   ttaCC‑AGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAG‑CCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCt                                                                   >  2:449380/3‑151 (MQ=255)                                                                                         atcatgCTGGGCGGCATGCCGGAAAAATTTT‑CCC‑TTAG‑CCGTAAAGAGTTTTCCCTCCCACCCCCTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTa                                                              <  2:567774/146‑1 (MQ=255)                                                                                           cAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAG‑CCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTaa                                                             >  1:385793/1‑150 (MQ=255)                                                                                                tggggggatttCCGGGGAATGTTG‑CCC‑TAAA‑CCGTGAAGAGTTTTCCCGACCACCCGGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACAcc                                                         <  1:330881/140‑1 (MQ=255)                                                                                                tgggcgacacccgggTGAATGTTT‑CCC‑TTGG‑CCGTGAGGAGTCTTCCCCCCCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAac                                                            <  2:242202/133‑1 (MQ=255)                                                                                                     ggcATTCCGGAAAAAGTTT‑CCC‑TTAC‑CCGTGAAAGGTTTTTCCCCCCCCCCCGTACCCCCCCCCCTTTTTCCAATGACGGGCAAATTGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCt                                                   <  2:117718/150‑1 (MQ=255)                                                                                                           ccGGTGAAAGTTT‑CCC‑TTAG‑CCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCa                                             >  1:136576/1‑150 (MQ=255)                                                                                                           ccGGTAAAATTTC‑CCC‑TTGG‑CCGGAAAGTGTTTTTCCCCGCACCCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCt                                                   <  1:413148/144‑1 (MQ=255)                                                                                                               tGAAAGTTT‑CCC‑TTAG‑CCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCAtctc                                         >  1:943210/1‑150 (MQ=255)                                                                                                                 aaaTTTC‑CCC‑TGAC‑CCTTGAGGTGTTTTCCCCCCCACCCCTTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCa                                        <  2:411002/149‑1 (MQ=255)                                                                                                                  aatttt‑CCCTTTGG‑CGGTGAAGTGTTTCCCCCCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCAtctc                                         <  1:292758/145‑1 (MQ=255)                                                                                                                               aG‑CCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCg                          >  2:842999/1‑151 (MQ=255)                                                                                                                                g‑CCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGc                         >  1:823450/1‑151 (MQ=255)                                                                                                                                       aaGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGc                   >  1:784316/1‑151 (MQ=255)                                                                                                                                                 ccTCCCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGaa          <  1:536585/150‑1 (MQ=255)                                                                                                                                                        aCCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCg  >  1:889508/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                       TGCGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACC‑AGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAG‑CCGTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCG  >  NC_000913/3001017‑3001314

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