Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,128,169	T→G	D76A (GAT→GCT)	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,169	0	T	G	100.0%	195.1 / NA	49	D76A (GAT→GCT)	metJ	transcriptional repressor, S‑adenosylmethionine‑binding
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (27/22);  total (27/22)

CCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAA  >  NC_000913/4128025‑4128315                                                                                                                                                 |                                                                                                                                                   ccTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGcat                                                                                                                                              >  2:527420/1‑151 (MQ=255)  cTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGcatc                                                                                                                                             >  1:297743/1‑151 (MQ=255)   tGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGcatc                                                                                                                                             >  2:355774/1‑150 (MQ=255)    gATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGcatca                                                                                                                                            <  1:555966/150‑1 (MQ=255)       aaGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCg                                                                                                                                        >  2:324696/1‑151 (MQ=255)        aGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCgg                                                                                                                                       >  2:16218/1‑151 (MQ=255)                gcATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGtt                                                                                                                               >  1:271505/1‑151 (MQ=255)                          aTTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTaaa                                                                                                                     <  1:311918/151‑1 (MQ=255)                            tCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAg                                                                                                                    <  2:234743/150‑1 (MQ=255)                             ccACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGc                                                                                                                  >  1:741472/1‑151 (MQ=255)                                   cgcgCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTGATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATg                                                                                                               >  2:599601/1‑148 (MQ=255)                                     cgccgcTCTTTTTTGCTTTAGTATTCCCACGTCTCAGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGaa                                                                                                          <  1:158073/151‑1 (MQ=255)                                          ctctTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCt                                                                                                     >  2:644388/1‑151 (MQ=255)                                              tttttGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCtt                                                                                                    <  1:475356/148‑1 (MQ=255)                                                tttGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGca                                                                                                <  2:438975/150‑1 (MQ=255)                                                tttGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGca                                                                                                <  2:640316/150‑1 (MQ=255)                                                        aGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTc                                                                                       <  2:117875/151‑1 (MQ=255)                                                              cccACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGtt                                                                                 >  1:5473/1‑151 (MQ=255)                                                                cACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTgg                                                                               <  1:453231/151‑1 (MQ=255)                                                                  cGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTa                                                                             <  1:313710/151‑1 (MQ=255)                                                                       ccGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTg                                                                        >  2:185582/1‑151 (MQ=255)                                                                               aTCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGgt                                                                >  2:436536/1‑151 (MQ=255)                                                                                 ccccATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGgtt                                                               <  1:495193/150‑1 (MQ=255)                                                                                  cccATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGgttgt                                                             >  2:761377/1‑151 (MQ=255)                                                                                        tcACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCt                                                       >  2:443917/1‑151 (MQ=255)                                                                                          aCGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTga                                                     <  2:607861/151‑1 (MQ=255)                                                                                           cGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTgac                                                    <  1:370106/151‑1 (MQ=255)                                                                                           cGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTgac                                                    >  2:293638/1‑151 (MQ=255)                                                                                                 aTCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGAcgcg                                              >  1:401590/1‑151 (MQ=255)                                                                                                  tctcTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACgcgt                                             >  1:252064/1‑151 (MQ=255)                                                                                                   ctctTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACgcgtg                                            >  2:312721/1‑151 (MQ=255)                                                                                                     ctTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACgcgt                                             <  1:795535/148‑1 (MQ=255)                                                                                                       tttGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACgcgtgc                                           <  2:482969/148‑1 (MQ=255)                                                                                                         tGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTc                                       >  2:775714/1‑150 (MQ=255)                                                                                                              cTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGt                                 >  1:276938/1‑151 (MQ=255)                                                                                                              cTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGt                                 >  2:505444/1‑151 (MQ=255)                                                                                                                   gggATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGCCGACGCGTGCGTTCATCGGTGAgg                             >  1:750648/1‑150 (MQ=255)                                                                                                                   gggATTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGa                            <  2:15396/151‑1 (MQ=255)                                                                                                                      aTTTCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGAttt                         <  1:260987/151‑1 (MQ=255)                                                                                                                       tttCGTCGCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGAtttt                        >  1:221645/1‑151 (MQ=255)                                                                                                                             cgCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAacac                  <  2:351958/151‑1 (MQ=255)                                                                                                                              gCTGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACAcc                 >  2:705626/1‑151 (MQ=255)                                                                                                                                tGCGCTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCtt               >  1:609954/1‑151 (MQ=255)                                                                                                                                  cgcTCTTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCtt               >  2:629968/1‑149 (MQ=255)                                                                                                                                     tctTTACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAag            >  2:827853/1‑149 (MQ=255)                                                                                                                                         tACGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGaa       <  2:463499/150‑1 (MQ=255)                                                                                                                                          aCGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATg     <  1:62926/151‑1 (MQ=255)                                                                                                                                           cGCAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTTACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATgg    <  1:84972/151‑1 (MQ=255)                                                                                                                                             cAGAGCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGaa  <  2:388078/151‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                   CCTGATAAGCGTAGCGCATCAGGCGATTCCACTCCGCGCCGCTCTTTTTTGCTTTAGTATTCCCACGTCTCCGGGTTAATCCCCATCTCACGCATGATCTCTTTTGCCGCTTCCGGGATTTCGTCGCTGCGCTCTTTACGCAGATCGGCATCATCCGGCAAAGGTTGCCCGGTAAAGGCATGCAGAAACGCTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAA  >  NC_000913/4128025‑4128315

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