Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,751,949	0	C	T	21.6%	104.0 / 12.5	52	A52V (GCC→GTC)	recN	recombination and repair protein
Reads supporting (aligned to +/- strand):  ref base C (13/27);  new base T (11/0);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 5.24e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCACAACTGACCATCAGCAACTTTGCTATCGTTCGTGAGCTTGAGATTGATTTTCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACG‑CCAGCGGCTCTGCGCTGGCTGGAAGAAAACCAGCTTGAAGACGGGCATGAATGTTTGCTTCGTCGCGTGATCAGCAGCGA  >  NC_000913/2751801‑2752093                                                                                                                                                     |                                                                                                                                                  gCACAACTGACCATCAGCAACTTTGCTATCGTTCGTGAGCTTGAGATTGATTTTCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGc                                                                                                                                                   <  1:864027/149‑1 (MQ=255) gCACAACTGACCATCAGCAACTTTGCTATCGTTCGTGAGCTTGAGATTGATTTTCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGc                                                                                                                                                   <  1:702715/149‑1 (MQ=255)   actaCTGACCATCAGCAACTTTGCTATCGTTCGTGAGCTTGAGATTGATTTTCATAGCGGCATGACCGTAATAACTGGCGAGTCCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCCg                                                                                                                                                 <  2:647390/146‑1 (MQ=255)        tGACCATCAGCAACTTTGCTATCGTTCGTGAGCTTGAGATTGATTTTCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCCGACATg                                                                                                                                            >  1:202734/1‑149 (MQ=255)          acaaTCAGCAACTTTGCTATCGTTCGTGAGCTTGAGATTGATTTTCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCCGACATGGt                                                                                                                                          <  1:117612/146‑1 (MQ=255)          aCCATCAGCAACTTTGCTATCGTTCGTGAGCTTGAGATTGATTTTCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCCGACATGGt                                                                                                                                          <  1:280786/149‑1 (MQ=255)            cATCAGCAACTTTGCTATCGTTCGTGAGCTTGAGATTGATTTTCATAGCGGCATGACCGTAATAACTGGCGCGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGc                                                                                                                                        >  1:157061/1‑149 (MQ=255)            cATCAGCAACTTTGCTATCGTTCGTGAGCTTGAGATTGATTTTCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGc                                                                                                                                        >  1:92587/1‑149 (MQ=255)             aTCAGCAACTTTGCTATCGTTCGTGAGCTTGAGATTGATTTTCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCg                                                                                                                                       <  1:692090/149‑1 (MQ=255)                   aaCCTTGCTATCGTACGTGAGCTTGAGATTGATTTCCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCgg                                                                                                                                 <  1:243261/149‑1 (MQ=255)                      tttGGTATCGTTCGTGATCTGGAGATTGATTTTCATGGGCGCATGACCGTAATAACTGGCGGGCCCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCCGAAGCCGACATGGTGCGTACCTgcgc                                                                                                                              <  2:699709/149‑1 (MQ=255)                             tcgttcgtGAGCTTGAGATTGATTTTCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGTCGACATGGTGCGTACCGGCGCTGCTcgc                                                                                                                       >  1:516149/1‑149 (MQ=255)                                       gCTTGAGATTGATTTTCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCtgt                                                                                                             <  1:256582/149‑1 (MQ=255)                                           gagaTTTATTTTGATAGCGGCATGATGGTAATAACTGGCGAGACCGGCGCGGGTGAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTgcgc                                                                                                         <  2:161813/149‑1 (MQ=255)                                                    tttCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTtc                                                                                                >  1:646932/1‑149 (MQ=255)                                                     ttCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCCGTATCATTAAAAAAAAAAGTACAGGCGCTGCTCGCGCTTAACTGTGCGCCCGTTTTtct                                                                                               >  2:65313/1‑149 (MQ=255)                                                       cATAGCGGCATGACCGTAATAACTTGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTctct                                                                                             <  1:329587/149‑1 (MQ=255)                                                       cATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTctct                                                                                             <  1:863853/149‑1 (MQ=255)                                                        aTAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGAc                                                                                                                 <  2:70813/128‑1 (MQ=255)                                                        aTAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGAc                                                                                                                 >  1:70813/1‑128 (MQ=255)                                                              gCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGAt                                                                                      <  1:647633/149‑1 (MQ=255)                                                                     cGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCCGTATCATTAAAAAAAAAAGTACTGGCGCTGCTCGCGCTGACCTGAGCGCCCGTTTTTCTCTGAAAGACACG‑CGAg                                                                              >  2:627904/1‑149 (MQ=255)                                                                     cGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCCGTATCATTAAAAAAAAAAGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTATGAAGGATACG‑ccca                                                                              >  2:613701/1‑147 (MQ=255)                                                                                     aCCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACG‑CCAGCGGCTCTGCgctggct                                                              <  1:572716/149‑1 (MQ=255)                                                                                     aCCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCCGTAACATCAAAAAAAAAAGTACCGGCCCAGCTCGCGCTGACCGTTGTGCCCGTTTTTCTATGTAAGAAACGACCGGCGGC‑CTgcggtagcc                                                              >  2:74823/1‑143 (MQ=255)                                                                                      ccGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCCGTAACATCAAAAAAAAAAGTACCGGCGCTGCTCGCGCTGAACCGTCCGCCCGTTTTTCACTGACAGAAACG‑CCACCGGCGCTgcgccacacg                                                             >  2:612072/1‑143 (MQ=255)                                                                                         gcgcGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTGGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACG‑CCAGCGGCTCTGCGCTGGCTGgaa                                                          <  1:741413/149‑1 (MQ=255)                                                                                             gggTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACG‑CCAGCGGCTCTGCGCTGGCTGGAAGaaa                                                      <  1:628358/149‑1 (MQ=255)                                                                                             gggTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCCCTATCGTCGAAAAAGTGCGTACCGGCGCTACTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGCAAGTTACC‑CCACCGGCTCTGCGCTGGCTGGGAGaaa                                                      >  2:626985/1‑149 (MQ=255)                                                                                               gaaaaTCTTTTGCACTAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCCGAAATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACG‑CCAGCGGCTCTGCGCTGGCTGGAAGAAAAc                                                    <  1:273117/147‑1 (MQ=255)                                                                                               gTAAATCTATGGCAATAGATGCCCGCGGTCTTTGTCTCGGTGGTAGTGATAATGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACG‑CCAGCGGCTCTGCGCTGGCTGGAAGAAAAc                                                    <  1:257440/149‑1 (MQ=255)                                                                                                tAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCCGTATCATTAAAAAAAAAAGTACCGGAGCTGCTCGCGCTGACCTGTTAGGCCGTTTTTCTCTGAAAGATACG‑CCAGCGACAAGGCGCTGGCTGGCAGAacaca                                                   >  2:605316/1‑148 (MQ=255)                                                                                                  aaTCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCCGTATCATCAAAAAAAAGAGTCCCCGCGCTGCTCGCGCTGCCCTGTGAGCCCGTTTTTCTCTGAGAGATACG‑CCAGCGGCTCTGCGCTGGCTGGAAGAAAACCAg                                                 >  2:313291/1‑149 (MQ=255)                                                                                                     cTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACG‑CCAGCGGCTCTGCGCTGGCTGGAAGAAAACCAGCtt                                              <  1:662779/149‑1 (MQ=255)                                                                                                      tATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCCGTATCATTAAAAAAAAAAGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACG‑CCAGCGGCTCTGCGCTGGCTGGAAGAAAACCAGCTTg                                             >  2:474480/1‑149 (MQ=255)                                                                                                        ttGGAATAGTAGCCCTCGGTCTTTGTTTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTTACCTGTGCGCCCGTTTTTCTCTGAAAGATACG‑CCAGCGGCTCTGCGCTGGCTGGAAGAAAACCAGCTTGaa                                           <  2:445863/149‑1 (MQ=255)                                                                                                           cAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGaa                                                                                          >  1:244893/1‑100 (MQ=255)                                                                                                           cAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGaa                                                                                          <  2:244893/100‑1 (MQ=255)                                                                                                           cAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTCAAATTTACG‑CTTGCGGCTCTGCGCTGGCTGGAAGAAAACCAGCTTGAAGAc                                        >  1:760336/1‑149 (MQ=255)                                                                                                           acaacgATGCCCTCGGTCTTTTTATCCGTGGGCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACG‑CCAGCGGCTCTGCGCTGGCTGGAAGAAAACCAGCTTGAAGAc                                        <  1:626985/144‑1 (MQ=255)                                                                                                             aTAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACG‑CCAGCGGCTCTGCGCTGGCTGGAAGAAAACCAGCTTGAAGACgg                                      <  2:289901/149‑1 (MQ=255)                                                                                                                     ccTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACG‑CCAGCGGCTCTGCGCTGGCTGGAAGAAAACCAGCTTGAAGACgg                                      <  2:343161/141‑1 (MQ=255)                                                                                                                     ccTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACG‑CCAGCGGCTCTGCGCTGGCTGGAAGAAAACCAGCTTGAAGACgg                                      >  1:343161/1‑141 (MQ=255)                                                                                                                     ccTCGGTCTTTGTCTCGGTGGTCGCCGTATCATTAACAAAAAGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACG‑CCAGAGGCTCTGCGCTGGCTGGAAGAAAACCAGCTGAACGCCGGGCAtccac                              >  2:541531/1‑145 (MQ=255)                                                                                                                       gcGGTTTTTGGCTCGGTGTCCGTGTTGAAGGCGACATGGTGCGTACCGGCGCTGCCCGCGCTGGCCTGTGCGCCCGTTTTTCTCTGAAAGATACG‑CCGGCGGCTGTGCGGTGGCTGGAAGAAAACCAGCTTGAAGACGGGCATGAATGt                            <  2:96118/148‑1 (MQ=255)                                                                                                                         ggTCTTTGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACG‑CCAGCGGCTCTGCGCTGGCTGGAAGAAAACCAGCTTGAAGACGGGCATGAATGttt                          <  1:473014/149‑1 (MQ=255)                                                                                                                           gcTTTGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACG‑CCAGCGGCTCTGCGCTGGCTGGAAGAAAACCAGCTTGAAGACGGGCATGAATGTTTGc                        <  2:563249/148‑1 (MQ=255)                                                                                                                              ttGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACG‑CCAGCGGCTCTGCGCTGGCTGGAAGAAAACCAGCTTGAAGACGGGCATGAATGTTTGCTtc                     <  1:3762/149‑1 (MQ=255)                                                                                                                                   tcGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACG‑CCAGCGGCTCTGCGCTGGCTGGAAGAAAACCAGCTTGAAGACGGGCATGAATGTTTGGTTCGTcgc                >  1:707150/1‑149 (MQ=255)                                                                                                                                     ggtggtCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACG‑CCAGCGGCTCTGCGCTGGCTGGAAGAAAACCAGCTTGAAGACGGGCATGAATGTTTGCTTCGTCGCGt              <  1:867745/149‑1 (MQ=255)                                                                                                                                         gtCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACG‑CCAGCGGCTCTGCGCTGGCTGGAAGAAAACCAGCTTGACGACGGGCATGAATGTTTGCTTCGTCGCGTGATc          >  1:500205/1‑149 (MQ=255)                                                                                                                                          tCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACG‑CCAGCGGCTCTGCGCTGGCTGGAAGAAAACCAGCTTGAAGACGGGCATGAATGTTTGCTTCGTCGCGTGATca         >  2:78476/1‑149 (MQ=255)                                                                                                                                                tGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACG‑CCAGCGGCTCTGCGCTGGCTGGAAGAAAACCAGCTTGAAGACGGGGATGAATGTTTGCTTCGTCGCGTGATCAGCAGCg   >  1:456213/1‑149 (MQ=255)                                                                                                                                                 gAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACG‑CCAGCGGCTCTGCGCTGGCTGGAAGAAAACCAGCTTGAAGACGGGCATGAATGTTTGCTTCGTCGCGTTATCAGCAGCGa  >  2:589024/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                  GCACAACTGACCATCAGCAACTTTGCTATCGTTCGTGAGCTTGAGATTGATTTTCATAGCGGCATGACCGTAATAACTGGCGAGACCGGCGCGGGTAAATCTATTGCAATAGATGCCCTCGGTCTTTGTCTCGGTGGTCGCGCTGAAGCCGACATGGTGCGTACCGGCGCTGCTCGCGCTGACCTGTGCGCCCGTTTTTCTCTGAAAGATACG‑CCAGCGGCTCTGCGCTGGCTGGAAGAAAACCAGCTTGAAGACGGGCATGAATGTTTGCTTCGTCGCGTGATCAGCAGCGA  >  NC_000913/2751801‑2752093

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