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

TACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAG‑TTT‑CCC‑TTAGCCGTGAAGAG‑TGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAA  >  NC_000913/3001021‑3001309                                                                                                                                                    |                                                                                                                                                   tACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAG‑TTT‑CCC‑TTAGCCGTGAAGAG‑TGTTTCCTCGCAACCCGTA‑c                                                                                                                                            >  1:894624/1‑151 (MQ=255)  aCGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGTAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAG‑GTT‑CCC‑TTAGCCGTGAAGAG‑TGTTTCCTCGCAACCCGTA‑cc                                                                                                                                           >  2:628526/1‑151 (MQ=255)   cGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAG‑TTT‑CCC‑TTAGCCGTGAAGAG‑TGTTTCCTCGCAACCCGTA‑ccc                                                                                                                                          <  1:73241/151‑1 (MQ=255)    gAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAG‑TTT‑CCC‑TTAGCCGTGAAGAG‑TGTTTCCTCGCAACCCGTA‑ccc                                                                                                                                          >  2:431037/1‑150 (MQ=255)            aaaCCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAG‑TTT‑CCC‑TTAGCCGTGAAGAG‑TGTTTCCTCGCAACCCGTA‑CCCGCCACGCtt                                                                                                                                 >  2:236729/1‑151 (MQ=255)               ccATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAG‑TTT‑CCC‑TTAGCCGGGAAGAG‑TGTTTCCTCGCAACCCGTA‑CCCGCCACGCtttt                                                                                                                               >  2:682323/1‑150 (MQ=255)                  tttGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAG‑TTT‑CCC‑TTAGCCGTGAAGAG‑TGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTAcc                                                                                                                            >  2:711570/1‑150 (MQ=255)                    tGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAG‑TTT‑CCC‑TTAGCCGTGAAGAG‑TGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCAtt                                                                                                                         >  2:256261/1‑151 (MQ=255)                             cGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAG‑TTT‑CCC‑TTAGCCGTGAAGAG‑TGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCa                                                                                                                 >  2:70112/1‑150 (MQ=255)                                       aaaaaaCGGGAGGTAAGGGAAGAGCCAATGGGGGCTTTCCTGACAAGAAGGCTGGGCGGTTTCCCGGTGAAAT‑TTT‑CCC‑TGAGCGGTGAGGGG‑TGTTCCCTCGCACCCCGTC‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGt                                                                                                        <  2:398208/146‑1 (MQ=255)                                            aCAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAG‑TTT‑CCC‑TTAGCCGTGAAGAG‑TGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACcgc                                                                                                 >  1:878815/1‑151 (MQ=255)                                             cAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAG‑TTT‑CCC‑TTAGCCGTGAAGAG‑TGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACcgcg                                                                                                >  2:924774/1‑151 (MQ=255)                                                 aTGTACTGGAAGAGCCAAGGGCGGCTTTCCCGCCCAGAAAGTTTGGCGGCATCCGGGTAAAAT‑TTC‑CCC‑TTCCCGGTAAAGGT‑TTTTCCCTCGCCCCCCTTA‑CCCGCCACGCTTTTACCATTGACGGGCAGATGGGCGTGAACCGCGACg                                                                                             <  2:951926/150‑1 (MQ=255)                                                    tACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAG‑TTT‑CCC‑TTAGCCGTGAAGAG‑TGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAAc                                                                                         >  1:899053/1‑151 (MQ=255)                                                     agtGGAAGAGCCAAGGGCGGCTTCCCTGACAAGCAAGCTGGGCGGCATCCCGGTAAAAG‑TTC‑CCC‑TAAGCCGAGAAGGG‑TTTTCCCCCGCACCCCTTC‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACg                                                                                             <  2:51024/144‑1 (MQ=255)                                                       tGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAG‑TTT‑CCC‑TTAGCCGTGAAGAG‑TGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACAtt                                                                                      >  2:244583/1‑151 (MQ=255)                                                         gAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAG‑TTT‑CCC‑TTAGCCGTGAAGAG‑TGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTg                                                                                     >  1:822711/1‑150 (MQ=255)                                                              gcaagtGGCGCCTTTCCCGACCAGCAGTCTGGGCGCCTTCCCGGAAAAAT‑TTC‑CCC‑TTAGCGGGAAAGAT‑TTTTCCCCCCCCCCCCGCC‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAgg                                                                                <  1:225169/145‑1 (MQ=255)                                                              gCCAAGGGCGGCTTTCCTAACCAGAAAG‑TTGGGGGGCTTCCGGGAAAAA‑TTTCCCC‑TTAGCGGTGAGGAG‑TTTTCCCCCCCACCCCGTC‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAgg                                                                                <  2:14362/150‑1 (MQ=255)                                                                 aaTGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAG‑TTT‑CCC‑TTAGCCGTGAAGAG‑TGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTta                                                                             >  1:329300/1‑150 (MQ=255)                                                                  ttGGCGCCTTCCCTGACCAGCAGGTTGGGCGCCATTCGGGTAAATT‑TTC‑CCC‑TTGGCCGTAAAGTG‑TTTTCCCTCGCACCCCTTC‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTtat                                                                            <  2:855247/149‑1 (MQ=255)                                                                           ttCCTGACCAGCAAGCTGGGCGGCATCCGGGTGAAAG‑TTC‑CCC‑TTAGCGGGAAAGAG‑TTTTCCCTCGCAACCCGTC‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGAt                                                                   <  2:9735/150‑1 (MQ=255)                                                                           tcccTGACCAGCAATGTTGGCGGCATTCGGGTAAATG‑TTT‑CCC‑TTAGCCGTGAAGTG‑TGTTCCCCCCCACCCCGAC‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTg                                                                      <  1:592922/145‑1 (MQ=255)                                                                             ccTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAG‑TTT‑CCC‑TTAGCCGTGAAGAG‑TGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGtt                                                                >  2:992998/1‑151 (MQ=255)                                                                                  ccAGCAAGCTGGGCGCCTTCCCGGT‑AAAT‑TTT‑CCCTTTCGCCGGGAAGGT‑TTTTCCCGCGCACCCCGTC‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCt                                                              <  1:161503/148‑1 (MQ=255)                                                                                   cAGCAAGCTTGGCGGCATTCCGGTGAAAG‑TTT‑CCC‑TTAGCCGTGAAGAG‑TGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCt                                                          >  1:496812/1‑151 (MQ=255)                                                                                    agaaaGCTGGGC‑GCATTTCGGGAAAAT‑TTT‑CCC‑TTTGGCGGTAAGGGTTGTTCCCCCGCACCCGTTC‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTa                                                         <  2:533265/148‑1 (MQ=255)                                                                                       aggCTGGGCGGCTTTCCGGTAAAAT‑TTC‑CCC‑TTAGCCGTAAAGAG‑TTTTCCCCCGCACCCCGTC‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTa                                                         <  2:57774/146‑1 (MQ=255)                                                                                        aggttggggggcttCCGGGTAAAA‑TTT‑CCCTTTGGCGGTGAAGGT‑TGTTCCCCCCCACCCCGTC‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAac                                                       <  2:329429/138‑1 (MQ=255)                                                                                             ggcggcATTCCGGTGAAAG‑TTT‑CCC‑TTAGCCGTGAAGAG‑TGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGgcg                                                >  2:406829/1‑151 (MQ=255)                                                                                             ggcggc‑TTCCGGGTAAAGTTTT‑CCT‑TTAGCCGTAAAGTT‑TTTTCCCCCGCACCCCGTC‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGgcg                                                <  1:711352/151‑1 (MQ=255)                                                                                                ggcATCCGGGTGAAAT‑TTC‑CCG‑TTAGCGGTGAGGGT‑TTTTTCCCCGCACCCCCTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGgc                                                 <  2:64261/147‑1 (MQ=255)                                                                                                         gTGAAAG‑TTT‑CCC‑TTAGCCGTGAAGAG‑TGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCAtctc                                    >  1:23537/1‑151 (MQ=255)                                                                                                          ggtaaattttcccc‑TTAGCGGGGAAGAG‑TGTTCCCCCCCCCCCCGCA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCAtc                                      <  1:70088/139‑1 (MQ=255)                                                                                                            aatg‑TTT‑CCC‑TTCGCCGTGAAGAG‑TTTTCCCCCGCACCCCGTC‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTa                                            <  1:236639/137‑1 (MQ=255)                                                                                                              aatttt‑CCCTTTAGCCGTAAAGAG‑TTTTCCCCCGCACCCCGTC‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCAtctc                                    <  1:992481/145‑1 (MQ=255)                                                                                                                  ttcccc‑TTAGCGGTAAAGGG‑TTTTCCCCCGCACCCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGc                              <  2:903388/147‑1 (MQ=255)                                                                                                                                tGAAGTG‑TGTTCCCCCGCACCCCTTC‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTc                 <  2:535634/150‑1 (MQ=255)                                                                                                                                tGAAGAG‑TGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCt                >  1:678755/1‑151 (MQ=255)                                                                                                                                  aaGAT‑TGTTCCCCCGCCCCCGATA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCAt                       <  1:183615/142‑1 (MQ=255)                                                                                                                                  aaGAG‑TGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGc              >  1:447987/1‑151 (MQ=255)                                                                                                                                    gag‑TGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCgg            >  2:799365/1‑151 (MQ=255)                                                                                                                                     ggtttttcccccGCACCCCGTC‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGAATCTCACGGGCACTCCATCGCTTCTGCgg            <  2:985255/140‑1 (MQ=255)                                                                                                                                               ccGCAACCCATCCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATa   <  1:670742/150‑1 (MQ=255)                                                                                                                                                 gCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATaa  >  2:974623/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                   TACGAGTCATCAAACCATTTGTCGGTGGCGGTTTTGGTAATAAACAGGATGTACTGGAAGAGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAG‑TTT‑CCC‑TTAGCCGTGAAGAG‑TGTTTCCTCGCAACCCGTA‑CCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTCTGTCTAACACCGGCGCTTATGCATCTCACGGGCACTCCATCGCTTCTGCGGGGGGGAATAA  >  NC_000913/3001021‑3001309

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