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%	86.7 / NA	29	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/20);  total (9/20)

AACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACA‑GTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATT  >  NC_000913/4223401‑4223875                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                              aaCCTTACGGACGAAATCCTGCTGTTGCATCGGTGTTAGTAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTCCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCTCTCCAGTAATTTCAGGCCACGCGTGAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCTGGTTTTCTGCCGCGTTTCGCGTGAATGGGTGCGACCATGACT‑GTCTCCTTGTTCTGTATCg                                                                                                                                                                                                                                  <  2:295494/251‑1 (MQ=255) aaCCTGACGCACGAAACCCTGCTGTTGCATCGTTGTTAGCAGGCGGTGTGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCTTTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCTGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACA‑GTCTCCTTTTTCTGTATCg                                                                                                                                                                                                                                  <  1:100483/251‑1 (MQ=255)               aaCCCTGCTGTTGCATCGGGGTTAGCAGGGGGTGGGTCGTGGAATTGGGTAACCCGGCGTGTTGCGCCAGTTCCGTGAGTGCCACACGGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACA‑GTCTCCTTTTTCTGTATCGTGGAAAtcattttc                                                                                                                                                                                                                    <  2:323779/251‑1 (MQ=255)                        gTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACA‑GTCTCCTTTTTCTGTATCGTGGAAATcattttcattttattg                                                                                                                                                                                                           <  1:341314/251‑5 (MQ=255)                        gTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACA‑GTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTa                                                                                                                                                                                                             <  1:491993/249‑1 (MQ=255)                        gTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACA‑GTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTAt                                                                                                                                                                                                            <  1:470745/250‑1 (MQ=255)                        gTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTGCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACA‑GTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTAtt                                                                                                                                                                                                           <  2:131043/251‑1 (MQ=255)                             aTCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGAGGCAACGGCGGGGTTTCTGACGCGGTTCGCGGGGAAGGGTGCGACCAGGAAA‑AGCTCCTTTTTTTGTAACGTGGAAAACATTTTCATTTTTAttttt                                                                                                                                                                                                        >  1:461464/1‑246 (MQ=255)                               cGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGATTGCCACACTGCCATTGGAGTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGTCGGGTTTTCTGCCGCGTTTCGCTGGAATGGGTGCGACCATGACA‑GTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGc                                                                                                                                                                                                     <  1:288012/250‑1 (MQ=255)                                  ggTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTTGCGGGGGTGGGTGGGGCCATGACA‑GTCTCCCTTTTTTGTGTTGTGGGAATTACTTTTATTTTTTTTGtttggtaa                                                                                                                                                                                                  >  1:282582/1‑244 (MQ=255)                                              gTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGTTTTGCGGGAATGGGGGCGAACATGACA‑GTCTCCTTTTTTTTTATCGTGGGAATCATTTTTATTTTTTTTGTTATCTTATGGAATAGTTaat                                                                                                                                                                                     >  1:222075/1‑249 (MQ=255)                                                                     cGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACA‑GTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGt                                                                                                                                                              <  1:163330/251‑1 (MQ=255)                                                                      ggCTTGTAGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACAGCGTGTTTTCTGCCGCGTTTCGCGGGAATTGGTGCGACCATGACA‑GTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGtt                                                                                                                                                             <  1:285783/251‑1 (MQ=255)                                                                                gcCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGTTTCGCGGGAATGGGGGCGACCCTGACA‑GTCTCCTTTTTCTGTATCGTGGGAATCCTTTTTATTTTTTTTGTTATCTAATGCAATAGTTACTTAAATGATCCggtaagtttaagttttacaaatct                                                                                                                                                   >  1:345936/1‑228 (MQ=255)                                                                                         gTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGCTTCGCGGGGATTGGTGCGAACATG‑CACGACTCCTTTTTCTTTATCGTTGGAATCATTTTCCTTTTTTTTTTTAGCTAATGCAATAGTTACTGTACTGATCCGATGAGTTAATGTTGAAAAAATCTCATGTTGCg                                                                                                                                          >  1:356504/1‑251 (MQ=255)                                                                                                                    gCAATCCACTCCAGTAATTGCAGGCCACGCGTTAAAGACTGAGCCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACA‑GTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGGCGCTTTTCCCCCAGATGCgc                                                                                                               <  2:292295/251‑2 (MQ=255)                                                                                                                              ccAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACA‑GTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGt                                                                                                     <  2:5166/251‑1 (MQ=255)                                                                                                                                                     aaaGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCGCCGTTCGCGGGGAAGGGGGCGGCCATCACA‑GGCTCCTTTTTTTGTATCGTGGGAAACCTTTTCCTTTTTTTTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGGGTTTATGTTGAACAAATCTCATGTTTCGGGGGGGGCGCTTTTTACACAGGTGCGTTTTTGGCCGTATTGTTTGTTTGATTTttttta                                                                               >  2:425303/1‑250 (MQ=255)                                                                                                                                                         aCTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACA‑GTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAAGTTTTATTAAATCt                                                                          <  2:264499/251‑1 (MQ=255)                                                                                                                                                                                  ggTGGCAACGGCGGGTTTTTTGCCGCGGTTCGCGGGAATGGGGGCGGCCCTGACA‑GGCTCCTTTTTCTGTATCGTGGGAATTATTTTTATTTTTTTTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTTCCACAGATGCGTTTATGCCAGGATGGTTTGTTGGATTTTTTTTAAATCTGGGTTgggggatcggaagagcaca                                                  >  1:417839/1‑232 (MQ=255)                                                                                                                                                                                  ggTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACA‑GTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAg                                                                  <  2:417839/234‑1 (MQ=255)                                                                                                                                                                                                gTTTTCTGCCGCGTTTCGCGGTAATGGGTGCGACCATGACA‑GTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGAGCCGATGAGTTAATGGTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATTCGTTTATGCCAGTATGGGTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGGAGCAAAGTg                                   <  2:79740/251‑1 (MQ=255)                                                                                                                                                                                                     cTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACA‑GTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGaac                               >  1:18585/1‑250 (MQ=255)                                                                                                                                                                                                       gtcgcgTTGCGCGGGAATCGGTGCGTCTATGACA‑GTTTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGGTAGCTAACGCAATATTTACTGAACTGATCCTATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAAACTGGGTTGAGCGTGTCGGGAGCGAGTGTGAGCAGCAAAGAGGaac                               <  2:8886/246‑1 (MQ=255)                                                                                                                                                                                                       gCCGCGTTTCGCGGTAATGGGTGCGACCATGACA‑GTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCTCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTTAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGaaca                              <  1:215399/249‑1 (MQ=255)                                                                                                                                                                                                                             gCGACCATGACA‑GTCTCCTTTTTCTGTAACGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTTGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAAc       <  2:172916/250‑1 (MQ=255)                                                                                                                                                                                                                                 ccATGACA‑GTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTAtt  <  2:213887/251‑1 (MQ=255)                                                                                                                                                                                                                                           tctcCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACtgcgtgc                     <  2:409364/223‑1 (MQ=255)                                                                                                                                                                                                                                           tctcCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACtgcgtgc                     >  1:409364/1‑223 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                              AACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACA‑GTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATT  >  NC_000913/4223401‑4223875

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