Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	43,457	Δ1 bp	9.4%	intergenic (+23/+70)	PP_0040 → / ← cadA‑I	hypothetical protein/cadmium translocating P‑type ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	43,457	0	T	.	9.4%	126.1 / 10.9	53	intergenic (+23/+70)	PP_0040/cadA‑I	hypothetical protein/cadmium translocating P‑type ATPase
Reads supporting (aligned to +/- strand):  ref base T (19/27);  new base . (3/2);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 6.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.88e-01

AAGTCGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCATGGATC  >  NC_002947/43313‑43599                                                                                                                                                 |                                                                                                                                                aaGTCGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATc                                                                                                                                             >  2:270643/1‑149 (MQ=255) aaGTCGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATc                                                                                                                                             >  3:255694/1‑149 (MQ=255)   gTCGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATc‑gc                                                                                                                                          >  1:331161/1‑149 (MQ=255)      ggTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCa                                                                                                                                       <  2:77277/149‑1 (MQ=255)       gTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGC‑GATCAGCGCAc                                                                                                                                      <  3:90326/149‑1 (MQ=255)           cGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGAc                                                                                                                                   <  1:339387/148‑1 (MQ=255)             tAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGGGGTGTTC‑GCGCACGGcccc                                                                                                                                >  2:203248/1‑148 (MQ=255)                          gAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGtacttctactt                                                                                                                     >  6:101106/1‑147 (MQ=255)                                 aaaGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCGGATA‑GCCCACGACCCGTACTTCTACTTCAGGTCAtt                                                                                                            >  2:116752/1‑149 (MQ=255)                                   aGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGAAC‑GCGCCCGACCCGTACGTGTACTTCAGGtacttta                                                                                                          >  2:86713/1‑143 (MQ=255)                                      tGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATA‑GCGCACGACCCGTACTTCTACTTCAGGCCATTTACCt                                                                                                       >  4:237236/1‑149 (MQ=255)                                        cATTGGCCTGCGGCTTACGGCATGGCTCACATGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTTAGAGAAAAGTACGGGTCGTGAGC‑GATCAGCGCACGACCCGtacttttatct                                                                                                                     >  1:358850/1‑130 (MQ=37)                                           tGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCt                                                                                                    >  5:103370/1‑147 (MQ=255)                                            ggCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTa                                                                                                  <  2:331161/148‑1 (MQ=255)                                                  cGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTc                                                                                             <  7:310133/147‑1 (MQ=255)                                                  cGGATTACGGCATGGCTCACCTGAGTGCTGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACTGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCAt                                                                                           <  7:313078/149‑1 (MQ=255)                                                      ttACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGc                                                                                                     <  2:263875/135‑1 (MQ=255)                                                      ttACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCACAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGc                                                                                                     >  1:263875/1‑135 (MQ=255)                                                        aCGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTg                                                                                      <  4:255694/148‑1 (MQ=255)                                                         cGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGATGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGaa                                                                                    <  1:116752/149‑1 (MQ=255)                                                         cGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGaa                                                                                    <  1:203248/149‑1 (MQ=255)                                                                    aCCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGAt                                                                          >  6:184924/1‑148 (MQ=255)                                                                      cTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGAGAGAGAAAAGTACGGGGCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCa                                                                                                                   <  7:173584/105‑1 (MQ=255)                                                                         aGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAg                                                                                                                  >  2:261523/1‑103 (MQ=255)                                                                         aGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAg                                                                                                                  <  1:261523/103‑1 (MQ=255)                                                                             cGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTc                                                                >  7:12529/1‑149 (MQ=255)                                                                                    gCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTc                                                                                                          <  7:66061/100‑1 (MQ=255)                                                                                    gCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTc                                                                                                          >  8:66061/1‑100 (MQ=255)                                                                                     cccTTGAAAGCCTCAAGCGCGTCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGGGCGGATA‑GCGCCCCCCCCCGACTTGTAATTAAGGGGATTTTCCAGCGTTTCTTGATTATTTAAGAGCAGGGATCAGGTGGTTGTGGGGTAg                                                        >  2:319886/1‑149 (MQ=255)                                                                                      ccTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGt                                                                     >  4:266965/1‑135 (MQ=255)                                                                                      ccTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGt                                                                     <  3:266965/135‑1 (MQ=255)                                                                                       cTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGc                                                      >  5:196100/1‑149 (MQ=255)                                                                                        ttGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATA‑GCGCACGACCCGTACTTATACTTCAGGTCATTTACCTGCTTAGCTTCATGATTGAAGCTGAGGGATAACGTGCTTCTGCGGTAGGc                                                      >  6:154025/1‑148 (MQ=255)                                                                                          gAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGcc                                                     <  8:12529/147‑1 (MQ=255)                                                                                               ccTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGgc                                              <  5:154025/149‑1 (MQ=255)                                                                                                cTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACCACCCGTACTTCTACTTCAGGTCATTTACATTCTTAGCTTCAAGATTGAAGCTCAGGGATAAGGTGCTTCTGCGGTAGGCCAAAAGgcg                                             >  7:81156/1‑149 (MQ=255)                                                                                                   aaGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCaa                                          <  2:159549/149‑1 (MQ=255)                                                                                                        cATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGt                                     <  1:86713/149‑1 (MQ=255)                                                                                                            gggATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGCCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGa                                  <  3:237236/148‑1 (MQ=255)                                                                                                              gATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGCATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAga                               <  1:270643/149‑1 (MQ=255)                                                                                                                   acgaGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGgcgc                                            <  2:46424/131‑1 (MQ=255)                                                                                                                   acgaGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGgcgc                                            >  1:46424/1‑131 (MQ=255)                                                                                                                          agagaAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAgtg                   <  6:196100/149‑1 (MQ=255)                                                                                                                            agaAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGtgtg                 <  5:184924/149‑1 (MQ=255)                                                                                                                               gaaGTACGGGTCGTGCGC‑GATCAGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAg              >  4:232058/2‑149 (MQ=255)                                                                                                                                aaGTACGGGTCGTGCGCTGATC‑GCGCACGACCCCTACTTCTACTTCAGGTCATTTACCTGCTTATCTTCAATATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACACGCGCAACGCCGTGAATCCCACCAACAGAGTGGAGc             >  6:95482/1‑149 (MQ=255)                                                                                                                                      gggTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCAt       <  5:95482/149‑1 (MQ=255)                                                                                                                                        gTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTg                                   <  5:8243/119‑1 (MQ=255)                                                                                                                                        gTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCATgg     <  2:159395/149‑1 (MQ=255)                                                                                                                                        gTCGTGCGC‑GATCAGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAgt                    <  7:151084/134‑1 (MQ=255)                                                                                                                                        gTCGTGCGC‑GATCAGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAgt                    >  8:151084/1‑134 (MQ=255)                                                                                                                                         tCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCATGGa    <  5:238093/149‑1 (MQ=255)                                                                                                                                           gTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCATGGATc  <  6:103370/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                AAGTCGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCATGGATC  >  NC_002947/43313‑43599

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑