Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 43,457 Δ1 bp 19.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_00294743,4570T.19.4% 79.1 / 15.6 31intergenic (+23/+70)PP_0040/cadA‑Ihypothetical protein/cadmium translocating P‑type ATPase
Reads supporting (aligned to +/- strand):  ref base T (14/11);  new base . (3/3);  total (17/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCAT  >  NC_002947/43317‑43594
                                                                                                                                            |                                                                                                                                          
cGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCa                                                                                                                                  <  1:448667/150‑1 (MQ=255)
                    tGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGGCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCa                                                                                                              >  1:373912/1‑150 (MQ=255)
                     gAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAg                                                                                                             >  2:250950/1‑150 (MQ=255)
                              aGCTGCATTGGCCTGCGGCTTACGCCATGGCTCACCTGAGTGCGGACCGGTCCGGGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTcc                                                                                                    <  1:30266/150‑1 (MQ=255)
                                cTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCt                                                                                                  >  1:290265/1‑150 (MQ=255)
                                  gCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCGATCAGCGCACGACCCGta                                                                                                                         >  2:52744/1‑125 (MQ=255)
                                  gCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCGATCAGCGCACGACCCGta                                                                                                                         <  1:52744/125‑1 (MQ=255)
                                     ttGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCGATCAGCGCACGACCCGta                                                                                                                         <  2:486163/122‑1 (MQ=255)
                                     ttGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCGATCAGCGCACGACCCGta                                                                                                                         >  1:486163/1‑122 (MQ=255)
                                          cTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCGATCAGCGCACGAc                                                                                                                              <  1:120877/112‑1 (MQ=255)
                                          cTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCGATCAGCGCACGAc                                                                                                                              >  2:120877/1‑112 (MQ=255)
                                                    cGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAg                                                                              <  2:290265/150‑1 (MQ=255)
                                                     ggCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTTAGAGAAAAGTACTGGACGTGCGCTGATCTCGCACGACCCGTACATCTACTACAGATCATTTCCCTGCTTCGCTTCATGATTGaatc                                                                             >  2:137337/1‑148 (MQ=255)
                                                       cATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCAttt                                                                                                      >  2:366096/1‑123 (MQ=255)
                                                       cATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCAttt                                                                                                      <  1:366096/123‑1 (MQ=255)
                                                             tCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGAt                                                                     <  1:250950/150‑1 (MQ=255)
                                                                         ggAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCAt                                                                                      <  2:377505/121‑1 (MQ=255)
                                                                         ggAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCAt                                                                                      >  1:377505/1‑121 (MQ=255)
                                                                                 ccTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGTCCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTa                                                    <  1:475244/147‑1 (MQ=255)
                                                                                 ccTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTa                                                    >  2:475244/1‑147 (MQ=255)
                                                                                  cTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGt                                                                <  2:179621/134‑1 (MQ=255)
                                                                                  cTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGt                                                                >  1:179621/1‑134 (MQ=255)
                                                                                     gAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGtac                                                                                                                        <  2:155095/75‑1 (MQ=255)
                                                                                     gAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGtac                                                                                                                        >  1:155095/1‑75 (MQ=255)
                                                                                            cccatcgcATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCa                                      >  1:142160/6‑150 (MQ=255)
                                                                                                           tCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTCATCGCGCACGTCCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAgacga                       <  1:137337/150‑1 (MQ=255)
                                                                                                             gacgaGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGAcc                     >  2:240260/1‑150 (MQ=255)
                                                                                                              acgaGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCa                    >  1:365115/1‑150 (MQ=255)
                                                                                                              acgaGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCa                    >  1:365021/1‑150 (MQ=255)
                                                                                                                           aaGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGcc       >  2:328585/1‑150 (MQ=255)
                                                                                                                                cGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCAt  <  1:328585/150‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                          
CGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATCGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGGAGCCCTCAT  >  NC_002947/43317‑43594

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