Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 906,219 T→C 7.5% V76A (GTG→GCG)  ampD → N‑acetyl‑anhydromuramyl‑L‑alanine‑amidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947906,2190TC 7.5% 151.0 / 2.0 53V76A (GTG→GCG) ampDN‑acetyl‑anhydromuramyl‑L‑alanine‑amidase
Reads supporting (aligned to +/- strand):  ref base T (26/23);  new base C (2/2);  total (28/25)
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.95e-01

GAAGGCGAATCGATTTCCCTGCTGGTGATCCACAACATCAGTCTGCCCCCGGCCTGTTTCGGCACGGGCAAGGTGCAGCAATTCTTCCAGAACCGCCTGGCCCCTGATGAACACCCGTATTTTGCCAGCATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGACTCAGTTCGTGTCGTTGCATGACCGTGCCTGGCATGCGGGTGTGTCGAGCTTTGACGGGCGCGAAGGCTGTAACGACTTTTCCATCGGCATCGAGCTGGAA  >  NC_002947/906074‑906358
                                                                                                                                                 |                                                                                                                                            
gAAGGCGAATCGATTTCCCTGCTGGTGATCCACAACCTCAGTCTGCCCCCGGCCTGTTTCGGCACAGGCAAGGTGCAGCAATTCTTCCAGAACCGCCTGGCCCCTGATGAACACCCGTATTTTGCCAGCATCAACCACCTGACCGTGTCg                                                                                                                                          >  1:63506/1‑150 (MQ=255)
      gAATCGATTTCCCTGCTGGTGATCCACAACATCAGTCTGCCCCCGGCCTGTTTCGGCACGGGCAAGGTGCAGCAATTCTTCCAGAACCGCCTGGCCCCTGATGAACACCCGTATTTTGCCAGCATCAACCACCTGACCGTGTCGGCGCAc                                                                                                                                    <  2:238185/150‑1 (MQ=255)
                   tgctgGTGATCCACAACATCAGTCTGCCCCCGGCCTGTTTCGGCACGGGCAAGGTGCAGCAATTCTTCCAGAACCGCCTGGCCCCTGATGAACACCCGTATTTTGCCAGCATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGAGc                                                                                                                       <  1:1233497/150‑1 (MQ=255)
                               acaacaTCAGTCTGCCCCCGGCCTGTTTCGGCACGGGCAAGGTGCAGCAATTCTTCCAGAACCGCCTGGCCCCTGATGAACACCCGTATTTTGCCAGCATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑Cgg                                                                                                          >  1:697613/1‑150 (MQ=255)
                               acaacaTCAGTCTGCCCCCGGCCTGTTTCGGCACGGGCAAGGTGCAGCAATTCTTCCAGAACCGCCTGGCCCCTGATGAACACCCGTATTTTGCCAGCATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑Cgg                                                                                                          <  1:242475/150‑1 (MQ=255)
                                  acaTCAGTCTGCCCCCGGCCTGTTTCGGCACGGGCAAGGTGCAGCAATTCTTCCAGAACCGCCTGGCCCCTGATGAACACCCGTATTTTGCCAGCATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGa                                                                                                       <  2:1062807/150‑1 (MQ=255)
                                         tCTGCCCCCGGCCTGTTTCGGCACGGGCAAGGTGCAGCAATTCTTCCAGAACCGCGTGGCCCCTGATGAACACCCGTATTTTGCCAGCATCAACCACCTGCCCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGACTCAGtt                                                                                                >  2:254208/1‑150 (MQ=255)
                                              ccccGGCCTGTTTCGGCACGGGCAAGGTGCAGCAATTCTTCCAGAACCGCCTGGCCCCTGATGAACACCCGTATTTTGCCAGCATCAACCACCTGACCTTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGACTCAGTTCgtgt                                                                                           >  1:976307/1‑150 (MQ=255)
                                              ccccGGCCTGTTTCGGCACGGGCAAGGTGCAGCAATTCTTCCAGAACCGCCTGGCCCCTGATGAACACCCGTATTTTGCCAGCATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGACTCAGTTCgtgt                                                                                           <  2:63506/150‑1 (MQ=255)
                                               cccGGCCTGTTTCGGCACGGGCAAGGTGCAGCAATTCTTCCAGAACCGCCTGGCCCCTGATGAACACCCGTATTTTGCCAGCATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGACTCAGTTCgtgt                                                                                           <  2:1284761/149‑1 (MQ=255)
                                               cccGGCCTGTTTCGGCACGGGCAAGGTGCAGCAATTCTTCCAGAACCGCCTGGCCCCTGATGAACACCCGTATTTTGCCAGCATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGACTCAGTTCgtgt                                                                                           >  1:1284761/1‑149 (MQ=255)
                                               cccGGCCTGTTTCGGCACGGGCAAGGTGCAGCAATTCTTCCAGAACCGCCTGGCCCCTGATGAACACCCGTATTTTGCCAGCATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGACTCAGTTCGTGTc                                                                                          >  1:444496/1‑150 (MQ=255)
                                                  ggCCTGTTTCGGCACGGGCAAGGTGCAGCAATTCTTCCAGAACCGCCTGGCCCCTGATGAACACCCGTATTTTGCCAGCATCAACCACCTGACCTTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGACTCAGTTCGTGTCGtc                                                                                       <  2:976307/150‑2 (MQ=255)
                                                  ggCCTGTTTCGGCACGGGCAAGGTGCAGCAATTCTTCCAGAACCGCCTGGCCCCTGATGAACACCCGTATTTTGCCAGCATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGACTCAGTTCGTGTCGtt                                                                                       <  1:846734/150‑1 (MQ=255)
                                                             gCACGGGCAAGGTGCAGCAATTCTTCCAGAACCGCCTGGCCCCTGATGAACACCCGTATTTTGCCAGCATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGACTCAGTTCGTGTCGTTGCATGACCGTg                                                                            >  1:1252082/1‑150 (MQ=255)
                                                                cGGGCAATGTGCAGCAATTCTTCCAGAACCGCCTGGCCCCTGATGAACACCCGTATTTTGCCAGCATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGACTCAGTTCGTGTCGTTGCATGACCGTGCCt                                                                         <  2:8235/150‑1 (MQ=255)
                                                                  ggCAAGGTGCAGCAATTCTTCCAGAACCGCCTGGCCCCTGATGATCACCCGTATTTTGCCAGCATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGACTCAGTTCGTGTCGTTGCATGACCGTGCCTgg                                                                       <  2:1028217/150‑1 (MQ=255)
                                                                   gCAAGGTGCAGCAATTCTTCCAGAACCGCCTGGCCCCTGATGCACACCCGAGTTTTGCCAGCATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGAGGGTGACGGCG‑CGGTGACTCAGTTCGTGTCGTTGCATGACCGTGCCTGGc                                                                      >  1:1327038/1‑150 (MQ=255)
                                                                   gCAAGGTGCAGCAATTCTTCCAGAACCGCCTGGCCCCTGATGAACACCCGTATTTTGCCAGCATCAACAACCTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGACTCAGTTCGTGTCGTTGCATGACCGTGCCTGGc                                                                      >  1:349082/1‑150 (MQ=255)
                                                                     aaGGTGCAGCAATTCTTCCAGAACCGCCTGGCCCCTGATGAACACCCGTATTTTGCCAGCATCAACCACCTGACCGCGTCGGCGCACCTGTTCGTCGAGCGTGa                                                                                                                   >  2:631966/1‑104 (MQ=255)
                                                                     aaGGTGCAGCAATTCTTCCAGAACCGCCTGGCCCCGGATGAACACCCGTATTTTGCCAGCATCAACCACCTGACCGCGTCGGCGCACCTGTTCGTCGAGCGTGa                                                                                                                   <  1:631966/104‑1 (MQ=255)
                                                                       ggTGCAGCAATTCTTCCAGAACCGCCTGGCCCCTGATGAACACCCGTATTTTGCCAGCATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGt                                                                                                         <  1:188537/111‑1 (MQ=255)
                                                                       ggTGCAGCAATTCTTCCAGAACCGCCTGGCCCCTGATGAACACCCGTATTTTGCCAGCATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGt                                                                                                         >  2:188537/1‑111 (MQ=255)
                                                                        gTGCAGCAATTCTTCCAGAACCGCCTGGCCCCTGATGAACACCCGTATTTTGCCAGCATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGACTCAGTTCGTGTCGTTGCATGa                                                                                 >  1:606946/1‑134 (MQ=255)
                                                                        gTGCAGCAATTCTTCCAGAACCGCCTGGCCCCTGATGAACACCCGTATTTTGCCAGCATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGACTCAGTTCGTGTCGTTGCATGa                                                                                 <  2:606946/134‑1 (MQ=255)
                                                                         tGCAGCAATTCTTCCAGAACCGCCTGGCCCCTGATGAACACCCGTATTTTGCCAGCATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGACTGAGTTCGTGTCGTTGCATGACCGTGCCTGGCATGCgg                                                                <  2:549847/150‑1 (MQ=255)
                                                                          gcagcaGTTCTTCCAGAACCGCCTGGCCCCTGATGAACACCCGTATTTTGCCAGCATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGACTCAGTTCGTGTCGTTGCATGAGCGTGCCTGGCATGCggg                                                               >  1:1190787/1‑150 (MQ=255)
                                                                           cagcaATTCTTCCAGAACCGCCTGGCCCCTGATGAACACCCGTATTTTGCCAGCATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGACTCAGTTCGTGTCGTTGCATGACCGTGCCTGGCATGCGGgt                                                              <  2:1387569/150‑1 (MQ=255)
                                                                             gcaATTCTTCCAGAACCGCCTGGCCCCTGATGAACACCCGTATTTTGCCAGCATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGa                                                                                                                         <  1:35337/90‑1 (MQ=255)
                                                                             gcaATTCTTCCAGAACCGCCTGGCCCCTGATGAACACCCGTATTTTGCCAGCATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGa                                                                                                                         >  2:35337/1‑90 (MQ=255)
                                                                                 ttcttcCAGAACCGCCTGGCCCCTGATGAACACCCGTATTTTGCCAGCATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGACTCAGTTCGTGTCGTTGCATGACCGTGCCTGGCATGCGGGTGTGTCg                                                        >  2:1240303/1‑150 (MQ=255)
                                                                                    ttcCAGAACCGCCTGGCCCCTGATGAACACCCGTATTTTGCCAGCATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGACTCAGTTCGTGTCGTTGCATGACCGTGCCTGGCATGCGGGTGTGTCGAGc                                                     <  1:1240303/150‑1 (MQ=255)
                                                                                     tcCAGAACCGCCTGGCCCCTGATGAACACCCGTATTTTGCCAGCATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGACTCAGTTCGTGTCGTTGCATGACCGTGCCTGGCATGCGGGTGTGTCGAGCt                                                    >  1:1371932/1‑150 (MQ=255)
                                                                                       cAGAACCGCCTGGCCCCTGATGAACACCCGTATTTTGCCAGCATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGACTCAATTCGTTTCGTTGCATGACCGTGCCTGGCATGCGGGTGTGTCGAGCttt                                                  <  1:1396155/150‑1 (MQ=255)
                                                                                         gAACCGCCTGGCCCCTGATGAACACCCGTATTTTGCCAGCATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGAGTCAGTTCGTGTTGTTGCATGACCGTGCCTGGCATGAGGGTGCGTCGAGCTTTGa                                                >  2:469062/1‑150 (MQ=255)
                                                                                                     cccTGATGAACACCCGTATTTTGCCAGCATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGACTCCGTTCGGGCCGTTGCATGACCGTGCCTGGCATGCGGGTGTGTCGAGCTTTGACGGGCGCGAAgg                                    <  2:126288/150‑1 (MQ=255)
                                                                                                     cacTGATGAACACCCGTATTTTGCCAGCATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGACTCAGTTCGTGTCGTTGCATGACCGTGCCTGGCATGCGGGTGTGTCGAGCTTTGACGGGCGCGAAgg                                    >  2:1347759/3‑150 (MQ=255)
                                                                                                      ccTGATGAACACCCGTATTTTGCCAGCATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGACTCAGTTCGTGTCGTTGCATGACCGTGCCTGGCATGCGGGTGTGTCGAGCTTTGACGGGCGCGAAGGc                                   >  2:16121/1‑150 (MQ=255)
                                                                                                        tgatgaACACCCGTATTTTGCCAGCATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGACTCAGTTCGTGTCGTTGCATGACCGTGCCTGGCATGCGGGTGTGTCGAGCTTTGACGGGCGCGAAGGCTg                                 >  2:685654/1‑150 (MQ=255)
                                                                                                               cacCCGTATTTTGCCAGCATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGACTCAGTTCGTGTCGTTGCATGACCGTGCCTGTCTTGCGGGTGTGTCGAGCTTTGACGGGCGCGAAGGCTGTAACGAc                          >  1:1224036/1‑150 (MQ=255)
                                                                                                                  ccGTATTTTGCCAGCATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGACTCAGTTCGTGTCGTTGCATGACCGTGCCTGGCATGCGGGTGTGTCGAGCTTTGACGGGCGCGAAGGCTGTAACGACttt                       >  2:643387/1‑150 (MQ=255)
                                                                                                                             cAGCATCAACCACCTGACCGTGTCGGCGCCACTGCAACTACAGCGTCCCGAAGACTATGGCCCAG‑CCCTGTCTTGCCATGACCGTGCCTGGCATGCGGGGGTGCCGGGCTTGGACGGGCCCGAAGGCTGTAACCCCCCTTCCCTCGGCAt            >  1:26586/1‑150 (MQ=255)
                                                                                                                               gCATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGACTCAGTTCGTGTCGTTGCATGACCGTGCCTGGCATGCgg                                                                <  1:920937/96‑1 (MQ=255)
                                                                                                                               gCATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGACTCAGTTCGTGTCGTTGCATGACCGTGCCTGGCATGCgg                                                                >  2:920937/1‑96 (MQ=255)
                                                                                                                                cATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGa                                                                                                       >  2:616848/1‑56 (MQ=255)
                                                                                                                                cATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGa                                                                                                       <  1:616848/56‑1 (MQ=255)
                                                                                                                                    aaCCACCTGACCGCGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGACTCAGTTCGTGTCGTTGCATGACCGTGCCTGGCATGCGGGTGTGTCGAGCTTTGACGGGCGCGAAGGCTGTAACGACTTTTCCATCGGCATCGAGCTg     >  2:321511/1‑150 (MQ=255)
                                                                                                                                     accaccTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGACTCAGTTCGTGTCGTTGCATGACCGTGCCTGGCATGCGGGTGTGTCGAGCTTTGACGGGCGCGAAGGCTGTAACGACTTTTCCATCGGCATCGAGCTgg    >  1:1326902/1‑150 (MQ=255)
                                                                                                                                       caccTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGACTCAGTTCGTGTCGTTGCATGACCGTGCCTGGCATGCGGGTGTGTCGAGCTTTGACGGGCGCGAAGGCTGTAACGACTTTTCCATCGGCATCGAGCTGGaa  <  1:685654/150‑1 (MQ=255)
                                                                                                                                       caccTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGACTCAGTTCGTGTCGTTGCATGACCGTGCCTGGCATGCGGGTGTGTCGAGCTTTGACGGGCGCGAAGGCTGTAACGACTTTTCCATCGGCATCGAGCTGGaa  <  1:1190788/150‑1 (MQ=255)
                                                                                                                                       caccTGACCGCGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGACTCAGTTCGTGTCGTTGCATGACCGTGCCTGGCATGCGGGTGTGTCGAGCTTTGACGGGCGCGAAGGCTGTAACGACTTTTCCATCGGCATCGAGCTGGaa  <  1:321511/150‑1 (MQ=255)
                                                                                                                                           tGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGACTCAGTTCGTGTCGTTGCATGACCGTGCCTGGCATGCGGGTGTGTCGAGCTTTGACGGGCGCGAAGGCTGTaa                              <  2:334788/118‑1 (MQ=255)
                                                                                                                                           tGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGACTCAGTTCGTGTCGTTGCATGACCGTGCCTGGCATGCGGGTGTGTCGAGCTTTGACGGGCGCGAAGGCTGTaa                              >  1:334788/1‑118 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
GAAGGCGAATCGATTTCCCTGCTGGTGATCCACAACATCAGTCTGCCCCCGGCCTGTTTCGGCACGGGCAAGGTGCAGCAATTCTTCCAGAACCGCCTGGCCCCTGATGAACACCCGTATTTTGCCAGCATCAACCACCTGACCGTGTCGGCGCACCTGTTCGTCGAGCGTGACGGCG‑CGGTGACTCAGTTCGTGTCGTTGCATGACCGTGCCTGGCATGCGGGTGTGTCGAGCTTTGACGGGCGCGAAGGCTGTAACGACTTTTCCATCGGCATCGAGCTGGAA  >  NC_002947/906074‑906358

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑