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

GTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAACCATGAAACCATGAACGGAATTAAAGCGCCGATTGGCACCGG  >  NC_000913/3613966‑3614220                                                                                                                               |                                                                                                                                 gTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAgc                                                                                                                      <  2:854811/139‑1 (MQ=255) gTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAgc                                                                                                                      <  2:133972/139‑1 (MQ=255)  ttCGATGCCGAGGCGGCGGCAGAACACTTCCGCGCACTACTCGATAACCGTCAACGTCACGCCTGGCAGGAACTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGAAAAACAGAGCGGCGACTCACCCTGAAAAgcg                                                                                                                     >  1:521892/1‑139 (MQ=255)    cGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGcc                                                                                                                   >  1:582189/1‑139 (MQ=255)     gATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCct                                                                                                                  >  2:1195071/1‑139 (MQ=255)                 gcggcAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCttt                                                                                                        >  1:1070883/1‑137 (MQ=255)                 gcggcAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCttt                                                                                                        <  2:1070883/137‑1 (MQ=255)                                   gcAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGc                                                                                    <  2:757264/139‑1 (MQ=255)                                   gcAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGc                                                                                    <  2:191550/139‑1 (MQ=255)                                      gTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCcta                                                                                                                 <  1:267241/107‑1 (MQ=255)                                      gTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCcta                                                                                                                 >  2:267241/1‑107 (MQ=255)                                       tGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGcc                                                                                        >  2:123274/1‑131 (MQ=255)                                       tGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGcc                                                                                        <  1:123274/131‑1 (MQ=255)                                                 ccGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAGTCACCCTGAAAAGCGCcta                                                                                                                 >  1:1056980/1‑96 (MQ=255)                                                 ccGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGa                                                                                                                           >  1:37274/1‑86 (MQ=255)                                                 ccGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGa                                                                                                                           <  2:37274/86‑1 (MQ=255)                                                 ccGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCcta                                                                                                                 <  2:1056980/96‑1 (MQ=255)                                                 ccGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGt                                                                      >  2:1168806/1‑139 (MQ=255)                                                   gTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGttt                                                                    >  2:537302/1‑139 (MQ=255)                                                      aaCGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAGAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTAt                                                                 >  2:607752/1‑139 (MQ=255)                                                           cACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGc                                                                                    <  2:656083/115‑1 (MQ=255)                                                           cACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGc                                                                                    >  1:656083/1‑115 (MQ=255)                                                                       gAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTa                                                <  1:475178/139‑1 (MQ=255)                                                                             gCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAg                                                                                                                       >  1:367049/1‑62 (MQ=255)                                                                             gCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAg                                                                                                                       <  2:367049/62‑1 (MQ=255)                                                                             gCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGcc                                                                                        >  1:1110504/1‑93 (MQ=255)                                                                             gCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGcc                                                                                        <  2:1110504/93‑1 (MQ=255)                                                                                          ttgATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCt                                                                            >  2:1229222/1‑92 (MQ=255)                                                                                          ttgATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCt                                                                            <  1:1229222/92‑1 (MQ=255)                                                                                          ttgATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAaccatgaaaccatgaa                             <  2:582189/139‑1 (MQ=255)                                                                                             aTGTTAAAGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAACCATGAAACCATGAACgg                          >  1:1105477/1‑139 (MQ=255)                                                                                                    aGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGcc                                                                                        >  1:1458917/1‑70 (MQ=255)                                                                                                    aGCACTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGcc                                                                                        <  2:1458917/70‑1 (MQ=255)                                                                                                       aCTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAACCATGAAACCATGAACGGAATTAAAgcg                <  1:537302/139‑1 (MQ=255)                                                                                                        cTCAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAACCATGAAACCATGAACGGAATTAAAgcgc               >  2:638157/1‑139 (MQ=255)                                                                                                          cAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTTCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTa                                                <  1:850426/104‑1 (MQ=255)                                                                                                          cAGTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTa                                                >  2:850426/1‑104 (MQ=255)                                                                                                            gTAAAACAGAGCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAACCATGAAACCATGAACGGAATTAAAGCGCCGAt           <  2:1377549/139‑1 (MQ=255)                                                                                                                   agagCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAACCATGAAACCATGAACGGAATTAAAGCGCCGATTGGCAcc    <  1:1195071/139‑1 (MQ=255)                                                                                                                     agCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCg                                                                                <  2:561893/61‑1 (MQ=255)                                                                                                                     agCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCg                                                                                >  1:561893/1‑61 (MQ=255)                                                                                                                     agCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAACCATGAAACCATGAACGGAATTAAAGCGCCGATTGGCACCgg  >  2:231608/1‑139 (MQ=255)                                                                                                                     agCTGCAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAACCATGAAACCATGAACGGAATTAAAGCGCCGATTGGCACCgg  >  2:216157/1‑139 (MQ=255)                                                                                                                          cAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCCGTTTAAAAACCATGAAACCATGAACGGAATTAAAGCGCCGAt           <  2:291887/125‑1 (MQ=255)                                                                                                                          cAAATCACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAACCATGAAACCATGAACGGAATTAAAGCGCCGAt           >  1:291887/1‑125 (MQ=255)                                                                                                                               |                                                                                                                                 GTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAGCGCCTACTATCCTTTCCTGCAAGAACTGGCCCTGCCCCGTCCTTTCCGTTTTATCGCTCCCTCGCAGTTTAAAAACCATGAAACCATGAACGGAATTAAAGCGCCGATTGGCACCGG  >  NC_000913/3613966‑3614220

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