Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 61,271 T→C 10.2% P147P (CCT→CCC) ‡ PP_0052 → beta‑lactamase domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294761,2710TC10.2% 100.7 / 3.3 41P147P (CCT→CCC) ‡PP_0052beta‑lactamase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (18/17);  new base C (2/2);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.10e-01

ACATTACCCAGGTGCAGAAAGTGTTTGGCGCCTTGTTCAACGCAGAGCCCGGCTTCGCCCGTGATGGCAGCCAGTTCGATGTGCTGCTCGAAGACGAGGAAGGTTTCCGGATTGGCAACCTGCAGGCCCGTGCGTTGCACACGCCTGGGCACACACCCGCGTGCATGAGTTTCATGATCGAAGATGCCGGCGAGATCGCGGTGTTCGTTGGAGACACCTTGTTCATGCCCGACTATGGCACCGCCCGCTGCGACTTCCCCGGGGCCGATGCCCGCACCTTGTACCGCTCAATC  >  NC_002947/61126‑61418
                                                                                                                                                  |                                                                                                                                                   
acaTTACCCAGGTGCAGAAAGTGTTTGGCGCCTTGTTCAACGCAGAGCCCGGCTTCGCCCGTGATGGCAGCCAGTTCGATGTGCTGCTCGAAGACGAGGAAGGTTTCCGGATTGGCAACCTGCAGGCCCGTGCGTTGCACACGCCTggg                                                                                                                                                  <  5:48811/149‑1 (MQ=255)
  aTTACCCAGGTGCAGAAAGTGTTTGGCGCCTTGTTCAACGCAGAGCCCGGCTTCGCCCGTGATGGCAGCCAGTTCGATGTGCTGCTCGAAGACGAGGAAGGTTTCCGGATTGGCAACCTGCAGGCCCGTGCGTTGCACACGCCTGGGca                                                                                                                                                >  6:96709/1‑149 (MQ=255)
   ggcatcggccccggggaagtcgcagcgggcggtgccatagtcgggcatgaacaaggtgtctccaacgaacaccgcgatctcgccggcatcttcgatcatgaaactcatgcacgcgggtgtgtGCCCGTGCGTTGCACACGCCTGGGcac                                                                                                                                               >  7:137451‑M1/123‑149 (MQ=255)
              cAGAAAGTGTTTGGCGCCTTGTTCAACGCAGAGCCCGGCTTCGCCCGTGATGGCAGCCAGTTCGATGTGCTGCTCGAAGACGAGGAAGGTTTCCGGATTGGCAACCTGCAGGCCCGTGCGTTGCACACGCCTGGGCACACACCCGCGt                                                                                                                                     <  5:88367/148‑1 (MQ=255)
              cAGAAAGTGTTTGGCGCCTTGTTCAACGCAGAGCCCGGCTTCGCCCGTGATGGCAGCCAGTTCGATGTGCTGCTCGAAGACGAGGAAGGTTTCCGGATTGGCAACCTGCAGGCCCGTGCGTTGCACACGCCTGGGCACACACCCGCGTg                                                                                                                                    >  5:51510/1‑149 (MQ=255)
              cAGAAAGTGTTTGGCGCCTTGTTCAACGCAGAGCCCGGCTTCGCCCGTGATGGCAGCCAGTTCGATGTGCTGCTCGAAGACGAGGAAGGTTTCCGGATTGGCAACCTGCAGGCCCGTGCGTTGCACACGCCTGGGCACACACCCGCGTg                                                                                                                                    >  5:46313/1‑149 (MQ=255)
                  aaGTGGTTGGCGCCTTGTTCAACGCAGAGCACGGCTTCGCTCTTGATGGCAGCCAGTTCGATGTGCTTCTCGAAGACGAGGAAGGTTTCATGATTTGGAAACTTCAGTCCCGTGCTTTGCAGTCGCCTTGGCACACACGCGGGTGCAt                                                                                                                                 <  4:191496/148‑1 (MQ=255)
                              ccTTGTTCAACGCAGAGCCCGGCTTCGCCCGTGATGGCAGCCAGTTCGATGTGCTGCTCGAAGACGAGGAAGGTTTCCGGATTGGCAACCTGCAGGCCCGTGCGTTGCACACGCCTGGGCACACACCCGCGTGCATGAGTTTCATGATc                                                                                                                    >  6:71645/1‑149 (MQ=255)
                                  gTTCAACGCAGAGCCCGGCTTCGCCCGTGATGGCAGCCAGTTCGATGTGCTGCTCGAAGACGAGGAAGGTTTCCGGATTGGCAACCTGCAGGCCCGTGCGTTGCACACGCCTGGGCACACACCCGCGTGCATGAGTTTCATGATCGaa                                                                                                                 <  6:46313/148‑1 (MQ=255)
                                  gTTCAACGCAGAGCCCGGCTTCGCCCGTGATGGCAGCCAGTTCGATGTGCTGCTCGAAGACGAGGAAGGTTTCCGGATTGGCAACCTGCAGGCCCGTGCGTTGCACACGCCTGGGCACACACCCGCGTGCATGAGTTTCATGATCGaa                                                                                                                 <  6:51510/148‑1 (MQ=255)
                                  gTTCAACGCAGAGCCCGGCTTCGCCCGTGATGGCAGCCAGTTCGATGTGCTGCTCGAAGACGAGGAAGGTTTCCGGATTGGCAACCTGCAGGCCCGTGCGTTGCACACGCCTGGGCACACACCCGCGTGCATGAGTTTCATGATCGaa                                                                                                                 <  5:96709/148‑1 (MQ=255)
                                   ttCAACGCAGAGCCCGGCTTCGCCCGTGATGGCAGCCAGTTCGATGTGCTGCTCGAAGACGAGGAAGGTTTCCGGATTGGCAACCTGCAGGCCCGTGCGTTGCACACGCCTGGGCACACACCCGCGTGCATGAGTTTCATGATc                                                                                                                    <  4:71174/144‑1 (MQ=255)
                                   ttCAACGCAGAGCCCGGCTTCGCCCGTGATGGCAGCCAGTTCGATGTGCTGCTCGAAGACGAGGAAGGTTTCCGGATTGGCAACCTGCAGGCCCGTGCGTTGCACACGCCTGGGCACACACCCGCGTGCATGAGTTTCATGATCGagat                                                                                                               >  3:71174/1‑146 (MQ=255)
                                     cAACGCAGAGCCCGGCTTCGCCCGTGATGGCAGCCAGTTCGATGTGCTGCTCGAAGACGAGGAAGGTTTCCGGATTGGCAACCTGCAGGCCCGTGCGTTGCACACGCCTGGGCACACACCCGCGTGCATGAGTTTCATGATCGAAGAt                                                                                                              <  3:171135/148‑1 (MQ=255)
                                                   gCTTCGCCCGTGATGGCAGCCAGTTCGATGTGCTGCTCGAAGACGAGGAAGGTTTCCGGATTGGCAACCTGCAGGCCCGTGCGTTGCACACGCCTGGGCACACACCCGCGTGCATGAGTTTCATGATCGAAGATGCCGGCGAGATcgcg                                                                                               <  2:310839/149‑1 (MQ=255)
                                                       cGCCCGTGATGGCAGCCAGTTCGATGTGCTGCTCGAAGACGAGGAAGGTTTCCGGATTGGCAACCTGCAGGCCCGTGCGTTGCACACGCCTGGGCACACACCCGCGTGCATGAGTTTCATGATCGAAGATGCCGGCGAGATCGCGgt                                                                                             <  4:275658/147‑1 (MQ=255)
                                                          ccGTGATGGCAGCCAGTTCGATGTGCTGCTCGAAGACGAGGAAGGTTTCCGGATTGGCAACCTGCAGGCCCGTGCGTTGCACACGCCTGGGCACACACCCGCGTGCATGAGTTTCATGATCGAAGATGCCGGCGAGATCGCGGTGTTCg                                                                                        >  5:226681/1‑149 (MQ=255)
                                                            gTGATGGCAGCCAGTTCGATGTGCTGCTCGAAGACGAGGAAGGTTTCCGGATTGGCAACCTGCAGGCCCGTGCGTTGCACACGCCTGGGCACACACCCGCGTGCATGAGTTTCATGATCGAAGATGCCGGCGAGATc                                                                                                  <  2:41507/137‑1 (MQ=255)
                                                            gTGATGGCAGCCAGTTCGATGTGCTGCTCGAAGACGAGGAAGGTTTCCGGATTGGCAACCTGCAGGCCCGTGCGTTGCACACGCCTGGGCACACACCCGCGTGCATGAGTTTCATGATCGAAGATGCCGGCGAGATc                                                                                                  >  1:41507/1‑137 (MQ=255)
                                                             tGATGGCAGCCAGTTCGATGTGCTGCTCGAAGACGAGGAAGGTTTCCGGATTGGCAACCTGCAGGCCCGTGCGTTGCACACGCCTGGGCACTCACCCGCGTGCATGAGTTTCATGATCGAAGATGCCGGCGAGATCGCGGTGTTCGTTg                                                                                     <  4:321431/149‑1 (MQ=255)
                                                               aTGGCAGCCAGTTCGATGTGCTGCTCGAAGACGAGGAAGGTTTCCGGATTGGCAACCTGCAGGCCCGTGCGTTGCACACGCCTGGGCACACACCCGCGTGCATGAGTTTCATGATCGAAGATGCCGGCGAGATCGCGGTGTTCGTTGga                                                                                   <  2:131860/149‑1 (MQ=255)
                                                                tGGCAGCCAGTTCGATGTGCTGCTCGAAGACGAGGAAGGTTTCCGGATTGGCAACCTGCAGGCCCGTGCGTTGCACACGCCGGGCACACACCCGCGTGCATGAGTTTCATGATCGAAGATGCCGGCGAGATAGCGGTGTTCGTTGgaga                                                                                 >  3:35350/1‑149 (MQ=255)
                                                                       cAGTTCGATGTGCTGCTCGAAGACGAGGAAGGTTTCCGGATTGGCAACCTGCAGGCCCGTGCGTTGCACACGCCTGGGCACACACCCGCGTGCATGAGTTTCATGATCGAAGATGCCGGCGAGATCGCGGTGTTCGTTGGAGACACCtt                                                                           >  8:14773/1‑149 (MQ=255)
                                                                            cGATGTGCTGCTCGAAGACGAGGAAGGTTTCCGGATTGGCAACCTGCAGGCCCGTGCGTTGCACACGCCTGGGCACACACCCGCGTGCATGAGTTTCATGATCGAAGATGCCGGCGAGATCGCGGTGTTCGTTGGAGACACCTTGTTc                                                                       <  6:160457/148‑1 (MQ=255)
                                                                                         gAAGACGAGGAAGGTTTCCGGATTGGCAACCTGCAGGCCCGTGCGTTGCACACGCCTGGGCACACACCCGCGTGCATGAGTTTCATGATCGAAGATGCCGGCGAGATCGCGGTGTTCGTTGGAGACACCTTGTTCATGCCCGACTAt                                                           <  5:71645/147‑1 (MQ=255)
                                                                                         gAAGACGAGGAAGGTTTCCGGATTGGCAACCTGCAGGCCCGTGCGTTGCACACGCCTGGGCACACACCCGCGTGCATGAGTTTCATGATCGAAGATGCCGGCGAGATCGCGGTGTTCGTTGGAGACACCTTGTTCATGCCCGACTATgg                                                         >  7:311041/1‑149 (MQ=255)
                                                                                                  gAAGGTTTCCGGATTGGCAACCTGCAGGCCCGTGCGTTGCACACGCCTGGGCACACACCCGCGTGCATGAGTTTCATGATCGAAGATGCCGGCGAGATCGCGGTGTTCGTTGGAGACACCTTGTTCATGCCCGACTATGGCAccgcccg                                                >  4:198833/1‑149 (MQ=255)
                                                                                                       tttCCGGATTGGCAACCTGCAGGCCCGTGCGTTGCACACGCCTGGGCACACACCCGCGTGCATGAGTTTCATGATCGAAGATGCCGGCGAGATCGCGGTGTTCGTTGGAGACACCTTGTTCATGCCCGACTATGGCACCGCCCGCTGCg                                           >  4:32935/1‑149 (MQ=255)
                                                                                                                  gCAACCTGCAGGCCCGTGCGTTGCACACGCCTGGGCACACACCCGCGTGCATGAGTTTCATGATCGAAGATGCCGGCGAGATCGCGGTGTTCGTTGGAGACACCTTGTTCATGCCCGACTATGGCACCGCCCGCTGCGACTTCCCCggg                                >  3:96665/1‑149 (MQ=255)
                                                                                                                  gCAACCTGCAGGCCCGTGCGTTGCACACGCCTGGGCACACACCCGCGTGCATGAGTTTCATGATCGAAGATGCCGGCGAGATCGCGGTGTTCGTTGGAGACACCTTGTTCATGCCCGACTATGGCACCGCCCGCTGCGACTTCCCCggg                                >  3:27011/1‑149 (MQ=255)
                                                                                                                   gcgggtgtgtGCCCAGGCGTGTGCAACGCACGGGCACACACCCGCGTGCATGAGTTTCATGATCGAAGATGCCGGCGAGATcgcg                                                                                               <  5:185936‑M1/75‑1 (MQ=255)
                                                                                                                   gcgggtgtgtGCCCAGGCGTGTGCAACGCACGGGCACACACCCGCGTGCATGAGTTTCATGATCGAAGATGCCGGCGAGATcgcg                                                                                               >  6:185936‑M1/11‑85 (MQ=255)
                                                                                                                   gcgggtgtgtGCCCAGGCGTGTGCAACGCACGGGCACACACCCGCGTGCATGAGTTTCATGATCGAAGATGCCGGCGAGATcgc                                                                                                <  7:210194‑M1/74‑1 (MQ=255)
                                                                                                                   gcgggtgtgtGCCCAGGCGTGTGCAACGCACGGGCACACACCCGCGTGCATGAGTTTCATGATCGAAGATGCCGGCGAGATcgc                                                                                                >  8:210194‑M1/11‑84 (MQ=255)
                                                                                                                   cAACCTGCAGGCCCGTGCGTTGCACACGCCTGGGCACACACCCGCGTGCATGAGTTTCATGATCGAAGATGCCGGCGAGATCGCGGTGTTCGTTGGAGACACCTTGTTCATGCCCGACTATGGCACCGCCCGCTGCGACTTCCCCgggg                               >  1:197529/1‑149 (MQ=255)
                                                                                                                                    cGTTGCACACGCCTGGGCACACACCCGCGTGCATGAGTTTCATGATCGAAGATGCCGGCGAGATCGCGGTGTTCGTTGGAGACACCTTGTTCATGCCCGACTATg                                                          >  8:137701/1‑105 (MQ=255)
                                                                                                                                      ttGCACACGCCGGGCACACACCCGCGTGCATGAGTTTCATGATCGAAGATGCCGGCGAGATAGCGGTGTTCGTTGGAGACACCTTGTTCATGCCCGACTATGGCACCGCCCGCTGCGACTTCCCCGGGGCCGATGCCCGCACCTTGTAc           <  4:35350/149‑1 (MQ=255)
                                                                                                                                           acacGCCTGGGCACACACCCGCGTGCATGAGTTTCATGATCGAAGATGCCGGCGAGATCGCGGTGTTCGTTGGAGACACCTTGTTCATGCCCGACTATGGCACCGCCCGCTGCGACTTCCCCGGGGCCGATGCCCGCACCTTGTACCGc        >  4:239371/1‑149 (MQ=255)
                                                                                                                                             acGCCTGGGCACACACCCGCGTGCATGAGTTTCATGATCGAAGATGCCGGCGAGATCGCGGTGTTGGTTGGAGACACCTTGTTCATGCCCGACTATGGCACCGCCCGCTGCGACTTCCCCGGGGCCGATGCCCGCACCTTGTACCGCt       >  2:111367/1‑148 (MQ=255)
                                                                                                                                                ccTGGGCACACACCCGCGTGCATGAGTTTCATGATCGAAGATGCCGGCGAGATCGCGGTGTTCGTTGGAGACACCTTGTTCATGCCCGACTATGGCACCGCCCGCTGCGACTTCCCCGGGGCCGATGCCCGCACCTTGTACCGCTCAAt   <  2:197529/149‑1 (MQ=255)
                                                                                                                                                 cTGGGCACACACCCGCGTGCATGAGTTTCATGTTCGAAGATGCCGGCGAGATCGCGGTGTTCGTTGGAGACACCTTGTTCATGCCCGACTATGGCACCGCCCGCTGCGACTTCCCCGGGGCCGATGCCCGCACCTTGTACCGCTCAATc  <  6:144971/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ACATTACCCAGGTGCAGAAAGTGTTTGGCGCCTTGTTCAACGCAGAGCCCGGCTTCGCCCGTGATGGCAGCCAGTTCGATGTGCTGCTCGAAGACGAGGAAGGTTTCCGGATTGGCAACCTGCAGGCCCGTGCGTTGCACACGCCTGGGCACACACCCGCGTGCATGAGTTTCATGATCGAAGATGCCGGCGAGATCGCGGTGTTCGTTGGAGACACCTTGTTCATGCCCGACTATGGCACCGCCCGCTGCGACTTCCCCGGGGCCGATGCCCGCACCTTGTACCGCTCAATC  >  NC_002947/61126‑61418

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