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%	52.8 / NA	20	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 (2/18);  total (2/18)

GTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTAT‑CGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGA  >  NC_000913/4223394‑4223868                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                        gTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTc                                                                                                                                                                                                                                   <  1:69752/251‑1 (MQ=255)              aCGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGTTAACCCGGCTTGTTTCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAAGTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTAT‑CGTGGaaa                                                                                                                                                                                                                     <  1:129746/251‑1 (MQ=255)                                gtGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATTACAGTCTCCTTTTTCTGTAT‑CGTGGAAATCATTTTCATTTTTATTg                                                                                                                                                                                                   <  2:15596/250‑1 (MQ=255)                                    aTCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTTTCCTTTTTCTGTAT‑CGTGGAAATCATGTTCATTTTTATTGTTAg                                                                                                                                                                                               <  1:46618/251‑1 (MQ=255)                                       gTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTAT‑CGTGGAAATCATTTTCATTTTTATTGTTAGCTa                                                                                                                                                                                            <  1:165262/251‑1 (MQ=255)                                                        ggTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTAT‑CGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTactgaa                                                                                                                                                                           <  2:1288/251‑1 (MQ=255)                                                                                 tgttgCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGGATGGGGGCGACCAAGACAGTCTCCTTTTTCTTTATTCTGGGAAAT‑ATTTTTATTTTTTTTGTTTGATCATGCAATAGTTACTGTACTGATACGATggggttatgtttaa                                                                                                                                                   >  1:137887/1‑237 (MQ=255)                                                                                          aGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTAT‑CGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCAt                                                                                                                                         <  1:4677/251‑1 (MQ=255)                                                                                                      gCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCTATAAAGACTGAGCCTAGCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTCTCTGCCGCGTTTGGCTGGAATGTGTGCAACCATGACAGTCGCCTTTTTGTGTAT‑CGTGGAAACCATTTTCATGTTAATTCTTAGCTACTGCAATAGGTACTGCACTTACCCGATGAGTCAATGGTGAGTAAAGCTCATGTTTCgtggtg                                                                                                                              <  2:251650/250‑1 (MQ=255)                                                                                                                cATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTAT‑CGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTAc                                                                                                                   <  1:245885/251‑1 (MQ=255)                                                                                                                    ggATTCGGCAATCCACGCCAGTAACTTCAGGCCACGCGTTAAAGACTGAACCTGTCCATTCGTTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTGCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTAT‑CGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACcac                                                                                                                <  1:13752/250‑1 (MQ=255)                                                                                                                           gCAATCCACTCCAGTAATCTCAGGCCACGCGTTAAAGACTGAACCTGTTCAGTCGCTGGTGCGGTGGCAACGGCTGTTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGGCCATGACATGCTCCTTTTTCTGTAT‑CGTGGAAATCATTTTCATTTTTATTGTTAGCTAAGGCAATAGTTACTGAACGGATCCGATGAGTTAATGTTGTACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGt                                                                                                        <  1:127813/251‑1 (MQ=255)                                                                                                                                  aCTCCAGTAATTTCAGGCCGCGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTAT‑CGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAAGCTCATGTTGCGTGGTTGTCGCTTTTACCACAGATGCGTTTATGcc                                                                                                 <  2:181991/251‑1 (MQ=255)                                                                                                                                                   ccACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTAT‑CGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAAt                                                                                <  1:42042/251‑1 (MQ=255)                                                                                                                                                                    aaCCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTAT‑CGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGt                                                               <  2:250753/251‑1 (MQ=255)                                                                                                                                                                            cAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTAT‑CGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTGTATTAAATCTGGGTTGAGCgtg                                                       <  1:97573/251‑1 (MQ=255)                                                                                                                                                                                   tGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCTCGGGAATGGGTGCGACCATGACAGTCTCCCTTTTTTTTAT‑CGTGGCCATAATTTTCATTTTTTTTTTTTGCTTATGCAATAGTTTATGTAATTTTCCGCTGAGGTAATGTTGCACAAATCTCATCCTGAGTGGTGGTCGCTTTTTCCACAAGTGCGTTTTTGCCCGTATGGTTTGTTCAATTTttattatattttggtttagcgagttagggg                                                >  1:2622/1‑227 (MQ=255)                                                                                                                                                                                                                 gcgTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTAT‑CGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACtgc                  <  1:170610/251‑1 (MQ=255)                                                                                                                                                                                                                            gAATGGGTGCGACCATGACAGTCTCCTGTTTCTGTAT‑CGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTGGAATTTTTAGTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACGGCGTGCGCAGTTa       <  2:108089/251‑1 (MQ=255)                                                                                                                                                                                                                                 ggTGCGACCATGACAGTCTCCTTTTTCTGTAT‑CGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTGTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACGGCGTGCGCAGTGAAATGa  <  2:231008/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                        GTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTAT‑CGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGA  >  NC_000913/4223394‑4223868

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