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

CATGCGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TT‑A‑G‑CCGTGAAGAGTGTTTCCTCGCAACCCGT‑A‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCG  >  NC_000913/3001015‑3001314                                                                                                                                                           |                                                                                                                                                          cATGCGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TT‑A‑G‑CCGTGAAGAGTGTTTCCTCGCaa                                                                                                                                                           <  2:329195/150‑1 (MQ=255) cATGCGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TT‑A‑G‑CCGTGAAGAGTGTTTCCTCGCAAc                                                                                                                                                          >  1:469560/1‑151 (MQ=255)    gCGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TT‑A‑G‑CCGTGAAGAGTGTTTCCTCGCAACCCg                                                                                                                                                       >  2:828409/1‑151 (MQ=255)      gTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TT‑A‑G‑CCGTGAAGAGTGTTTCCTCGCAACCCGT‑a                                                                                                                                                    <  1:1026461/151‑1 (MQ=255)            gTCATCAAACCATTTGTCGGTGGCGGTTTTGGAAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCCTTCCGGTGAAAGTTT‑TCC‑TT‑A‑G‑CCCTGAAGAGTGTTTCCTCGCAACCCGT‑A‑CCCGcc                                                                                                                                             >  1:607565/1‑151 (MQ=255)                              gggggggtttttGGAAAAAAACAGGTTGTACGGAAAGAGCCAGGGGCGCCTTCCCGGACCAGCAAGCTTGGCGCCTTCCCGGTAAAATTTT‑CCC‑TT‑G‑C‑CCGGGAAGAGTGTTCCCTCGCACCCCGT‑C‑CCCGCCACGCTTTTACCATTGACg                                                                                                                           <  1:381175/143‑1 (MQ=255)                                  gggttttgGGAAAAAAACGGGTTGTAGGGAAAGAGCCAAGGGCGGCTTCCCTGCCCAGCAGGCTGGGCGG‑ATTCCGGGAAAAATTT‑CCCTTT‑A‑G‑CCGTGAAGGGGTTTCCCGCGCACCCGTT‑A‑CCCGCCACGCTTTTACCATTGACg                                                                                                                           <  2:379613/143‑1 (MQ=255)                                  gCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TT‑A‑G‑CCGTGAAGAGTGTTTCCTCGCAACCCGT‑A‑CCCGCCACGCTTTTACCATTGACGGGc                                                                                                                        >  2:134138/1‑150 (MQ=255)                                          ggTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TT‑A‑G‑CCGTGAAGAGTGTTTCCTCGCAACCCGT‑A‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCg                                                                                                               >  2:451592/1‑151 (MQ=255)                                           gTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TT‑A‑G‑CCGTGAAGAGTGTTTCCTCGCAACCCGT‑A‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑t                                                                                                             >  2:628025/1‑151 (MQ=255)                                            taataaACATGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TT‑A‑G‑CCGTGAAGAGTGTTTCCTCGCAACCCGT‑A‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑Tg                                                                                                            >  1:659201/1‑151 (MQ=255)                                             aaaaaaCGGGATGTAGGGGAAGAGCAAGTGGGGCCTTTCCGGACAAGCAAGTTGGGCGGCATTCGGGTAAAATTTC‑CCT‑TT‑G‑G‑CGGTGAAGGGTGTTCCCCCCCACCCCGT‑C‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑t                                                                                                             <  1:904689/146‑1 (MQ=255)                                              ataaACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TT‑A‑G‑CCGTGAAGAGTGTTTCCTCGCAACCCGT‑A‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGa                                                                                                           >  2:1009490/1‑150 (MQ=255)                                                aaaCGGGATGTAGTGGAGGGGCAAAGGGCGGAATCCCTGCCCAG‑AAGGTTGGGGGCCTTTCGGTTAAAATTTCCCC‑TC‑A‑G‑CCGGAAAGAGT‑TTTCCCCCCCACCCGTCC‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAAcc                                                                                                        <  2:944626/151‑1 (MQ=255)                                                aaaCAGGTTGTACTAAAAGAGCAAAGGGGGGCATTCCAGACAAGGAGGTTGGGCAGCA‑TACGGTAAAGATTT‑CCC‑TT‑A‑GCCCAAGATGTTTTTTTCCCCACACCCCGT‑C‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAAc                                                                                                         <  1:106355/150‑1 (MQ=255)                                                aaaCAGAATGTTGTGAAAGCCCCAAGGGCGCCTTCCCTGACCACAATGCTGGGGGGC‑TTCCGGTGAAAGTTT‑CCC‑CT‑AGG‑CGGGGAGGATTTTTTCCTCCCCCCCCGC‑C‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAAcc                                                                                                        <  2:740673/151‑1 (MQ=255)                                                     ggATGTACGGAAAGAGCAAAGGGCGGCTTCCCTGCCCAGAAAGTTGGGGGGCTTCCGGGTAAAATTTC‑CCC‑TG‑A‑C‑CCGGAAGGAGTGTTCCCCCCCACCCCGT‑C‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACcg                                                                                                       <  1:139603/147‑1 (MQ=255)                                                          tACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TT‑A‑G‑CCGTGAAGAGTGTTTCCTCGCAACCCGT‑A‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAAc                                                                                              >  1:802875/1‑151 (MQ=255)                                                             gggAAGAGCAAAGGGCGCCTTCCCTGCCCAGAAGGCTTGGCGGCTTCCCGGTAAAATTTC‑CCC‑TT‑A‑G‑CGGTGAGGTGTTTTCCCCCGCACCCCGT‑C‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAACAt                                                                                            <  1:348363/149‑1 (MQ=255)                                                                aagaccaaggggcgccTTTCCTGCCCAGCAGGTTGGGCGGCATCCGGGTAAAAGTTT‑CCC‑TT‑G‑G‑CCGTAAAGGGTTTTCCCCCGCACCCCGT‑C‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAACATTGa                                                                                         <  1:126415/140‑1 (MQ=255)                                                                aaGAGCAAAGGGCGGCATCCCTGACCAGCAAGTTTGGCGGCATTCGGGTAAATTTTC‑CCC‑TT‑A‑G‑CGTGGAGGGGTGTTCCCCCGCCCCCCGC‑C‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAACATTg                                                                                          <  2:697490/149‑1 (MQ=255)                                                                  gagcaaggggcggcTTCCCTGCCCAGAAAGC‑TGGCGGACCTCCCGTTAAAGCTT‑CCCGTT‑G‑G‑CGGAGAGGAGTGTTTCCTGGCCCCCGCT‑C‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAACATTgt                                                                                         <  2:463711/140‑2 (MQ=255)                                                                  gagCCAAGGGCGCCATTCCTGACAAGAAAGCTGGGCGGCATTCCGGTAAATGTTC‑CCC‑TT‑A‑G‑CCGAGAAGGTTTTTCCCCCGCCCCCCGT‑C‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAACATTGaaa                                                                                       <  2:967554/150‑1 (MQ=255)                                                                     ccAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCAGTGAAAGTTT‑CCC‑TT‑A‑G‑CCGTGAAGAGTGTTTCCTCGCAACCCGT‑A‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAACATTGAAAGGtt                                                                                   >  2:820360/1‑151 (MQ=255)                                                                       aaTGGCGTCATCCGTGACCAGAAAGTTGGGCGGCATGCGGGTGAAATTTC‑CCC‑TG‑A‑G‑CCGTGAAGTGTTTTCCCCCGCACCCCGT‑C‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAACATTGAAAGGTta                                                                                  <  1:971738/150‑1 (MQ=255)                                                                       aaTGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TT‑A‑G‑CCGTGAAGAGTGTTTCCTCGCAACCCGT‑A‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAACATTGAAAGGTtat                                                                                 >  1:1022201/1‑151 (MQ=255)                                                                       aaTGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TT‑A‑G‑CCGTGAAGAGTGTTTCCTCGCAACCCGT‑A‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAACATTGAAAGGTtat                                                                                 >  1:807941/1‑151 (MQ=255)                                                                         tGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TT‑A‑G‑CCGTGAAGAGTGTTTCCTCGCAACCCGT‑A‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAACATTGAAAGGTtat                                                                                 >  1:495470/1‑149 (MQ=255)                                                                             ggcATCCCTACCCAGCAGGTTGGGCGCCTTTCCGGTAAAATTTC‑CCT‑TT‑G‑G‑CCGAGAAGTGTGTTCCCGCCCAACCCGC‑A‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAACATTGAAAGGTTATAg                                                                               <  2:769526/147‑1 (MQ=255)                                                                                   ccTGACCAGAAAGCTGGGCGCCATCCCGGTGAAATTTC‑CCC‑TT‑G‑C‑CCGGGAAGAGTGTTCCCCCCCACCCCGT‑C‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGt                                                                      <  1:480790/150‑1 (MQ=255)                                                                                             aggttgggcgccATGCGGGTGAAAGTTC‑CCC‑TT‑A‑G‑GCGAAAGGATTTTTCCCTCGCACCCCGT‑C‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCt                                                                   <  2:694897/137‑1 (MQ=255)                                                                                                cTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TT‑A‑G‑CCGTGAAGAGTGTTTCCTCGCAACCCGT‑A‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCg                                                        >  2:937919/1‑151 (MQ=255)                                                                                                          tgccGTTGAAAGTCT‑CCC‑TT‑A‑G‑CCGAGAAGGGTGTTCCCCCCCCCCCTGA‑C‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTa                                                 <  2:628086/146‑1 (MQ=255)                                                                                                           tCCGGTGAAAGTTT‑CCC‑TT‑A‑G‑CCGTGAAGAGTGTTTCCTCGCAACCCGT‑A‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGc                                              >  2:116668/1‑150 (MQ=255)                                                                                                            tcgggaaaaCTTT‑CCC‑TTGA‑G‑CCGGGAAGATTTTTCCCGCCCACCCCGC‑C‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCg                                                        <  2:3344/134‑1 (MQ=255)                                                                                                             ccggtaaaaTTT‑CCC‑TT‑T‑GCCCGTAAAGGGTGTTTCCCCGCACCCCGT‑C‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCAt                                            <  1:165628/146‑1 (MQ=255)                                                                                                                  aatgTTC‑CCT‑TT‑A‑G‑CGGAAAAGAGTTTTCCCCCGCCCCCCGT‑C‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCAtctc                                         <  2:197019/145‑1 (MQ=255)                                                                                                                  aaaTTTT‑CCC‑TT‑A‑G‑CCGGGAAGGGTTTTTCCCCCCCCCCCGG‑T‑CCCCGCACCGTTTTTCCAATGGCCGGGAAAAGGGCGTTGACCCGCG‑CCGAACATTTGAAAGTTTTAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCAc                                       <  1:1023529/150‑1 (MQ=255)                                                                                                                       tt‑CCC‑TT‑A‑G‑CCGTGAAGAGTGTTTCCTCGCAACCCGT‑A‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCAc                                 >  2:386659/1‑151 (MQ=255)                                                                                                                            c‑TT‑A‑G‑CCGTGAAGAGTGTTTCCTCGCAACCCGT‑A‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTcc                              >  2:163414/1‑150 (MQ=255)                                                                                                                                        tGAAGAGTGTTTCCTCGCAACCCGT‑A‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCt                     >  1:692343/1‑151 (MQ=255)                                                                                                                                                 ttcccccGCACCCCGT‑C‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCggggg              <  1:232846/143‑1 (MQ=255)                                                                                                                                                  tccccgcgcaCCCGT‑ACCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAAAACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCgggggg             <  1:387088/141‑1 (MQ=255)                                                                                                                                                    ccccGCACCCCGT‑C‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGaa          <  1:268115/147‑1 (MQ=255)                                                                                                                                                       cGCAACCCGT‑A‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATaa       >  1:371551/1‑150 (MQ=255)                                                                                                                                                           aCCCGT‑A‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCg  >  1:385500/1‑151 (MQ=255)                                                                                                                                                           |                                                                                                                                                          CATGCGTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TT‑A‑G‑CCGTGAAGAGTGTTTCCTCGCAACCCGT‑A‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCG  >  NC_000913/3001015‑3001314

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