Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,128,361	T→C	Y12C (TAC→TGC)	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,361	0	T	C	100.0%	220.6 / NA	57	Y12C (TAC→TGC)	metJ	transcriptional repressor, S‑adenosylmethionine‑binding
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (34/23);  total (34/23)

CTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAA  >  NC_000913/4128215‑4128508                                                                                                                                                   |                                                                                                                                                    cTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATgg                                                                                                                                                 >  2:219388/1‑151 (MQ=255)  ttCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATggg                                                                                                                                                >  1:395669/1‑151 (MQ=255)   tCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCGTCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATgg                                                                                                                                                 <  1:113091/149‑1 (MQ=255)   tCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACTCCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATggg                                                                                                                                                <  2:59146/150‑1 (MQ=255)      cacaGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTACGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTg                                                                                                                                             >  1:281037/1‑150 (MQ=255)          gcagcTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGatat                                                                                                                                         >  2:495335/1‑150 (MQ=255)                cGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTcgcc                                                                                                                                  <  1:270577/151‑1 (MQ=255)                cGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTcgc                                                                                                                                   <  2:288644/150‑1 (MQ=255)                            cGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGc                                                                                                                      <  2:610002/151‑1 (MQ=255)                              tGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGcc                                                                                                                     >  2:680410/1‑150 (MQ=255)                                      gttgttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATAc                                                                                                            <  1:99045/151‑1 (MQ=255)                                         gttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTa                                                                                                         <  2:744797/151‑1 (MQ=255)                                          ttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTa                                                                                                         >  2:439551/1‑150 (MQ=255)                                           tCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTaa                                                                                                        <  2:612776/150‑1 (MQ=255)                                             aCCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATcc                                                                                                     >  1:235139/1‑151 (MQ=255)                                              ccTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCct                                                                                                    <  2:191909/151‑1 (MQ=255)                                              ccTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCct                                                                                                    <  1:715499/151‑1 (MQ=255)                                               cTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCctc                                                                                                   >  1:331736/1‑151 (MQ=255)                                                tGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCctct                                                                                                  >  2:782605/1‑151 (MQ=255)                                                 gacgacGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCctc                                                                                                   <  2:531673/149‑1 (MQ=255)                                                  acgacgCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTc                                                                                                <  1:494183/151‑1 (MQ=255)                                                    gacgCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGt                                                                                              >  2:536675/1‑151 (MQ=255)                                                       gcgtgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCaa                                                                                           <  2:426768/151‑1 (MQ=255)                                                         gtgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAAt                                                                                          >  1:555552/1‑150 (MQ=255)                                                              ttCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTg                                                                                    >  1:763582/1‑151 (MQ=255)                                                                cATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTga                                                                                   <  1:457144/150‑1 (MQ=255)                                                                cATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTg                                                                                    >  1:549035/1‑149 (MQ=255)                                                                 aTCGGGGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCTCTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTg                                                                                    >  1:757429/1‑148 (MQ=255)                                                                       gAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccagac                                                                           <  1:520032/151‑1 (MQ=255)                                                                           aTTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACcac                                                                        >  2:106945/1‑150 (MQ=255)                                                                            tttttAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACcacat                                                                      >  1:714184/1‑150 (MQ=255)                                                                               ttAACACCTTAAGAGGAATGGAAAACGTAATCTTTTTGACGTGTTCCCTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTACATTCAGCCATGAGATACTTAATCCCCTGCGTCAATAAATTGAGACCAGACtacacaa                                                                    <  2:647048/150‑2 (MQ=255)                                                                                 aaCACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGtt                                                                  <  1:1987/150‑1 (MQ=255)                                                                                 aaCACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTg                                                                 >  1:133106/1‑151 (MQ=255)                                                                                    acCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGAtg                                                              <  1:219408/151‑1 (MQ=255)                                                                                     ccTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGAtgt                                                             >  2:574883/1‑151 (MQ=255)                                                                                         aaGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGt                                                         >  2:468570/1‑151 (MQ=255)                                                                                            agGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACt                                                      >  2:680163/1‑151 (MQ=255)                                                                                             ggAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTg                                                     >  2:686159/1‑151 (MQ=255)                                                                                               aaTGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGAc                                                   >  1:283068/1‑151 (MQ=255)                                                                                                 tGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCg                                                 >  1:790061/1‑151 (MQ=255)                                                                                                      aCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTaaa                                             <  2:705296/150‑1 (MQ=255)                                                                                                         gTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCg                                         >  1:199185/1‑151 (MQ=255)                                                                                                          tAATCTTTTTGACTTGTTCTCTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCg                                         <  2:618853/150‑1 (MQ=255)                                                                                                          tAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGc                                        >  1:294821/1‑151 (MQ=255)                                                                                                               tttttGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGt                                   >  1:60862/1‑151 (MQ=255)                                                                                                                ttttGACTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAg                                    <  2:567292/149‑1 (MQ=255)                                                                                                                      cTTGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGt                            >  1:632143/1‑151 (MQ=255)                                                                                                                       ttGTTCACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGt                            >  2:814629/1‑150 (MQ=255)                                                                                                                            cACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAgg                       <  2:499253/150‑1 (MQ=255)                                                                                                                            cACTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGc                      >  1:604835/1‑151 (MQ=255)                                                                                                                             aCTCTTCTTGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCg                     >  2:159628/1‑151 (MQ=255)                                                                                                                                cttcttACCGTGCTCAGAGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTAGGTACTGACCGTAACCCCGCATAGTTTACCGAACAGGCGtt                   >  2:644030/1‑150 (MQ=255)                                                                                                                                 ttcttGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTAc                 >  1:286561/1‑151 (MQ=255)                                                                                                                                  tcttGCCGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTa                  <  1:519906/149‑1 (MQ=255)                                                                                                                                       ccGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGAc           >  1:323584/1‑151 (MQ=255)                                                                                                                                       ccGTGCTCAGCGCATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGAc           >  2:673634/1‑151 (MQ=255)                                                                                                                                                gcgcATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTaa  <  2:780478/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CTTCGCACAGCAGCTCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAA  >  NC_000913/4128215‑4128508

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: ‑