Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,261,842	0	T	G	20.0%	112.3 / 12.2	45	T12P (ACC→CCC)	prs	ribose‑phosphate diphosphokinase
Reads supporting (aligned to +/- strand):  ref base T (6/30);  new base G (9/0);  total (15/30)
Fisher's exact test for biased strand distribution p-value = 5.65e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGG‑TTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAAC  >  NC_000913/1261704‑1261980                                                                                                                                             |                                                                                                                                           aCTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGG‑TTTGGAATACGTTGTGCTAG‑TTCCGGGGTg                                                                                                                                           >  1:506083/1‑140 (MQ=255)     gatgatGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGG‑TTGGCGATACGGTTTGGTTGTTTCCGgggggggt                                                                                                                                        >  2:573162/1‑138 (MQ=255)     gatgatGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGG‑TTGGCAATACGTTGTGCTAG‑TTCCGGGGGGGGGt                                                                                                                                       >  1:1207456/1‑140 (MQ=255)     gatgatGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGG‑TTGGCAATACGGTGTGCTAG‑TTCCGGGGGGGCgg                                                                                                                                       >  1:281877/1‑139 (MQ=255)       tgatgaAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGG‑TTGGCAATACGTTGTGCTAG‑TTCCGGGGGGGCGTTa                                                                                                                                     >  1:659055/1‑140 (MQ=255)               atatCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGG‑TTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCaaa                                                                                                                             <  2:625755/140‑1 (MQ=255)               atatCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGG‑TTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCaaa                                                                                                                             <  2:682695/140‑1 (MQ=255)                   caccacCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGG‑TTGGCAATACGTTGTGCTAG‑TTCCGGGGGGGCGTTACCAGCAAAAAGc                                                                                                                         >  1:972345/1‑140 (MQ=255)                     ccaccGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGG‑TTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCtt                                                                                                                       <  1:23389/140‑1 (MQ=255)                         cgcgTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGGTTTGGCAA‑ACCTTGTGTTAG‑TTTCCGGGGGGCGTTGCCAACAAAAAGCTTcttt                                                                                                                   >  2:1225550/1‑137 (MQ=255)                           cgTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGG‑TTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATc                                                                                                                 <  2:466722/140‑1 (MQ=255)                           cgTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGG‑TTGGCAATACGTTGTGCTAG‑TTCCGGGGGGGCGTTACCAGCAAAAAGCTTCATATc                                                                                                                 >  1:1084816/1‑140 (MQ=255)                           cgTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGG‑TTGGCAATACGTTGTGCTAG‑TTCCGGGGGGGCGTTACCAGCAAAAAGCTTCATATc                                                                                                                 >  1:502215/1‑140 (MQ=255)                           cgTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGG‑TTGGCAATACGGTGTGCTAG‑TTCCGGGGGGGCGGTACCAGCAAAAAGCTTCATATc                                                                                                                 >  2:11747/1‑140 (MQ=255)                             tACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGG‑TTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAg                                                                                                               >  1:436107/1‑140 (MQ=255)                               cATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGG‑TTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGc                                                                                                             <  2:192172/140‑1 (MQ=255)                                aTTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGG‑TTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCa                                                                                                            <  2:22792/140‑1 (MQ=255)                                aTTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGG‑TTGGCAATACGTTGTGCTAG‑TTCCGGGGGGGCGGTACCAGCAAAAAACTTCATATCAGGcc                                                                                                            >  2:985284/1‑139 (MQ=255)                                  tttCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGG‑TTGGCAATACGTTGTGCTAG‑TTCCGGGGGGGGGTTACCAGCAAAAAGCTTCATATCAGGCACg                                                                                                          >  2:1036112/1‑140 (MQ=255)                                             tGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGG‑TTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTc                                                                                               <  1:627073/140‑1 (MQ=255)                                                aCGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGG‑TTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAgg                                                                                            <  1:573162/140‑1 (MQ=255)                                                      aCTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGG‑TTGGCAATACGTTGTGCTAG‑TTCCGGGGGGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCg                                                                                      >  1:1063434/1‑140 (MQ=255)                                                          cGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGG‑TTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCa                                                                                  <  1:908243/140‑1 (MQ=255)                                                          cGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGG‑TTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCa                                                                                  <  2:502215/140‑1 (MQ=255)                                                                        cGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGG‑TTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAAt                                                                    <  1:168645/140‑1 (MQ=255)                                                                        cGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGG‑TTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAAt                                                                    <  1:1225550/140‑1 (MQ=255)                                                                                     cgTCGCCGAGTGAAGTGTACAGGCGG‑TTGGCAAAACGTTGTGCTAG‑TTCCGGGGGGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAAACCCAGGCATGCGTCCCTTGGTGGAAAGAAACTGCCGAAAACtg                                                       >  2:134824/1‑140 (MQ=255)                                                                                       tCGCCGAGTGAAGTGTACAGGCGG‑TTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACtgtg                                                     <  2:659055/140‑1 (MQ=255)                                                                                       tCGCCGAGTGAAGTGTACAGGCGG‑TTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACtgtg                                                     <  2:1063434/140‑1 (MQ=255)                                                                                       tCGCCGAGTGAAGTGTACAGGCGG‑TTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACtgtg                                                     <  1:11747/140‑1 (MQ=255)                                                                                             aGTGAAGTGTACAGGCGG‑TTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGcc                                               <  2:34464/140‑1 (MQ=255)                                                                                             aGTGAAGTGTACAGGCGG‑TTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGcc                                               <  2:1195747/140‑1 (MQ=255)                                                                                             aGTGAAGTGTACAGGCGG‑TTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGcc                                               <  1:815757/140‑1 (MQ=255)                                                                                                  aGTGTACAGGCGG‑TTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCa                                          <  2:807197/140‑1 (MQ=255)                                                                                                        cAGGCGG‑TTGGCAATACGTTGTGGTAG‑TTCCGGGGGGGGGTTACCAGCAAAAAACTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGGCCATTGGGGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATcc                                    >  2:395275/1‑140 (MQ=255)                                                                                                              g‑TTTGCAATACGTTGTGCTAG‑TTCCGGGGGGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTccc                              >  2:387656/1‑139 (MQ=255)                                                                                                              g‑TTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCAGGCCAGGCATGATCCTTTCCa                              <  1:302828/140‑1 (MQ=255)                                                                                                                 tGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAgg                            <  1:1043022/140‑1 (MQ=255)                                                                                                                 tGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAgg                            <  2:436107/140‑1 (MQ=255)                                                                                                                 tGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAgg                            <  2:281877/140‑1 (MQ=255)                                                                                                                   gCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCg                          <  2:1084816/140‑1 (MQ=255)                                                                                                                   gCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCg                          <  2:209867/140‑1 (MQ=255)                                                                                                                   gCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCg                          <  1:774663/140‑1 (MQ=255)                                                                                                                               gCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAaga              >  2:1257577/1‑140 (MQ=255)                                                                                                                               gCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAaga              <  1:985284/140‑1 (MQ=255)                                                                                                                                cTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCCTCATATCAGGCACGAGAAGAGCCTCAGGCACGCGTCCATTGGTGGAAAGGATCTGCTGAAAACTGTGCGGGCCAGTCATGATCCTTTCCAGGCTGTGTATTAAagag             <  2:1228487/140‑1 (MQ=255)                                                                                                                                 tAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAgc            >  2:21491/1‑140 (MQ=255)                                                                                                                                  aG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAgcg           <  1:21491/140‑1 (MQ=255)                                                                                                                                   g‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAgcgc          >  1:297099/1‑140 (MQ=255)                                                                                                                                     ttCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGcgcg         >  2:861838/1‑140 (MQ=255)                                                                                                                                         ggggTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGc     >  2:631535/1‑140 (MQ=255)                                                                                                                                         ggggTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGc     >  1:355009/1‑140 (MQ=255)                                                                                                                                          gggTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCa    <  1:1052078/140‑1 (MQ=255)                                                                                                                                            gTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAAc  <  2:1207456/140‑1 (MQ=255)                                                                                                                                             |                                                                                                                                           ACTGGATGATGAAAATATCACCACCGCGTACATTTTCATTAATTTGTACGCTGACTTCGCCATCGCTAAAGCGACCTACAGCGGCGTCGCCGAGTGAAGTGTACAGGCGG‑TTGGCAATACGTTGTGCTAG‑TTCCGGGGTGGCGTTACCAGCAAAAAGCTTCATATCAGGCACGAGAAGAACCTCAGGCATGCGTCCATTGGTGGAAAGAATCTGCCGAAAACTGTGCGGGCCAGGCATGATCCTTTCCAGGCGGTGTATTAAAGAGCGCGATGCAAC  >  NC_000913/1261704‑1261980

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