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.6% 123.3 / ‑7.2 41*771* (TAG→TAApbpCpenicillin‑insensitive murein repair transglycosylase, inactive transpeptidase domain protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (22/18);  total (23/18)
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

TCTGTTGTTGTCTGTACCAGCGTACGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTAACGGTTCGCCATTCAGAAACCAC  >  NC_000913/2644884‑2645148
                                                                                                                                 |                                                                                                                                       
tCTGTTGTTGTCTGTACCAGCGTACGAGATGGCGCGGATTATAATGAGCAACAGGGCCGATGACTGTTGGTGAAGGTAACATTTTGTTACAATAAAATGAGATTTAGCAACGATCAAGGACAAAACAGACtctgtctca                                                                                                                                >  1:130003/1‑131 (MQ=255)
  tgttgttgtCTGTACCAGCGTACGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGAc                                                                                                                              <  2:1850292/139‑1 (MQ=255)
            tGTACCAGCGTACGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACt                                                                                                                    >  2:458434/1‑139 (MQ=255)
                    cGTACGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGAt                                                                                                            >  2:1913715/1‑139 (MQ=255)
                    cGTACGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGAt                                                                                                            >  1:720288/1‑139 (MQ=255)
                      tACGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGAttt                                                                                                          >  2:498542/1‑139 (MQ=255)
                        cGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCa                                                                                                                                  <  2:955505/113‑1 (MQ=255)
                        cGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCa                                                                                                                                  >  1:955505/1‑113 (MQ=255)
                        cGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGc                                                                                                        >  1:1101432/1‑139 (MQ=255)
                          agaTGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGccc                                                                                                      <  2:1101432/139‑1 (MQ=255)
                              ggCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACg                                                                                                  <  1:1913715/139‑1 (MQ=255)
                                         taATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCAc                                                                                       <  2:720288/139‑1 (MQ=255)
                                               gCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCaa                                                                                 >  1:236342/1‑139 (MQ=255)
                                                        ccGTTGACTATTGATGAAGGTAACATTTTTTTAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCa                                                                       >  1:908565/1‑139 (MQ=255)
                                                                     atgaAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTc                                                           <  1:458434/139‑1 (MQ=255)
                                                                          ggTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATg                                                      <  1:498542/139‑1 (MQ=255)
                                                                               cATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAggg                                                 <  2:236342/139‑1 (MQ=255)
                                                                                    ttttAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGcc                                                                                                       <  1:1397136/80‑1 (MQ=255)
                                                                                    ttttAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGcc                                                                                                       >  2:1397136/1‑80 (MQ=255)
                                                                                    ttttAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACg                                            >  2:1957808/1‑139 (MQ=255)
                                                                                     tttAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATg                                                      <  1:1467049/128‑1 (MQ=255)
                                                                                     tttAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATg                                                      >  2:1467049/1‑128 (MQ=255)
                                                                                         aaataaaatGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTgcg                                       <  1:318293/139‑1 (MQ=255)
                                                                                               aatGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCcgcg                                 <  2:908565/139‑1 (MQ=255)
                                                                                                      tttAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGtt                          >  2:118689/1‑139 (MQ=255)
                                                                                                        tAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGt                                <  1:622277/131‑1 (MQ=255)
                                                                                                        tAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGt                                >  2:622277/1‑131 (MQ=255)
                                                                                                         aGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTAAc                       >  1:1917665/1‑139 (MQ=255)
                                                                                                          gCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTAACg                      <  1:118689/139‑1 (MQ=255)
                                                                                                               aaTCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTAACGGTTCg                 >  2:560228/1‑139 (MQ=255)
                                                                                                               aaTCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTAACGGTTCg                 >  2:33341/1‑139 (MQ=255)
                                                                                                                aTCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTAACGGTTCGc                >  1:989623/1‑139 (MQ=255)
                                                                                                                   aGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTAACGGTTCGCCAt             >  1:1017479/1‑139 (MQ=255)
                                                                                                                         aaaaCAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGt                                >  2:962180/1‑114 (MQ=255)
                                                                                                                         aaaaCAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGt                                <  1:962180/114‑1 (MQ=255)
                                                                                                                           aaCAGATTATTGCATGACAAATTTCACTGTCGCGATTTGccc                                                                                                      <  2:1088012/42‑1 (MQ=255)
                                                                                                                           aaCAGATTATTGCATGACAAATTTCACTGTCGCGATTTGccc                                                                                                      >  1:1088012/1‑42 (MQ=255)
                                                                                                                           aaCAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATaa                                                                          >  2:1147640/1‑70 (MQ=255)
                                                                                                                           aaCAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATaa                                                                          <  1:1147640/70‑1 (MQ=255)
                                                                                                                           aaCAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTAACGGTTCGCCATTCAGAAac     >  2:1754048/1‑139 (MQ=255)
                                                                                                                              aGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTAACGGTTCGCCATTCAGAAaccac  <  1:33341/139‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                       
TCTGTTGTTGTCTGTACCAGCGTACGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTAACGGTTCGCCATTCAGAAACCAC  >  NC_000913/2644884‑2645148

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