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%	59.7 / NA	18	I142I (ATT→ATC)	nikA	nickel/heme ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (9/9);  total (9/9)

TTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAACCATGAAACCATGAACGGAATTAAAGCGCCGATTGGCACCGGACC  >  NC_000913/3613967‑3614223                                                                                                                              |                                                                                                                                    ttCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAgcg                                                                                                                        >  2:293253/1‑139 (MQ=255)                gcggcAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTcc                                                                                                         >  1:139229/1‑139 (MQ=255)                             tCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGc                                                                                            <  1:141993/139‑1 (MQ=255)                                      tGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCg                                                                                   >  1:798492/1‑139 (MQ=255)                                                   tCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGtttt                                                                      >  1:28448/1‑139 (MQ=255)                                                                tggctggAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGc                                                         <  1:293253/139‑1 (MQ=255)                                                                tggctggAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGc                                                         <  1:672504/139‑1 (MQ=255)                                                                                      ttgttgATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAgc                                                                                                                         >  2:577210/1‑53 (MQ=255)                                                                                      ttgttgATGTTAAAGCACTCAGTAAAACAGAGCTGAAAATCACCCTGAAAAgc                                                                                                                         <  1:577210/53‑1 (MQ=39)                                                                                                       cTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAACCATGAAACCATGAACGGAATTAAAgcgc                  >  2:644180/1‑139 (MQ=255)                                                                                                           gTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAACCATGAAACCATGAACGGAATTAAAGCGCCGAt              <  2:28448/139‑1 (MQ=255)                                                                                                                    agCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAaaaccatg                                          >  1:379768/1‑102 (MQ=255)                                                                                                                    agCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAaaaccatg                                          <  2:379768/102‑1 (MQ=255)                                                                                                                      cTGCAAATCACCCTGAAAAGCGCCTACTATTCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAaccatgaaaccatgaa                                <  1:482979/110‑1 (MQ=255)                                                                                                                      cTGCAAATCACCCTGAAAAGCGCCTACTATTCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAaccatgaaaccatgaa                                >  2:482979/1‑110 (MQ=255)                                                                                                                      cTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAaccatgaaaccatgaa                                >  2:867772/1‑110 (MQ=255)                                                                                                                      cTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAaccatgaaaccatgaa                                <  1:867772/110‑1 (MQ=255)                                                                                                                       tGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCGCCCTCGCAGTTTAAAAACCATGAAACCATGAACGGAATTAAAGCGCCGATTGGCACCGGAcc  <  2:798492/139‑1 (MQ=255)                                                                                                                              |                                                                                                                                    TTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAACCATGAAACCATGAACGGAATTAAAGCGCCGATTGGCACCGGACC  >  NC_000913/3613967‑3614223

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