Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	646,747	0	T	G	42.3%	83.4 / 11.8	53	T255P (ACC→CCC)	citG	2‑(5''‑triphosphoribosyl)‑3'‑dephosphocoenzyme‑A synthase
Reads supporting (aligned to +/- strand):  ref base T (4/26);  new base G (22/0);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 1.11e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTA  >  NC_000913/646598‑646893                                                                                                                                                       |                                                                                                                                                     tAATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGt                                                                                                                                                      <  2:671197/150‑1 (MQ=255)         ccAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCG‑GATg                                                                                                                                              >  2:577236/1‑149 (MQ=255)             aTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AATG‑CCCCCt                                                                                                                                       <  1:364255/151‑1 (MQ=255)              tttATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGG‑ATTG‑cccc                                                                                                                                         >  2:668942/1‑149 (MQ=255)              tttATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCG‑AATG‑ccccc                                                                                                                                        >  2:569621/1‑149 (MQ=255)                tATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCG‑GATG‑CCCCCcttt                                                                                                                                    >  2:364486/1‑151 (MQ=255)                tATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCG‑AATG‑CCCCCcttt                                                                                                                                    >  2:225981/1‑151 (MQ=255)                tATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGGCGTAAAG‑CCCCCtt                                                                                                                                      >  2:50026/1‑150 (MQ=255)                   cAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AATG‑CCCCCTTTTTg                                                                                                                                  <  1:665107/150‑1 (MQ=255)                   cAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AATG‑CCCCCTTTTTGc                                                                                                                                 <  1:68566/151‑1 (MQ=255)                    aaGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTTG‑AATG‑CCCCCctttt                                                                                                                                   >  2:750288/1‑148 (MQ=255)                     aGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTTG‑GATGTCCCCCttttt                                                                                                                                   >  2:192266/1‑148 (MQ=255)                       tGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGGCG‑GAAG‑CCCCCctttttgta                                                                                                                               >  2:263149/1‑145 (MQ=255)                           tAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AATG‑CCCCCTTTTTGCAATAATg                                                                                                                          <  1:512547/150‑1 (MQ=255)                             aataatTAAATCTGTGCTACAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGTGGGCGG‑AACG‑ACCCCCTTTTTTATTAATGtt                                                                                                                        >  1:731014/1‑150 (MQ=255)                                     aaTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTcg                                                                                                               <  1:577198/151‑1 (MQ=255)                                         tgtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCG‑GAAG‑CCCCCCTTTTTCAATAAAATTTTCGCCCCgcgc                                                                                                            >  2:578927/1‑150 (MQ=255)                                           tgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTTG‑AATGACCCCCTTTTTTCAAAAAT‑TTTGCGccgcccgccg                                                                                                          >  1:545180/1‑148 (MQ=255)                                            gCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGt                                                                                                         <  1:284905/150‑1 (MQ=255)                                            gCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTa                                                                                                        <  1:662048/151‑1 (MQ=255)                                                  aaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTTG‑AAAG‑CCCCCCTTTTGTTAAAATGTTTTTGCCTTgcgccggtggca                                                                                                    >  2:734184/1‑141 (MQ=255)                                                  aaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCG‑AATG‑CCCCCCTTTTTCTATAATGTTTtcggcgccc                                                                                                              >  1:521877/1‑133 (MQ=255)                                                    ccAGGTAAGGATCAGTATGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGACCGGC‑GGGGGTTCG‑AATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAgcgc                                                                                                <  1:58410/151‑1 (MQ=255)                                                             gATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGAGGGTCG‑AATT‑CCCCCTTTTTTGAATAATGGTTTCGCCTCGCGCTGTAGCCATCGCCGGcccccc                                                                                       >  2:632339/1‑151 (MQ=255)                                                              aTCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCG‑AATG‑CCCCCTTTTTTTAATAATTTTTTTGCCCCCCGCCGGAGACAGCGCACGgccccc                                                                                       >  2:54456/1‑150 (MQ=255)                                                                 aGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGATCG‑ACTG‑TCCCCCTTTTCTGCCAATGGTTTCGCCTCCCGCCGTAGTTAACGCCGGGCCCCCTCGc                                                                                   >  2:344595/1‑151 (MQ=255)                                                                   tAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCG‑GATG‑CCCCCCTTTTTCAATAATTTTTTCGCCTCCCGCTGTAGCCAGCGCAGGccccccccccc                                                                                  >  2:447682/1‑145 (MQ=255)                                                                   tAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCG‑GATG‑CCCCCCTTTTTCAATAATGTTTTCGCCCCCCCCCGTAGCCCGCGCAGGcccccccccccc                                                                                 >  2:252989/1‑145 (MQ=255)                                                                    aGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCAc                                                                                <  2:575094/151‑1 (MQ=255)                                                                               gccgccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGc                                                                      <  1:306708/150‑1 (MQ=255)                                                                               gccgccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCa                                                                     <  2:753991/151‑1 (MQ=255)                                                                                   ccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGTCCCCCCTCGCCGCCGCGCGATGCAACGt                                                                 <  2:394156/151‑1 (MQ=255)                                                                                    cTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTTG‑CATG‑CCCCCTTTTTTCAATTAAATTTTCGCCCCGCGCCGGAGCCAGCGCCAGCCCCCCCCGCCACcccgc                                                                           >  1:856396/1‑140 (MQ=255)                                                                                     tGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTTG‑AATG‑CCCCCCTTTTTTAATAATTTTTTCGCCCCCCGCCGTTGCCAGCGCCGGCCCCCCCCCCCCCCCCCCGGCGCAACgct                                                                >  1:335997/1‑148 (MQ=255)                                                                                           gagaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATc                                                          <  2:621507/150‑1 (MQ=255)                                                                                              aTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCg                                                      <  2:342281/151‑1 (MQ=255)                                                                                                ttcgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGt                                                     <  2:280375/150‑1 (MQ=255)                                                                                                 tcgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGtt                                                    <  2:242449/150‑1 (MQ=255)                                                                                                  cgttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTg                                                   <  2:140253/150‑1 (MQ=255)                                                                                                   gttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGGTG‑CAAG‑GCCCCCCTTTGCAAAAATTTTTTCGCCCCGCGCTGTGGCCAGCCCAGGCCCCCCCCCCCCCCCCGCGCTGCCACGATGGTATCCCCGt                                                     >  1:598098/1‑147 (MQ=255)                                                                                                     tcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTCG‑AATG‑CCCCCTTTTTTCAATAATGTTTTCGCCCCGCGCTGTAGCCAGCGCAGGCCCCCCCCCCCCCCCCCCGCTGCAACGTTGGTATCGCCGTTGAt                                                 >  2:429800/1‑149 (MQ=255)                                                                                                      cgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCATCGCGCGATGCAACGTTGGTATCGCCGTTGATCg                                               <  2:196508/150‑1 (MQ=255)                                                                                                           cacTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCa                                         <  1:515667/151‑1 (MQ=255)                                                                                                                ccTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTTG‑GATG‑CCCCCCTTTTTTAATAATTTTTTCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCCCCCCCCCCCCCCATGCAACGTTGGTATCGCCGGTGGTCGCCATCAGTAgg                                    >  1:66020/1‑151 (MQ=255)                                                                                                                   gTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAg                                  <  2:786934/150‑1 (MQ=255)                                                                                                                               gagaTAATCGAGATCGGC‑GGGGGTTCG‑AATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGc                     <  2:224826/151‑1 (MQ=255)                                                                                                                                agaTAATCGAGATCGGC‑GGGGGTTCG‑AATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCa                    <  2:211849/151‑1 (MQ=255)                                                                                                                                    aaTCGAGATCGGC‑GGGGGTTCG‑AATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGTGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGc                <  2:416638/151‑1 (MQ=255)                                                                                                                                       cGAGATCGGC‑GGGGGTTCG‑AATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCa              <  1:653863/150‑1 (MQ=255)                                                                                                                                          gaTCGGC‑GGGGGTTCG‑AATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGtt           <  1:499568/150‑1 (MQ=255)                                                                                                                                          gaTCGGC‑GGGGGTTCG‑AATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGGATCGAGCAATGCCAGtt           >  2:399476/1‑150 (MQ=255)                                                                                                                                            tCGGC‑GGGGGTTCG‑AATG‑CCCCCTTTTTGCAATAATGTTTTCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCCCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGGATCGAGCAATGCCAGTTc          >  1:502278/1‑149 (MQ=255)                                                                                                                                              ggc‑ggGGGTTCG‑AATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGa      <  2:452136/151‑1 (MQ=255)                                                                                                                                                    gggTTCG‑AATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTa  >  1:737396/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                     TAATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTCG‑AATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTA  >  NC_000913/646598‑646893

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