Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,614,979	A→G	H438R (CAC→CGC)	nikA →	Ni(2(+)) 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,979	0	A	G	100.0%	181.1 / NA	51	H438R (CAC→CGC)	nikA	Ni(2(+)) ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (32/19);  total (32/19)

GAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATCGCGACATTCTGACCCGTCTGCATGACGA  >  NC_000913/3614861‑3615114                                                                                                                       |                                                                                                                                        gaagaGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCa                                                                                                                    <  1:5685/140‑1 (MQ=255)      gagCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGAt                                                                                                               >  2:701461/1‑140 (MQ=255)       agCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGAtt                                                                                                              >  2:217897/1‑140 (MQ=255)       agCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGAtt                                                                                                              >  2:410247/1‑140 (MQ=255)            tATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCg                                                                                                         >  2:16941/1‑140 (MQ=255)                    cTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCg                                                                                                 <  1:841696/140‑1 (MQ=255)                      cGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCt                                                                                               >  1:10451/1‑140 (MQ=255)                      cGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCt                                                                                               >  2:755723/1‑140 (MQ=255)                      cGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCt                                                                                               >  2:713563/1‑140 (MQ=255)                      cGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCt                                                                                               <  1:567625/140‑1 (MQ=255)                        tCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCg                                                                                                 >  1:457317/1‑136 (MQ=255)                        tCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCg                                                                                                 <  2:457317/136‑1 (MQ=255)                        tCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCtga                                                                                             <  1:621526/140‑1 (MQ=255)                        tCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCtga                                                                                             >  2:224181/1‑140 (MQ=255)                        tCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCtga                                                                                             >  1:833766/1‑140 (MQ=255)                        tCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCtga                                                                                             >  1:592531/1‑140 (MQ=255)                        tCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCtga                                                                                             <  1:217897/140‑1 (MQ=255)                        tCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCtga                                                                                             <  1:410247/140‑1 (MQ=255)                           gcgcGGCCGTCGTTTTGGCCTGAATTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCtgattg                                                                                          <  2:489144/140‑1 (MQ=255)                            cgcgACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCtgattga                                                                                         >  2:647974/1‑140 (MQ=255)                                        tttGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCgg                                                                             <  1:701461/140‑1 (MQ=255)                                                           gCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCAt                                                          >  1:322363/1‑140 (MQ=255)                                                            cACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATg                                                         >  1:737603/1‑140 (MQ=255)                                                               cTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACg                                                      <  2:388904/140‑1 (MQ=255)                                                                tGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGa                                                     <  1:488251/140‑1 (MQ=255)                                                                       cgcCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCaa                                              >  2:325070/1‑140 (MQ=255)                                                                       cgcCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCaa                                              >  2:447556/1‑140 (MQ=255)                                                                         ccATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACg                                            >  2:497079/1‑140 (MQ=255)                                                                            tatGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCa                                         >  2:364714/1‑140 (MQ=255)                                                                                    acacGCCTTCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGt                                 >  1:766830/1‑140 (MQ=255)                                                                                    acacGCCTTCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGt                                 >  1:75025/1‑140 (MQ=255)                                                                                      acGCCTTCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTAt                               >  1:752337/1‑140 (MQ=255)                                                                                      acGCCTTCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTAt                               >  2:833709/1‑140 (MQ=255)                                                                                      acGCCTTCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTAt                               >  1:616172/1‑140 (MQ=255)                                                                                       cGCCTTCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATc                              >  1:105801/1‑140 (MQ=255)                                                                                         ccTTCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACG‑AACGCAACGTCAGGCGCTGTATcgcg                           <  2:833766/140‑1 (MQ=255)                                                                                          cTTCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATcgcg                           <  2:105801/140‑1 (MQ=255)                                                                                          cTTCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATcgcg                           <  1:364714/140‑1 (MQ=255)                                                                                            tCCTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATCGCGAc                         >  1:699014/1‑140 (MQ=255)                                                                                             ccTCAGTTCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATCGCGACa                        >  1:781178/1‑140 (MQ=255)                                                                                                    tCAATGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATCGCGACATTCTGAc                 <  2:766830/140‑1 (MQ=255)                                                                                                      aaTGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATCGCGACATTCTGAccc               <  1:521174/140‑1 (MQ=255)                                                                                                       aTGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATCGCGACATTCTGACCCg              <  1:713563/140‑1 (MQ=255)                                                                                                        tGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATCGCGACATTCTGACCCGt             >  1:621727/1‑140 (MQ=255)                                                                                                        tGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATCGCGACATTCTGACCCGt             >  1:609536/1‑140 (MQ=255)                                                                                                        tGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATCGCGACATTCTGACCCGt             <  1:447556/140‑1 (MQ=255)                                                                                                        tGCGCGTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATCGCGACATTCTGACCCGt             >  1:191434/1‑140 (MQ=255)                                                                                                             gTACCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATCGCGACATTCTGACCCGTCTGCa        >  1:4207/1‑140 (MQ=255)                                                                                                               aCCGTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATCGCGACATTCTGACCCGTCTGCATg      >  2:719882/1‑140 (MQ=255)                                                                                                                  gTCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATCGCGACATTCTGACCCGTCTGCATGACg   <  2:322363/140‑1 (MQ=255)                                                                                                                   tCACGCGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATCGCGACATTCTGACCCGTCTGCATGACGa  >  2:680052/1‑140 (MQ=255)                                                                                                                       |                                                                                                                                        GAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTCACACGCTGACTTCCAGGCACAGCAAGGATTAGCCGACAAACCGCTGATTGATAAAGAGATCGGCGAAGTGCTGGCGACCCATGACGAAACGCAACGTCAGGCGCTGTATCGCGACATTCTGACCCGTCTGCATGACGA  >  NC_000913/3614861‑3615114

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