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.2 / 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 (8/19);  total (8/19)

ACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTT‑GAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGT  >  NC_000913/4223402‑4223872                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                            aCCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCAGCGGTTTGCGGGGATGGGTGCGGCCCTGAAAGTCTCCCTTTTTTGTATCGt                                                                                                                                                                                                                               >  1:82266/1‑251 (MQ=255)          aCGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTTCC‑CGGTTCGCGGGAATGGGGGGGGCCCTGGCAGTCTCCTTTTTCTGTTTCGTggggaacat                                                                                                                                                                                                                      >  1:35919/1‑244 (MQ=255)                                  gTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATg                                                                                                                                                                                              <  2:97854/251‑1 (MQ=255)                                                   cGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTactgaactg                                                                                                                                                                             <  1:13012/251‑1 (MQ=255)                                                    gTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAtcatttt                                                                                                                                                                                                                   <  2:112278/212‑1 (MQ=255)                                                      ggAATTGGGTTACCCGGCTTGGTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACGCCAGTGATGTCAGGCCACGCGGTAAAGTCTGGACCTGTGCAGTCGCTGGTGCGGTGGCAACGGCGGG‑GTTTTTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTAGCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGAt                                                                                                                                                                           <  2:132605/250‑1 (MQ=255)                                                       gAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATcc                                                                                                                                                                         <  1:133619/251‑1 (MQ=255)                                                                aaCCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTTAGCCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTTCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGCAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTAAACGGATCCGATGAGtt                                                                                                                                                                 <  1:226792/250‑1 (MQ=255)                                                                                                cacaCTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑GTTTTTGCCGCGGTTTGCGGGAAAGGGGGGGGCCATGACAGTCTCCTTTTTTTTTTTTGTGGAAATCATTTTTATTTTTTTTTTTTGCTAATGCAATAGTTACTGAACTGTTCCGATGAGTTAATTTT‑GAACAAATTTCATTTTGCGtggtgggg                                                                                                                               >  1:147511/1‑249 (MQ=255)                                                                                                                tCGGCAATCCACTCCAGGAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTT‑GAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGAt                                                                                                                <  1:227323/250‑1 (MQ=255)                                                                                                                 cGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑GTTTCTGCCGCGTTTTGTGGGAATGGGGGCGGCCCCGAAAGGCTCCTTTTTTTGTATCGGGGAAATCATTTTTATTTTTTTTTTTTGCTAATGCAATAATTAATGAACTTATCCGCAGAGTTAATGTTTGAA‑AAATATCATGTTGCGTGGGGGTCGGTTTTatcccagaatg                                                                                                              >  1:369610/1‑241 (MQ=255)                                                                                                                                 tAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑GTTTCTGCCGCGTTTTGCGGGAAAGGGTGCGACCATGACAGTCTCCTTTTTTTGTATCGTGGGAATCATTTTTATTTTTTTTTTTAGCTAATGCAATAGTTACTGAACTGATCCGGTGAGTTAATGTT‑TAACAAATCTCATGTTGCGTGGTGGGCGCTTTTACCACAGATGCGTTTTTGCCCGTATgg                                                                                              >  1:108982/1‑251 (MQ=255)                                                                                                                                   aTTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGGCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGGTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTT‑GAACAAATCTCATGTTGCGTGGTGGTCTCTTTTACCACAGATGCGTTTATGCCAGTAGGGtt                                                                                            <  1:355031/251‑1 (MQ=255)                                                                                                                                             ccgcgTTAAAGACTGAACCTGTCCAGGCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGTCCATGACAGTCTCCTTTTTCTGTAGCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTT‑GAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAAttt                                                                                  <  2:197630/250‑1 (MQ=255)                                                                                                                                              cgcgTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑GTTTTTGCCGCGGTTCGCGGGAATGGGTGCGACCATGGCAGTCTCCTTTTTTTGTTTCGTGGGAATAATTTTCATTTTTTTTTTTTGATAATGCAATAGGTACTGAAATGATCCCGTGAGGTTATGTT‑TAACAAAACTTATGTTTCGTGGGGGGCGCTTTTTCCACAGATGCGGTTTTGACAGTATGGTTTTTTGAAtttt                                                                                 >  1:303606/1‑251 (MQ=255)                                                                                                                                                    aaaGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTATTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTT‑GAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTttatta                                                                            <  2:82266/250‑1 (MQ=255)                                                                                                                                                        aCTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTT‑GAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCt                                                                       <  1:240978/251‑1 (MQ=255)                                                                                                                                                                 gTCCTGTCGCTGGGGCGTTGGGAACGGCGCG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTTCGACCATGACAGTTTCCGTGTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTT‑GAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGTTGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGa                                                                <  2:132519/249‑1 (MQ=255)                                                                                                                                                                      ggcGCTGGTGCGTGGGCAACGTCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTT‑GAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTc                                                         <  2:181851/249‑1 (MQ=255)                                                                                                                                                                                               g‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTT‑GAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTg                                >  2:346635/1‑251 (MQ=255)                                                                                                                                                                                                                cgcgGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTT‑GAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGGGAGCAGCAAAGTGGAACAACTGCGTgcgc                >  2:61395/1‑251 (MQ=255)                                                                                                                                                                                                                  cgGGAATGGGTGCGACCATGACAGTCTCCTTATTATGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAAGAGTTACTGAACTGAACCGATGAGTTAATGAT‑GAACAAAGCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATTCGTGTAGGCCAGTATGGTTTGTTGAATGTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGGGTGAGCAGCAAAGTGGAACAACTGCGTACGCAg              <  2:175435/251‑1 (MQ=255)                                                                                                                                                                                                                     gAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTT‑GAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTa           <  2:147511/251‑1 (MQ=255)                                                                                                                                                                                                                     gAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTT‑GAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTa           <  2:358583/251‑1 (MQ=255)                                                                                                                                                                                                                           gTGCGACCATGACAGTCTTCTTTTTCTGTATCGTGGAAATCATTTTCAGTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTT‑GAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGaa     <  1:249913/251‑1 (MQ=255)                                                                                                                                                                                                                             gCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTGTATTGTTAGCTAATGCCATAGTTACTGAACTGATACGATGAGTTAATGTT‑GAACAAATCTCATGTTGCTTGGTGGTCGCTTTTGCCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGTAGCAAGTGTTAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACg   <  2:324351/251‑1 (MQ=255)                                                                                                                                                                                                                              cGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTT‑GAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGt  <  1:342851/251‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                            ACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGG‑TTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTT‑GAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGT  >  NC_000913/4223402‑4223872

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