Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,311,405	T→C	T44A (ACC→GCC)	apbE ←	putative thiamine‑synthetic flavin transferase lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,311,405	0	T	C	100.0%	191.2 / NA	54	T44A (ACC→GCC)	apbE	putative thiamine‑synthetic flavin transferase lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/33);  total (21/33)

CGAGTCGTTAAAGCGCATCAGCGCGGAATCTTTTTTATAGGTCGAAAGCAGCTGATCGTCGGCGTCCAGTTGGGTCTGAATCTTTTCTTTAAGTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGTACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGC  >  NC_000913/2311262‑2311552                                                                                                                                                |                                                                                                                                                    cGAGTCGTTAAAGCGCATCAGCGCGGAATCTTTTTTATAGGTCGAAAGCAGCTGATCGTCGGCGTCCAGTTGGGTCTGAATCTTTTCTTTAAGTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATa                                                                                                                                              <  1:177689/151‑1 (MQ=255)  gAGTCGTTAAAGCGCATCAGCGCGGAATCTTTTTTATAGGTCGAAAGCAGCTGATCGTCGGCGTCCAGTTGGGTCTGAATCTTTTCTTTAAGTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAg                                                                                                                                             <  1:124936/151‑1 (MQ=255)      cGTTAAAGCGCATCAGCGCGGAATCTTTTTTATAGGTCGAAAGCAGCTGATCGTCGGCGTCCAGTTGGGTCTGAATCTTTTCTTTAAGTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGtttt                                                                                                                                         <  1:452535/151‑1 (MQ=255)        ttAAAGCGCATCAGCGCGGAATCTTTTTTATAGGTCGAAAGCAGCTGATCGTCGGCGTCCAGTTGGGTCTGAATCTTGTCCTTAAGTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCGCCCATAGttttt                                                                                                                                        <  2:108264/150‑2 (MQ=255)               gcATCAGCGCGGAATCTTTTTTATAGGTCGAAAGCAGCTGATCGTCGGCGTCCAGTTGGGTCTGAATCTTTTCTTTAAGTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAg                                                                                                                                <  1:690844/151‑1 (MQ=255)               gcATCAGCGCGGAATCTTTTTTATAGGTCGAAAGCAGCTGATCGTCGGCGTCCAGTTGGGTCTGAATCTTTTCTTTAAGTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAg                                                                                                                                <  1:588308/151‑1 (MQ=255)                cATCAGCGCGGAATCTTTTTTATAGGTCGAAAGCAGCTGATCGTCGGCGTCCAGTTGGGTCTGAATCTTTTCTTTAAGTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCa                                                                                                                                  >  1:444377/1‑148 (MQ=255)                    aGCGCGGAATCTTTTTTATAGGTCGAAAGCAGCTGATCGTCGGCGTCCAGTTGGGTCTGAATCTTTTCTTTAAGTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACgg                                                                                                                           >  2:261095/1‑151 (MQ=255)                    aGCGCGGAATCTTTTTTATAGGTCGAAAGCAGCTGATCGTCGGCGTCCAGTTGGGTCTGAATCTTTTCTTTAAGTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACgg                                                                                                                           >  1:380535/1‑151 (MQ=255)                         ggAATCTTTTTTATAGGTCGAAAGCAGCTGATCGTCGGCGTCCAGTTGGGTCTGAATCTTTTCTTTAAGTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACt                                                                                                                      >  1:848915/1‑151 (MQ=255)                            aTCTTTTTTATAGGTCGAAAGCAGCTGATCGTCGGCGTCCAGTTGGGTCTGAATCTTTTCTTTAAGTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTc                                                                                                                    >  2:44337/1‑150 (MQ=255)                             tCTTTTTTATAGGTCGAAAGCAGCTGATCGTCGGCGTCCAGTTGGGTCTGAATCTTTTCTTTAAGTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTc                                                                                                                    >  2:464611/1‑149 (MQ=255)                                tttttATAGGTCGAAAGCAGCTGATCGTCGGCGTCCAGTTGGGTCTGAATCTTTTCTTTAAGTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAg                                                                                                                  <  2:348991/148‑1 (MQ=255)                                tttttATAGGTCGAAAGCAGCTGATCGTCGGCGTCCAGTTGGGTCTGAATCTTTTCTTTAAGTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGt                                                                                                                 <  2:21141/149‑1 (MQ=255)                                tttttATAGGTCGAAAGCAGCTGATCGTCGGCGTCCAGTTGGGTCTGAATCTTTTCTTTAAGTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTa                                                                                                                <  1:216806/150‑1 (MQ=255)                                       aGGTCGAAAGCAGCTGATCGTCGGCGTCCAGTTGGGTCTGAATCTTTTCTTTAAGTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTggg                                                                                                         <  2:38285/150‑1 (MQ=255)                                             aaaGCAGCTGATCGTCGGCGTCCAGTTGGGTCTGAATCTTTTCTTTAAGTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGc                                                                                                  <  1:363241/151‑1 (MQ=255)                                                    cTGATCGTCGGCGTCCAGTTGGGTCTGAATCTTTTCTTTAAGTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTtgt                                                                                           >  1:479485/1‑151 (MQ=255)                                                      gATCGTCGGCGTCCAGTTGGGTCTGAATCTTTTCTTTAAGTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTgg                                                                                         <  2:405285/151‑1 (MQ=255)                                                         cgtcgGCGTCCAGTTGGGTCTGAATCTTTTCTTTAAGTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGttt                                                                                      <  1:474587/151‑1 (MQ=255)                                                          gtcgGCGTCCAGTTGGGTCTGAATCTTTTCTTTAAGTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGtttt                                                                                     >  1:901981/1‑151 (MQ=255)                                                            cgGCGTCCAGTTGGGTCTGAATCTTTTCTTTAAGTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTg                                                                                   <  1:872449/151‑1 (MQ=255)                                                                gTCCAGTTGGGTCTGAATCTTTTCTTTAAGTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCg                                                                               <  2:344839/151‑1 (MQ=255)                                                                 tCCAGTTGGGTCTGAATCTTTTCTTTAAGTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGc                                                                              >  1:355053/1‑151 (MQ=255)                                                                       tGGGTCTGAATCTTTTCTTTAAGTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGcaacca                                                                         >  1:18002/1‑150 (MQ=255)                                                                            cTGAATCTTTTCTTTAAGTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACcaacaa                                                                     >  2:276403/1‑149 (MQ=255)                                                                              gAATCTTTTCTTTAAGTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAaagaa                                                                 <  2:887066/151‑1 (MQ=255)                                                                              gAATCTTTTCTTTAAGTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAaaga                                                                  >  2:791537/1‑150 (MQ=255)                                                                                    tttctttAAGTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAagaaga                                                               <  1:337846/147‑1 (MQ=255)                                                                                    tttctttAAGTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAagaaga                                                               <  1:449187/147‑1 (MQ=255)                                                                                     ttctttAAGTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGcgcg                                                          >  1:623400/1‑151 (MQ=255)                                                                                      tctttAAGTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGcgcg                                                          <  1:464423/150‑1 (MQ=255)                                                                                       ctttAAGTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCgg                                                         >  2:425281/1‑150 (MQ=255)                                                                                             gTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGcagccag                                                  <  1:791588/151‑1 (MQ=255)                                                                                               tCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGcagccagc                                                 <  2:683898/150‑1 (MQ=255)                                                                                                    gcTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGcc                                           >  2:877062/1‑151 (MQ=255)                                                                                                      tGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCAc                                         >  2:225167/1‑151 (MQ=255)                                                                                                       gCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCTCGGCAGCCAGCAGTGCCa                                          <  1:640008/149‑1 (MQ=255)                                                                                                        cGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCAcgc                                       <  1:367636/151‑1 (MQ=255)                                                                                                         gTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCAcgcg                                      <  2:265959/151‑1 (MQ=255)                                                                                                             ggCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGt                                   >  2:437351/1‑150 (MQ=255)                                                                                                              gCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTaa                                 <  1:57629/151‑1 (MQ=255)                                                                                                                  cAATGCCCGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCt                             <  1:461333/151‑1 (MQ=255)                                                                                                                       cccGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTAttt                        >  2:850153/1‑151 (MQ=255)                                                                                                                        ccGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTc                       <  2:462259/151‑1 (MQ=255)                                                                                                                         cGGGATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTcc                      >  1:746031/1‑151 (MQ=255)                                                                                                                           ggATGCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCAt                    <  1:51449/151‑1 (MQ=255)                                                                                                                               gCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACg                 <  1:112481/150‑1 (MQ=255)                                                                                                                               gCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACg                 <  1:291438/150‑1 (MQ=255)                                                                                                                               gCTGGCACGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGt                >  1:139399/1‑151 (MQ=255)                                                                                                                                      cGCCAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTGCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGACACGCGGGTAAAGCTTATTTCCATACGTGATTAt          <  2:101149/150‑1 (MQ=255)                                                                                                                                         cAGAAGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTc      >  1:661577/1‑151 (MQ=255)                                                                                                                                            aaGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATg   <  2:300219/151‑1 (MQ=255)                                                                                                                                             aGGCACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTAACCTCATGc  <  1:336355/151‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                    CGAGTCGTTAAAGCGCATCAGCGCGGAATCTTTTTTATAGGTCGAAAGCAGCTGATCGTCGGCGTCCAGTTGGGTCTGAATCTTTTCTTTAAGTTCGGCGCTGCGTTTGGCGTCAATGCCCGGGATGCTGGCACGCCAGAAGGTACCCATAGTTTTGCCTTCAAGAACGGTAACTTCAGTAGCGTGGGTTTTGGCGGGTTGTGGTTTTTGATCGCAACCAACAAAGAAGAGCGCGGCAGCCAGCAGTGCCACGCGGGTAAAGCTTATTTCCATACGTGATTATCCTCATGC  >  NC_000913/2311262‑2311552

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