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

TTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAACCATGAAACCATGAACGGAATTAAAGCGCCGATTGGCACCGGACCGT  >  NC_000913/3613967‑3614225                                                                                                                              |                                                                                                                                      ttCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACATAGCTGCAAATCACCCTGAAAAgcg                                                                                                                          >  2:474860/1‑139 (MQ=255)   cGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGcc                                                                                                                        <  1:295190/139‑1 (MQ=255)   cGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGcc                                                                                                                        >  1:668530/1‑139 (MQ=255)           aGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATc                                                                                                                 <  1:234659/138‑1 (MQ=255)           aGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATc                                                                                                                 >  2:234659/1‑138 (MQ=255)                gcggcAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCttt                                                                                                             <  2:392577/137‑1 (MQ=255)                gcggcAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCttt                                                                                                             >  1:392577/1‑137 (MQ=255)                        aaaCTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGaa                                                                                                   <  1:150806/139‑1 (MQ=255)                        aaaCTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGaa                                                                                                   <  2:295971/139‑1 (MQ=255)                              ccGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCcta                                                                                                                      >  2:656642/1‑114 (MQ=255)                              ccGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCcta                                                                                                                      <  1:656642/114‑1 (MQ=255)                              ccGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGcc                                                                                             <  1:1058602/139‑1 (MQ=255)                              ccGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGcc                                                                                             <  2:668530/139‑1 (MQ=255)                              ccGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGcc                                                                                             <  1:1303210/139‑1 (MQ=255)                                                ccGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGcc                                                                                                                        <  2:978991/94‑1 (MQ=255)                                                ccGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGcc                                                                                                                        >  1:978991/1‑94 (MQ=255)                                                ccGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAAc                                                                                                  <  1:1057736/116‑1 (MQ=255)                                                ccGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAAc                                                                                                  >  2:1057736/1‑116 (MQ=255)                                                                         cTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCttt                                                                               <  2:1223783/110‑1 (MQ=255)                                                                         cTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCttt                                                                               >  1:1223783/1‑110 (MQ=255)                                                                                  ataagagacagtGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAaaaccatgaaa                                         >  1:1092105/12‑139 (MQ=255)                                                                                         ttgATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAaccatgaaaccatgaa                                  <  2:1189603/139‑1 (MQ=255)                                                                                         ttgATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAaccatgaaaccatgaa                                  >  1:1189603/1‑139 (MQ=255)                                                                                                   aacACTCAGTAAAACAGAGCTGCAAATCACCCTAAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGct                                                                                        <  1:827368/73‑2 (MQ=255)                                                                                                   aGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGcc                                                                                        >  2:827368/1‑75 (MQ=255)                                                                                                       cTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCt                                                                                 >  1:797612/1‑78 (MQ=255)                                                                                                       cTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCt                                                                                 <  2:797612/78‑1 (MQ=255)                                                                                                       cTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAaccatgaaaccatgaa                                  <  2:789684/125‑1 (MQ=255)                                                                                                       cTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAaccatgaaaccatgaa                                  >  1:789684/1‑125 (MQ=255)                                                                                                        tCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGATCTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAaccatgaaaccatgaa                                  <  1:746878/124‑1 (MQ=255)                                                                                                        tCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAaccatgaaaccatgaa                                  >  2:746878/1‑124 (MQ=255)                                                                                                          aGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAACCATGAAACCATGAACGGAATTAAAGCGCCGa                 >  1:742036/1‑139 (MQ=255)                                                                                                          aGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAACCATGAAACCATGAACGGAATTAAAGCGCCGa                 >  1:723970/1‑139 (MQ=255)                                                                                                                    agCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGcc                                                                                             >  1:347381/1‑53 (MQ=255)                                                                                                                    agCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGcc                                                                                             <  2:347381/53‑1 (MQ=255)                                                                                                                    agCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAACCATGAAACCATGAACGGAATTAAAGCGCCGAt                >  1:224999/1‑130 (MQ=255)                                                                                                                    agCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAACCATGAAACCATGAACGGAATTAAAGCGCCGAt                <  2:224999/130‑1 (MQ=255)                                                                                                                      cTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTaaa                                                   >  1:778379/1‑93 (MQ=255)                                                                                                                      cTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAACCATGAAACCATGAACGGAATTAAAGCGCCGATTGGCACCGGAc     >  2:1024508/1‑139 (MQ=255)                                                                                                                      cTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTAGCCCTGCCCCGTCCTGTCCGTTTTATCGCTCCCTCGCAGTTTaaa                                                   <  2:778379/93‑1 (MQ=255)                                                                                                                         cAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAACCATGAAACCATGAACGGAATTAAAGCGCCGATTGGCACCGGACCGt  <  2:742036/139‑1 (MQ=255)                                                                                                                              |                                                                                                                                      TTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAACCATGAAACCATGAACGGAATTAAAGCGCCGATTGGCACCGGACCGT  >  NC_000913/3613967‑3614225

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