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%	188.8 / NA	48	W4G (TGG→GGG)	metJ	transcriptional repressor, S‑adenosylmethionine‑binding
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (19/29);  total (19/29)

AGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCAT  >  NC_000913/4128240‑4128533                                                                                                                                                   |                                                                                                                                                    aGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCa                                                                                                                                                 <  1:296217/151‑1 (MQ=255)  gCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAg                                                                                                                                                <  1:368114/151‑1 (MQ=255)    gTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGcc                                                                                                                                              >  2:229161/1‑151 (MQ=255)      gACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCa                                                                                                                                             <  2:642801/150‑1 (MQ=255)             gttgttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATAc                                                                                                                                     <  1:206871/151‑1 (MQ=255)                   cACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAAt                                                                                                                                <  2:606159/150‑1 (MQ=255)                      cTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCctc                                                                                                                            >  1:764889/1‑151 (MQ=255)                              gcgtgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCaa                                                                                                                    <  1:368446/151‑1 (MQ=255)                                gtgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAAt                                                                                                                   >  1:234473/1‑150 (MQ=255)                                gtgcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAAt                                                                                                                   >  2:457169/1‑150 (MQ=255)                                  gcgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATaa                                                                                                                 >  1:100824/1‑150 (MQ=255)                                  gcgtTCATCGGTGAGGATTTTAAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATaaa                                                                                                                <  2:602369/151‑1 (MQ=255)                                   cgtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAAt                                                                                                               >  2:599027/1‑151 (MQ=255)                                    gtTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAAtt                                                                                                              <  2:252810/151‑1 (MQ=255)                                       cATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTga                                                                                                            <  2:476442/150‑1 (MQ=255)                                          cGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagacc                                                                                                        >  1:376278/1‑151 (MQ=255)                                          cGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagacc                                                                                                        <  2:656302/151‑1 (MQ=255)                                              gAGGATTTTTAACCCCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTgagacgagac                                                                                                    <  1:260548/151‑1 (MQ=255)                                                ggATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAgaccagacca                                                                                                  <  1:795738/151‑1 (MQ=255)                                                 gATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACcac                                                                                                 >  2:662313/1‑151 (MQ=255)                                                  aTTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACcac                                                                                                 >  1:604469/1‑150 (MQ=255)                                                     tttAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACcacaca                                                                                              <  1:269972/150‑1 (MQ=255)                                                       tAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGtt                                                                                           >  2:817576/1‑151 (MQ=255)                                                        aaCACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTg                                                                                          <  2:470861/151‑1 (MQ=255)                                                            ccTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGAtg                                                                                       <  1:8341/150‑1 (MQ=255)                                                                 agagGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGt                                                                                  >  2:286113/1‑150 (MQ=255)                                                                  gagGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTAc                                                                                >  2:143474/1‑151 (MQ=255)                                                                   agGAATTGAAACCGTAATCTGTTTGACTTGTTCACTCTTCTTGCCGTACTCAGCGTATGGGCTGATATATTCGCCGCTCCCTACAGCCATGAGATACTTAATCCTCGTCGTCAATAAATTGAGACCAGACCACACAGTTTATGTGGGTAcg                                                                               <  1:256502/151‑2 (MQ=255)                                                                   agGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACt                                                                               <  1:636934/151‑1 (MQ=255)                                                                       aTGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGAcc                                                                           <  2:444230/151‑1 (MQ=255)                                                                                     cTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCAt                                                               >  1:366303/1‑149 (MQ=255)                                                                                       ttttGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATa                                                              <  2:382666/148‑1 (MQ=255)                                                                                         ttGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGttt                                                          >  1:247716/1‑150 (MQ=255)                                                                                             cTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGt                                                     >  1:254048/1‑151 (MQ=255)                                                                                                  tCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAgg                                                >  1:647370/1‑151 (MQ=255)                                                                                                         tcttGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGtt                                            <  1:74133/148‑1 (MQ=255)                                                                                                            tGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCg                                        <  1:669439/149‑1 (MQ=255)                                                                                                             gCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTg                                      >  2:854975/1‑150 (MQ=255)                                                                                                             gCCGTGCTCAGCGTATGGGCGGATATATTCTCCGCTCCCTTCAGCCATGAGATGCTTAATCCTCGTCGTCAATAAATGGAGACCAGGCCACACAGTTGATGTGGGTACTGACCGTAAACCGGCATAGGTGACCGTACAGGCGTTACCGGGa                                     <  1:705711/151‑1 (MQ=255)                                                                                                                     cAGCGTATGGGCTGATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCgt                               <  1:275174/149‑1 (MQ=255)                                                                                                                                  gATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGt                  <  1:40659/149‑1 (MQ=255)                                                                                                                                  gATATATTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTc                 <  2:764966/150‑1 (MQ=255)                                                                                                                                    tatatTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGc              >  1:155309/1‑151 (MQ=255)                                                                                                                                     atatTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACCCAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGc              <  1:859841/150‑1 (MQ=255)                                                                                                                                       atTCTCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAGTGCACCTGTCGGCGt            <  1:292565/150‑1 (MQ=255)                                                                                                                                       atTCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGt            <  1:598961/150‑1 (MQ=255)                                                                                                                                         tCGCCGCTCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGa          <  2:273592/150‑1 (MQ=255)                                                                                                                                                tCCCTTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCat  >  1:478572/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                    AGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGACCGTAAACCCGCATAGTTTACCGTACAGGCGTTACCGTGACATCGTGTAATGCACCTGTCGGCGTGATAATGCAT  >  NC_000913/4128240‑4128533

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