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%	101.3 / NA	31	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/22);  total (9/22)

GAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCG‑ACCA‑TGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCAT‑GTTGCGTGGTGGTCGCTTTTACCACA‑GATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGT  >  NC_000913/4223413‑4223880                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     gAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGTGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCG‑ACCA‑TGACAGTCTCCTTTTTCTGTATCGTGGAAAtcattt                                                                                                                                                                                                                             <  1:124623/251‑1 (MQ=255)   aaCCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCG‑ACCA‑TGACAGTCTCCTTTTTCTGTATCGTGGAAAtcattttc                                                                                                                                                                                                                           <  2:63824/251‑1 (MQ=255)     cccTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAAGTGGGTACCCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCG‑ACCA‑TGACAGTCTCCTTTTTCTGTATCGTGGAAAtcattttcat                                                                                                                                                                                                                         <  1:204508/251‑1 (MQ=255)           tgTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGTTTTGCGGGGATGGGGGCG‑GCCA‑TGGCAGGCCCCTTTTTTTTTATCGTGGGAATCatttttattttttt                                                                                                                                                                                                                   >  1:186983/1‑249 (MQ=255)               gCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCGTGTTGCGCCAGTTCCGTGAGTGCATCACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTTGCTGGTGCGGTGGGAACGGCGGGTTGTCTGCGGCGTTGCGCGGGAATGGGTGCT‑ACAA‑TTACAGTCTCCTTTTTCTGTACCGTGGAAATGATTTTCATTTTTATTGt                                                                                                                                                                                                                <  2:14445/250‑1 (MQ=255)                 aTCGTGGTTAGCCGGCGGTGGGTCGTGGAAGTGGGTAACCCGGCGTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCG‑ACCA‑TGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAg                                                                                                                                                                                                             <  2:2111/251‑1 (MQ=255)                   cGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGGCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCG‑ACCA‑TGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCt                                                                                                                                                                                                           <  2:139491/251‑1 (MQ=255)                          aGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGCTTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTTCCGCGATTCGCGGGGATAGGTGGG‑GCCA‑TGACACTCTCCTTTTTCTGTTTTGGGGGAATCATTTTTCTTTTTTTTTTTAGcatattct                                                                                                                                                                                                     >  1:17495/1‑243 (MQ=255)                           gCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAAGTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCG‑ACCA‑TGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAAt                                                                                                                                                                                                   <  1:219235/251‑1 (MQ=255)                                  gTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGGTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCACCGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCG‑ACCA‑TGACAGTCTCCTTTTTCTGTAGCGTGGAAATCATGTTCATTTTTATTGTTAGCTAATTCAATAGTTact                                                                                                                                                                                            <  2:186816/251‑1 (MQ=255)                                     ggTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCG‑ACCA‑TGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTactgaa                                                                                                                                                                                         <  1:9324/251‑1 (MQ=255)                                      gTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCG‑ACCA‑TGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTAATTTTTATTGTTAGCTAATGCAATAGTTactgaac                                                                                                                                                                                        <  2:186983/251‑1 (MQ=255)                                        cGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCGGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCG‑ACCA‑TGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTactgaactg                                                                                                                                                                                      <  2:126513/251‑1 (MQ=255)                                          gggAATTGGGTAACCCGTCTTTTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCG‑ACCA‑TGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTactgaactgaa                                                                                                                                                                                    <  2:190663/250‑2 (MQ=255)                                           ggAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGGGTGGGGGGGGCCCT‑TGA‑AGTCTCCTTTTTTTGTTTAGTGGAAAATATTTTTATTTTTTTTGTTAAGTAATAGAAAAGTTAATGAACTGAtt                                                                                                                                                                                   >  1:131634/1‑250 (MQ=255)                                             aaTTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGTTTCGCGGGAATGGGTGCG‑GCCA‑TGACAGTCTCCTTTTTCTGTATCGTGGGAATCATTTTCATTTTTTTTTTTTGCTAATGCAATAGTTAATGAACTGATCCg                                                                                                                                                                                 >  1:144062/1‑251 (MQ=255)                                                         cGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGGTTCGCGGGGATGGGGGCG‑ACCA‑TGACAGTCTCCTTTTTCTGTATCGTGGGAATCATTTTTATTTTTTTTTTTTGGTAATTGAATAGGTACTGAACTGTTTCGCTGAGGTAATGt                                                                                                                                                                     >  1:198714/1‑251 (MQ=255)                                                                  gcgcCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGTAATGGTTGCG‑ACCA‑TGACATTCTCCTTTTTCTGTATCGTGGAAAGCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACGGAACTGATCCGATGAGTTAATGTTGAACAAAt                                                                                                                                                            <  1:155948/251‑1 (MQ=255)                                                                           ccGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCG‑ACCA‑TGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCAT‑Gtt                                                                                                                                                   <  1:99284/250‑1 (MQ=255)                                                                                 gTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTGAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTTCCCCGTTTCGCGGGGGTGGGTGCG‑ACCA‑TGACAGTCTCCTTTTTCTGTTTCTTGGAGATACTTTTCATTTTTTTTTTTTTCTTATGCAATAGTTTCTTCACTTATTCGATGTGTTATTTTTTAACAAATCTCAT‑Gttttgatcg                                                                                                                                            >  2:52601/1‑244 (MQ=255)                                                                                 gTGCCACACTGCCATTGGATGCGGCAAGCCACTCCAGTAATTTCAGGCCACGCGTTGAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCG‑ACCA‑TGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCAT‑GTTg                                                                                                                                                  <  1:52601/245‑1 (MQ=255)                                                                                                 ggATTCGGTAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCG‑ACCA‑TGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCAT‑GTTGCGTGGTGGTCTCTTTTACcaca                                                                                                                            <  2:26570/251‑1 (MQ=255)                                                                                                                    aGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCG‑ACCA‑TGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCAT‑GTTGCGTGGTGGTCGCTTTTACCACA‑GATGCGTTTATGCCAGTAt                                                                                                        <  1:73955/251‑1 (MQ=255)                                                                                                                            cAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCG‑ACCA‑TGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCAT‑GTTTCGTGGTGGTCGCTTTTACCACA‑GATGCGTTTATGCCAGTATGGTttgtt                                                                                                <  1:177246/251‑1 (MQ=255)                                                                                                                                ccACGCGTTAACGACTGAACCTGTACAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCAGCGGTACGAGGGAAGGGGTGCG‑‑CCACTCAAAGGCTCCTTTTTTTCTATCGCGGGAATTATTTTTATTTTTTTTGTTTGATAATGCAATCGTTATTGAAAAGATCCGATGCGTTAATATTGAAAAAAT‑TCAAAGTTGCGTTGGGGGCGCTTTTTCCACAAGA‑GCGTTTATTTCAATATGGttgtttgtgt                                                                                            >  2:226326/1‑243 (MQ=255)                                                                                                                                   cgcgTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCTGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCG‑ACCA‑TGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCAT‑GTTGCGTGGTGGTCGCTTTTACCACA‑GATGCGTTTATGCCAGTATGGTTTGTTGAAtttt                                                                                         <  1:132831/251‑1 (MQ=255)                                                                                                                                                                                      tgtCTGCCGCGTTTCGCGGGAATGGGTGCG‑ACCA‑GGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATTTTGAACAAATCTCAT‑GTTGCGTGGTGGTCGCTTTTACCACA‑GATGCGTTTATTCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTgg                                       <  1:254977/248‑1 (MQ=255)                                                                                                                                                                                              gcgTTTCGCGGGAATGGGTGCG‑ACCA‑TGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCAT‑GTTGCGTGGTGGTCGCTTTTACCACA‑GATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACtgc                              <  2:144062/251‑1 (MQ=255)                                                                                                                                                                                                          aattggTGCG‑ATCA‑TGATAGTCTCCTTGTTCTGTATCGTGGAAATCAGTTTCAGTTTTATTGTTCGCTAATGCAATAGTTACTGATCTGAGCCGATGAGTTAATGTTGAACAAATCTCAT‑GTTGCGTGGTGGTCGCTTTTACCACA‑GATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTAGTAAATCTTGGTTGTGCGTGTCTGTAGCAAGAGTGATCAGCAAAGTGGAACAACTGCGTGAGCatttta                  <  2:155518/247‑6 (MQ=255)                                                                                                                                                                                                               gagCG‑ATCA‑TGACAGTCTCCTTTGTCTGAATCGTGGAAGTCATTTTCCTTTTTATTGTTATCTACTGCAATAGTTACTGATCTGATCCGATGAGTTAATGTTGAACGAGTCTAAT‑GCTGCGTGGTGGTCGCTTTTACCACA‑GATGCGTTTATGCCAGTATGGTTTGTTGACTTTTTATTAAATCTGGGTTGAGCGTGTCGGGTGCACGTGTGAGCAGCAAAGGGGAACAACTGCGTGCGCCGTTAAATGa              >  2:96802/3‑250 (MQ=255)                                                                                                                                                                                                                            aCAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCAT‑GTTGCGTGGTGGTCGCTTTTACCACA‑GATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTTCGTGCGCAGTTAAATGAACGTATTCTGGt  >  2:139114/1‑251 (MQ=255)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     GAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCG‑ACCA‑TGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCAT‑GTTGCGTGGTGGTCGCTTTTACCACA‑GATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGT  >  NC_000913/4223413‑4223880

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: ‑