Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,223,638	T→C	intergenic (‑10/‑190)	iclR ← / → metH	DNA‑binding transcriptional repressor IclR/cobalamin‑dependent methionine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,223,638	0	T	C	100.0%	109.6 / NA	33	intergenic (‑10/‑190)	iclR/metH	DNA‑binding transcriptional repressor IclR/cobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (10/23);  total (10/23)

AGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCT  >  NC_000913/4223393‑4223877                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                cgTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGTGTAACCCGGCTTGTTGCGCCAGTTCCGTTAGTGCCACACTGCCTTTGGATTCGGCAATCCCCTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCTCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGGGTCGCGGGAATGGGTGCGACCATGACAGTCTCcttttt                                                                                                                                                                                                                                            <  2:502740/250‑1 (MQ=255)              gacccccGAACCCCGGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTGGTGGAAGTGGGTACCCCGGCTTGTTGCGCCATGTCCGGGAGTGCCCCACCGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTTTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTgg                                                                                                                                                                                                                                 <  2:369443/243‑1 (MQ=255)                cGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAt                                                                                                                                                                                                                             <  1:219689/251‑1 (MQ=255)                  cACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAtca                                                                                                                                                                                                                           <  1:123245/251‑1 (MQ=255)                                gTTGCATGGTGATTAGCGGGCGGGGGGTCGTGGAATTGG‑TAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTAtt                                                                                                                                                                                                             <  2:14831/250‑1 (MQ=255)                                                     ggtggtTCGTGGAATTGTGTATCCCTGCTGGTTGCGCCAGTTCCGAGAGTGCCACACTCCCATTGGATTCGGCAATCCACTCCATTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGTTGGCAACGGCGGGTTGTCTTCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCCTTTTTATTGTTATCGAATGCAATTGTTac                                                                                                                                                                                        <  2:433757/251‑1 (MQ=255)                                                               ggAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGAt                                                                                                                                                                               <  1:137777/250‑1 (MQ=255)                                                                gAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTg                                                                                                                                                                                                            <  2:142204/220‑1 (MQ=255)                                                                     gggTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGCTTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTTCCGCCTTTCGCGGGGATGGGGGCGGCCCTGAACGTCTCCCTTTTTTTTATCGGGGAAATCATTTTTATTTTTTTTGTTTGCTAATGCAATaatttattgactggttcgatgg                                                                                                                                                                        >  2:411807/1‑230 (MQ=255)                                                                                        gcCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGGTTTGCGGGGAAGGGGGCGACCATGCAAGACTCCTTTTTTTGTGTCGTGGGAATCATTTTTATTTTTTTTTTTTGGTTATGCAATAATTAATGAACTGATCCGATGGGTTAATGtttaaacaattt                                                                                                                                                     >  1:487549/1‑241 (MQ=255)                                                                                                               gCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGGATGGGGGCGACCCTGACAGTCTCCTTTTTCTTTTTCGTGGGAATCATTTTTATTTTTATTTTTTGCTTATGCAATAGTTACTGAAATGATCCGATGAGTTTATGTTGAACAAATATCATGTTTCGTGGGGGTCGCtttt                                                                                                                              >  1:339106/1‑251 (MQ=255)                                                                                                                ccATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTa                                                                                                                             <  1:84768/251‑1 (MQ=255)                                                                                                                           ggCAATCCCCGCCAGTAATTGCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTGTTCTGCCGCGTGTCGGTGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTGATCTTTAGCTAATGCAATAGTTGCTGAACTGATCCTATGAGTTAAGTTTGAGCTAATCGCATGTTGCTTGGTGGTCGCGGTTACCATAGATGCg                                                                                                                  <  2:320456/251‑1 (MQ=255)                                                                                                                           ggCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGGTTCGCGGGGATGGGTGCGACCATGACAGTCTCCTTTTTTTGTATCGTGGGAAACATTTTTATTTTTTTTGTTAGATAATGCAATAGTTACTGAACTGATCCGATGAGGTAATTTTGAACAAAACACATGTTGGGTGGTGGTCGCTTTTATCACAGATGCg                                                                                                                  >  1:291747/1‑251 (MQ=255)                                                                                                                           ggCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACATGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGTTTCGCGGGGATGGGGGCGACCATAACAGTCTCCTTTTTTTGTATTGTGGGAATCATTTTTCTTTTTTTTTTTTGCTAATGCAATAGTTACTGAACTGTTCCGATGAGTTAATGTTGAACAAATCTCCTGTTGCGTGGTGGTCGCTTTTACCACAGATGCg                                                                                                                  >  1:490576/1‑251 (MQ=255)                                                                                                                              aaTCCACGCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTTGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATTACAGTCGCCTTTTTCTGTATCGTGGAAATCATTTTCATTGTTATTGTTAGCGAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTGGAACAAATCTCATGTTGCGTGTTGGTCGCGTTTACCTCAGATGCGttt                                                                                                               <  2:320088/251‑1 (MQ=255)                                                                                                                                       cAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTa                                                                                                      <  1:389891/251‑1 (MQ=255)                                                                                                                                              ttCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGttt                                                                                                <  1:411807/250‑1 (MQ=255)                                                                                                                                                     cACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGa                                                                                                                                                             <  1:181910/182‑1 (MQ=255)                                                                                                                                                     cACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGGTTCGCGGGGATGGGTGCGGCCCCGACAGGCTCCCTTTTTTGTATCGTGGGAATACTTTTTTTTTTTTTTTTTTGCTAATGCCATTGTTTTTGGACTGATCCGTTGAGTTAATTTTgt                                                                                                                                                             >  2:181910/1‑181 (MQ=255)                                                                                                                                                          gTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGGTTTGCGGGGATGGGTGCGACCCTGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTTTTAGCTAATGCAATAGTTACTTAACTGATCCGATGTGTTAATGTTGTACAAATCTCCTGTTTCGTGGTGGTCGCTTTTACCACAGATGCGGTTATGACCGTATGGGTTGTTTGATTTttat                                                                                   >  1:145483/1‑251 (MQ=255)                                                                                                                                                            tgaagtcTGAACCGGTCCAGTCGCTGGGGCGGTGGCAACGGCGGGTGTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTTCTTTTTCTGTATCGTGGAAATCATTTTCATTTGTATTGTTAGCTAATGCAATAGTTACTGATCTGACCCGATGAGTTAATGTTGTACAAATCTCATGTTGCGTGGTGGTCGCTTTTGCCACAGATGCGGTTATGCCAGTAGGGTTTGTTGAATTTttatt                                                                                  <  1:45134/244‑1 (MQ=255)                                                                                                                                                                            ccAGGCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCgt                                                                 <  2:407625/251‑1 (MQ=255)                                                                                                                                                                                 cGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCg                                                             <  1:191127/250‑1 (MQ=255)                                                                                                                                                                                         cGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACTGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAgt                                                    <  1:10325/251‑1 (MQ=255)                                                                                                                                                                                                         ttttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGACAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTgg                                    >  1:233278/1‑251 (MQ=255)                                                                                                                                                                                                               gCCGCGTTTCGCGGGAATGGCTGCTACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATTCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACATAGGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTTTGTTGAGCGTGTCGGGCGCAAGTGTGAGCAGCAAATTGGaac                                 <  2:31775/248‑1 (MQ=255)                                                                                                                                                                                                                           gggAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAACTCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAAATGCTCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTAAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGt                  >  2:239941/1‑251 (MQ=255)                                                                                                                                                                                                                           gggAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGCGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTTAGCAGCAAAGTGGAACAACTGCGTTCGCAGt                  <  2:29356/251‑1 (MQ=255)                                                                                                                                                                                                                             gAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTa                <  1:47812/251‑1 (MQ=255)                                                                                                                                                                                                                                   gTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGATCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGaa          <  1:260036/251‑1 (MQ=255)                                                                                                                                                                                                                                      cGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCAGTTTTATTGTTAGCTAATGCAATAGGTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGGTGAAGTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGGGAGCAGCAAAGTGGAACAACGGCGTGCGCAGTTAAATGAACGt       <  2:59090/251‑1 (MQ=255)                                                                                                                                                                                                                                         ccATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTAtt    >  1:211051/1‑251 (MQ=255)                                                                                                                                                                                                                                           aTGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTTCCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTGAATCTGTGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCt  <  2:24448/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                AGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCT  >  NC_000913/4223393‑4223877

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