Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,645,013 C→T *771* (TAG→TAA 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,645,0130CT97.4% 122.6 / ‑6.6 39*771* (TAG→TAApbpCpenicillin‑insensitive murein repair transglycosylase, inactive transpeptidase domain protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (18/20);  minor base G (0/1);  total (18/21)
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

TGTACCAGCGTACGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTAACGGTTCGCCATTCAGAAACCAC  >  NC_000913/2644896‑2645148
                                                                                                                     |                                                                                                                                       
tGTACCAGCGTACGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACt                                                                                                                    >  1:468560/1‑139 (MQ=255)
  tACCAGCGTACGAGATGGCGCGGATTATACTGAGCACCAGGGCCGTTGAACATTAATGAAGGTAACTTTTTTTTAAAATAATATGAGTTTTATCACCAATCAAGGACAAAACATATAATTGCATGTCAAATTTCAgtgt                                                                                                                  >  1:320773/1‑139 (MQ=255)
        cGTACGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGAt                                                                                                            >  1:214460/1‑139 (MQ=255)
        cGTACGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGAt                                                                                                            <  1:680795/139‑1 (MQ=255)
         gTACGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGAtt                                                                                                           <  2:630182/139‑1 (MQ=255)
             gagaTGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGcc                                                                                                       <  2:600451/139‑1 (MQ=255)
             gagaTGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGcc                                                                                                       <  2:214460/139‑1 (MQ=255)
             gagaTGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGcc                                                                                                       <  2:141392/139‑1 (MQ=255)
               gaTGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGcccc                                                                                                     >  2:630633/1‑139 (MQ=255)
                  ggCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACg                                                                                                  >  1:449812/1‑139 (MQ=255)
                  ggCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACg                                                                                                  >  2:651757/1‑139 (MQ=255)
                   gcgcGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGt                                                                                                 <  2:572314/139‑1 (MQ=255)
                    cgcgGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATCGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTc                                                                                                >  1:108891/1‑139 (MQ=255)
                         agacagaTGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTAAGCAACAATCAGGGACAAAACAGAGTATTGCATGACAAATTTCACTGTCGCGATTGGCCCCACGTcatcca                                                                                           <  2:786412/133‑1 (MQ=255)
                                tGAGCAACAGTGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAg                                                                                    <  1:211802/139‑1 (MQ=255)
                                   gCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCaa                                                                                 >  2:293397/1‑139 (MQ=255)
                                    cAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAAt                                                                                <  2:670582/139‑1 (MQ=255)
                                    cAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAAt                                                                                <  2:449812/139‑1 (MQ=255)
                                          ggCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATaa                                                                          <  2:325567/139‑1 (MQ=255)
                                          cgCCGTTGACTAATGATGAAGGTAACCTTTTTTTAAAATAAAATGAGTTTTAGCAGCAATCAGGGCCAAAACAGATTATGGTATGAGAAATTTCACTTTCCCGACTTGCCCCACGTCATCCATCACGAGCAttttgtct                                                                          >  1:786412/2‑134 (MQ=255)
                                           gCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAAt                                                                         <  2:775121/139‑1 (MQ=255)
                                                         atgaAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTc                                                           >  1:45469/1‑139 (MQ=255)
                                                            aaGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCaaa                                                        >  2:627686/1‑139 (MQ=255)
                                                             aGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAAt                                                       <  2:100948/139‑1 (MQ=255)
                                                                tAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCa                                                    >  2:65604/1‑139 (MQ=255)
                                                                tAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCa                                                    <  1:627686/139‑1 (MQ=255)
                                                                        ttttAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACg                                            >  2:228062/1‑139 (MQ=255)
                                                                           taaaataaaaTGAGTTTTAGCAACAATCAGGGACAAAACAGATTATCGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTg                                         <  2:108891/139‑1 (MQ=255)
                                                                               ataaaatGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGcc                                     <  2:468560/139‑1 (MQ=255)
                                                                                              gCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTAACg                      <  1:532601/139‑1 (MQ=255)
                                                                                                 acaaTCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTAACGGtt                   >  2:433436/1‑139 (MQ=255)
                                                                                                             aaaaCAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAAt                                                                                >  2:581123/1‑66 (MQ=255)
                                                                                                             aaaaCAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAAt                                                                                <  1:741245/66‑1 (MQ=255)
                                                                                                             aaaaCAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAAt                                                                                <  1:581123/66‑1 (MQ=255)
                                                                                                             aaaaCAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAAt                                                                                >  2:741245/1‑66 (MQ=255)
                                                                                                             aaaaCAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCCCCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTAACGGTTCGCCATTCa          >  1:56098/1‑136 (MQ=255)
                                                                                                             aaaaCAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTAACGGTTCGCCATTCa          <  2:56098/136‑1 (MQ=255)
                                                                                                               aaCAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTAACGGTTCGCCATTCAGAAac     >  2:796116/1‑139 (MQ=255)
                                                                                                                  aGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTAACGGTTCGCCATTCAGAAaccac  <  1:651757/139‑1 (MQ=255)
                                                                                                                     |                                                                                                                                       
TGTACCAGCGTACGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTAACGGTTCGCCATTCAGAAACCAC  >  NC_000913/2644896‑2645148

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