Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,128,197	T→C	T67A (ACC→GCC)	metJ ←	transcriptional repressor, S‑adenosylmethionine‑binding

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,128,197	0	T	C	100.0%	194.0 / NA	52	T67A (ACC→GCC)	metJ	transcriptional repressor, S‑adenosylmethionine‑binding
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (26/26);  total (26/26)

TCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTC  >  NC_000913/4128052‑4128342                                                                                                                                                  |                                                                                                                                                  tCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAgg                                                                                                                                              >  2:331884/1‑151 (MQ=255) tCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAg                                                                                                                                               <  1:462494/150‑1 (MQ=255)    aCTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGc                                                                                                                                             <  2:89038/149‑1 (MQ=255)           gccgcTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGaa                                                                                                                                     >  1:102800/1‑150 (MQ=255)             cgcTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACg                                                                                                                                  >  1:537223/1‑151 (MQ=255)                 ctTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCt                                                                                                                                <  1:314226/149‑1 (MQ=255)                   tttttGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGc                                                                                                                            >  1:489951/1‑151 (MQ=255)                           ttAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACagcagc                                                                                                                    >  1:400069/1‑151 (MQ=255)                             aGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTc                                                                                                                  >  1:592419/1‑151 (MQ=255)                                aTTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCt                                                                                                               >  1:726193/1‑151 (MQ=255)                                  tCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTg                                                                                                              <  2:341979/150‑1 (MQ=255)                                  tCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGt                                                                                                             >  2:241494/1‑151 (MQ=255)                                   cccACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGtt                                                                                                            >  1:241830/1‑151 (MQ=255)                                          ctcCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCg                                                                                                     >  1:316000/1‑151 (MQ=255)                                          ctcCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCg                                                                                                     >  2:295476/1‑151 (MQ=255)                                           tcCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGt                                                                                                    >  2:351434/1‑151 (MQ=255)                                                  tAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCa                                                                                              <  1:450038/150‑1 (MQ=255)                                                   aaTCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAgg                                                                                            <  1:513902/151‑1 (MQ=255)                                                   aaTCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAg                                                                                             <  1:518525/150‑1 (MQ=255)                                                           tctcACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCAc                                                                                    >  1:758901/1‑151 (MQ=255)                                                             tcACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCt                                                                                  >  2:208006/1‑151 (MQ=255)                                                                          ctTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGAcgc                                                                          <  1:498509/146‑1 (MQ=255)                                                                            tttGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTGGCGTGACGCAGGTTGTTCACCTGACGACgcgtgc                                                                      <  2:522261/148‑1 (MQ=255)                                                                                ccGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATc                                                               <  2:686595/151‑1 (MQ=255)                                                                                 cGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCg                                                              <  2:762381/151‑1 (MQ=255)                                                                                 cGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCg                                                              >  1:406539/1‑151 (MQ=255)                                                                                    ttCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTg                                                           >  2:356897/1‑151 (MQ=255)                                                                                      ccGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGa                                                          <  2:570719/150‑1 (MQ=255)                                                                                      ccGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGa                                                          <  2:411403/150‑1 (MQ=255)                                                                                       cGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAgg                                                        >  2:266440/1‑151 (MQ=255)                                                                                           aTTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGAttt                                                    <  2:489928/151‑1 (MQ=255)                                                                                             ttCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGAtttt                                                   <  2:452363/150‑1 (MQ=255)                                                                                                  cgCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAaca                                              <  1:170722/150‑1 (MQ=255)                                                                                                   gCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACAcc                                            >  2:694037/1‑151 (MQ=255)                                                                                                      gcgcTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCtt                                          >  2:97551/1‑150 (MQ=255)                                                                                                       cgcTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTaa                                        <  1:151563/151‑1 (MQ=255)                                                                                                        gcTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAag                                       <  1:69336/151‑1 (MQ=255)                                                                                                         ctctTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTa                                         >  2:412269/1‑148 (MQ=255)                                                                                                          tctTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAag                                       >  1:9825/1‑149 (MQ=255)                                                                                                          tctTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAag                                       >  2:743620/1‑149 (MQ=255)                                                                                                             ttACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAgg                                    <  2:165341/149‑1 (MQ=255)                                                                                                             ttACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGa                                   <  2:515120/150‑1 (MQ=255)                                                                                                             ttACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGa                                   <  1:331898/150‑1 (MQ=255)                                                                                                               aCGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATg                                >  1:242250/1‑151 (MQ=255)                                                                                                                        ggCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGt                        >  2:135550/1‑150 (MQ=255)                                                                                                                           atcatcCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATc                    <  2:667130/151‑1 (MQ=255)                                                                                                                             catcCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCt                   <  1:102499/150‑1 (MQ=255)                                                                                                                                  ggCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCGCAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGa             <  2:453823/151‑1 (MQ=255)                                                                                                                                      aaGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTg         <  1:589118/151‑1 (MQ=255)                                                                                                                                       aGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGt        <  2:689371/151‑1 (MQ=255)                                                                                                                                       aGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGt        >  1:685844/1‑151 (MQ=255)                                                                                                                                             cccGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCActc  >  2:28812/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                  TCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTC  >  NC_000913/4128052‑4128342

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