Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	43,457	Δ1 bp	9.8%	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.8%	95.7 / 7.1	41	intergenic (+23/+70)	PP_0040/cadA‑I	hypothetical protein/cadmium translocating P‑type ATPase
Reads supporting (aligned to +/- strand):  ref base T (13/21);  new base . (2/2);  total (18/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

CAAGTCGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTC‑ATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCATGGATCGCAA  >  NC_002947/43312‑43603                                                                                                                                                  |                                                                                                                                                     cAAGTCGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCt                                                                                                                                                      >  8:181634/1‑146 (MQ=255) cAAGTCGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCt                                                                                                                                                      <  7:181634/146‑1 (MQ=255)   aGTCGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATc‑g                                                                                                                                                >  7:81094/1‑149 (MQ=255)              tAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCg                                                                                                                                     <  8:22441/149‑1 (MQ=255)                    gtgtCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGtacttc                                                                                                                               >  6:307842/1‑149 (MQ=255)                       tCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGtacttctac                                                                                                                            <  8:81094/149‑1 (MQ=255)                         aTGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCGGGTA‑GCGCCCGACCCGgacatgtact                                                                                                                           >  8:117175/1‑148 (MQ=255)                              ccTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGC‑GATCAGCGCACGa                                                                                                                                         >  7:94004/1‑129 (MQ=255)                              ccTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGC‑GATCAGCGCACGa                                                                                                                                         <  8:94004/129‑1 (MQ=255)                                tCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTAGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGCGTCGTGCGCTAAGC‑GCGCACGCACCGAACTTCTAAATCAGGTCa                                                                                                                   >  4:257607/1‑149 (MQ=255)                                 cAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCTACGAGTGAGAGAAAAGTACGGGTCGTGCCCTGAAA‑GCGCACGACCCCGACATCTAATTAAGGTAAt                                                                                                                  >  1:320615/1‑149 (MQ=255)                                 cAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCGGATC‑GCGCACGACCCGTACTTCTAGTTGAGGTCAt                                                                                                                  >  6:114676/1‑149 (MQ=255)                                     gCTGCATTGGCCTGCGGCTTACGGGATGTCTCACCTGAGTGCGGAAATGACCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑TCGCACGACCCGTACTTCTACTTCAGGTCATTTcc                                                                                                              <  8:243962/149‑1 (MQ=255)                                               ccTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGtacttctactt                                                                                                                          >  2:265320/1‑127 (MQ=255)                                               ccTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGtacttctactt                                                                                                                          <  1:265320/127‑1 (MQ=255)                                                           ggCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCCCGACCCGTACATCTACTTCAGGTCATTTCCCTGCTTAGCTTC‑ATGATTGAAg                                                                                       >  5:258255/1‑149 (MQ=255)                                                            gCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTC‑ATGATTGAAGc                                                                                      <  5:307842/149‑1 (MQ=255)                                                                    cACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTC‑ATGATTGAAGCTCAGGGa                                                                               >  8:83486/1‑148 (MQ=255)                                                                      ccTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTC‑ATGATTGAAGCTCAGGGATCa                                                                            <  7:117175/149‑1 (MQ=255)                                                                       cTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTC‑ATGATTGAAGCTCAGGGATCAg                                                                           <  7:83486/149‑1 (MQ=255)                                                                                gAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATA‑GCGCACGACCCGTACTTATACTTCAGGTCATTTCCCTGCTTAGCTTA‑ATGATTTAAGCTCAGGGATCAGGTGCTTATg                                                                  >  6:257568/1‑149 (MQ=255)                                                                                       ccTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTC‑ATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAgg                                                           <  3:205567/149‑1 (MQ=255)                                                                                       ccTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATA‑GCGCACGCCCCGTACATATAATTCAGGTGAATTACCTGCTTATCTTA‑ATGATTTAAGCTAAGGGATGAGGTGCGTCTTCGGTatg                                                           >  6:298838/1‑147 (MQ=255)                                                                                            gaaGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTC‑ATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAAc                                                      <  4:126493/148‑1 (MQ=255)                                                                                               gCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTC‑ATGATTGAAGCTCAGGGATCa                                                                            >  5:319410/1‑124 (MQ=255)                                                                                               gCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTC‑ATGATTGAAGCTCAGGGATCa                                                                            <  6:319410/124‑1 (MQ=255)                                                                                                  tCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTC‑ATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGgcg                                                 >  2:234338/1‑148 (MQ=255)                                                                                                    aaGCGCATCGGGAGCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTC‑ATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCaa                                              <  3:34745/149‑1 (MQ=255)                                                                                                      gcgcATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATA‑GCGCACGAACCGTACTTCTAATTCAGGTGATTTACCTGCTTAGCTTG‑ATGATTGAAGCTAAAGGATCAGGTGCTTCTGCGGTATGCCAACAAGCGCaatg                                            >  2:223877/1‑147 (MQ=255)                                                                                                          aTCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTC‑ATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGtt                                        <  5:298838/149‑1 (MQ=255)                                                                                                           tCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTC‑ATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTg                                       <  1:164217/149‑1 (MQ=255)                                                                                                             gggATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTC‑ATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGaa                                     <  1:234338/149‑1 (MQ=255)                                                                                                                    acgaGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTC‑ATGATTGAAGCTCAGGGATCAGGTGCTTCTGc                                                                 >  6:340125/1‑114 (MQ=255)                                                                                                                    acgaGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTC‑ATGATTGAAGCTCAGGGATCAGGTGCTTCTGc                                                                 <  5:340125/114‑1 (MQ=255)                                                                                                                         tGAGAGAAAAGTACGGGTCGTGGGCGGATA‑GCGCCCGACCCCGAAGTATACTTGAGGTGATTTACCAGCTTAGCTTG‑ATGATTGAAGATCAGGGAGAAGGGGGGTGTGGGGGAGGGCAAACGGCGCAAGGCGGTGAAGACGAccca                            >  8:5903/1‑144 (MQ=255)                                                                                                                         tGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTC‑ATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACa                          <  7:5903/148‑1 (MQ=255)                                                                                                                           agagaAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTC‑ATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAgt                        <  6:258255/148‑1 (MQ=255)                                                                                                                              gaAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCGGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTC‑ATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCATCAGtgtg                     <  3:257607/148‑1 (MQ=255)                                                                                                                                      cGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTC‑ATGAATGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCa            <  4:241737/149‑1 (MQ=255)                                                                                                                                        ggTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCTCTGCTTAGCTTC‑ATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCATg          <  7:257517/149‑1 (MQ=255)                                                                                                                                          tCGTGCGC‑GATCAGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTC‑ATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGa                   >  4:300662/1‑138 (MQ=255)                                                                                                                                          tCGTGCGC‑GATCAGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTC‑ATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGa                   <  3:300662/138‑1 (MQ=255)                                                                                                                                             tGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTC‑ATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCATGGATCg     <  3:286773/149‑1 (MQ=255)                                                                                                                                                  tGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCGATGATTGAAGCTCAGGGATAAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACAAACAGTGTGGAGCCCTCATGGATCGaaa  >  1:85992/1‑148 (MQ=255)                                                                                                                                                  tGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTC‑ATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGtt                                        <  4:124074/109‑1 (MQ=255)                                                                                                                                                  tGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTC‑ATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGtt                                        >  3:124074/1‑109 (MQ=255)                                                                                                                                                  tGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTC‑ATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCATGGATCGCaa  >  6:183941/1‑147 (MQ=255)                                                                                                                                                  |                                                                                                                                                     CAAGTCGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTC‑ATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCATGGATCGCAA  >  NC_002947/43312‑43603

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