Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,170,272 G→C 22.2% P35A (CCA→GCA)  ldtC ← L,D‑transpeptidase linking Lpp to murein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,170,2720GC22.2% 35.2 / 8.1 18P35A (CCA→GCA) ldtCL,D‑transpeptidase linking Lpp to murein
Reads supporting (aligned to +/- strand):  ref base G (7/7);  new base C (2/2);  total (9/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.51e-01

GAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAA  >  NC_000913/1170167‑1170380
                                                                                                         |                                                                                                            
gagCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTgg                                                                                                              >  4:272212/1‑106 (MQ=255)
gagCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTgg                                                                                                              <  3:272212/106‑1 (MQ=255)
    gagaAAGCCGACGTTGTATTTTTTGGCGATGGCTTACAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCa                                                                                                  >  4:1078807/1‑114 (MQ=255)
    gagaAAGCCGACGTTGTATTTTTTGGCGATAGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCa                                                                                                  <  3:1078807/114‑1 (MQ=255)
                      ttttttGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGcc                                                                                                   <  4:425644/95‑1 (MQ=255)
                      ttttttGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGcc                                                                                                   >  3:425644/1‑95 (MQ=255)
                              cGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGcc                                                                                  <  4:462695/104‑1 (MQ=255)
                              cGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGcc                                                                                  >  3:462695/1‑104 (MQ=255)
                                aTGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCa                                                             <  3:32366/123‑1 (MQ=255)
                                aTGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCa                                                             >  4:32366/1‑123 (MQ=255)
                                                caccGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTggcg                                                                                                            >  3:404903/1‑60 (MQ=255)
                                                caccGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTggcg                                                                                                            <  4:404903/60‑1 (MQ=255)
                                                 accGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGc                                                                                                    <  3:565889/67‑1 (MQ=255)
                                                 accGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGc                                                                                                    >  4:565889/1‑67 (MQ=255)
                                                                                              ccTGGCCGCTGCCGCCAGCGGGCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGaa                                                     <  3:333472/68‑1 (MQ=25)
                                                                                              ccTGGCCGCTGCCGCCAGCGGGCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGaa                                                     >  4:333472/2‑69 (MQ=25)
                                                                                              acTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGt                    >  1:717697/1‑102 (MQ=255)
                                                                                              acTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGt                    <  2:717697/102‑1 (MQ=255)
                                                                                                  gccGCTGCCGCCAGCGGGCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCgag                              >  2:957125/2‑88 (MQ=255)
                                                                                                  gccGCTGCCGCCAGCGGGCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCgag                              <  1:957125/87‑1 (MQ=255)
                                                                                                  cccGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTg                >  2:124565/1‑102 (MQ=255)
                                                                                                  cccGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTg                <  1:124565/102‑1 (MQ=255)
                                                                                                        ggcggcAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCaaaaa  <  3:632266/110‑1 (MQ=255)
                                                                                                        ggcggcAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCaaaaa  >  4:632266/1‑110 (MQ=255)
                                                                                                         |                                                                                                            
GAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAA  >  NC_000913/1170167‑1170380

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: