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

GTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TT‑A‑G‑CCGTGAAGAGTGTTTCCTCGCAACCCGT‑A‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACC  >  NC_000913/3001020‑3001320                                                                                                                                                             |                                                                                                                                                         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)                                                                                                                                             tcccccccaccccGT‑C‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCg                        <  2:469429/134‑1 (MQ=255)                                                                                                                                               ccccGCACCCCGT‑C‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGaa                <  1:268115/147‑1 (MQ=255)                                                                                                                                                  cGCAACCCGT‑A‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATaa             >  1:371551/1‑150 (MQ=255)                                                                                                                                                    ccccccGT‑A‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAg           <  2:619738/147‑1 (MQ=255)                                                                                                                                                      aCCCGT‑A‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCg        >  1:385500/1‑151 (MQ=255)                                                                                                                                                       cccGT‑A‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCGc       >  1:753346/1‑151 (MQ=255)                                                                                                                                                             a‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCGCTTAcc  >  1:973317/1‑151 (MQ=255)                                                                                                                                                             |                                                                                                                                                         GTACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TT‑A‑G‑CCGTGAAGAGTGTTTCCTCGCAACCCGT‑A‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCG‑TGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAAAGTCGCTTACC  >  NC_000913/3001020‑3001320

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