Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	436,516	0	C	A	41.5%	90.2 / 18.1	53	Q295K (CAA→AAA)	thiL	thiamine monophosphate kinase
Reads supporting (aligned to +/- strand):  ref base C (24/7);  new base A (0/22);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 3.50e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTC  >  NC_000913/436366‑436666                                                                                                                                                        |                                                                                                                                                        cTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGc                                                                                                                                                         >  1:146487/1‑151 (MQ=255)    tctcGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTAGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGC‑aa                                                                                                                                                      <  2:1057989/150‑1 (MQ=255)      tcGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGC‑AAATg                                                                                                                                                   >  2:1146419/1‑151 (MQ=255)         ccATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGC‑AAATGAcc                                                                                                                                                >  2:762164/1‑151 (MQ=255)          cATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGC‑AAATGAccg                                                                                                                                               >  1:133796/1‑151 (MQ=255)          cATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGC‑AAATGAccg                                                                                                                                               >  1:871734/1‑151 (MQ=255)                     ggAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGC‑AAATGACCGCCGATATCGa                                                                                                                                     >  1:953892/1‑150 (MQ=255)                         cAGGCGCTGCGTTGGGCGCTCTGGGGGGTTGAATATAACGTGTTGTTTTTCCCTGTG‑CGGAAATGGACCGGTGGGGCGTGGGTGGTGTCCCTCGCACTCTTGGGCGTCCCTTTTCCCTGTTCGGGGA‑AAATGACCGCCGATATCGAAggg                                                                                                                                 <  2:636704/150‑1 (MQ=255)                                  cgctgGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGtttt                                                                                                                        >  1:1100202/1‑150 (MQ=255)                                      gggCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCg                                                                                                                   >  1:520596/1‑151 (MQ=255)                                      gggCGCTCTCGGGGGTTGAAGATTAGGGGTTGTTTTCCCCTGTGCGGGAACTGACCC‑GGGGCGCGCTGGGTGTGGCTCCCGGACCCCGGGGCGGCCTTTTTCCCTGTAGGGGGA‑AAATGACCGCCGATATCGAAGGGCTTTGt                                                                                                                           <  2:388908/143‑1 (MQ=255)                                        gcgcTTTGTGGGGGGAAAGTTTACGAGTTGTTTTCCGCTGTGCCGTAACAAACCG‑GTGGCGCGCTGGATGTGGCTCTGGTTCTTTTGGGGGTCCCGTTTCCCTATGAGGGAC‑AAATTCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTg                                                                                                                 <  1:1001202/151‑1 (MQ=255)                                             ctctGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGc                                                                                                            >  2:956118/1‑151 (MQ=255)                                             ctctGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGc                                                                                                            >  1:533290/1‑151 (MQ=255)                                                 ggggttGAAGATTACGAGTTGTTTTCCCCTGTGCCGAAACTGACCC‑GGGGCGCGTTGGTTGTGGCTCTCGGACACCGGGGCGTCCCTTTTACCTGTACGGGGA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCg                                                                                                           <  1:924020/143‑1 (MQ=255)                                                 gggggTGAATATTACGTGTTGTTTTTCACTGTGCGGGAACTGACCC‑GGGGCGCGCTGGATGGGGCTCTCGGACACCGGGGCGTCCCTTTTCCCTGTATGGGGA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACgg                                                                                                             <  1:1191507/143‑1 (MQ=255)                                                 gggggTGAAGATTGCGTGTGGTTTTTCACGGAACGGGAACGGATCG‑GGGGGGCGTGGGTTGGGGGTCCCGAACACCGGGGGGTACCTTTTACGTGTGGCGAGC‑AAATGGCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGa                                                                                                          <  2:743651/146‑1 (MQ=255)                                                 gggggTGAAGATTACGAGTTGTGTTTCCCTGTGCCGGAACTGACCC‑GTGGCGCGCTGGATGGGGCTCTCGGACACCGGGGCGTCCCGTTTCCCTGTATCGGGA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACg                                                                                                              <  1:237444/142‑1 (MQ=255)                                                   cGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCt                                                                                                      >  2:994031/1‑151 (MQ=255)                                                      tGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGtt                                                                                                   >  2:197951/1‑151 (MQ=255)                                                        aaGATTACGAGTTTTTTTTCCCTGTGCCGGAACTGACCC‑GTGGCGCGTTGTTTGTGGCTCTCGCCCCCGTGGGCGT‑CCTTTTTCCTTTTTCGGGCAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTAc                                                                                                 <  2:586158/151‑1 (MQ=255)                                                        aaGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGTCGGCGAACCTGtt                                                                                                   >  2:830680/1‑149 (MQ=255)                                                            ttACGAGTTTTTTTTCCCGGTGCCGGAATGGACCC‑GTGGGGCGTTGGTTGTGGCTCTCGGACACCTGGGCTTACCGTTTCCCTGTATCGGGA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAtt                                                                                              <  1:546618/150‑1 (MQ=255)                                                            ttACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACAtt                                                                                              >  1:515659/1‑150 (MQ=255)                                                               cGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGAc                                                                                          >  2:30567/1‑151 (MQ=255)                                                                gAGTTGTTTTTCACTGTGCCGGAACTGACCC‑GTGGCGCGCTGGATGTGGCTCTCGGACCCCGGGGCGTCCCTTTTCCCTTTATCGGAA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTa                                                                                             <  1:989156/147‑1 (MQ=255)                                                                 aGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACt                                                                                         >  1:898334/1‑150 (MQ=255)                                                                   ttttttttCATTGTGCCGAAACTAACCG‑GGGGCGCGTGGGTTGTGGCTCTCGGCCCCTTGGGCGTCCTTTTTCCTTGTATCGGGA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTg                                                                                        <  2:172621/144‑1 (MQ=255)                                                                   ttGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTgg                                                                                       >  2:470945/1‑150 (MQ=255)                                                                    tgttTTTCATGGTCCGGGATCAGAACC‑GTGGCCCGCTGGATGTGGCTCGCGGCCACCGGGGCGTCCCTTTTCCTTTTTTCGGGA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAg                                                                                            <  1:932177/144‑1 (MQ=255)                                                                    tgtgTTCCACTGTCCCGGAACTGAACC‑GGGGGGCGTTGGATGTGGCTCTCGGACACCGGGGCGTACCTTTTCCTTGTAGCGGGA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGa                                                                                      <  1:786766/150‑1 (MQ=255)                                                                      tttttCACTGTGCGGGAACTGACCG‑GTGGCGCGGTGGTTGTGGCTCTGGAACACCGGGGCGTCCCTTTTCCCTGTTGGGGGA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGaaa                                                                                    <  1:942051/148‑1 (MQ=255)                                                                      tgtTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGaa                                                                                     >  1:803419/1‑149 (MQ=255)                                                                        tttCACTGTGCCGGAACTGACCC‑GGGGCGCGCTGGATGTGGCTCTCGGACCCCTGGGCGTCCCGTTTCCCTGTATGGGGA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAg                                                                                   <  1:610056/149‑1 (MQ=255)                                                                           cACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGatat                                                                              >  2:438528/1‑151 (MQ=255)                                                                             cTGTCCCGGCACTAAACC‑GGGGCAGGGTGGAGGTTGCTCTCGAACCCCGGGGCGTCCTTTTTGCCTTTATCGGGC‑AAAGGCCCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGata                                                                               <  1:853244/148‑1 (MQ=255)                                                                                 gCCGGACCTGAACG‑GTGGCGCGTTGGATGTGTCTCTCGGACCCCTGGGCTTCCCTTTTACCTGTATGGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCa                                                                         <  1:800073/150‑1 (MQ=255)                                                                                    ggAACTGACCC‑GGGGCGCGTGGGAGGTGGCTCCCGGACACCGGGGCGTCCCTTTTACCTGTAGCGGGA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGAcc                                                                          <  2:404142/146‑1 (MQ=255)                                                                                    ggAACTGACCC‑GGGGCGCGCTGGATGTGGCTCTCGGACCCGGGGGCGTACCTTTTCCTTTTGTCGGAA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtttt                                                                     <  1:1061579/151‑1 (MQ=255)                                                                                    ggAACTGAACC‑GGGGCGCGGTGGATGTGGCTCCGGGACCCGTGGGCGTCCTGTTTACCTGTAGCGGGA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCAtttt                                                                     <  2:79127/151‑1 (MQ=255)                                                                                    gaaccTGACCG‑GTGGCGCGCTGGATGTGGCTCTCGGCCCCGGGGGCTTCCCGTTTACCTGTATGGGAA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGAcc                                                                          <  1:416955/142‑1 (MQ=255)                                                                                      aaCTGACCC‑GTGGCGCGCTGGATGTGGCTCTCGGACCCCGGGGCGTCCCTTTTCCCTGTATCGGAA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTg                                                                    <  1:490376/150‑1 (MQ=255)                                                                                      aaCTGACCC‑GTGGCGCGCTGGAGGTGCCTCTCGGACACCGGGGCGTCCCTTTTCCCTGTACGGGGA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGc                                                                   <  1:197953/151‑1 (MQ=255)                                                                                         tGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCAc                                                                >  2:510584/1‑151 (MQ=255)                                                                                            accggTGGCGCGCTGGTTGTGGCTCGCGCCCACGTGGGCGTCCCTTTTACCTGTAGGGGAA‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCa                                                                         <  2:559775/136‑1 (MQ=255)                                                                                                           gATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTcgcg                                               >  1:616895/1‑151 (MQ=255)                                                                                                                                 gCGTACCGTTTACCTGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGt                          >  2:129989/1‑150 (MQ=255)                                                                                                                                 gCGTACCGTTTACCTGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGt                          >  1:456003/1‑150 (MQ=255)                                                                                                                                      cccttttccttttgcGGGCAAAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAAt                       <  1:852476/135‑1 (MQ=255)                                                                                                                                       cGTTTACCTGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGt                   >  2:652803/1‑151 (MQ=255)                                                                                                                                         tttACCTTTATGGGGA‑AATTAACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATc                      <  1:873521/146‑1 (MQ=255)                                                                                                                                            aCCTGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCAt              >  1:699303/1‑151 (MQ=255)                                                                                                                                                  acggggA‑AAATAACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCa               <  1:470959/141‑1 (MQ=255)                                                                                                                                                     gggA‑AAAGGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACt         <  1:328747/147‑1 (MQ=255)                                                                                                                                                        c‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTc  >  2:725805/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                        CTTTCTCGCCATGTTGAACCGGAACAGGCGCTGCGCTGGGCGCTCTCTGGCGGTGAAGATTACGAGTTGTGTTTCACTGTGCCGGAACTGAACC‑GTGGCGCGCTGGATGTGGCTCTCGGACACCTGGGCGTACCGTTTACCTGTATCGGGC‑AAATGACCGCCGATATCGAAGGGCTTTGTTTTATTCGTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTC  >  NC_000913/436366‑436666

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