Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	445,558	0	A	C	25.6%	82.8 / 15.4	39	V370G (GTA→GGA)	yajR	putative major facilitator superfamily transport protein YajR
Reads supporting (aligned to +/- strand):  ref base A (16/13);  new base C (0/10);  total (16/23)
Fisher's exact test for biased strand distribution p-value = 2.23e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGCAAACGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTG‑CGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGA  >  NC_000913/445420‑445694                                                                                                                                            |                                                                                                                                         aGCAAACGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTG‑CGCCAGCGAGAAATAc                                                                                                                                         >  1:366289/1‑140 (MQ=255)     aaCGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTG‑CGCCAGCGAGAAATAccccc                                                                                                                                     <  2:1218129/140‑1 (MQ=255)      aCGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTG‑CGCCAGCGAGAAATACCCCCt                                                                                                                                    <  2:649892/140‑1 (MQ=255)      aCGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTG‑CGCCAGCGAGAAATACCCCCt                                                                                                                                    <  1:1570743/140‑1 (MQ=255)       cGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTG‑CGCCAGCGAGAAATACCCCCTg                                                                                                                                   <  2:366289/140‑1 (MQ=255)        gCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTG‑CGCCAGCGAGAAATACCCCCTGa                                                                                                                                  <  2:717597/140‑1 (MQ=255)                  cGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTG‑CGCCAGCGAGAAATACCCCCTGACCGTCAAACa                                                                                                                        <  2:401981/140‑1 (MQ=255)                  cGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTG‑CGCCAGCGAGAAATACCCCCTGACCGTCAAACa                                                                                                                        <  1:819700/140‑1 (MQ=255)                   gCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTG‑CGCCAGCGAGAAATACCCCCTGACCGTCAAACAt                                                                                                                       <  2:1238612/140‑1 (MQ=255)                   gCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTG‑CGCCAGCGAGAAATACCCCCTGACCGTCAAACAt                                                                                                                       <  1:294479/140‑1 (MQ=255)                     cTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGGCTGTCAGCCACACTGCGGCCAGCATTG‑CGCCAGCGAGAAAAACCCCCTGACCGTCAAACATGc                                                                                                                     >  2:624229/1‑140 (MQ=255)                      tCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGGACTGGGGGCTGGCAGCCACACCGCGGCCCGCCTTG‑CGCCAGCGGGAAAAAACCCCCCACCGGCAAAAAAGcc                                                                                                                    >  2:434972/1‑140 (MQ=255)                                cAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTG‑CGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAg                                                                                                          <  2:1441528/140‑1 (MQ=255)                                  aTGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTG‑CGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAgcc                                                                                                        <  1:371593/140‑1 (MQ=255)                                     ttCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGCCTGTCAGCCACCCTGCGGCCAGCATTG‑CGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAgccgcc                                                                                                     <  2:155915/140‑1 (MQ=255)                                                         cAACCTGCTGACATCCGGCGTTTTTTTCATGGTCCGGGCGACTGCCAGCCACACTGCGGCCAGCATTG‑CGCCAGCGAAAAATCCCCCCTGCCCGTAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATcgcca                                                                                 <  2:42591/140‑1 (MQ=255)                                                         cAAACTGCTGACATCCGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTG‑CGCCAGCGAGAAATCCCCCCTGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATcgcca                                                                                 <  1:762611/140‑1 (MQ=255)                                                             ctgctgACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTG‑CGCCCGCGAGAAAAACCCCCCGACCCGCAAACATGCCGTTAATCCCGCCCCCCCGCGAACCGCCCAAccccacgcc                                                                             >  2:265518/1‑140 (MQ=255)                                                                   aCATACGGCGGTTCTTTCAGGGTACGGGCGACTGTCAGCCACACGGCGGCCAGCATGG‑CGCCAGCGAAAAATCCCCCCTGACCGTAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGaaa                                                                       <  2:856853/140‑1 (MQ=255)                                                                        gggcggTTCTTTCATGGTACTGGCGCCTGTCAGCCCCCCTGCGGCCAGCATTG‑CGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAActggc                                                                  <  2:640229/139‑1 (MQ=255)                                                                        cggcggTTCTTTCATGGTACTGGCGACTGTCAGCCACACCGCGGCCCGCATTG‑CGCCAGCGAGAAAAACCCCCCGACCCGCCAACAAGCCGTTAATCCAGCCCCCCCGCGAACCGCCAAACCCCCCCCCCAGAAAccggc                                                                  >  2:1071537/1‑140 (MQ=255)                                                                             tttCTTTCATGGTCTTGGCGCCTGCCAGCCCCACTGCGGCAACCTTTG‑CGCCCGGGGGAAATCCCCCCTGACCGTAAACCATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCtggtg                                                             <  2:1281630/139‑1 (MQ=255)                                                                             gTTCTTTCATGGTACGGGCGCCTGTCAGCCACACTGCGGCCAGCTTTG‑CGCCGGCGAGAAATCCCCCCTGACCGCCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCtggtg                                                             <  2:120404/140‑1 (MQ=255)                                                                                        gTACTGGCGCCTGTCGCCAACTTGGCGGCCAGCATTG‑CCCCAGCGAAAATTCCCCCCTGACCGACAACCTTGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACAcc                                                  <  2:821269/140‑1 (MQ=255)                                                                                        gTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTG‑CGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACAcc                                                  >  2:1553190/1‑140 (MQ=255)                                                                                          aCTGGCGACTGTCAGCCACACTGCGGCCAGCATTG‑CGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCa                                                >  1:1209469/1‑140 (MQ=255)                                                                                             ggCGACTGTCAGCCACACTGCGGCCAGCATTG‑CGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCg                                             >  1:885051/1‑140 (MQ=255)                                                                                              gCGACTGTCAGCCACACTGCGGCCAGCATTG‑CGCCAGCGAGAAATa                                                                                                                                          >  1:930570/1‑46 (MQ=255)                                                                                              gCGACTGTCAGCCACACTGCGGCCAGCATTG‑CGCCAGCGAGAAATa                                                                                                                                          <  2:930570/46‑1 (MQ=255)                                                                                              gCGACTGTCAGCCACACTGCGGCCAGCATTG‑CGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGc                                            >  1:635853/1‑140 (MQ=255)                                                                                               ggccTGTCAGCCACCCTGCGGCCAGCTTGG‑CGCCAGCGAGAAACCCCCCCGGCCCTTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGcc                                           <  2:243794/137‑1 (MQ=255)                                                                                                 cccgccAGCAACACTGCGCCCAGCATTG‑CGCCAGCGAGAAATCCCCCCTGACCGTAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGt                                         <  1:1401289/135‑1 (MQ=255)                                                                                                 aCTGTCAGCCACACTGCGGCCAGCATTG‑CGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGt                                         >  1:55206/1‑140 (MQ=255)                                                                                                 aCTGTCAGCCACACTGCGGCCAGCATTG‑CGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAAAAACTGGCTGGTGGGGGAAACACCCATCGCCGt                                         >  2:1207217/1‑140 (MQ=255)                                                                                                  cTGTCAGCCACACTGCGGCCAGCATTG‑CGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTa                                        <  1:1493790/140‑1 (MQ=255)                                                                                                  cTGGCACCCACACGGCGGCCAG‑ATTGCCGCCAGCAAGAAACCCCCCCTGCCCGACACACATGCCGTAAATCCGGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTa                                        <  2:363740/140‑1 (MQ=255)                                                                                                   tgccacccACACGGCGGCCAGCTTTG‑CCCCAGCAAAAAATCCCCCCTGCCCGTCAAACATGCCGTTAACCCAGCCCCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTAc                                       <  2:188687/134‑1 (MQ=255)                                                                                                   tGTCAGCCACACTGCGGCCAGCATTG‑CGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTAc                                       >  2:704511/1‑140 (MQ=255)                                                                                                   tGTCAGCCACACTGCGGCCAGCATTG‑CGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTAc                                       >  1:22301/1‑140 (MQ=255)                                                                                                   tGTCAGCCACACTGCGGCCAGCATTG‑CGCCAGCGAGAAATACCCCCTGACCGCCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCGGGCGGAGCAAACACCCATCGCCGAAc                                       >  2:624047/1‑140 (MQ=255)                                                                                                          caccCTGCGGCCAGCATTG‑CGCCGGCGAAAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTa                                <  1:265518/140‑1 (MQ=255)                                                                                                            cacGGCGGCCAGCATTG‑CGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAAc                              <  1:547649/140‑1 (MQ=255)                                                                                                             acTGCGGCCAGCATTG‑CGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAAcc                             >  1:922525/1‑140 (MQ=255)                                                                                                             acTGCGGCCAGCATTG‑CGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAAcc                             >  1:125841/1‑140 (MQ=255)                                                                                                               tGCGGCCAGCATTG‑CGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACctg                           >  2:668690/1‑140 (MQ=255)                                                                                                               tGCGGCCAGCATTG‑CGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACctg                           >  2:294610/1‑140 (MQ=255)                                                                                                                 cGGCCAGCATTG‑CGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAACCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACctgct                         >  2:1113590/1‑140 (MQ=255)                                                                                                                     cAGCATTG‑CGCCAGCGAGAAATCCCCCCTGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCg                     <  1:153637/140‑1 (MQ=255)                                                                                                                      aGCTTTG‑CGCCAGCGAAAAATCCCCCCTGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGa                    <  1:1071537/140‑1 (MQ=255)                                                                                                                      aGCATTG‑CGCCAGCGAGAAATCCCCCCTGCCCGTCAACCATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGa                    <  1:434972/140‑1 (MQ=255)                                                                                                                      aGCATTG‑CGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGa                    >  1:648297/1‑140 (MQ=255)                                                                                                                      aGCATGG‑CGCCAGCGAAAAATCCCCCCTGCCCGTAAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGa                    <  1:1518581/140‑1 (MQ=255)                                                                                                                        cATTG‑CGCCGGCGAGAAATCCCCCCTGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGAct                  <  2:408538/140‑1 (MQ=255)                                                                                                                        cATTG‑CCCCAGCGAGAAATCCCCCCTGCCCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGAct                  <  1:343562/140‑1 (MQ=255)                                                                                                                          ttG‑CGCCAGCGAGAAATCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGActct                <  1:1148156/140‑1 (MQ=255)                                                                                                                            g‑cgcCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCttt              <  1:668690/140‑1 (MQ=255)                                                                                                                            g‑cgcCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCttt              <  1:1553190/140‑1 (MQ=255)                                                                                                                              cgcCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTa             >  2:13442/1‑140 (MQ=255)                                                                                                                                         aaTCCCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGa  <  1:1340756/140‑1 (MQ=255)                                                                                                                                            |                                                                                                                                         AGCAAACGCACTTTTAACGCCTCGTTTGCGGCAATGTTCGCCGGGATTTCAATGCGCAAACTGCTGACATACGGCGGTTCTTTCATGGTACTGGCGACTGTCAGCCACACTGCGGCCAGCATTG‑CGCCAGCGAGAAATACCCCCTGACCGTCAAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGA  >  NC_000913/445420‑445694

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