Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,163,371	0	C	T	45.5%	51.7 / 26.3	44	L38F (CTC→TTC)	thiK	thiamine kinase
Reads supporting (aligned to +/- strand):  ref base C (1/23);  new base T (20/0);  total (21/23)
Fisher's exact test for biased strand distribution p-value = 1.19e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GA‑GTTTTCTCATTGAACATC‑AGGGCCAGC‑GTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACAACTACCCGCATGCATTGCACCGAAGCCGCATTTATATCTCCGTG  >  NC_000913/1163230‑1163512                                                                                                                                                 |                                                                                                                                                aGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GG‑GGTTTTTCtttt                                                                                                                                           >  1:864339/1‑147 (MQ=255)     gggCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GG‑GGTTTTTCCTTgt                                                                                                                                          >  1:829485/1‑144 (MQ=255)      ggCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCTGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GG‑GTGTTTTCTTTTAACATC‑a                                                                                                                                   >  1:446533/1‑150 (MQ=255)      ggCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GA‑GTTTTCTCATTGAACATC‑Ag                                                                                                                                  <  1:716504/151‑1 (MQ=255)       gCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GA‑GTTTTCTCATTGAACATC‑Agg                                                                                                                                 <  1:336619/151‑1 (MQ=255)       gCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GA‑GTTTTCTCATTGAACATC‑Agg                                                                                                                                 <  1:770037/151‑1 (MQ=255)           aaTTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GA‑GTTTTCTCATTGAACATC‑AGGGCCa                                                                                                                             <  1:956690/151‑1 (MQ=255)             ttATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GA‑GGTTTTTTCTTTAACATC‑AAGGc                                                                                                                               >  1:266889/1‑147 (MQ=255)              tATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GA‑GTTTTCTCATTGAACATC‑AGGGCCAGc                                                                                                                           <  2:82535/150‑1 (MQ=255)                   ggTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GA‑GTTTTCTCACTGAACATC‑AGGGCCAGC‑Gtttt                                                                                                                     <  2:848156/150‑1 (MQ=255)                      cAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GG‑GGTTTTTCTTTTAAAAAC‑CAGGCCCGC‑GGtttttt                                                                                                                  >  1:168607/1‑150 (MQ=255)                        ggTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GA‑GTTTTCTCATTGAACATC‑AGGGCCAGC‑GTTTTGTtgtg                                                                                                               <  1:124506/151‑1 (MQ=255)                           aaaGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GA‑GTTTTCTCATTGAACATC‑AGGGCCAGC‑GTTTTGTTGTGCGt                                                                                                            <  2:495583/151‑1 (MQ=255)                             aGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GA‑GTTTTCTCATTGAACATC‑AGGGCCAGC‑GTTTTGTTGTGCGTCa                                                                                                          <  2:298251/151‑1 (MQ=255)                                 gCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GG‑GGTTTTTTTTTTTAAATC‑AAGGCCCGCGGTTttgtttgcggtcgcc                                                                                                        >  2:594544/1‑140 (MQ=255)                                 gCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GAGGTTTTTTCTTTTAAAAAA‑AAGCCCA‑C‑GTTTTTTTGTGGGTCTGcc                                                                                                       >  2:732357/1‑150 (MQ=255)                                   cGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGCCACGACGTTTAATAGTGGGCTTAGTGGC‑GGCGG‑GTTTTTTCATTAACCATT‑TGGGCCAAC‑‑CTTtgttgttggggc                                                                                                          >  1:533198/1‑139 (MQ=255)                                     tttCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GG‑GGTTTTTCTTTTTAAAAC‑AGGGGCCGC‑GGTTTTTtttg                                                                                                               >  1:848459/1‑138 (MQ=255)                                          cagcaATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GA‑GTTTTCTCATTGAACATC‑AGGGCCAGC‑GTTTTGTTGTGCGTCAGCCGCACGATCCt                                                                                             <  1:347384/151‑1 (MQ=255)                                           agcaATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCTGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GA‑GGTTTTTTCTTTTACATC‑AGAGGCAGC‑GGGTTTTTTTGTGTCAGCCGCAccggcct                                                                                             >  2:26914/1‑144 (MQ=255)                                               ataatCTCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GA‑GTTTTCTCATTGAACATC‑AGGGCCAGC‑GTTTTGTTGTGCGTCAGCCGCACGATCCTGATgc                                                                                        <  1:102998/151‑1 (MQ=255)                                                    cccATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GA‑GTTTTCTCATTGAACATC‑AGGGCCAGC‑GTTTTGTTGTGCGTCAGCCGCACGATCCTGATGcgccg                                                                                    <  2:117868/150‑1 (MQ=255)                                                     ccATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GG‑GGTTTTTCTTTGTACATC‑ACGGCCCGC‑GGTTTTTTTTTTGTGAGGCGCACGCACCCGATGcgccgc                                                                                   >  2:568952/1‑150 (MQ=255)                                                           cgcgcgACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GA‑GTTTTCTCATTGAACATC‑AGGGCCAGC‑GTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCcgc                                                                            <  2:50978/151‑1 (MQ=255)                                                             cgcgACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGCGGG‑GG‑GTTTTTTTTTTTAAATCC‑AGGGCC‑GC‑GGTTTTTTTTGGGTCAGGCGCCAGGCCCTCATGCGCCCCCGTCcg                                                                             >  1:543269/1‑148 (MQ=255)                                                                 acgaATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GA‑GTTTTCTCATTGAACATC‑AGGGCCAGC‑GTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCtt                                                                      <  1:544482/151‑1 (MQ=255)                                                                 acgaATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GA‑GTTTTCTCATTGAACATC‑AGGGCCAGC‑GTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCtt                                                                      <  2:232063/151‑1 (MQ=255)                                                                  cgaATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GA‑GTTTTCTCATTGAACATC‑AGGGCCAGC‑GTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTg                                                                     <  2:414537/151‑1 (MQ=255)                                                                   gaATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GG‑GGTTTTTCTTTGAAAAAC‑AAGGCCACC‑GGTTTTTTGTGCGGCAGCCGCACGAACCCGATGCGCCGCAGGCCCCGGTCtt                                                                      >  2:982938/1‑149 (MQ=255)                                                                    aaTTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GG‑GGTTTTTCTTTTAACATA‑AATGACAGC‑GTTTTGTTTTGGGGCAGCCGCCCGATCCCTGTGCGCCGCAGGcccccgt                                                                         >  1:175091/1‑140 (MQ=255)                                                                              gcgcTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GA‑GTTTTCTCATTGAACATC‑AGGGCCAGC‑GTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCgg                                                         <  2:795316/151‑1 (MQ=255)                                                                                gcTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GG‑GTTTTTTTCTTTAACATC‑AGGGGCCGG‑GGTTTTTTTTGCGTCAGCCGCCCCATCCTGCTGCGCCGCCGCCCGCGGTCTTGCGCCCGTATCggg                                                        >  1:394599/1‑150 (MQ=255)                                                                                gcTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GG‑GGTTTTTCCTTTTAAAAC‑AAGGGCAGC‑GGTTTTTTTTGGGTCAGCCGCACCAACCCGATGCGCCGCAGGGCCCGGTCTTTCGCCCGTATCggg                                                        >  2:727652/1‑150 (MQ=255)                                                                                 cTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GC‑GCTTTCTCTTTCTACCTC‑CAGGCCCCC‑GGTTTTTTGTGGGGCAGCCGCCCCCACCTAAGGCGCCCCCGCCCCCGGTCTTGGGCCCGGCTCGGGc                                                       >  1:750554/1‑150 (MQ=255)                                                                                           cAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GG‑GGTTTTTTTTTTAAAATCAAGGGCCGCG‑GTTTTTTTTTGGGTC‑GCCCCACCAACCTGTTTCGCCCCAGGCCCCCGTCTTTCGCCCGTGTCGGGGttttttacca                                             >  2:860552/1‑143 (MQ=255)                                                                                            aGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GA‑GTTTTCTCATTGAACATC‑AGGGCCAGC‑GTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACAACt                                           <  2:624071/151‑1 (MQ=255)                                                                                                     ccGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GA‑GTTTTCTCATTGAACATC‑AGGGCCAGC‑GTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACAACTACCCgcatg                                  <  2:477749/151‑1 (MQ=255)                                                                                                       gTCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GA‑GTTTTCTCATTGAACATC‑AGGGCCAGC‑GTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACAACTACCCgcatgca                                <  2:618371/151‑1 (MQ=255)                                                                                                        tCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GG‑GGTTTTTCTTTTAAAAAC‑CAGGCCAGCGGTTTTTTTTTGGGT‑AGCCGCACGGACCTGATGCGCCCCCGGCCCCGGTCTTTCGCCCGCATCCGGGTTTTTTACAACTAaccccccgc                                 >  2:660419/1‑143 (MQ=255)                                                                                                          acgacgTTTAATAGTGGGCTTAGTGGC‑GG‑GG‑GGTTTTTTTTTTTAAAAC‑AAGGCCAGC‑GGTTTTTTTTGGGTCAGCCCCACCATCCCGATTCCCCCCCGGCCGCCGTCTTTCGCCAGGATCCGGCCTTTTTACAAACACCCCCCAGCAtt                              >  2:971018/1‑150 (MQ=255)                                                                                                            gacgTTTAATAGTGGGCTTAGTGGC‑GG‑GG‑GGTTTTTCATTTAAAATC‑AAGGCCCGC‑GGTTTTTTTTTCGTCAGCCCCACGATCCCGATGCGCCCCAGTCCGCGGTCTTTCGCCCGTATCGGGGGTTTTTACAACTACCCCCCTGCATTGcc                           >  2:978943/1‑150 (MQ=255)                                                                                                                              tAGTGGC‑GG‑GA‑GTTTTCTCATTGAACATC‑AGGGCCAGC‑GTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACAACTACCCGCATGCATTGCACCGAAGCCGCATTtata          <  2:36502/150‑1 (MQ=255)                                                                                                                               aGTGGC‑GG‑GA‑GTTTTCTCATTGAACATC‑AGGGCCAGC‑GTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACAACTACCCGCATGCATTGCACCGAAGCCGCATTTatat         <  2:480179/150‑1 (MQ=255)                                                                                                                                      gg‑gA‑GTTTTCTCATTGAACATC‑AGGGCCAGC‑GTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACAACTACCCGCATGCATTGCACCGAAGCCGCATTTATATCTCCGTg  <  2:788214/151‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                AGACGGGCGAAATTATCTGGTCAGGTAAAGGTGCCGTTTCGCAGCAATAATCCCATCACGCGCGACGAATTGCTGTCGCGCTTTTTCCCGCAGTATCATCCCGTCACGACGTTTAATAGTGGGCTTAGTGGC‑GG‑GA‑GTTTTCTCATTGAACATC‑AGGGCCAGC‑GTTTTGTTGTGCGTCAGCCGCACGATCCTGATGCGCCGCAGTCCGCGTTCTTGCGCCAGTATCGGGCTTTATCACAACTACCCGCATGCATTGCACCGAAGCCGCATTTATATCTCCGTG  >  NC_000913/1163230‑1163512

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