Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,462,261 G→A 7.2% D401N (GAT→AAT) ‡ PP_4798 → membrane‑bound lytic murein transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,462,2610GA 7.2% 162.5 / 2.2 55D401N (GAT→AAT) ‡PP_4798membrane‑bound lytic murein transglycosylase
Reads supporting (aligned to +/- strand):  ref base G (22/29);  new base A (2/2);  total (24/31)
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.55e-01

GGTGAGCCTGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGT‑TGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGG  >  NC_002947/5462119‑5462405
                                                                                                                                               |                                                                                                                                                
ggTGAGCCTGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGT‑TGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCa                                                                                                                                            <  3:240705/149‑1 (MQ=255)
     gCCTGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGT‑TGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGATAATGCCATCGGCATTacc                                                                                                                                       <  3:262026/149‑3 (MQ=255)
          cTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGT‑TGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCggg                                                                                                                                    <  7:197912/147‑1 (MQ=255)
          cTCGGGGATCGCTTAAGCGGCTGGGGCTTGATTGCCGGGAGCTGGCCGAAAGAGGACCTGCCGATCAGCCGCACCGAGCTCATGGAGT‑TTCAGAACTTGCTGAACAGCCGCGGGCATGATGCGTGTAATGCCGATGGCATTATCggg                                                                                                                                    <  8:250379/147‑1 (MQ=255)
               ggATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGT‑TGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATa                                                                                                                             <  1:38449/149‑1 (MQ=255)
                   cGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGT‑TGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCg                                                                                                                         >  4:288216/1‑149 (MQ=255)
                      ttCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGT‑TGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGc                                                                                                                        >  4:232912/1‑147 (MQ=255)
                           cGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGT‑TGCAGAACTTGCTGAACAGCCGCTGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCaa                                                                                                                 <  1:270731/149‑1 (MQ=255)
                           cGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGT‑TGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCaa                                                                                                                 <  7:199363/149‑1 (MQ=255)
                           cGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGT‑TGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCaa                                                                                                                 <  4:282829/149‑1 (MQ=255)
                            ggCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGT‑TGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAAt                                                                                                                <  8:58129/149‑1 (MQ=255)
                                  ggCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGT‑TGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCaa                                                                                                          >  6:386027/1‑149 (MQ=255)
                                       cATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGT‑TGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGcc                                                                                                      >  2:325984/1‑148 (MQ=255)
                                       cATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGT‑TGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGcc                                                                                                      <  1:325984/148‑1 (MQ=255)
                                       cATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGT‑TGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGcc                                                                                                      >  7:22635/1‑148 (MQ=255)
                                         ttGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGT‑TGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCTGCAAGGCAATCCGCAATGCCCAg                                                                                                   >  7:301189/1‑149 (MQ=255)
                                            ccGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGT‑TGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCggg                                                                                                                                    <  3:219311/113‑1 (MQ=255)
                                            ccGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGT‑TGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCggg                                                                                                                                    >  4:219311/1‑113 (MQ=255)
                                             cGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGACGCACCGAGCGCATTGAGT‑TGGAGAACTTGCTGAACAGCCGCGGGCATGAAGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGTCAATCCGCAATGCCCAGCAAg                                                                                               <  2:193287/149‑1 (MQ=255)
                                                   cTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGT‑TGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTgg                                                                                         <  5:386027/149‑1 (MQ=255)
                                                      gCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGT‑TGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGAt                                                                                      >  8:319326/1‑149 (MQ=255)
                                                       ccGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGT‑TGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATg                                                                                     <  3:232912/149‑1 (MQ=255)
                                                        cGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGT‑TGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATgg                                                                                    >  4:309007/1‑149 (MQ=255)
                                                                  ccTGCCGCTCAGCCGCACCGAGCGCATGGAGT‑TGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGc                                                                                                                                             <  7:302974/82‑1 (MQ=255)
                                                                  ccTGCCGCTCAGCCGCACCGAGCGCATGGAGT‑TGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGc                                                                                                                                             >  8:302974/1‑82 (MQ=255)
                                                                       cGCTCAGCCGCACCGAGCGCATGGAGTGTGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTa                                                                      >  1:17722/1‑149 (MQ=255)
                                                                       cGCTCAGCCGCACCGAGCGCATGGAGT‑TGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTAc                                                                     <  7:319326/149‑1 (MQ=255)
                                                                         cTCAGCCGCACCGAGCGCATGGAGT‑TGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTAccc                                                                   >  5:90762/1‑149 (MQ=255)
                                                                         cTCAGCCGCACCGAGCGCATGGAGT‑TGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTAc                                                                     >  3:50313/1‑147 (MQ=255)
                                                                         cTCAGCCGCACCGAGCGCATGGAGT‑TGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTAc                                                                     <  4:50313/147‑1 (MQ=255)
                                                                            aGCCGCACCGAGCGCATGGAGT‑TGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGAc                                                                <  3:229560/149‑1 (MQ=255)
                                                                                  aCCGAGCGCATGGAGT‑TGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCcac                                                            >  1:341023/1‑147 (MQ=255)
                                                                                   ccGAGCGCATGGAGT‑TGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAg                                                         >  1:145487/1‑149 (MQ=255)
                                                                                         gcATGGAGT‑TGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCtt                                                   <  1:253736/149‑1 (MQ=255)
                                                                                           aTGGAGT‑TGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTga                                                 <  5:40190/149‑1 (MQ=255)
                                                                                              gAGT‑TGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTgaga                                               <  6:90762/148‑1 (MQ=255)
                                                                                                     cAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTgc                                         <  2:145487/148‑1 (MQ=255)
                                                                                                      aGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTgcgc                                       >  8:170731/1‑149 (MQ=255)
                                                                                                        aaCTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGAc                                                                >  3:158071/1‑122 (MQ=255)
                                                                                                        aaCTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGAc                                                                <  4:158071/122‑1 (MQ=255)
                                                                                                        aaCTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCca                                     <  8:22635/149‑1 (MQ=255)
                                                                                                              cTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGt                                <  2:17722/148‑1 (MQ=255)
                                                                                                              cTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTa                               <  3:309007/149‑1 (MQ=255)
                                                                                                               tGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGt                                >  4:120384/1‑147 (MQ=255)
                                                                                                                 aaCAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAg                             <  7:170731/148‑1 (MQ=255)
                                                                                                                      ccGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTcc                       >  5:161037/1‑149 (MQ=255)
                                                                                                                        gcgGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCgg                     <  2:341023/149‑1 (MQ=255)
                                                                                                                              atgatgCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGc               >  8:13789/1‑149 (MQ=255)
                                                                                                                                   gCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTgcgc          >  3:314374/1‑149 (MQ=255)
                                                                                                                                    ccgatAATGCCATCGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGcgcg         >  1:382175/5‑149 (MQ=255)
                                                                                                                                     gggTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCgg        <  3:120384/149‑1 (MQ=255)
                                                                                                                                      gatAATGCCATCGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTgc                                         <  8:136573/113‑1 (MQ=255)
                                                                                                                                      gatAATGCCATCGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTgc                                         >  7:136573/3‑115 (MQ=255)
                                                                                                                                        tAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGACCGCTGCGCGGccc     <  8:210962/149‑1 (MQ=255)
                                                                                                                                            gCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCgg  >  5:241321/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                
GGTGAGCCTGCTCGGGGATCGCTTCAGCGGCTGGGGCTTCATTGCCGGCAGCTGGCCGAAAGAGGACCTGCCGCTCAGCCGCACCGAGCGCATGGAGT‑TGCAGAACTTGCTGAACAGCCGCGGGCATGATGCGGGTAATGCCGATGGCATTATCGGGGCCAATACCCGCAAGGCAATCCGCAATGCCCAGCAAGGGCTGGGATGGCCGGCGGATGGCTACCCGACCCACAAGCTGCTTGAGAGCCTGCGCCAGCAGTAAGCATTCCGGGGCCGCTGCGCGGCCCCGG  >  NC_002947/5462119‑5462405

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: ‑