Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,186,666	(C)5→6	coding (750/825 nt)	rcnA →	membrane protein conferring nickel and cobalt resistance

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,186,661	1	.	C	100.0%	158.5 / NA	44	A249P (GCC→CCC)	rcnA	membrane protein conferring nickel and cobalt resistance
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (18/26);  total (18/26)

GACACTGGGCGCAACACTGGTCGTCAGTTTCAGCATTGGCCTGGCGTTAACGCTTGTCACCGTAGGCGTTGGCGCAGCAATCAGCGTTCAGCAGGTCGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCG‑CCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGTCGGTGTGTATATGGGCGTACATGGCTTCATGGGCATAATGCGATAACTCAATCCTTTAGGCTTGCATGTTATCTTCAGCCCAGGA  >  NC_000913/2186525‑2186780                                                                                                                                          |                                                                                                                        gACACTGGGCGCAACACTGGTCGTCAGTTTCAGCATTGGCCTGGCGTTAACGCTTGTCACCGTAGGCGTTGGCGCAGCAATCAGCGTTCAGTAGGTCGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCG‑cc                                                                                                                       <  2:769808/139‑1 (MQ=255)   cacTGGGCGCAACACTGGTCGTCAGTTTCAGCATTGGCCTGGCGTTAACGCTTGTCACCGTAGGCGTTGGCGCAGCAATCAGCGTTCAGCAGGTCGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCG‑cccc                                                                                                                     <  2:470418/139‑1 (MQ=255)       gggCGCAACACTGGTCGTCAGTTTCAGCATTGGCCTGGCGTTAACGCTTGTCACCGTAGGCGTTGGCGCAGCAATCAGCGTTCAGCAGGTCGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCG‑ccccccta                                                                                                                 <  1:214380/139‑4 (MQ=255)        ggCGCAACACTGGTCGTCAGTTTCAGCATTGGCCTGGCGTTAACGCTTGTCACCGTAGGCGTTGGCGCAGCAATCAGCGTTCAGCAGGTCGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCG‑cccccctat                                                                                                                >  1:726249/1‑135 (MQ=255)         gcgcAACACTGGTCGTCAGTTTCAGCATTGGCCTGGCGTTAACGCTTGTCACCGTAGGCGTTGGCGCAGCAATCAGCGTTCAGCAGGTCGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCG‑cccccctatt                                                                                                               >  2:783220/1‑134 (MQ=255)          cgcAACACTGGTCGTCAGTTTCAGCATTGGCCTGGCGTTAACGCTTGTCACCGTAGGCGTTGGCGCAGCAATCAGCGTTCAGCAGGTCGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCGCCCCCCTAttt                                                                                                               >  1:485762/1‑139 (MQ=255)                   ggTCGTCAGTTTCAGCATTGGCCTGGCGTTAACGCTTGTCACCGTAGGCGTTGGCGCAGCAATCAGCGTTCAGCAGGTCGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCGCCCCCCTATTTTTCCAGTCt                                                                                                      >  2:56512/1‑139 (MQ=255)                       gtcAGTTTCAGCATTGGCCTGGCGTTAACGCTTGTCACCGTAGGCGTTGGCGCAGCAATCAGCGTTCAGCAGGTCGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCGCCCCCCTATTTTTCCAGTCTGTTg                                                                                                  >  2:519429/1‑139 (MQ=255)                         cAGTTTCAGCATTGGCCTGGCGTTAACGCTTGTCACCGTAGGCGTTGGCGCAGCAATCAGCGTTCAGCAGGTCGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCGCCCCCCTATTTTTCCAGTCTGTTGAt                                                                                                >  2:887121/1‑139 (MQ=255)                                aGCATTGGCCTGGCGTTAACGCTTGTCACCGTAGGCGTTGGCGCAGCAATCAGCGTTCAGCAGGTCGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTa                                                                                         <  1:887121/139‑1 (MQ=255)                                 gCATTGGCCTGGCGTTAACGCTTGTCACCGTAGGCGTTGGCGCAGCAATCAGCGTTCAGCAGGTCGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTAg                                                                                        <  2:169470/139‑1 (MQ=255)                                  cATTGGCCTGGCGTTAACGCTTGTCACCGTAGGCGTTGGCGCAGCAATCAGCGTTCAGCAGGTCGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGt                                                                                       <  1:519429/139‑1 (MQ=255)                                  cATTGGCCTGGCGTTAACGCTTGTCACCGTAGGCGTTGGCGCAGCAATCAGCGTTCAGCAGGTCGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGt                                                                                       <  2:485762/139‑1 (MQ=255)                                    ttGGCCTGGCGTTAACGCTTGTCACCGTAGGCGTTGGCGCAGCAATCAGCGTTCAGCAGGTCGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGTCg                                                                                     <  1:575051/139‑1 (MQ=255)                                    ttGGCCTGGCGTTAACGCTTGTCACCGTAGGCGTTGGCGCAGCAATCAGCGTTCAGCAGGTCGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGTCg                                                                                     <  2:484532/139‑1 (MQ=255)                                    ttGGCCTGGCGTTAACGCTTGTCACCGTAGGCGTTGGCGCAGCAATCAGCGTTCAGCAGGTCGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGTCg                                                                                     <  1:150262/139‑1 (MQ=255)                                      ggCCTGGCGTTAACGCTTGTCACCGTAGGCGTTGGCGCAGCAATCAGCGTTCAGCAGGTCGCAAAACGCTGGAGCGGATTTAACACTCTCGGTAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGTCGgt                                                                                   >  2:575903/1‑139 (MQ=255)                                      ggCCTGGCGTTAACGCTTGTCACCGTAGGCGTTGGCGCAGCAATCAGCGTTCAGCAGGTCGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGTCGgt                                                                                   <  2:898937/139‑1 (MQ=255)                                      ggCCTGGCGTTAACGCTTGTCACCGTAGGCGTTGGCGCAGCAATCAGCGTTCAGCAGGTCGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGTCGgt                                                                                   <  1:275750/139‑1 (MQ=255)                                      ggCCTGGCGTTAACGCTTGTCACCGTAGGCGTTGGCGCAGCAATCAGCGTTCAGCAGGTCGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGTCGgt                                                                                   <  1:850881/139‑1 (MQ=255)                                      ggCCTGGCGTTAACGCTTGTCACCGTAGGCGTTGGCGCAGCAATCAGCGTTCAGCAGGTCGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGTCGgt                                                                                   <  2:364508/139‑1 (MQ=255)                                             cGTTAACGCTTGTCACCGTAGGCGTTGGCGCAGCAATCAGCGTTCAGCAGGTCGCAAAACGCTGGAGCGGATTTAACACTCTCGGTAAACGCGCCCCCCTATTTTTCCAGACTGTTGATTGGCTTAGTCGGTGTGtata                                                                            <  1:575903/139‑1 (MQ=255)                                             cGTTAACGCTTGTCACCGTAGGCGTTGGCGCAGCAATCAGCGTTCAGCAGGTCGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGTCGGTGTGtata                                                                            <  1:56512/139‑1 (MQ=255)                                                         tCACCGTAGGCGTTGGCGCAGCAATCAGCGTTCAGCAGGTCGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGTCGGTGTGTATATGGGCGTACATg                                                                <  2:726249/139‑1 (MQ=255)                                                                aGGCGTTGGCGCAGCAATCAGCGTTCAGCAGGTCGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGTCGgt                                                                                   <  2:382068/113‑1 (MQ=255)                                                                aGGCGTTGGCGCAGCAATCAGCGTTCAGCAGGTCGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGTCGgt                                                                                   >  1:382068/1‑113 (MQ=255)                                                                  gCGTTGGCGCAGCAATCAGCGTTCAGCAGGTCGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGTCGGTGTGTATATGGGCGTACATGGCTTCATgg                                                       >  2:871189/1‑139 (MQ=255)                                                                   cGTTGGCGCAGCAATCAGCGTTCAGCAGGTCGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGTCGGTGTGTATATGGGCGTACATGGCTTCATggg                                                      <  1:393974/139‑1 (MQ=255)                                                                    gTTGGCGCAGCAATCAGCGTTCAGCAGGTCGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGt                                                                                       >  2:24985/1‑105 (MQ=255)                                                                    gTTGGCGCAGCAATCAGCGTTCAGCAGGTCGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGt                                                                                       <  1:24985/105‑1 (MQ=255)                                                                             gcaATCAGCGTTCAGCAGGTCGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGTCGGTGTGTATATGGGCGTACATGGCTTCATGGGCATAATGCGa                                            <  1:871189/139‑1 (MQ=255)                                                                                       ttCAGCAGGTCGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGTCGGTGTGTATATGGGCGTACATGGCTTCATGGGCATAATGCGATAACTCAATc                                  >  1:833836/1‑139 (MQ=255)                                                                                        tCAGCAGGTCGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGTCGGTGTGTATATGGGCGTACATGGCTTCATGGGCATAATGCGATAACTCAATcc                                 >  2:609267/1‑139 (MQ=255)                                                                                          agcagGTCGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGTCGGTGTGTATATGGGCGTACATGGCTTCATGGGCATAATGCGATAACTCaa                                    >  1:632404/1‑134 (MQ=255)                                                                                          agcagGTCGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGTCGGTGTGTATATGGGCGTACATGGCTTCATGGGCATAATGCGATAACTCaa                                    <  2:632404/134‑1 (MQ=255)                                                                                               gTCGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGTCGGTGTGTATATGGGCGTAc                                                                   <  2:355681/98‑1 (MQ=255)                                                                                               gTCGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGTCGGTGTGTATATGGGCGTAc                                                                   >  1:355681/1‑98 (MQ=255)                                                                                               gTCGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGTCGGTGTGTATATGGGCGTACATGGCTTCATGGGCATAATGCGATAACTCAATCCTTTAGGc                          >  2:508434/1‑139 (MQ=255)                                                                                                 cGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGt                                                                                       <  1:445679/76‑1 (MQ=255)                                                                                                 cGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGt                                                                                       >  2:445679/1‑76 (MQ=255)                                                                                                 cGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGTCGgt                                                                                   <  2:505442/80‑1 (MQ=255)                                                                                                 cGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGTCGgt                                                                                   >  1:505442/1‑80 (MQ=255)                                                                                                                       aaCACTCTCGCTAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGTCGGTGTGTATATGGGCGTACATGGCTTCATGGGCATAATGCGATAACTCAATCCTTTAGGCTTGCATGTTATCTTCAGCCCAGga  <  1:783220/139‑1 (MQ=255)                                                                                                                             ctcGCTAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGTCGGTGTGTATATGGGCGTACATGGCTTCATGGGCATAATGCGATAACTCAATCCTTTAGGCTTGCATGTTATCTTCAGCCCAGga  <  2:298306/133‑1 (MQ=255)                                                                                                                             ctcGCTAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGTCGGTGTGTATATGGGCGTACATGGCTTCATGGGCATAATGCGATAACTCAATCCTTTAGGCTTGCATGTTATCTTCAGCCCAGga  >  1:298306/1‑133 (MQ=255)                                                                                                                              tcGCTAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGt                                                                                       >  1:283619/1‑47 (MQ=255)                                                                                                                              tcGCTAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGt                                                                                       <  2:283619/47‑1 (MQ=255)                                                                                                                                  tAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGt                                                                                       <  1:93285/43‑1 (MQ=255)                                                                                                                                  tAAACGCGCCCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGt                                                                                       >  2:93285/1‑43 (MQ=255)                                                                                                                                          |                                                                                                                        GACACTGGGCGCAACACTGGTCGTCAGTTTCAGCATTGGCCTGGCGTTAACGCTTGTCACCGTAGGCGTTGGCGCAGCAATCAGCGTTCAGCAGGTCGCAAAACGCTGGAGCGGATTTAACACTCTCGCTAAACGCG‑CCCCCTATTTTTCCAGTCTGTTGATTGGCTTAGTCGGTGTGTATATGGGCGTACATGGCTTCATGGGCATAATGCGATAACTCAATCCTTTAGGCTTGCATGTTATCTTCAGCCCAGGA  >  NC_000913/2186525‑2186780

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: ‑