Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,646,583 (G)6→7 coding (743/2313 nt) pbpC ← penicillin‑insensitive murein repair transglycosylase, inactive transpeptidase domain protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,646,5771.G90.0% 64.8 / 1.3 20coding (749/2313 nt)pbpCpenicillin‑insensitive murein repair transglycosylase, inactive transpeptidase domain protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (7/11);  total (8/12)
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.05e-01

GCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACTTTATTACGCGCGGCTTCGGC  >  NC_000913/2646440‑2646683
                                                                                                                                          |                                                                                                          
gCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCg                                                                                                           <  1:188812/139‑1 (MQ=255)
                   ttGCGCCAGTTCTTCCAGACGTCGTTGAAGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATggg                                                                                         <  2:677205/139‑1 (MQ=255)
                     gcgcCAGTTCTTCCAGACGTCGTTGAAGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGtt                                                                                       <  1:32153/139‑1 (MQ=255)
                     gcgcCAGTTCTTCCAGACGTCGTTGAAGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGtt                                                                                       <  1:363443/139‑1 (MQ=255)
                             tcttcCAGACGTCGTTGAAGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCtt                                                                               <  2:375419/139‑1 (MQ=255)
                                        tcgtTGAAGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTAc                                                                    <  2:86792/139‑1 (MQ=255)
                                          gtTGAAGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGttc                                                                                      >  1:432415/1‑119 (MQ=255)
                                          gtTGAAGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGttc                                                                                      <  2:432415/119‑1 (MQ=255)
                                                   cGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTAc                                                                    >  1:517006/1‑128 (MQ=255)
                                                   cGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTAc                                                                    <  2:517006/128‑1 (MQ=255)
                                                     gCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGAc                                                       <  1:710136/139‑1 (MQ=255)
                                                        tCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACcac                                                    <  2:237791/139‑1 (MQ=255)
                                                            aCGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGtt                                                                                       >  1:627324/1‑100 (MQ=255)
                                                            aCGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGtt                                                                                       <  2:627324/100‑1 (MQ=255)
                                                                  tagtGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACg                                                           >  1:687550/1‑121 (MQ=255)
                                                                  tagtGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACg                                                           <  2:687550/121‑1 (MQ=255)
                                                                            gTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCg                                >  2:195727/1‑139 (MQ=255)
                                                                                  tttgcttttACCGAGCATCATGCGCGAACACACCAGTTCCCGTCGCGGCCTTTGTCGGGGGGGCCAGCCAGATGGGCTGCTCCCTTGACTCTTTTACCTGCTCACCGCACCACACACCTCGCACTGCCACCCGTCCGAg                          >  2:385551/1‑139 (MQ=255)
                                                                                                 catcatGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGttct                                                                                     >  1:77380/1‑65 (MQ=255)
                                                                                                 catcatGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGttct                                                                                     <  2:77380/65‑1 (MQ=255)
                                                                                                      tGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACTTTATTACGCGCGGCtt      >  1:106915/1‑139 (MQ=255)
                                                                                                          cgAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACTTTATTACGCGCGGCTTCGGc  >  1:404946/1‑139 (MQ=255)
                                                                                                                                          |                                                                                                          
GCAACCGCCCTTTCCAGTTTTGCGCCAGTTCTTCCAGACGTCGTTGAAGACCGGCATCCAACGTAGTAGTGATTTTGTCGCTTTTGCTTTTACCGAGCATCATGCGCGAAAACAGCGGTGCCAGTTGCGGCATTTGTCGGGGGGCCAGCCAGATGGGTTCTTCCCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACTTTATTACGCGCGGCTTCGGC  >  NC_000913/2646440‑2646683

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