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,0130CT95.8% 70.9 / NA 24*771* (TAG→TAApbpCpenicillin‑insensitive murein repair transglycosylase, inactive transpeptidase domain protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (11/12);  minor base G (0/1);  total (11/13)
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

TTGTTCTGTTGTTGTCTGTACCAGCGTACGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTAACG  >  NC_000913/2644880‑2645128
                                                                                                                                     |                                                                                                                   
ttGTTCTGTTGTTGTCTGTACCAGCGTACGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTg                                                                                                                <  2:723511/139‑1 (MQ=255)
    tCTGTTGTTGTCTGTACCAGCGTACGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATg                                                                                                            >  2:742354/1‑139 (MQ=255)
    tCTGTTGTTGTCTGTACCAGCGTACGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATg                                                                                                            >  2:401121/1‑139 (MQ=255)
                  tACCAGCGTACGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGt                                                                                              >  2:653192/1‑139 (MQ=255)
                          tACGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATCAATTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGAttt                                                                                      >  1:189068/1‑139 (MQ=255)
                            cGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGc                                                                                    <  1:653192/139‑1 (MQ=255)
                                    cgcgGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTc                                                                            >  2:648407/1‑139 (MQ=255)
                                             taATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGGCGCGATTTGCCCCACGTCATCCATCAc                                                                   <  1:401121/139‑1 (MQ=255)
                                                  aGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTg                                                                                     >  2:657998/1‑116 (MQ=255)
                                                  aGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTg                                                                                     <  1:657998/116‑1 (MQ=255)
                                                   gCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCaa                                                             >  1:272115/1‑139 (MQ=255)
                                                     aaCAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAg                                                                <  1:340885/134‑1 (MQ=255)
                                                     aaCAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAg                                                                >  2:340885/1‑134 (MQ=255)
                                                           gCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGAGTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAAt                                                     <  2:538878/139‑1 (MQ=255)
                                                              gTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCAc                                                  <  1:347137/139‑1 (MQ=255)
                                                                         atgaAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTc                                       <  1:648407/139‑1 (MQ=255)
                                                                            aaGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTCTCGCGATTTGCCCCAc                                                                               <  2:265332/96‑1 (MQ=255)
                                                                            aaGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTCTCGCGATTTGCCCCAc                                                                               >  1:265332/1‑96 (MQ=255)
                                                                                 aaCATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAg                               >  2:494561/1‑139 (MQ=255)
                                                                                        ttttAAAATTAAATGAGTTTTAGCAACAATCAGGCACAAATCAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACg                        <  1:494561/139‑1 (MQ=255)
                                                                                                            tAGCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTaa    <  2:272115/139‑1 (MQ=255)
                                                                                                              gCAACAATCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTAACg  <  1:368146/139‑1 (MQ=255)
                                                                                                                    aTCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGccc                                                                                  <  2:54534/53‑1 (MQ=255)
                                                                                                                    aTCAGGGACAAAACAGATTATTGCATGACAAATTTCACTGTCGCGATTTGccc                                                                                  >  1:54534/1‑53 (MQ=255)
                                                                                                                                     |                                                                                                                   
TTGTTCTGTTGTTGTCTGTACCAGCGTACGAGATGGCGCGGATTATAATGAGCAACAGGGCCGTTGACTATTGATGAAGGTAACATTTTTTTAAAATAAAATGAGTTTTAGCAACAATCAGGGACAAAACAGACTATTGCATGACAAATTTCACTGTCGCGATTTGCCCCACGTCATCCATCACCAGCAATTGATAATCACCTTTATCCGTCAAATGCAGGGTCACGTTGCGCCCGCGTTCAGTTAACG  >  NC_000913/2644880‑2645128

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