Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	43,457	Δ1 bp	8.3%	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	.	8.3%	146.7 / 7.5	48	intergenic (+23/+70)	PP_0040/cadA‑I	hypothetical protein/cadmium translocating P‑type ATPase
Reads supporting (aligned to +/- strand):  ref base T (23/21);  new base . (2/2);  total (25/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.10e-01

AGTCGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGT‑GCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCATG  >  NC_002947/43314‑43595                                                                                                                                                 |                                                                                                                                            aGTCGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGT‑GCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATc‑gc                                                                                                                                      >  1:203253/1‑150 (MQ=255) aGTCGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGT‑GCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATc‑gc                                                                                                                                      >  2:492048/1‑150 (MQ=255)   tCGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTCAGGCATCGCTGCCCTAGGTCGCGTACGGGCTATGCTTCGGC‑CTGGCTGTTCGGGATTTATGAGTGATAGAAAAGTACGGGTCGTGCGCTGATC‑gcgc                                                                                                                                    <  2:689355/150‑1 (MQ=255)    cGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGT‑GCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCa                                                                                                                                   <  2:203253/150‑1 (MQ=255)                  gtgtCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGT‑GCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATc‑gc                                                                                                                                      >  2:749561/1‑133 (MQ=255)                  gtgtCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGT‑GCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATc‑gc                                                                                                                                      <  1:749561/133‑1 (MQ=255)                       aTGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGT‑GCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGtacttctacttc                                                                                                                >  1:493584/1‑150 (MQ=255)                       aTGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGT‑GCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGtacttctacttc                                                                                                                >  1:425267/1‑150 (MQ=255)                          aGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGT‑GCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGa                                                                                                                                <  1:272794/131‑1 (MQ=255)                          aGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGT‑GCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGa                                                                                                                                >  2:272794/1‑131 (MQ=255)                                      gCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGT‑GCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGc                                                                                                 >  1:559024/1‑150 (MQ=255)                                        aTTGGCCTGCGGCTTACGGCATGGCTCACCTGAGT‑GCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCtt                                                                                               >  1:249365/1‑150 (MQ=255)                                           ggCCTGCGGCTTACGGCATGGCTCACCTGAGT‑GCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGc                                                                                            <  1:417335/150‑1 (MQ=255)                                                    cTTACGGCATGGCTCACCTGAGT‑GCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGAtt                                                                                   >  1:86204/1‑150 (MQ=255)                                                                 tCACCTGAGT‑GCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGAt                                                                      <  2:493584/150‑1 (MQ=255)                                                                 tCACCTGAGT‑GCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGC‑GATCAGCGCACGAc                                                                                                                               <  1:381265/93‑1 (MQ=38)                                                                 tCACCTGAGT‑GCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGC‑GATCAGCGCACGAc                                                                                                                               >  2:381265/1‑93 (MQ=38)                                                                                    gCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAg                                                    <  2:249365/150‑1 (MQ=255)                                                                                    gCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTCCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAg                                                    >  1:566183/1‑150 (MQ=255)                                                                                       cTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGcc                                                 >  1:350172/1‑150 (MQ=255)                                                                                          gAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTccc                                                                                                    >  1:239936/1‑96 (MQ=255)                                                                                          gAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTccc                                                                                                    <  2:239936/96‑1 (MQ=255)                                                                                           aaaGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGtacttctactt                                                                                                                 >  1:733621/1‑82 (MQ=255)                                                                                           aaaGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGtacttctactt                                                                                                                 <  2:733621/82‑1 (MQ=255)                                                                                              gCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGgc                                          <  2:559024/150‑1 (MQ=255)                                                                                               ccTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTccc                                                                                                    <  1:18215/91‑1 (MQ=255)                                                                                               ccTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTccc                                                                                                    >  2:18215/1‑91 (MQ=255)                                                                                                 tCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCa                                       >  2:225368/1‑150 (MQ=255)                                                                                                 tCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCGGCGTTGGGCCAAAAGGCGCa                                       >  2:781750/1‑150 (MQ=255)                                                                                                         aTCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGAtt                                                                                   <  2:787329/98‑1 (MQ=255)                                                                                                         aTCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGAtt                                                                                   >  1:787329/1‑98 (MQ=255)                                                                                                          tCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGa                              <  2:350172/150‑1 (MQ=255)                                                                                                             ggATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAga                           >  1:274975/1‑150 (MQ=255)                                                                                                             ggATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTCCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAga                           <  2:566183/150‑1 (MQ=255)                                                                                                                 cgacgaGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAgacgac                       <  2:274975/150‑1 (MQ=255)                                                                                                                    cgaGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCaa                    <  1:284023/150‑1 (MQ=255)                                                                                                                              aaaaGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAg                                                                          <  1:759042/86‑1 (MQ=255)                                                                                                                              aaaaGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAg                                                                          >  2:759042/1‑86 (MQ=255)                                                                                                                                aaGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAgg                                                                         >  2:165309/1‑85 (MQ=255)                                                                                                                                aaGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAgg                                                                         <  1:165309/85‑1 (MQ=255)                                                                                                                                  gTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCt      >  1:315325/1‑150 (MQ=255)                                                                                                                                      gggTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGAGCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCATg  <  1:225368/150‑1 (MQ=255)                                                                                                                                       ggTCGTGCGC‑GATCAGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGa           <  2:464062/140‑1 (MQ=255)                                                                                                                                       ggTCGTGCGC‑GATCAGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGa           >  1:464062/1‑140 (MQ=255)                                                                                                                                           gTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAgg                                                   >  1:190196/1‑96 (MQ=255)                                                                                                                                           gTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAgg                                                   <  2:190196/96‑1 (MQ=255)                                                                                                                                               gcTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGtgtg             <  1:781923/130‑1 (MQ=255)                                                                                                                                               gcTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGtgtg             >  2:781923/1‑130 (MQ=255)                                                                                                                                                 |                                                                                                                                            AGTCGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGT‑GCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCATG  >  NC_002947/43314‑43595

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