Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,250,457 C→T 6.9% A76A (GCG→GCA cidB ← holin CidA anti‑holin regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,250,4570CT 6.9% 174.9 / 2.2 58A76A (GCG→GCAcidBholin CidA anti‑holin regulator
Reads supporting (aligned to +/- strand):  ref base C (29/25);  new base T (2/2);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

ACAGGCTGGTTGCGCTGTCCAGCGCCAACAGGTGCGCCAGCCCCCAGGAGCTGGCGATCGAGGCACTGCTGCCGGCGACCATGCCCACCATCAGGGCCGGCCAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGG  >  NC_002947/5250315‑5250592
                                                                                                                                              |                                                                                                                                       
aCAGGCTGGTTGCGCGGTCCAGCGCCAACAGGTGCGCCAGCCCCCAGGAGCTGGCGATCGAGGCACTGCTGCCGGCGACCATGCCCACCATCAGGGCCGGCCAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAg                                                                                                                                   <  5:244948/149‑1 (MQ=255)
    gCTGGTTGCGCTGTCCAGCGCCAACAGGTGCGCCAGCCCCCAGGAGCTGGCGATCGAGGCACTGCTGCCGGCGACCATGCCCACCATCAGGGCCGGCCAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGa                                                                                                                               >  1:255631/1‑149 (MQ=255)
        gTTGCGCTGTCCAGCGCCAACAGGTGCGCCAGCCCCCAGGAGCTGGCGATCGAGGCACTGCTGCCGGCGACCATGCCCACCATCAGGGCCGGCCAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACa                                                                                                                             <  1:11765/147‑1 (MQ=255)
        gTTGCGCTGTCCAGCGCCAACAGGTGCGCCAGCCCCCAGGAGCTGGCGATCGAGGCACTGCTGCCGGCGACCATGCCCACCATCAGGGCCGGCCAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACa                                                                                                                             <  2:38429/147‑1 (MQ=255)
        gTTGCGCTGTCCAGCGCCAACAGGTGCGCCAGCCCCCAGGAGCTGGCGATCGAGGCACTGCTGCCGGCGACCATGCCCACCATCAGGGCCGGCCAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGt                                                                                                                           >  4:173136/1‑149 (MQ=255)
             gcTGTCCAGCGCCAACAGGTGCGCCAGCCCCCAGGAGCTGGCGATCGAGGCACTGCTGCCGGCGACCATGCCCACCATCAGGGCCGGCCAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCg                                                                                                                      >  6:162927/1‑149 (MQ=255)
                      cgcCAACAGGTGCGCCAGCCCCCAGGAGCTGGCGATCGAGGCACTGCTGCCGGCGACCATGCCCACCATCAGGGCCGGCCAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCa                                                                                                             >  3:214722/1‑149 (MQ=255)
                          aaCAGGTGCGCCAGCCCCCAGGAGCTGGCGATCGAGGCACTGCTGCCGGCGCCCATGCCCACCATCAGGGCCGGCCAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGCGCACGcc                                                                                                         >  6:369328/1‑149 (MQ=255)
                               gTGCGCCAGCCCCCAGGAGCTGGCGATCGAGGCACTGCTGCCGGCGACCATGCCCACCATCAGGGCCGGCCAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATc                                                                                                     >  5:261871/1‑148 (MQ=255)
                                     cAGCCCCCAGGAGCTGGCGATCGAGGCACTGCTGCCGGCGACCATGCCCACCATCAGGGCCGGCCAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAAc                                                                                                  >  6:84334/1‑145 (MQ=255)
                                     cAGCCCCCAGGAGCTGGCGATCGAGGCACTGCTGCCGGCGACCATGCCCACCATCAGGGCCGGCCAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAAc                                                                                                  <  5:84334/145‑1 (MQ=255)
                                       gCCCCCAGGAGCTGGCGATCGAGGCACTGCTGCCGGCGACCATGCCCACCATCAGGGCCGGCCAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTgg                                                                                            <  6:261871/149‑1 (MQ=255)
                                       gCCCCCAGGAGCTGGCGATCGAGGCACTGCTGCCGGCGACCATGCCCACCATCAGGGCCGGCCAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTgg                                                                                            >  6:137141/1‑149 (MQ=255)
                                       gCCCCCAGGAGCTGGCGATCGAGGCACTGCTGCCGGCGACCATGCCCACCATCAGGGCCGGCCAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTgg                                                                                            <  2:255631/149‑1 (MQ=255)
                                        cccccAGGAGCTGGCGATCGAGGCACTGCTGCCGGCGACCATGCCCACCATCAGGGCCGGCCAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTggg                                                                                           >  3:347336/1‑149 (MQ=255)
                                         ccccAGGAGCTGGCGATCGAGGCACTGCTGCCGGCGACCATGCCCACCATCAGGGCCGGCCAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGgt                                                                                          <  7:14298/149‑1 (MQ=255)
                                              ggAGCTGGCGATCGAGGCACTGCTGCCGGCGACCATGCCCACCATCAGGGCCGGCCAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTgc                                                                                     <  5:137141/149‑1 (MQ=255)
                                                aGCTGGCGATCGAGGCACTGCTGCCGGCGACCATGCCCACCATCAGGGCCGGCCAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTgcgc                                                                                   <  6:193654/149‑1 (MQ=255)
                                                     gCGATCGAGGCACTGCTGCCGGCGACCATGCCCACCATCAGGGCCGGCCAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACTGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATa                                                                              >  8:38191/1‑149 (MQ=255)
                                                              gCACTGCTGCCGGCGACCATGCCCACCATCAGGGCCGGCCAATGACGCGCCAGCATGGCCCGCTGCTCCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCa                                                                       >  2:165244/1‑147 (MQ=255)
                                                                aCTGCTGCCGGCGACCATGCCCACCATCAGGGCCGGCCAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAgg                                                                   >  1:225357/1‑149 (MQ=255)
                                                                    ctgcCGGCGACCATGCCCACCATCAGGGCCGGCCAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACTGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGt                                                               <  7:38191/149‑1 (MQ=255)
                                                                     tgcCGGCGACCATGCCCACCATCAGGGCCGGCCAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTa                                                              >  4:154966/1‑149 (MQ=255)
                                                                       ccGGCGACCATGCCCACCATCAGGGCCGGCCAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGNCGTa                                                              <  5:162927/147‑1 (MQ=255)
                                                                       ccGGCGACCATGCCCACCATCAGGGCCGGCCAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATg                                                            <  4:214722/149‑1 (MQ=255)
                                                                          gCGACCATGCCCACCATCAGGGCCGGCCAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACTGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAg                                                         >  5:298951/1‑149 (MQ=255)
                                                                           cGACCATGCCCACCATCAGGGCCGGCCAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGc                                                           >  6:211754/1‑146 (MQ=255)
                                                                           cGACCATGCCCACCATCAGGGCCGGCCAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGc                                                           <  5:211754/146‑1 (MQ=255)
                                                                               cATGCCCACCATCAGGGCCGGCCAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTTACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCa                                                    >  5:55569/1‑149 (MQ=255)
                                                                                     caccaTCAGGGCCGGCCAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACt                                                                            >  6:349775/1‑119 (MQ=255)
                                                                                     caccaTCAGGGCCGGCCAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACt                                                                            <  5:349775/119‑1 (MQ=255)
                                                                                     caccaTCAGGGCCGGCCAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAg                                                                    <  6:242243/127‑1 (MQ=255)
                                                                                     caccaTCAGGGCCGGCCAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAg                                                                    >  5:242243/1‑127 (MQ=255)
                                                                                       ccaTCAGGGCCGGCCAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCcaca                                            <  3:173136/149‑1 (MQ=255)
                                                                                       ccaTCAGGGCCGGCCAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCcaca                                            <  4:347336/149‑1 (MQ=255)
                                                                                             gggCCGGCCAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACTGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCAcagccagc                                      <  6:298951/149‑1 (MQ=255)
                                                                                               gccggccAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCAcagcca                                        <  5:286742/145‑1 (MQ=255)
                                                                                               gccggccAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCAcagcca                                        >  6:286742/1‑145 (MQ=255)
                                                                                                      aaTGACGCGCCAGCATGGCCCGCTGCTCCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACAcc                             <  1:165244/149‑1 (MQ=255)
                                                                                                            gcgcCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCAc                        >  1:138286/1‑148 (MQ=255)
                                                                                                             cgcCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCAc                        >  1:190596/1‑147 (MQ=255)
                                                                                                              gcCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCagcagc                                   <  7:164404/135‑1 (MQ=255)
                                                                                                              gcCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCagcagc                                   >  8:164404/1‑135 (MQ=255)
                                                                                                               ccAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTgg                                                                                            >  4:197027/1‑77 (MQ=255)
                                                                                                               ccAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTgg                                                                                            <  3:197027/77‑1 (MQ=255)
                                                                                                               ccAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAg                     >  2:89200/1‑148 (MQ=255)
                                                                                                                cAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGa                   >  5:99349/1‑149 (MQ=255)
                                                                                                                  gCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCAcag                                           >  7:222613/1‑123 (MQ=255)
                                                                                                                  gCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCAcag                                           <  8:222613/123‑1 (MQ=255)
                                                                                                                         ccGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAgg                    >  3:228708/1‑139 (MQ=255)
                                                                                                                         ccGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAgg                    <  4:228708/139‑1 (MQ=255)
                                                                                                                            ctgctgCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGAca       <  1:89200/149‑1 (MQ=255)
                                                                                                                                 gcCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAgg  >  4:61065/1‑149 (MQ=255)
                                                                                                                                  ccAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGc                                                <  6:342385/102‑1 (MQ=255)
                                                                                                                                  ccAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGc                                                >  5:342385/1‑102 (MQ=255)
                                                                                                                                   cAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAgg  >  6:230557/1‑147 (MQ=255)
                                                                                                                                            aCCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGt                                                                                              >  2:247062/1‑46 (MQ=255)
                                                                                                                                            aCCGCGAAGACGACAGTCACCGGGCCGAGCACGCCCATCAACCAGt                                                                                              <  1:247062/46‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                       
ACAGGCTGGTTGCGCTGTCCAGCGCCAACAGGTGCGCCAGCCCCCAGGAGCTGGCGATCGAGGCACTGCTGCCGGCGACCATGCCCACCATCAGGGCCGGCCAATGACGCGCCAGCATGGCCCGCTGCTGCCAGATCGGCACCGCGAAGGCGACAGTCACCGGGCCGAGCACGCCCATCAACCAGTGGGTGTTGCGCGCATACTCGGCATAGGCCGTATGCAGCGGCACGGCCACAGCCAGCAGCAACACCGGCACCAGGATCAACGGCGACAGCAGG  >  NC_002947/5250315‑5250592

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