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%	87.0 / NA	27	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 (9/18);  total (9/18)

GCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGGAATGGGTGCGACCAT‑GACAGTCTCTTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGAT‑CCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTT  >  NC_000913/4223398‑4223862                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                   gCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGTTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGACTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGGAATGGGTGCGACCAT‑GACAGTCTCCTTTTTCTGTa                                                                                                                                                                                                                          <  1:150418/251‑1 (MQ=255)          aCGCACGAAACCCTTCTGTTGAATCGTGGTTAGCAGGCGGGGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCTCACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGGAATGGGTGCGACCAT‑GACAGTCTCCTTTTTCTGTAGCGTGG‑aaa                                                                                                                                                                                                                <  2:222281/251‑1 (MQ=255)              aCGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGGAATGGGTGCGACCAT‑GACAGTCTCCTTTTTCTGTATCGTGG‑GAAtcat                                                                                                                                                                                                            <  1:71236/251‑1 (MQ=255)                 aagtcctgaggttgactCGTGGGTAGGAGGTGGTGAGTCGTGGAATTGGTGATCCCGGCGTGTTGCGCTAGTTCCGTTAGTGACACCCTGCCATTGGATTCAGCAATCCACTCCAGGAATATCGGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGGGGCATCGTCGGGTTGTCTGCCGC‑GTTTCGCGGGAATGGGTGCGGCCAT‑GACAGTCTCCTGTTTCTTTATCGTGG‑AAAgcatttt                                                                                                                                                                                                         <  2:195814/235‑1 (MQ=255)                 aaaCCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGTCGC‑GTTTCGCGGGGATGGGGGCGGGCAT‑GACAGGCTCCTTTTTCttttttgtggggataatttt                                                                                                                                                                                                          >  2:102316/1‑232 (MQ=255)                     ccTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAAGCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGGAATGGGTGCGACCAT‑GACAGTCTCCTTTTTCTGTATCGTGG‑AAAtcattttcatt                                                                                                                                                                                                     <  1:225533/251‑1 (MQ=255)                                agcGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGGTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCCACGGCGGGTTTTCTGCCGC‑GTTTCGCGGGAATGGGTGCGACCAT‑GACAGTCTCCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAg                                                                                                                                                                                          <  1:93913/249‑1 (MQ=255)                                  cGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGGAATGGGTGCGACCAT‑GACAGTCTCCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCt                                                                                                                                                                                        <  1:184732/251‑1 (MQ=255)                                         aGCAGGGGGTGGGGCGTGGAATTGAGTAACCCGGTGTGTTGCGCCAGTTCCGTGAGGGCCACACTGCTATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCGGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGGAATGGGTGCGACCAT‑GACAGTCTCCTTTTTCTGTATCGTGG‑AAAGCAGTTTCATTTCTATTGATAGCTAATGCaa                                                                                                                                                                                 <  2:117581/251‑1 (MQ=255)                                         aGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAAGTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGGAATGGGTGCGACCAT‑GACAGTCGCCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTAGTGTTAGCTAATGCaa                                                                                                                                                                                 <  1:234620/251‑1 (MQ=255)                                            aGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCCC‑GGTTTGCGGGGAGGGGGGGGACCAA‑GACAGACTCCTTTTTTTGTATCGTGG‑AAAACATTTTTTTTTTTTTTGTTTGCTCATGAAAttg                                                                                                                                                                              >  1:83780/1‑249 (MQ=255)                                                   gggTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGCGTGCCACACTGCCATTGGCTTCGGCAACCCACGCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGCCCAGTCGCCGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGGAATGGGTGCGACCAT‑GACAGTCTCCTTTTTCTGTATCGTTG‑AAACCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTactga                                                                                                                                                                       <  1:80681/251‑1 (MQ=255)                                                                                                gTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCGGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTGCGCGGGAATGGGTGCGACCAT‑GACAGTCTCCTTTTTCTGTATCGTGG‑a                                                                                                                                                                                                                  <  1:13377/163‑1 (MQ=255)                                                                                                  gCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGGAATGGGTGCGACCAT‑GACAGTCTCCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGAT‑CCGATGAGTTAATGTTGAACAAATCTCATGTTGCgtggtg                                                                                                                        <  1:102316/250‑1 (MQ=255)                                                                                                      cacTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCTGTTGC‑CGGGGAGGGGGGGGACCATAGA‑AATCTCCTTTTTCTTTATCGTGG‑GAATCATTTTCATTTTTTTTGTTTTCTAATTCCATAGTTTCTGAACTAAT‑ACGATGAGTTAATGTTGTACAAATATCATGtttcttggagggcgc                                                                                                                   >  1:182405/1‑238 (MQ=255)                                                                                                       acTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGGAATGGGTGCGACCAT‑GACAGTCTCCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGAT‑CCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCt                                                                                                                  <  1:245863/251‑1 (MQ=255)                                                                                                                    tCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGGAATGGGTGCGACCAT‑GACAGTCTCCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGAT‑CCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGAt                                                                                                      <  1:167995/250‑1 (MQ=255)                                                                                                                             cACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGGAATGGGTGCGACCAT‑GACAGTCTCCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGAT‑CCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGc                                                                                            <  2:75291/251‑1 (MQ=255)                                                                                                                                  cAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCCC‑GTTTCGCGGGGGTGGGTGCGACCCT‑GGCAGTCTCCTTTTTTTGTTTCGGGG‑GAAACATTTTTATTTTTATTTTTTGCTAATTCAATAGTTATTGAACTGAT‑CCCATGAGTTAATGTTGTACAAAACTCATGTTGCGTGGGGGTCGCCTTTTCCCCAGATTCGTTTATGCCAGt                                                                                        >  2:17284/1‑250 (MQ=255)                                                                                                                                               ccACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGATTTTGCCGC‑GTTTCGCGGGAATGGGAGCGACCCA‑GAAAGTCTCCTTTTTCTGTATTGTGGGAAAAC‑TTTTCATTTTTATTTTTAGCTAATGCAAAAGTTACTGAAATGATTCCG‑TGAGTTTATGTTGAAAAAATCTCATTTTGCCGGGGGGTCGCTTTTACCACCAATGCGTTTATGCCAGTATGGGTTGGTGAAt                                                                          >  1:105321/1‑251 (MQ=255)                                                                                                                                                cACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGATTTTTGCCGC‑GGTTCGCGGGGATGGGGGGGCCCAT‑GACAGACTCCTTTTTTTGTTTCGTGG‑GAATCATTTTTATTTTTTTTTTTAGCTAATGCAATAGTTAATGAACTGAT‑ACGATGTGTTAATGTTGAAAAAATCTCATGTTGCGTGGGGGGCGCCTTTACCACCGATGCGTTTTTGCCCGTATGGGTTGTTGAAtt                                                                         >  1:231314/1‑251 (MQ=255)                                                                                                                                                                  ccTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGGAATGGGTGCGACCAT‑GACAGTCTCCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGAT‑CCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTg                                                       <  2:6592/251‑1 (MQ=255)                                                                                                                                                                              cTGGGGCGGTGGCAACGGCGGGTTTGCTGCCGC‑GTTTCGCGGGAATGGGTGCGACCAT‑GACAGTCTCCTGTTTCTGTATCGTGT‑AAAGCATGTTCATGTTTATGGATAGCTAATGCAATAGTTACGGAACTGAT‑CCGATGAGTTAATGGTGTGCAAAGCGCAGGGTGCGCGGTGGTCGCTTTTGCGACAGACGCGTTTATGCCCGTAAGTTGTGTTAAATTGTTATTAAATCTGGGTAGTGCGTGTCgg                                             <  2:269703/249‑1 (MQ=255)                                                                                                                                                                               tGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GGTTCGCGGGGATGGGGGCGGCCAT‑GAAAGTCTCCTTTTTCTGTATCGTGG‑GAATCATTTTCCTTTTTTTTGTTTGCTAATGCAATAATTAATGAAATGAT‑CCGTTGAGTTAATGTTGAACAAATCTCTTGTTTCGTGGTGGTCGCTTTTATCACCGATGCGGTTATGCCCGGATGGTTTGTTGAATTTTTTTTTAATTTGGGGTGAGCGGGTCggggg                                          >  1:267481/1‑249 (MQ=255)                                                                                                                                                                                                                   tCGCGGGAATGGGTGCGACCAT‑GACAGTCTCCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGAT‑CCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACagcgtgc        <  1:73079/250‑1 (MQ=255)                                                                                                                                                                                                                       gggAATGGGTGCGACCAT‑GACAGTCTCCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGAT‑CCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGt   >  1:39839/1‑251 (MQ=255)                                                                                                                                                                                                                         gAATGGGTGCGACCAT‑GACAGTCTCCTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGAT‑CCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGtt  >  1:101569/1‑250 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                   GCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGC‑GTTTCGCGGGAATGGGTGCGACCAT‑GACAGTCTCTTTTTTCTGTATCGTGG‑AAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGAT‑CCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTT  >  NC_000913/4223398‑4223862

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