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.3%	53.2 / 19.2	41	T274P (ACC→CCC)	xdhA	xanthine dehydrogenase, molybdenum binding subunit
Reads supporting (aligned to +/- strand):  ref base A (22/0);  new base C (0/19);  total (22/19)
Fisher's exact test for biased strand distribution p-value = 4.09e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATGCGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGT‑AATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGA  >  NC_000913/3001015‑3001305                                                                                                                                                          |                                                                                                                                              cATGCGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGT‑AATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAGCCGTGAA‑GAGTGTTTCCTCGCaa                                                                                                                                               <  2:100156/150‑1 (MQ=255)  aTGCGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGT‑AATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAGCCGTGAA‑GAGTGTTTCCTCGCAAcc                                                                                                                                             >  1:908865/1‑151 (MQ=255)           aGTCATCAAACCATTTGTCGGTGGCGGTTTTGGT‑AATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGc                                                                                                                                    >  2:195915/1‑151 (MQ=255)                     ccATTTGTCGGTGGCGGTTTTGGT‑AATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCtttt                                                                                                                           >  1:458136/1‑150 (MQ=255)                            tCGGTGGCGGTTTTGGT‑AATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCAtt                                                                                                                     >  1:684501/1‑149 (MQ=255)                               gggggtgTTTTGGTAAAAAAACAGGATGTACGGGAAGAGCAAGTGGCGCCTTCCCTGCCCAGAAAGTTGGGCGCCTTCCCGGTAAAAGTTC‑CCC‑TTATCCGTAAA‑GAGTTTTCCCCCCCACCCCTTCCCCGCCACGCTTTTACCATTg                                                                                                                    <  1:326595/142‑1 (MQ=255)                                      ttttGAA‑AAAAAACAGGATTTATGGGAAGAGCAAAGGGCGGCTTTCTTGCCAAGAAGGTTTGGC‑GTTTTCCGGGTAAATTTTCCCC‑TCAGCGGTAAA‑GGGTTTTCCCTCCCACCCCGCACCCGCCACGCTTTTACCATTGACGGGCaaa                                                                                                           <  2:564590/149‑1 (MQ=255)                                              aataaACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAGCCGTGAA‑GAGGGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTg                                                                                                   >  1:848168/1‑150 (MQ=255)                                              aaaaaaCAGGATGTACGGGAAGAGCAAAGGGCGGCATCCCTGCCCAGAAAGCTTGGCGGCATCCGGGTGAAATTTC‑CCC‑TTAGCGGTGAA‑GAGTTTTCCCCCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGa                                                                                                  <  2:100846/148‑1 (MQ=255)                                                 aagcAGGATGTAATGGAGGAGCAAATGCCGGCTTTCCTAACCAGAAAGTTTGGCGCCATTCGGGGAAAATTTC‑CCC‑TGAGCGGGAAG‑GTGTTTTCCCCCGCACCCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAAcc                                                                                               <  2:669758/148‑1 (MQ=255)                                                      ggATGTAGGGGAGGACCCAAGGGCGGCTTCCCTGCCCAGAAAGCTGGGCGCCTTTCGGGT‑AAAGTTT‑CCCTTTAGCCGAAAA‑GGGTGTTCCCTCGCACCCCCTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGa                                                                                           <  1:62332/150‑1 (MQ=255)                                                      ggATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGAc                                                                                          >  1:630635/1‑151 (MQ=255)                                                           tccgggAAGACCCATTGGCGCCTTTCCTGCCCAGCAAGCTTGGCGCCTTCCCGGTAAAATTTT‑CCC‑TTAGCGGTGAA‑GTTTTTTCCCCCCCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGa                                                                                       <  2:908880/145‑1 (MQ=255)                                                             cTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGAGAAAGTTT‑CCC‑TTAGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGaa                                                                                      >  1:417795/1‑148 (MQ=255)                                                              tGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACAtt                                                                                  >  2:32699/1‑151 (MQ=255)                                                                gAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTg                                                                                 >  1:648239/1‑150 (MQ=255)                                                                  agagCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGaa                                                                               >  2:561951/1‑150 (MQ=255)                                                                    agCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAgg                                                                            >  1:351580/1‑151 (MQ=255)                                                                      ccatgggcggcATCCCGCACCAACAGGCTTGGCGGCATTCGGGTAAATTTTC‑CCC‑TTAGCGGGAAG‑GGGTGTTCCCTCGCACCCCGCACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGtt                                                                          <  1:27803/146‑1 (MQ=255)                                                                        aaTGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTtat                                                                        >  2:298856/1‑151 (MQ=255)                                                                                   tCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGt                                                             >  1:633104/1‑151 (MQ=255)                                                                                         ccAGCAGGCTGGGCGGCTTCCCGGTAAAAGTTT‑CCC‑TTAGCGGTGAG‑GAGTGTTCCCCCCCACCCCGCCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTg                                                         <  1:198967/149‑1 (MQ=255)                                                                                         ccAGAAGGCTGGGGGGCCTCCCGGT‑AAATTTT‑CCC‑CTTGCCCGGAAGGGGTTTTCCCTCCCACCCCTCCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTg                                                         <  1:536853/149‑1 (MQ=255)                                                                                             cAAGTTGGGCGGCTTTCGGGTGAAATTTC‑CCC‑TTGGCGGTGAA‑GAGTTTTCCCCCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTa                                                     <  2:684517/149‑1 (MQ=255)                                                                                              atgtttGGCGGCATCCGGGTAAATTTTT‑CCC‑TTAGCGGTGAG‑GTTTTTTCCCCCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTa                                                     <  2:624504/144‑1 (MQ=255)                                                                                              aaGCTTGGCGGCTTTCGGGAAAAGTTTT‑CCT‑TTAGCCGTGAA‑GTGTTTTCCCTCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTa                                                     <  1:767512/148‑1 (MQ=255)                                                                                                      ctgctttgcgttaaaaTTTT‑CCT‑TTAGCCGAAAA‑GGTTTTTTCCTCGCACCCCGTTCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACAcc                                                <  1:452666/132‑1 (MQ=255)                                                                                                         cATTCCGGTGAAAGTTT‑CCC‑TTAGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTAt                                       >  1:866899/1‑151 (MQ=255)                                                                                                          aTTCCGGTGAAAGTTT‑CCC‑TTAGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATg                                      >  2:317728/1‑151 (MQ=255)                                                                                                           ttCGGGTAAAAGTTC‑CCC‑TTAGCCGAAAA‑GGGTTTTCCCCCGCACCCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTAt                                       <  2:474263/149‑1 (MQ=255)                                                                                                           ttCCGGTGAAAGTTT‑CCC‑TTAGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGc                                     >  1:51966/1‑151 (MQ=255)                                                                                                             ccGGTGAAATTTT‑CCT‑TTGGCCGTAAG‑GATTGTTTCCTCCCCCCCCGCCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCAt                                   <  1:161807/151‑1 (MQ=255)                                                                                                             ccGGTGAAAGTTT‑CCC‑TTAGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCAt                                   >  2:332766/1‑151 (MQ=255)                                                                                                                 tGAAAGTTT‑CCC‑TTAGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCAtctc                                >  1:181400/1‑150 (MQ=255)                                                                                                                  gAAAGTTT‑CCC‑TTAGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCAc                              >  2:65758/1‑151 (MQ=255)                                                                                                                        ttcccc‑TTAGCCGTGAA‑GAGTTTTCCCCCGCACCCCTTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGc                          <  2:144977/147‑1 (MQ=255)                                                                                                                                      tGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCt            >  1:637202/1‑151 (MQ=255)                                                                                                                                           gcgtctccccTCACACCCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATc                  <  1:132301/133‑1 (MQ=255)                                                                                                                                              ttttcccccGCACCCCGACCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGc          <  2:74097/138‑1 (MQ=255)                                                                                                                                              tgtTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCggggg     >  2:777381/1‑151 (MQ=255)                                                                                                                                                ttcccccGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCgggggg    <  1:298851/144‑1 (MQ=255)                                                                                                                                                  tCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGa  >  2:763045/1‑150 (MQ=255)                                                                                                                                                          |                                                                                                                                              CATGCGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGT‑AATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAGCCGTGAA‑GAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGA  >  NC_000913/3001015‑3001305

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: ‑