Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,128,386	A→C	W4G (TGG→GGG)	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,386	0	A	C	100.0%	229.1 / NA	56	W4G (TGG→GGG)	metJ	transcriptional repressor, S‑adenosylmethionine‑binding
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (28/28);  total (28/28)

GTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCAT  >  NC_000913/4128238‑4128533                                                                                                                                                     |                                                                                                                                                    gTAGCGTGACGCAGGTTGTTCACCTGACGACGTGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCtt                                                                                                                                                   <  1:516052/151‑1 (MQ=255)    gCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAg                                                                                                                                                <  2:92088/151‑1 (MQ=255)                      aCCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATcc                                                                                                                              >  1:635556/1‑151 (MQ=255)                      aCCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATcc                                                                                                                              <  1:86845/151‑1 (MQ=255)                           acgacgCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTc                                                                                                                         <  2:475118/151‑1 (MQ=255)                           acgacgCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTc                                                                                                                         <  1:544560/151‑1 (MQ=255)                                   tgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATa                                                                                                                  >  1:189443/1‑150 (MQ=255)                                      gtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAAtt                                                                                                              <  2:62679/151‑1 (MQ=255)                                         cATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTgag                                                                                                           <  2:99753/151‑1 (MQ=255)                                         cATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTga                                                                                                            <  2:19034/150‑1 (MQ=255)                                           tCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagac                                                                                                         >  1:3938/1‑151 (MQ=255)                                            cGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagac                                                                                                         <  2:19259/150‑1 (MQ=255)                                             ggTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagacca                                                                                                       <  2:261959/151‑1 (MQ=255)                                             ggTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagacc                                                                                                        >  2:26553/1‑150 (MQ=255)                                                      ttttAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACcacac                                                                                               >  2:181953/1‑150 (MQ=255)                                                      ttttAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACcacac                                                                                               >  2:281394/1‑150 (MQ=255)                                                           acacCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGa                                                                                         <  1:367922/151‑1 (MQ=255)                                                                ttAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTgg                                                                                    >  2:298937/1‑151 (MQ=255)                                                                 tAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTgg                                                                                    <  1:623398/150‑1 (MQ=255)                                                                     agGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCctc                                                                                                                            >  2:126827/1‑106 (MQ=255)                                                                     agGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCctc                                                                                                                            <  1:126808/106‑1 (MQ=255)                                                                     agGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACt                                                                               >  1:316946/1‑151 (MQ=255)                                                                        aaTGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGAc                                                                            >  1:678262/1‑151 (MQ=255)                                                                          tGGAAACCGTAATCTTTTTGACTTGTTCACTCTTATTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGAcc                                                                           >  1:405384/1‑150 (MQ=255)                                                                           ggAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGt                                                                         <  1:434324/151‑1 (MQ=255)                                                                           ggAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGt                                                                         >  2:289702/1‑151 (MQ=255)                                                                             aaaCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTaa                                                                       <  1:98314/151‑1 (MQ=255)                                                                                   tAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGc                                                                 >  1:255357/1‑151 (MQ=255)                                                                                        tttttGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGt                                                            >  1:397107/1‑151 (MQ=255)                                                                                         ttttGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCAt                                                               <  1:333679/147‑1 (MQ=255)                                                                                           ttGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGttt                                                          <  1:436947/150‑1 (MQ=255)                                                                                            tGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTAc                                                        >  1:139407/1‑151 (MQ=255)                                                                                            tGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTAc                                                        >  1:369112/1‑151 (MQ=255)                                                                                                 tGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTa                                                    <  2:378521/150‑1 (MQ=255)                                                                                                   ttCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTAc                                                   >  2:338244/1‑149 (MQ=255)                                                                                                    tCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAgg                                                >  1:344475/1‑151 (MQ=255)                                                                                                     cACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAgg                                                <  1:3081/150‑1 (MQ=255)                                                                                                           tcttGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTAc                                          <  2:410578/150‑1 (MQ=255)                                                                                                              tGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTAcc                                         <  2:283281/148‑1 (MQ=255)                                                                                                              tGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTAcc                                         <  1:466304/148‑1 (MQ=255)                                                                                                              tGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGt                                       <  1:237402/150‑1 (MQ=255)                                                                                                                ccGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGAc                                    <  1:53872/151‑1 (MQ=255)                                                                                                                    gCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCg                                >  2:89636/1‑151 (MQ=255)                                                                                                                         gCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTa                            >  2:304730/1‑150 (MQ=255)                                                                                                                           gTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATg                         >  1:66677/1‑151 (MQ=255)                                                                                                                           gTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATg                         >  2:476472/1‑151 (MQ=255)                                                                                                                              tGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGTATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCAc                      >  2:145414/1‑151 (MQ=255)                                                                                                                               gggCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCAc                      <  1:446773/150‑1 (MQ=255)                                                                                                                                ggCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCt                    >  2:497586/1‑151 (MQ=255)                                                                                                                                    gATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGt                  <  2:349946/149‑1 (MQ=255)                                                                                                                                    gATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTc                 <  1:548196/150‑1 (MQ=255)                                                                                                                                    gATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCg                <  1:387743/151‑1 (MQ=255)                                                                                                                                      tatatTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGc              >  1:588492/1‑151 (MQ=255)                                                                                                                                              ccgcTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAAt      >  2:415187/1‑151 (MQ=255)                                                                                                                                               cgcTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATg     >  1:60809/1‑151 (MQ=255)                                                                                                                                                  tCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCat  >  2:47432/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCAT  >  NC_000913/4128238‑4128533

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