Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,617,376	C→A	A121D (GCC→GAC)	nikD →	Ni(2(+)) ABC transporter ATP binding subunit NikD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,617,376	0	C	A	100.0%	170.3 / NA	49	A121D (GCC→GAC)	nikD	Ni(2(+)) ABC transporter ATP binding subunit NikD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (32/17);  total (32/17)

CATCAAAATTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCT‑GGCGTTAGGGAAA‑CCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTATCATCGCCG  >  NC_000913/3617248‑3617507                                                                                                                                   |                                                                                                                                    aataaaaaTTCCCACCACCAGCCAGACCCCGCGCAGCCCTTTAAACCCCCCCCCCCCCACCCCCCCCCCCGCGGGGAAACCCTGCGG‑GGTGTTAGGAAAC‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCgg                                                                                                                          <  2:380463/136‑1 (MQ=255)  ataaaaaTCCCCCCCACCACGCAGCCCCCGCGCGGCCCTTTTAACCCACGCCCCCCCGCCCCCCCCCCCGCGGGGAAACCCGGCGG‑GGTGTGAGGAAAC‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGt                                                                                                                         <  2:541360/137‑1 (MQ=255)      aaaTTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCT‑GGCGTTAGGGAAA‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTgggg                                                                                                                     >  1:380463/1‑140 (MQ=255)              accaTCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCT‑GGCGTTAGGGAAA‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGaaac                                                                                                             >  1:541360/1‑139 (MQ=255)              aaaatcatGCAAACCCCGCGCAGCCCTTTTAACCCCCTCCACACCATCCCCCCCCACGCG‑GGGAAACCCGGCTGGGGGTTAGGAAAC‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGaaaa                                                                                                             <  2:477314/137‑1 (MQ=255)                     tGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCT‑GGCGTTAGGGAAA‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCcgcg                                                                                                      >  1:137758/1‑140 (MQ=255)                      gcagcccccGCGCAGCGCTTTTAACCCCCCGCCCCCCACCCCCCCCCCCGGGCGGAAACCCTGCCG‑GGGGTGAGGAAAC‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCcgcgc                                                                                                     <  2:785414/134‑1 (MQ=255)                         gaccccGCGCGGCCCTTTTAACCCCCTCCACACCATCCCCCCCCCCGCGGGGAAACCCTGCCG‑GGGGTTGGGGAAC‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCcgcgcgcg                                                                                                  <  2:202759/137‑1 (MQ=255)                            cccGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCT‑GGCGTTAGGGAAA‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCcgcgcgcgctg                                                                                               >  1:477314/1‑137 (MQ=255)                            cccGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCT‑GGCGTTAGGGAAA‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGc                                                                                               >  1:958998/1‑140 (MQ=255)                            cccGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCT‑GGCGTTAGGGAAA‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGc                                                                                               >  1:1044281/1‑140 (MQ=255)                            cccGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCT‑GGCGTTAGGGAAA‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGCGGGGCTGGAAAACGCCGCGCGCGTGc                                                                                               >  1:390359/1‑140 (MQ=255)                               gcgcAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCT‑GGCGTTAGGGAAA‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGa                                                                                            >  2:296588/1‑140 (MQ=255)                                cgcAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCT‑GGCGTTAGGGAAA‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGaa                                                                                           >  2:289482/1‑140 (MQ=255)                                cgcAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCT‑GGCGTTAGGGAAA‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGaa                                                                                           >  1:625983/1‑140 (MQ=255)                                    gcgcCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCT‑GGCGTTAGGGAAA‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTg                                                                                       >  1:202759/1‑140 (MQ=255)                                     cgcCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCT‑GGCGTTAGGGAAA‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGTGTGCTGAAGCTGt                                                                                      >  1:564302/1‑140 (MQ=255)                                      cccTTTACCCCCCTCCCCACCACGCCCCCCCGCGCGC‑GGAAACCCGGCG‑GGGGGTGGGGAAACCCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTa                                                                                     <  2:605006/139‑1 (MQ=255)                                                  cTGCACACCATGCACACCCACGCGCGGGAAACCTGCCT‑GGCGTTAGGGAAA‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGAt                                                                         >  1:604790/1‑140 (MQ=255)                                                     cacaCCATGCACACCCACGCGCGGGAAACCTGCCT‑GGCGTTAGGGAAA‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAg                                                                      >  2:323009/1‑140 (MQ=255)                                                     cacaCCATGCACACCCACGCGCGGGAAACCTGCCT‑GGCGTTAGGGAAA‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAg                                                                      >  1:602866/1‑140 (MQ=255)                                                     cacaCCATGCACACCCACGCGCGGGAAACCTGCCT‑GGCGTTAGGGAAA‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAg                                                                      >  1:644355/1‑140 (MQ=255)                                                            tGCACACCCACGCGCGGGAAACCTGCCT‑GGCGTTAGGGAAA‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCggcggc                                                               >  1:579035/1‑140 (MQ=255)                                                            tGCACACCCACGCGCGGGAAACCTGCCT‑GGCGTTAGGGAAA‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCggcggc                                                               >  1:548129/1‑140 (MQ=255)                                                            tGCACACCCACGCGCGGGAAACCTGCCT‑GGCGTTAGGGAAA‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCggcggc                                                               >  1:135995/1‑140 (MQ=255)                                                            tGCACACCCACGCGCGGGAAACCTGCCT‑GGCGTTAGGGAAA‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGT‑GGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGGACCCGGTCGAGATGAGcggcggga                                                              >  1:42551/1‑138 (MQ=255)                                                             gCACACCCACGCGCGGGAAACCTGCCT‑GGCGTTAGGGAAA‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCa                                                              >  1:880929/1‑140 (MQ=255)                                                             gCACACCCACGCGCGGGAAACCTGCCT‑GGCGTTAGGGAAA‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCa                                                              >  1:643944/1‑140 (MQ=255)                                                             gCACACCCACGCGCGGGAAACCTGCCT‑GGCGTTAGGGAAA‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCa                                                              >  1:567945/1‑140 (MQ=255)                                                             cccccccccCGCGGGGGAACCCTGCTG‑GGCGTGAGGAAAC‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCa                                                              <  2:401504/135‑1 (MQ=255)                                                                 acCCACGCGCGGGAAACCTGCCT‑GGCGTTAGGGAAA‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGtt                                                          >  2:510878/1‑140 (MQ=255)                                                                  cccGCGCGGGGGAATCTTGCCT‑GGGTTAAGGGAAA‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTg                                                         <  1:4025/140‑1 (MQ=255)                                                                    cACGCGCGGGAAACCTGCCT‑GGCGTTAGGGAAA‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCa                                                       >  2:1081602/1‑140 (MQ=255)                                                                    cACGCGCGGGAAACCTGCCT‑GGCGTTAGGGAAA‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCa                                                       >  2:120689/1‑140 (MQ=255)                                                                    cACGCGCGGGAAACCTGCCT‑GGCGTTAGGGAAA‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCa                                                       >  1:173141/1‑140 (MQ=255)                                                                      cgcgcgGGAAACCTGCCT‑GGCGTTAGGGAAA‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAgc                                                     >  1:709058/1‑140 (MQ=255)                                                                      cgcgcgGGAAACCTGCCT‑GGCGTTAGGGAAA‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAgc                                                     >  1:266124/1‑140 (MQ=255)                                                                       gcgggGGAAGCCTGCTT‑GGGGGTAGGGAAC‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAgcg                                                    <  2:604790/140‑1 (MQ=255)                                                                       gcgggGGAACCCTGCC‑‑GGGGTTAAGGAAACCCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAgcg                                                    <  1:933473/140‑1 (MQ=255)                                                                         gcgGGAAACCTGCCT‑GGCGTTAGGGAAA‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCa                                                  >  1:925047/1‑140 (MQ=255)                                                                                     cct‑gGCGTTAGGGAAA‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAgcg                                                    >  2:119916/1‑126 (MQ=255)                                                                                     cct‑gGCGTTAGGGAAA‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAgcg                                                    <  1:119916/126‑1 (MQ=255)                                                                                         gggtatAGGGAAACCCCGCCGATGACGCTACGCTTACCGCTGAAATAGAAGCGGT‑GGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATgg                                   <  2:42551/135‑1 (MQ=255)                                                                                           ggtgaGGGAAC‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCg                                 <  2:548129/136‑1 (MQ=255)                                                                                              gagccata‑CCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAGCGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTg                              <  2:135995/133‑1 (MQ=255)                                                                                               ggggAAACCCCGCCGATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTg                              <  2:255063/139‑1 (MQ=255)                                                                                                           ccgATGACGCTACGCTTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCa                  >  1:660122/1‑140 (MQ=255)                                                                                                                      acgctTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTatcat       <  2:643944/140‑1 (MQ=255)                                                                                                                      acgctTACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTatcat       <  2:660122/140‑1 (MQ=255)                                                                                                                          ttACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTATCATCGcc   <  2:173141/140‑1 (MQ=255)                                                                                                                           tACCGCTGACATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTATCATCGCCg  <  1:120689/140‑1 (MQ=255)                                                                                                                                   |                                                                                                                                    CATCAAAATTGCCACCATCATGCAGAACCCGCGCAGCGCCTTTAATCCACTGCACACCATGCACACCCACGCGCGGGAAACCTGCCT‑GGCGTTAGGGAAA‑CCCGCCGATGACGCTACGCTTACCGCTGCCATAGAAGCGGTGGGGCTGGAAAACGCCGCGCGCGTGCTGAAGCTGTACCCGTTCGAGATGAGCGGCGGCATGTTGCAGCGCATGATGATTGCGATGGCGGTGCTGTGTGAATCACCGTTTATCATCGCCG  >  NC_000913/3617248‑3617507

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