Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 88,049 G→T 13.4% P50P (CCG→CCT ftsI → transpeptidase involved in septal peptidoglycan synthesis (penicillin‑binding protein 3)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00927388,0490GT13.4% 72.3 / 2.2 30P50P (CCG→CCTftsItranspeptidase involved in septal peptidoglycan synthesis (penicillin‑binding protein 3)
Reads supporting (aligned to +/- strand):  ref base G (14/12);  new base T (2/2);  total (16/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.77e-01

AATGCGCTCGGCGACCATAGCCGGGTGGAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATCCGTCACAAGAAAATATCGTAGTGCAAAAATAAGGATAAACGCGACGCATGAAAGCAGCGGCGAAAACGCAGAAACCAAAACGTCAGGAAGAACATGCCAACTTTATCAGTTGGCGTTTTGCGTTGTTATGCGGCTGTATTCTCCTGGCGCTGGCTTTTCTGCTCGGACGCGTAGCGTGGTTACAAGTTATCTCCCCGGATATGCTGGTGAAAGAGGGCGACATGCGTTCTCTTCGCGTTCAGCAAGTTTCCACCTCCCGCGGCATGATTACTGACCGTTCTGGTCGCCCGTTAGCGGTGAGCGTGCCGGTAAAAGCGATTTGGGCTGACCCGAAAGAAGTGCATGACGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCT  >  CP009273/87783‑88315
                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                          
aaTGCGCTCGGCGACCATAGCCGGGTGGAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATCCGTCACAAGAAAATATCGTAGTGCAAAAATAATGATAAACGCGACGCATGAAAGCAGCGGCGAAAACGCAGAAACCAAAACGTCAGGAAGAACATGCCAACTTTATCAGTTGGNGTTTTGCGTTGTTATGCGGCTGTATTCTCCTGGCGCTGGCTTTTCTGCTCGGACGCGTAGCGTGGTTACAAGTTATCTCCCCGGATATGCt                                                                                                                                                                                                                                                                    >  1:385999/1‑275 (MQ=255)
                                                         cagcaTGTTGATCCGTCACAAGAAAATATCGTAGTGCAAAAATAAGGATAAACGCGACGCATGAAAGCAGCGGCGAAAACGCAGAAACCAAAACGTCAGGAAGAACATGCCAACTTTATCAGTTGGCGTTTTGCGTTGTTATGCGGCTGTATTCTCCTGGCGCTGGCTTTTCTGCTCTGACGCGTAGCGTGGTTACAAGTTATCTCCCCGGATATGCTGGTGAAAGAGGGCGACATGCGTTCTCTTCGCGTTCAGCAAGTTTCCACATCACGCgg                                                                                                                                                                                                           >  1:275270/1‑275 (MQ=255)
                                                                  gATCCGTCACAAGAAAATATCGTAGTGCAAAAATAAGGATAAACGCGACGCATGAAAGCAGCGGCGAAAACGCAGAAACCAAAACGTCAGGAAGAACATGCCAACTTTATCAGTTGGCGTTTTGCGTTGTTATGCGGCTGTATTCTACTGGCGCTGGCTTTTCTGCTCGGACGCGTAGCGTGGTTACAAGTTATCTCCCCGGATATGCTGGTGAAAGAgg                                                                                                                                                                                                                                                         >  1:429002/1‑220 (MQ=255)
                                                                      cGTCACAAGAAAATATCGTAGTGCAAAAATAAGGATAAACGCGACGCATGAAAGCAGCGGCGAAAACGCAGAAACCAAAACGTCAGGAAGAACATGCCAACTTTATCAGTTGGCGTTTTGCGTTGTTATGCGGCTGTATTCTCCTGGCGCTGGCTTTTCTGCTCGGACGCGTAGCGTGGTTACAAGTTATCTCCCCGGATATGCTGGTGAAAGAGGGCGa                                                                                                                                                                                                                                                     <  2:428989/220‑1 (MQ=255)
                                                                                         aGTGAAAAAATAAGGATAAACGCGACGCATGAAAGCAGCGGCGAAAACGCAGAAACCAAAACGTCAGGAAGAACATGCCAACTTTATCAGTTGGCGTCTTGCGTTGTTATGCGGCTGTATTCTCCTGGCGCTGGCTTTTCTGCTCGGACGCGTAGCGTGGTTACAAGTCATCTCCCCGGATATGCTGGTGAAAGAGGGCGACATGCGTTCTCTTCGCGTTCAGCAAGTTTCCACCTCCCGCGGCATGATTACTGACCGTTCTGGTCGCCCGTTAg                                                                                                                                                                           <  1:259977/275‑1 (MQ=255)
                                                                                                                               cggcgAAAACGCAGAAACCAAAACGTCAGGAAGAACATGCCAACTTTATCAGTTGGCGTTTTGCGTTGTTATGCGGCTGTATTCTCCTGGCGCTGGCTTTTCTGCTCGGACGCGTAGCGTGGTTACAAGTTATCTCCCCGGATATGCTGGTg                                                                                                                                                                                                                                                                >  2:170828/1‑152 (MQ=255)
                                                                                                                                   gAAAACGCAGAAACCAAATCGTCAGGAAGAACATGCCAACTTTATCATTTGGCGTTTTGCGTTGTTATGCGGCTGTATTCTCCTGGCGCTGNCTTTTCTGCTCGGACGCGTAGCGTGGTTACAAGTTATCTCCCCGGATATGCTGGTGAAAGAGGGCGACATGCGTTCTCTTCGCGTTCAGCAAGTTTCCACCTCCCGCGGCATGATTACTGACCGTTCTGGTCGCCCGTTAGCGGTGAGCGTGCCTGTAAAAGCGATTTGGGCTGACCCGAaag                                                                                                                                 <  1:367824/275‑1 (MQ=255)
                                                                                                                                   gAAAACGCAGAAACCAAAACGTCAGGAAGAACATGCCAACTTTATCAGTTGGCGTTTTGCGTTGTTATGCGGCTGTATTCTCCTGGCGCTGGCTTTTCTGCTCGGACGCGTAGCGTGGTTACAAGTTATCTCCCCGGATATGCTGGTGAAag                                                                                                                                                                                                                                                            <  1:170834/152‑1 (MQ=255)
                                                                                                                                                                aaCATGCCAACTTTATCAGCTGGCGTTTTGCGTTGTTATGCGGCTGTATTCTCCTGGCGCTGGCTTTTCTGCTCGGACGCGTAGCGTGGTTACAAGTTATCTCCCCGGATATGCTGGTGAAAGAGGGCGACATGCGTTCTCTTCGCGTTCAGCAAGTTTCCACCTCCCGCGGCATGATTACTGACCGTTCTGGTCGCCCGTTAGCGGTGAGCGTGCCGGTAAAAGCGATTTGGGCTGACCCGAAAGAAGTGCATGACGCTGGCGGTATCAGCGTc                                                                                                    >  1:401132/1‑275 (MQ=255)
                                                                                                                                                                                    tGGCGTTTTGCGTTGTTATGCGGCTGTATTCTCCTGGCGCTGGCTTTTCTGCTCGGACGCGTAGCGTGGTTACAAGTTATCTCCCCGGATATGCTGGTGAAAGAGGGCGACATGCGTTCTCTTCGCGTTCAGCAAGTTTCCACCTCCCGCGGCATGATTACTGACCGTTCTGGTCGCTCGTTAGCGGTGAGCGTGCCGGTACAAGCGATTTGGGCTGACCCGAAAGAAGTGCATGACGCTGGCGGTATCAGCGTCGGAGACCGCTGGaactgtc                                                                                 >  1:262458/1‑269 (MQ=255)
                                                                                                                                                                                     ggCGTTTTGCGTTGTTATGCGGCTGTATTCTCCTGGCGCTGGCTTTTCTGCTCGGACGCGTAGCGTGGTTACAAGTTATCTCCCCGGATATGCTGGTGAAAGAGGGCGACATGCGTTCTCTTCGCGTTCAGCAAGTTTACCCCTCCCGCGGCCTGAt                                                                                                                                                                                                     >  2:372289/1‑157 (MQ=255)
                                                                                                                                                                                         ttttGCGTTGTTATGCGGCTGTATGCTCCTGGCGCTGGCTTTTCTGCTCGGACGCGTAGCGTGGTTACAAGTTATCTCCCCGGATATGCTGGTGAAAGAGGGCGACATGCGTTCTCTTCGCGTTCAGCTAGTTTCCACCTCCCGCGGCATGATTACt                                                                                                                                                                                                 <  1:372302/157‑1 (MQ=255)
                                                                                                                                                                                                        cGGCTGTATTCTCCTGGCGCTGGCTTTTCTGCTCGGACGCGTAGCGTGGTTACAAGTTATCTCCCCgg                                                                                                                                                                                                                                                                           >  2:267674/1‑68 (MQ=255)
                                                                                                                                                                                                            tGTATTCTCCTGGCGCTGGCTTTTCTGCTCGGACGCGTAGCGTGGTTACAAGTTATCTCCCCGGatat                                                                                                                                                                                                                                                                       <  1:267682/68‑1 (MQ=255)
                                                                                                                                                                                                             ggaTTCTCCTGGCGCTGGCTTTTCTGCTCGGACGCGTAGCGTGGTTACAAGTTATCTCCCCGGATATg                                                                                                                                                                                                                                                                      >  2:355380/3‑68 (MQ=255)
                                                                                                                                                                                                             gTATTCTCCTGGCGCTGGCTTTTCTGCTCGGACGCGTAGCGTGGTTACAAGTTATCTCCCCTGATATGCTGGTGAAAGAGGGCGACATGCGTTCTCTTCGCGTTCAGCAAGTTTCCACCTc                                                                                                                                                                                                                 >  1:284981/1‑121 (MQ=255)
                                                                                                                                                                                                                 tctcCTGGCGCTGGCTTTTCTGCTCGGACGCGTAGCGTGGTTACAAGTTATCTCCCCTGATATGCTGGTGAAAGAGGGCGACATGCGTTCTCTTCGCGTTCAGCAAGTTTCCACCTCCcgc                                                                                                                                                                                                             <  2:284973/121‑1 (MQ=255)
                                                                                                                                                                                                                 tctcCTGGCGCTGGCTTTTCTGCTCGGACGCGTAGCGTGGTTACAAGTTATCTCCCCGGATCTGCTgg                                                                                                                                                                                                                                                                  <  1:355392/68‑1 (MQ=255)
                                                                                                                                                                                                                 gctcCTGGCTCTTGCTTTTCTGCTCGGACGCGTAGCGTGGTTACAAGTTATCTCCCCGGATATGCTGGTGAAAGAGGGCGACATGCGTTCTCTTCGCGTTCAGCAAGTTTCCACCTCCCGCGGCATGATTACTGACCGTTCTGGTCGCCCGTTAGCGGTGAGCGTGCCGGTAAAAGCGATTTGGGCTGACCCGAAAGAAGTGCATGACGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATc                                                   <  1:39645/274‑1 (MQ=255)
                                                                                                                                                                                                                                        tCGGACGCGTAGCGTGGTTACAAGTTATCTCCCCGGATATGCTGGTGAAAGAGGGCGACATGCGTTCTCTTCGCGTTCAg                                                                                                                                                                                                                               >  1:246105/1‑80 (MQ=255)
                                                                                                                                                                                                                                            aCGCGTAGCGTGGTTACAAGTTATCTCCCCGGATATGCTGGTGAAAGAGGGCGACATGCGTTCTCTTCGCGTTCAGCAAg                                                                                                                                                                                                                           <  2:246098/80‑1 (MQ=255)
                                                                                                                                                                                                                                               cgTAGCGTGGTTACAAGTTATCTCCCCGGATATGCTGGTGAAAGAGGGCGACATGCGTTCTCTTCGCGTTCAGCAAGTTTCCACCTCCCGCGGCATGATTACTGACCGTTCTGGTCGCCCGtt                                                                                                                                                                             >  2:285181/1‑123 (MQ=255)
                                                                                                                                                                                                                                                   gCGTGGTTACAAGTTATCTCCCCGGATATGCTGGTGAAAGAGGGCGACATGCGTTCTCTTCGCGTTCAGCAAGTTTCCACCTCCCGCGGCATGATTACTGACCGTTCTGGTCGCCCGTTAGCg                                                                                                                                                                         <  1:285189/123‑1 (MQ=255)
                                                                                                                                                                                                                                                        gTTACAAGTTATCTCCCCTGATATGCTGGTGAAAGAGGGCGACATGCGTTCTCTTCGCGTTCAGCAAGTTTCCACCTCCcgc                                                                                                                                                                                                             >  2:235696/1‑82 (MQ=255)
                                                                                                                                                                                                                                                        gTTACAAGTTATCTCCCCGGATATGCTGGTGAAAGAGGGCGACATGCGTTCTCTTCGCGTTCAGCAAGTTTCCACCTCCCGCGGCATGATTACTGACCGTTCTGGTCGCCCGTTAGCGGTGAGCGTGCCGGTAAAAGCGATTTGGGCTGACCCGAAAGAAGTGCATGACGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGgcgc             >  1:309363/1‑274 (MQ=255)
                                                                                                                                                                                                                                                            cAAGTTATCTCCCCTGATATGCTGGTGAAAGAGGGCGACATGCGTTCTCTTCGCGTTCAGCAAGTTTCCACCTCCCGCGGCa                                                                                                                                                                                                         <  1:235703/82‑1 (MQ=255)
                                                                                                                                                                                                                                                                  aTCTCCCCGGATATGCTGGTGAAAGAGGGCGACATGCGTTCTCTTCGCGTTCAGCAAGTTTCCACCTCCcgc                                                                                                                                                                                                             >  2:31238/1‑72 (MQ=255)
                                                                                                                                                                                                                                                                  aTCTCCCCGGATATGCTGGTGAAAGAGGGCGACATGCGTTCTCTTCGCGTTCAGCAAGTTTCCACCTCCCGCGGCATGATTACTGACCGTTCTGGTCGCCCGt                                                                                                                                                                              >  1:370072/1‑103 (MQ=255)
                                                                                                                                                                                                                                                                  aTCTCCCCGGATATGCTGGTGAAAGAGAGCGACATGCGTTCTCTTCGCGTTCAGCAAGTTTCCACCTCCCGCGGCATGATTACTGACCGTTCTGGTCGCCCGTTAGCGGTGAGCGTGCCGGTAAAAGCGATTTGGGCTGACCCGAAAGAAGTGCATGACGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGGCTAAAGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGAGCTTTATTTATCt  >  1:292275/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                      ccccGGATATGCTGGTGAAAGAGGGCGACATGCGTTCTCTTCGCGTTCAGCAAGTTTCCACCTCCCGCGGCa                                                                                                                                                                                                         <  1:31239/72‑1 (MQ=255)
                                                                                                                                                                                                                                                                      ccccGGATATGCTGGTGAAAGAGGGCGACATGCGTTCTCTTCGCGTTCAGCAAGTTTCCACCTCCCGCGGCATGATTACTGACCGTTCTGGTCGCCCGTTAGc                                                                                                                                                                          <  2:370059/103‑1 (MQ=255)
                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                          
AATGCGCTCGGCGACCATAGCCGGGTGGAAAGGATCGCCACGGAAAAGCTGCAAATGCAGCATGTTGATCCGTCACAAGAAAATATCGTAGTGCAAAAATAAGGATAAACGCGACGCATGAAAGCAGCGGCGAAAACGCAGAAACCAAAACGTCAGGAAGAACATGCCAACTTTATCAGTTGGCGTTTTGCGTTGTTATGCGGCTGTATTCTCCTGGCGCTGGCTTTTCTGCTCGGACGCGTAGCGTGGTTACAAGTTATCTCCCCGGATATGCTGGTGAAAGAGGGCGACATGCGTTCTCTTCGCGTTCAGCAAGTTTCCACCTCCCGCGGCATGATTACTGACCGTTCTGGTCGCCCGTTAGCGGTGAGCGTGCCGGTAAAAGCGATTTGGGCTGACCCGAAAGAAGTGCATGACGCTGGCGGTATCAGCGTCGGTGACCGCTGGAAGGCGCTGGCTAACGCGCTCAATATTCCGCTGGATCAGCTTTCAGCCCGCATTAACGCCAACCCGAAAGGGCGCTTTATTTATCT  >  CP009273/87783‑88315

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