Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,614,092	T→C	I142I (ATT→ATC)	nikA →	nickel/heme ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,614,092	0	T	C	100.0%	100.4 / NA	30	I142I (ATT→ATC)	nikA	nickel/heme ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/16);  total (14/16)

CGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAACCATGAAACCATGAACGGAATTAAAGCGCCGATTGGCACCGG  >  NC_000913/3613965‑3614220                                                                                                                                |                                                                                                                                 cGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAg                                                                                                                       <  1:448990/139‑1 (MQ=255)   ttCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAgcg                                                                                                                     >  2:614515/1‑139 (MQ=255)    tCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGaaa                                                                                                                         >  2:25167/1‑134 (MQ=255)    tCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGaaa                                                                                                                         <  1:25167/134‑1 (MQ=255)                     gcAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGc                                                                                                   <  2:634886/139‑1 (MQ=255)                     gcAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGc                                                                                                   <  1:590680/139‑1 (MQ=255)                     gcAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCAGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAACGCGCCTACTATCCTTTCCttc                                                                                                   <  1:822293/139‑3 (MQ=255)                          aaaCTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGaa                                                                                              <  1:614515/139‑1 (MQ=255)                           aaCTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAAc                                                                                             <  1:1290520/139‑1 (MQ=255)                           aaCTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAAc                                                                                             <  1:65290/139‑1 (MQ=255)                                     cAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGTCCTGcc                                                                                   <  2:379742/139‑1 (MQ=255)                                     cAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGTCCTGcc                                                                                   >  1:379742/1‑139 (MQ=255)                                          cTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCGAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTc                                                                              <  1:1344675/139‑1 (MQ=255)                                                  ccGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGt                                                                      >  2:1227766/1‑139 (MQ=255)                                                     tCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGtttt                                                                   >  2:1011993/1‑139 (MQ=255)                                                         cGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCg                                                               >  2:1031064/1‑139 (MQ=255)                                                             aCGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTcc                                                           >  2:809352/1‑139 (MQ=255)                                                                                          gttgATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAaaccatgaaaccatga                              >  2:1015827/1‑139 (MQ=255)                                                                                           ttgATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTTGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAaccatgaaaccatgaa                             <  1:809352/139‑1 (MQ=255)                                                                                           ttgATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAaccatgaaaccatgaa                             <  1:1011993/139‑1 (MQ=255)                                                                                                     aGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAACCATGAAACCATGAACGGAATTaaa                   >  2:683865/1‑139 (MQ=255)                                                                                                          tCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAACCATGAAACCATGAACGGAATTAAAGCGcc              >  2:322830/1‑139 (MQ=255)                                                                                                            aGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGccc                                                                                  >  2:881009/1‑69 (MQ=255)                                                                                                            aGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGccc                                                                                  <  1:881009/69‑1 (MQ=255)                                                                                                            aGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAACCATGAAACCATGAACGGAATTAAAGCGCCGa            >  2:223155/1‑139 (MQ=255)                                                                                                             gTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAACCATGAAACCATGAACGGAATTAAAGCGCCGAt           <  1:1015827/139‑1 (MQ=255)                                                                                                                    agagCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAACCATGAAACCATGAACGGAATTAAAGCGCCGATTGGCAcc    <  1:766067/139‑1 (MQ=255)                                                                                                                      agCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTAt                                                                 <  2:639598/76‑1 (MQ=255)                                                                                                                      agCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTAt                                                                 >  1:639598/1‑76 (MQ=255)                                                                                                                      agCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAACCATGAAACCATGAACGGAATTAAAGCGCCGATTGGCACCgg  >  1:342258/1‑139 (MQ=255)                                                                                                                                |                                                                                                                                 CGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAACCATGAAACCATGAACGGAATTAAAGCGCCGATTGGCACCGG  >  NC_000913/3613965‑3614220

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