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%	108.0 / NA	34	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 (12/22);  total (12/22)

AGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTAT  >  NC_000913/4223393‑4223874                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              aGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCcttttt                                                                                                                                                                                                                                          <  1:288473/251‑1 (MQ=255)        cAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATc                                                                                                                                                                                                                                   <  1:68324/251‑1 (MQ=255)              gACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCTACCATGACAGTCTCCTTTTTCTGTATCGTGGa                                                                                                                                                                                                                              <  1:303341/250‑1 (MQ=255)                 gCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGA‑AAtc                                                                                                                                                                                                                         <  1:62986/251‑1 (MQ=255)                       accccTGCTTGTGCGTCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCTCCAGTGCTGTGATTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAATCTGTCGAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGA‑AAtcattttc                                                                                                                                                                                                                   <  2:280699/249‑1 (MQ=255)                          ccTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGA‑AAtcattttcatt                                                                                                                                                                                                                <  1:191156/251‑1 (MQ=255)                             gctgTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGGTTCGCGGGAATGGGGGCGCCCATGACAGTCTCCTTTTTTTTTATCGTGGA‑AAACATTTTCAttttt                                                                                                                                                                                                             >  1:286714/1‑251 (MQ=255)                                   gCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGA‑AATCATTTTCATTTTTATTGt                                                                                                                                                                                                        <  1:19541/250‑1 (MQ=255)                                                        gggTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTactg                                                                                                                                                                                   <  1:410940/250‑1 (MQ=255)                                                             gTGGACTTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTACTTTCCGGCCACGCGTTAAAGACTGCACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGTTTTGCGGGGATGGGTGCGAACATGACAATCTCCTTTTTTTGTGTCGTGGA‑GATCAATTTCATTTTTATTGTTTGATACTGCCATAGTTacttaactg                                                                                                                                                                              >  2:221876/1‑250 (MQ=255)                                                                         aaCCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGtt                                                                                                                                                                  <  1:128695/250‑1 (MQ=255)                                                                              ggCTTGTTGCGCCAGTGCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTTAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGATTTAATGtt                                                                                                                                                            <  1:172853/251‑1 (MQ=255)                                                                                                    aGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCgtgg                                                                                                                                      <  1:277047/251‑1 (MQ=255)                                                                                                       gCCACACGTCCAGTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCTACCATGACAGTCTCCTTTTTCTGTATCGTGGA‑AATCATGTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCgtggtg                                                                                                                                    <  1:403089/250‑1 (MQ=255)                                                                                                           cacTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGTTTTTCGGGGATGGGGGTGACCATGACCGTCTCCTTTTTCTGTATAGTGGA‑AAACATTTTCTTTTTTTTTGTTTGCTCATGGAAAAGTTATTGTAATGATCCAGTGAGTTATTTTTTAACAAATTTCATGTTGCGTGGTGGTCgg                                                                                                                               >  1:183372/1‑250 (MQ=255)                                                                                                                    tGGATTCGGCAAGCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACca                                                                                                                       <  2:289180/250‑1 (MQ=255)                                                                                                                             cAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGGCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGtt                                                                                                             <  1:8975/251‑1 (MQ=255)                                                                                                                              aaTCCACTCCAGTAATTTCAGTCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGA‑AATCATTTTCATTTTTATGGTGAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGttt                                                                                                            <  1:58614/251‑1 (MQ=255)                                                                                                                                       caggtaTTTCAGGACACTCGCTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTCTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGGCTCCTGTTTCTGTAGCGTGGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCAATAGGTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTa                                                                                                   <  2:9580/246‑1 (MQ=255)                                                                                                                                            aTTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCGCGGTTTGCGGGAATGGGTGGGACCATGACAGGCTCCTTTTTTTTTATCGTGGGGAAT‑ATTTTTATTTTTTTTTTTAGCTTATGCCATAGTTTCTGGACCGATCCGATGAGTTTAAGTTTACCCAATCTCCTTGTGCGTGGGGGTCGCTTTTACCACAGATGCGTTTATGTCCGTATgggt                                                                                              >  2:99114/1‑249 (MQ=255)                                                                                                                                                    ccgcgcgTTAAAGTCTGAACCTTGCCAGTCGCTGGTGCGGTGGCGACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTACTTTGTCTGTATCGCGGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATTGTTTGTTGAAt                                                                                      <  1:371429/248‑1 (MQ=255)                                                                                                                                                    ccACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAAt                                                                                      <  2:118642/251‑1 (MQ=255)                                                                                                                                                             aaaGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGGTTTGCGGGAATGGGGGGGACCATGACAGTCTCCTTTTTCTTTATAGTGGA‑AATCCTTTTTCTTTTTATTTTTATCTTATGCCATAGTTAATGAACTGATCCGGTGAGTTAATGTTTGACAAAACTCCTGTTGCGTGGGGGTCGCTTTTACCACAAAAGCGTTTTTGCCCGTATGGTTTGTTGACttttttttt                                                                              >  1:142834/1‑246 (MQ=255)                                                                                                                                                                                         cGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAgt                                                 <  1:427180/251‑1 (MQ=255)                                                                                                                                                                                                   gacggGTTTTCTGCCGCGTTTCGCGGGAATGGGAGCCACCATGACAGTATCCTTTTTCTGTATCGTGGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGagaagc                                        >  2:305214/3‑250 (MQ=255)                                                                                                                                                                                                       ggTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGt                                   >  2:354843/1‑251 (MQ=255)                                                                                                                                                                                                             cTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGaac                              >  2:314754/1‑250 (MQ=255)                                                                                                                                                                                                                  gcgTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTttgttg                                                                                         >  1:281382/1‑186 (MQ=255)                                                                                                                                                                                                                  gcgTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGA‑AATCATTTTCATGTTTATTGTTAGCTAATGCAATAGTTACTGAACTGAGCCGATGAGTTAATGTTGAACAAATCTCATGTTGTGTGGTGGGCGGTTTTACCACAGAGGCGTTTATGCCAGTATggttggttt                                                                                         <  2:281382/186‑2 (MQ=255)                                                                                                                                                                                                                       tCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAAt                                                                                                                                                    >  2:183648/1‑122 (MQ=255)                                                                                                                                                                                                                       tCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAAt                                                                                                                                                    <  1:183648/122‑1 (MQ=255)                                                                                                                                                                                                                               aTGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGA‑AATCATTTTCATTTTTATTGTTAGCTCATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTttg                                                                                            >  1:288958/1‑170 (MQ=255)                                                                                                                                                                                                                               aTGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTttg                                                                                            <  2:288958/170‑1 (MQ=255)                                                                                                                                                                                                                                 gggTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTGAAATg         <  2:61240/251‑1 (MQ=255)                                                                                                                                                                                                                                        aCCATGACAGTCTCCTTTTTCTGTATCGTGGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGCGTTAATGTTGAACCAATCTCATGTTGCGTGGAGGTCGCTTTTACCACAGATGCGTTTCTGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCACGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTAt  >  1:430861/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                              AGTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGA‑AATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTAT  >  NC_000913/4223393‑4223874

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