Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	646,747	0	T	G	42.5%	57.1 / 13.2	40	T255P (ACC→CCC)	citG	2‑(5''‑triphosphoribosyl)‑3'‑dephosphocoenzyme‑A synthase
Reads supporting (aligned to +/- strand):  ref base T (2/21);  new base G (17/0);  total (19/21)
Fisher's exact test for biased strand distribution p-value = 1.93e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑GC‑CCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTA  >  NC_000913/646598‑646893                                                                                                                                                       |                                                                                                                                                      tAATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGt                                                                                                                                                       <  2:323216/150‑1 (MQ=255)  aaTATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGttt                                                                                                                                                     >  2:760259/1‑150 (MQ=255)            aaTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑Gc‑cccc                                                                                                                                        <  1:862191/151‑1 (MQ=255)            aaTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑GAAT‑Gc‑cccc                                                                                                                                        >  2:732440/1‑151 (MQ=255)                 aTCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑GC‑CCCCttttt                                                                                                                                   <  1:825083/151‑1 (MQ=255)                    aaGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAAT‑GC‑CCCCttttt                                                                                                                                   >  2:599018/1‑148 (MQ=255)                       tGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAAT‑GC‑CCCCcttttttca                                                                                                                               >  1:410676/1‑145 (MQ=255)                           tAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑GC‑CCCCTTTTTGCaataat                                                                                                                           <  2:776587/149‑1 (MQ=255)                           tAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTGT‑GAAT‑GA‑CCCCCCTTTTCactaattt                                                                                                                         >  1:797662/1‑149 (MQ=255)                             aataatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑GC‑CCCCTTTTTGCAATAATGttt                                                                                                                       <  2:124757/151‑1 (MQ=255)                                  ttAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAATGGC‑CCCCTTTTTGAAAAAATtttttgcc                                                                                                                   >  1:142904/1‑147 (MQ=255)                                   tAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTA‑GGAT‑GC‑CCCCCTTTTTCAATAAAGTTTGCGcc                                                                                                                  >  1:561293/1‑150 (MQ=255)                                    aaaTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑GC‑CCCCTTTTTGCAATAATGTTTGCGCCt                                                                                                                 <  2:798556/150‑1 (MQ=255)                                       tCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGCC‑CAAT‑GC‑CCCC‑CTTTGCAATAAACGTTGCGCCTc                                                                                                                >  2:69800/1‑148 (MQ=255)                                        cTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAAT‑GC‑CCCCCTTTTTCAATAATGTTTTTGCCCcgc                                                                                                              >  1:364653/1‑149 (MQ=255)                                          gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑GC‑CCCCTTTTTGCAATAATGTTTGCGCCTcgcg                                                                                                             <  1:166385/148‑1 (MQ=255)                                          gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑GC‑CCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTg                                                                                                          <  1:221663/151‑1 (MQ=255)                                             cTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTG‑GCAT‑GC‑CCCCCTTTTTCCATAATGTTTTTGCCTcgcg                                                                                                             >  1:641496/1‑145 (MQ=255)                                                   aCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑GC‑CCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAgcg                                                                                                 <  1:189358/151‑1 (MQ=255)                                                     cAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GGAT‑GC‑CCCCCTTTTTCAAAAATTTTTTTGCCTCGCGCCGTAggc                                                                                                     >  2:254134/1‑143 (MQ=255)                                                      aGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GGAT‑TCGCCCCTTTTTTC‑ATAATTGTTTCGCCCCgcgccggaggc                                                                                                     >  1:552916/1‑137 (MQ=255)                                                       ggTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑GGAT‑GC‑CCCCCTTTTTTAATAAAGTTTTCGCCTCGCGCCGtgggc                                                                                                     >  1:252831/1‑139 (MQ=255)                                                                     ggTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCGGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCGGAGG‑GC‑CCCCTTTTTTATTCAAAGTTTTGTCCT‑GCGGTGGAGCCAGAGCAGGcccccccctcc                                                                                  >  1:38645/1‑143 (MQ=255)                                                                         aGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GTAT‑GC‑CCCCCTTTTGTAATAATGTTTTCGCCTCGCGCTGTAGCCCGCGCAGGCCCCCCTCCCCACcccgc                                                                           >  1:132827/1‑151 (MQ=255)                                                                            aCTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGCGGTTC‑GAAT‑GC‑CCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGAt                                                                        <  1:427747/151‑1 (MQ=255)                                                                                 cgccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAAT‑GC‑CCCCCTTTTTTAATAATTTTTTCGCCCCCCGCTGTAGGCCGCGCAGGCCCCCCCCCCCAcccc                                                                             >  2:98236/1‑139 (MQ=255)                                                                                  gccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GGAT‑GC‑CCCCCTTTTTCAATAATGTTTTTGCCTCGCGCTGTAGCCAGCGCCGGCCCCCCCCCCCACCCCGCGATGCaa                                                                    >  2:702957/1‑149 (MQ=255)                                                                                           gagaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑GC‑CCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATc                                                          <  2:804065/150‑1 (MQ=255)                                                                                                     tcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑GC‑CCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATc                                                <  2:538830/150‑1 (MQ=255)                                                                                                          acacTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑GC‑CCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCAt                                           <  1:279783/150‑1 (MQ=255)                                                                                                            acTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑GC‑CCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAg                                        <  2:319809/151‑1 (MQ=255)                                                                                                            acTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGGT‑GGAT‑GC‑CCCCCTTTTGTAAAAATATTTTCGCCCCGCGCTGTAGGCAGCGCCGGGCCCCCCCCCCACCCCGCGGGGGCACGTTGGTATCGCCGGTGGTCGCCAt                                           >  1:111257/1‑148 (MQ=255)                                                                                                               cccTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑GC‑CCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGt                                       <  1:858866/149‑1 (MQ=255)                                                                                                                         cTGCCGGAGATAATCGAGATCGGCGGGGGGTTG‑AAAT‑GC‑CCCCTTTTTCCAA‑AATTTTTGCGCCCCGCGCCGTAGCCCGCGCAGGCCCCCCCCCCCCCCCCCCCATGCCACGTTGGTATCGCCGTTGGTCGCCCTCCGGAGGAGGAAgg                             >  1:818477/1‑149 (MQ=255)                                                                                                                              ggAGATAATCGAGATCGGC‑GGGGGGTC‑GGGT‑GC‑CCCCTTTTTTTCCTAATGTTTTGGCCGCGCCGTGGGGGCAGGGGCGGCCCCCCCCCCCACCCCGCGCCGGCACGTTGGTGTCGCCGGTGTTCGGCACCAGTAGGGGCAAg                              >  1:4807/1‑143 (MQ=255)                                                                                                                                 gaTAATCGAGATCGGC‑GGGAGTTC‑GAAT‑GC‑CCCCTTTTTTCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGACACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGc                     <  1:257032/149‑1 (MQ=255)                                                                                                                                       cGAGATCGGC‑GGGGGTTC‑GAAT‑GC‑CCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCa              <  1:784822/150‑1 (MQ=255)                                                                                                                                             cggc‑ggGGGTTC‑GAAT‑GC‑CCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAgg       <  1:353894/151‑1 (MQ=255)                                                                                                                                               gc‑ggGGGTTC‑GAAT‑GC‑CCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGAt     <  2:568440/151‑1 (MQ=255)                                                                                                                                                  gggggTTC‑GAAT‑GC‑CCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCt   <  1:212421/151‑1 (MQ=255)                                                                                                                                                   ggggTTC‑GAAT‑GC‑CCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTa  <  2:195655/151‑1 (MQ=255)                                                                                                                                                   ggggTTC‑GAAT‑GC‑CCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTa  <  1:858603/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                      TAATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAT‑GC‑CCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTA  >  NC_000913/646598‑646893

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