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%	159.4 / NA	42	T67A (ACC→GCC)	metJ	transcriptional repressor, S‑adenosylmethionine‑binding
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/20);  total (22/20)

ATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACT  >  NC_000913/4128050‑4128341                                                                                                                                                    |                                                                                                                                                 aTTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCaa                                                                                                                                                >  1:470800/1‑150 (MQ=255) aTTCCACCCCGCGCCGCTCTTTGTTGCTTTAGTATTCCCTCGTCTCCGGGTTAACCCCCAGCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTAACGCAGATCGGCATCATCCGGCGAAGGTTGCCCGGCaaa                                                                                                                                               <  1:380354/151‑1 (MQ=255)    ccACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGc                                                                                                                                            <  2:93392/151‑1 (MQ=255)    ccACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGc                                                                                                                                            >  2:569999/1‑151 (MQ=255)     cACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGc                                                                                                                                            >  2:612350/1‑150 (MQ=255)        tCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATATCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCCGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGc                                                                                                                                        >  2:376822/1‑151 (MQ=255)         ccGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCa                                                                                                                                       <  2:255570/151‑1 (MQ=255)               cgcTCTTTTTTGCTTTAGTATTCCCTCGTCTCCGGGTTAATCCCCATCTCACGCATGATCGCTGTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCAGCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACg                                                                                                                                 <  1:223006/151‑1 (MQ=255)                gcTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGc                                                                                                                                >  1:295331/1‑151 (MQ=255)                     tttttGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCg                                                                                                                            <  2:556880/150‑1 (MQ=255)                                 tATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGc                                                                                                               >  2:122031/1‑151 (MQ=255)                                  aTTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCt                                                                                                              <  2:348912/151‑1 (MQ=255)                                   ttCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTg                                                                                                             >  1:129516/1‑151 (MQ=255)                                         cGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGt                                                                                                        >  1:113091/1‑150 (MQ=255)                                           tctcCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGc                                                                                                     >  1:338757/1‑151 (MQ=255)                                            ctcCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCg                                                                                                    >  2:493766/1‑151 (MQ=255)                                              ccGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGc                                                                                                     <  1:377230/148‑1 (MQ=255)                                                    tAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGc                                                                                              >  1:164574/1‑149 (MQ=255)                                                    tAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCa                                                                                             <  1:248637/150‑1 (MQ=255)                                                                  cGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTgac                                                                              <  2:540346/151‑1 (MQ=255)                                                                    cATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTgacg                                                                             >  2:191652/1‑150 (MQ=255)                                                                       gagctctTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGAcgc                                                                         <  1:336102/148‑1 (MQ=255)                                                                           tctTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACgcgtgc                                                                     >  2:106/1‑151 (MQ=255)                                                                                 gCCGCTGCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCAt                                                               <  1:453426/151‑1 (MQ=255)                                                                                  ccGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATc                                                              <  1:44766/151‑1 (MQ=255)                                                                                         cGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAgg                                                       >  2:294294/1‑151 (MQ=255)                                                                                             aTTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGAttt                                                   <  1:505685/151‑1 (MQ=255)                                                                                               ttCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGAttttt                                                 >  1:234797/1‑151 (MQ=255)                                                                                                 cgtcgCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTa                                                <  2:125724/150‑1 (MQ=255)                                                                                                  gtcgCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAac                                              >  2:98520/1‑151 (MQ=255)                                                                                                   tcgCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAac                                              >  2:107542/1‑150 (MQ=255)                                                                                                          gcTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTa                                        >  1:451536/1‑149 (MQ=255)                                                                                                            tctTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAag                                      >  1:451463/1‑149 (MQ=255)                                                                                                                    cAGATCGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGaa                            <  1:254346/151‑1 (MQ=255)                                                                                                                         cGGCATCATCCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGt                       <  1:612320/151‑1 (MQ=255)                                                                                                                             atcatcCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATc                   >  1:46408/1‑151 (MQ=255)                                                                                                                               catcCGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCtt                 <  1:449609/151‑1 (MQ=255)                                                                                                                                  ccGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCttttt              <  2:429254/151‑1 (MQ=255)                                                                                                                                  ccGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCtttt               <  2:173045/150‑1 (MQ=255)                                                                                                                                   cGGCAAAGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTg             <  1:493744/151‑1 (MQ=255)                                                                                                                                         aGGTTGCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGt       >  1:151054/1‑151 (MQ=255)                                                                                                                                              gCCCGGCAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCAct  >  2:23292/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                 ATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACT  >  NC_000913/4128050‑4128341

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