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%	79.0 / NA	25	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/16);  total (9/16)

GCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTG  >  NC_000913/4223421‑4223884                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                       gctgTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAtcattttcattt                                                                                                                                                                                                                         <  2:460417/249‑1 (MQ=255)   ggTTGCAGCGGGGTTAGCAGGAGGTGGGGCGTGTAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGAGTCGGCAATCCACGGCAGTAATTTGAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTGTCGCGGGAATGGGTGCGGCCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATGTTTAt                                                                                                                                                                                                                     <  2:138265/250‑1 (MQ=255)     ttGCATCGTGGTTAGCAGGCGGTGGGTCGTGGCATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGACACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTTCCCCGGTTCGCGGGGATGGGGGGGGCCCCGAAAGTGTCCTTTTTCTGTATCGTGGAACACATTtttatttttattt                                                                                                                                                                                                                   >  1:179419/1‑250 (MQ=255)         aTCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGTTTCGCGGGGATGGGGGCGGCCATGGCCGGCTCCTTTTTCTGTATCGTGGGAATCATTTTCATTTTTTTTGtttg                                                                                                                                                                                                               >  1:262159/1‑249 (MQ=255)                                        tGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGAt                                                                                                                                                                                 <  1:345169/250‑1 (MQ=255)                                                ccGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATg                                                                                                                                                                        <  2:92502/251‑1 (MQ=255)                                                                                     cATTGGATTCGGCAATCCACTCCAGGAATTTCAGGCCACGCGTTAAAGTCTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGGTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTTGTCGCTTTTAc                                                                                                                                   <  2:270490/251‑1 (MQ=255)                                                                                       ttGGATTCTGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTAcc                                                                                                                                  <  1:330211/250‑1 (MQ=255)                                                                                                  aatcccagccAGTAATGTCAGGCCACGCGGTAAAGTCTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGTGTTTTCTGCCGCGTTTCGCTGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCTATTAGCTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCCCAGATTCGttt                                                                                                                      <  2:128509/243‑1 (MQ=255)                                                                                                     ccACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTGTCTGCCGTGTTGCGCGGGAATGGGTGCTACCATGACAGTCTCCTTTGTCGGTATCGTGGAAATCATTTTCATTTGTATTGTGAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTTATCTAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTAt                                                                                                                    <  1:230553/250‑1 (MQ=255)                                                                                                     ccACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGGCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTAt                                                                                                                    <  2:245048/250‑1 (MQ=255)                                                                                                        cTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTgg                                                                                                                                                                                                                                        <  2:339636/131‑1 (MQ=255)                                                                                                            aGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGGATGGGGGGGGGCCTCGCCGTCCCCTTTTTCTTTTTAGTGGGAATCATTTTTATTTTTTTTTTTTGCTAATGCAATAGTTTATTAACTTATCCGGTGAGTTTATGTTGAACAAAACTCATGTTGCGGGGTGGTCGCTTTTTCCACAGATGCGTTTATTCCAGTAt                                                                                                            >  1:378946/1‑251 (MQ=255)                                                                                                                     aGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGAGCGGTGGCAACGGCGGGTTTTCTGCCGCGCTTCGCGGGGATGGGGGCGAACATGACAGACTCCTTTTTCTGTATCGTGGAAATCATTTTCTTTTTTATTTTTAGCTAATTCAATAGTTACTGGACTGATCCGATGGGTTAATGTTGTACAAATCTCATGTTTCGGCGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTTTGGTttgttg                                                                                                   >  1:207910/1‑251 (MQ=255)                                                                                                                                          aCCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGtt                                                                              <  1:331182/251‑1 (MQ=255)                                                                                                                                                               gTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAAGTTTTAGTAAATCTGGGTTGAGCGTGTCGGGAGCAAgtgt                                                         <  1:68821/251‑1 (MQ=255)                                                                                                                                                                         cggGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCaaa                                               <  1:154273/251‑1 (MQ=255)                                                                                                                                                                            gTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTTGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGTGTGTCGGGAGCAAGTGTGAGCAGCAAAGTg                                            <  1:233075/251‑1 (MQ=255)                                                                                                                                                                                                   aTGGGAGCGACCATGACAGTCTCCTTTTTCTGTCTCGTGGAAATCCTTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTaa                      >  2:299156/1‑250 (MQ=255)                                                                                                                                                                                                     gggTGCGACAATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATg                   >  2:466948/1‑251 (MQ=255)                                                                                                                                                                                                        tGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTggg                                                                                <  2:109396/187‑1 (MQ=255)                                                                                                                                                                                                        tGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTggg                                                                                >  1:109396/1‑187 (MQ=255)                                                                                                                                                                                                            aCCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACGGGTGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTTGGTGGAGCGTGTTGTGAGTAAGTGTGAGCATCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTAt            <  1:185486/251‑1 (MQ=255)                                                                                                                                                                                                                    aGTCTCCTTTTTCTGTATCGTAGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATACGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTATGGTGc    >  2:390265/1‑251 (MQ=255)                                                                                                                                                                                                                      tctcCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATACGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTACATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGAAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTg  >  2:91180/1‑251 (MQ=255)                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                       GCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCTG  >  NC_000913/4223421‑4223884

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