Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 43,461:1 +A 25.0% intergenic (+27/+66) 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,4611.A25.0% 33.3 / 13.4 16intergenic (+27/+66)PP_0040/cadA‑Ihypothetical protein/cadmium translocating P‑type ATPase
Reads supporting (aligned to +/- strand):  ref base . (6/6);  new base A (2/2);  total (8/8)
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.59e-01

CAAGTCGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGG  >  NC_002947/43312‑43585
                                                                                                                                                      |                                                                                                                            
cAAGTCGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATc                                                                                                                               <  2:491193/150‑1 (MQ=255)
    tCGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑gcgc                                                                                                                          >  1:370699/1‑150 (MQ=255)
                                                                         aGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGGGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGc                                                     <  2:370699/150‑1 (MQ=255)
                                                                               gAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGAcc                                                                                                                    <  2:433698/81‑1 (MQ=255)
                                                                               gAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGAca                                                                                                                    >  1:433698/1‑80 (MQ=255)
                                                                                      ccTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAgg                                                                                                   <  2:389002/91‑1 (MQ=255)
                                                                                      ccTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAgg                                                                                                   >  1:389002/1‑91 (MQ=255)
                                                                                                cTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGgcgc                              >  2:273383/1‑150 (MQ=255)
                                                                                                             ggATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCAt                                                                             <  2:156597/90‑1 (MQ=255)
                                                                                                             ggATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCAt                                                                             >  1:156597/1‑90 (MQ=255)
                                                                                                                          agagaAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAgtgt    <  1:273383/150‑1 (MQ=255)
                                                                                                                            agaAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTgg  >  2:476264/1‑150 (MQ=255)
                                                                                                                                   tACGGGTCGTGCGC‑GATCAGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAg                                                         >  1:433504/1‑88 (MQ=37)
                                                                                                                                   tACGGGTCGTGCGC‑GATCAGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAg                                                         <  2:433504/88‑1 (MQ=37)
                                                                                                                                        gTCGTGCGC‑GATCAGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCtt                                                   >  1:140028/1‑89 (MQ=38)
                                                                                                                                        gTCGTGCGC‑GATCAGCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCtt                                                   <  2:140028/89‑1 (MQ=38)
                                                                                                                                          cGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTc                                                                                            >  2:309363/1‑46 (MQ=255)
                                                                                                                                          cGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTc                                                                                            <  1:309363/46‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                            
CAAGTCGGTTTCGTAGCCGGTGTCATGAGCCTCAAAGCTGCATTGGCCTGCGGCTTACGGCATGGCTCACCTGAGTGCGGAAAGGCCCTTGAAAGCCTCAAGCGCATCGGGATCGACGAGTGAGAGAAAAGTACGGGTCGTGCGCTGATC‑GCGCACGACCCGTACTTCTACTTCAGGTCATTTCCCTGCTTAGCTTCATGATTGAAGCTCAGGGATCAGGTGCTTCTGCGGTAGGCCAACAGGCGCAAGGCGTTGAAGACGACCAACAGTGTGG  >  NC_002947/43312‑43585

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