Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 906,437 G→A 7.5% A150T (GCC→ACC)  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,4370GA 7.5% 154.1 / 2.0 53A150T (GCC→ACC) ampDN‑acetyl‑anhydromuramyl‑L‑alanine‑amidase
Reads supporting (aligned to +/- strand):  ref base G (25/24);  new base A (2/2);  total (27/26)
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.54e-01

CGGGTGTGTCGAGCTTTGACGGGCGCGAAGGCTGTAACGACTTTTCCATCGGCATCGAGCTGGAAGGCACCGACGACCTGCCCTATACCGATGCCCAGTATGCGGTGCTGGAGCAACTCACCCGGCACATCCGCAGTGCCTGGCCAGCCATCGACCTGAACCGCATTCAGGGCCACAGCGACATTGCCCCGCAGCGCAAGACCGACCCCGGCCCGGCCTTCGACTGGTCGCGCTACCGCCAGGCGCTGCAGCACAACGAGGACAAGGCATGAGTTTTCTGGTGTTGTTGCTG  >  NC_002947/906291‑906582
                                                                                                                                                  |                                                                                                                                                 
cGGGTGTGTCGAGCTTTGACGGGCGCGAAGGCTGTAACGACTTTTCCATCGGCATCGAGCTGGAAGGCACCGACGACCTGCCCTATACCGATGCCCAGTATGCGGTGCTGGAGCAACTCACCCGGCACATCCGCAGTGCCTGgccagcc                                                                                                                                                 <  3:225990/149‑1 (MQ=255)
           aGCTTTGACGGGCGCGAAGGCTGTAACGACTTTTCCATCGGCATCGAGCTGGAAGGCACCGACGACCTGCCCTATACCGATGCCCAGTATGCGGTGCTGGAGCAACTCACCCGGCACATCCGCAGTGCCTGGCCAGCCATCGACCTGaa                                                                                                                                      <  7:53297/149‑1 (MQ=255)
             cTTTGACGGGCGCGAAGGCTGTAACGACTTTTCCATCGGCATCGAGCTGGAAGGCTCCGACGACCTGCCCTATACCGATGCCCAGTATGCGGTGCTGGAGCAACTCACCCGGCACATCCGCAGTGCCTGGCCAGCCATCGACCTGAAc                                                                                                                                     <  6:82474/148‑1 (MQ=255)
                  aCGGGCGCGAAGGCTGTAACGACTTTTCCATCGGCATCGAGCTGGAAGGCACCGACGACCTGCCCTATACCGATGCCCAGTATGCGGTGCTGGAGCAACTCACCCGGCACATCCGAAGTGCCTGGCCAGCCATCGACCTGAACCGCAt                                                                                                                                >  1:147277/1‑148 (MQ=255)
                      gcgcGAAGGCTGTAACGACTTTTCCATCGGCATCGAGCTGGAAGGCACCGACGACCTGCCCTATACCGATGCCCAGTATGCGGTGCTGGAGCAACTCACCCGGCACATCCGCAGTGCCTGGCCTGCCATCGAcc                                                                                                                                          <  3:198413/134‑1 (MQ=255)
                      gcgcGAAGGCTGTAACGACTTTTCCATCGGCATCGAGCTGGAAGGCACCGACGACCTGCCCTATACCGATGCCCAGTATGCGGTGCTGGAGCAACTCACCCGGCACATCCGCAGTGCCTGGCCCGCCATCGAcc                                                                                                                                          >  4:198413/1‑134 (MQ=255)
                          gAAGGCTGTAACGACTTTTCCATCGGCATCGAGCTGGAAGGCACCGACGACCTGCCCTATACCGATGCCCAGTATGCGGTGCTGGAGCAACTCACCCGGCACATCCGCAGTGCCTGGCTAGCCATCGACCTGACCCGCATTAAGGGAca                                                                                                                       >  8:313509/1‑149 (MQ=255)
                           aaGGCTGTAACGACTTTTCCATCGGCATCGAGCTGGAAGGCACCGACGACCTGCCCTATACCGATGCCCAGTATGCGGTGCTGGAGCAACTCACCCGGCACATCCGCAGTGCCTGGCCAGCCATCGACCTGAACCGCATTCAGGGCcac                                                                                                                      <  6:169158/149‑1 (MQ=255)
                           aaGGCTGTAACGACTTTTCCATCGGCATCGAGCTGGAAGGCACCGACGACCTGCCCTATACCGATGCCCAGTATGCGGTGCTGGAGCAACTCACCCGGCACATCCGCAGTGCCTGGCCAGCCATCGACCTGAACCGCATTCAGGGCca                                                                                                                       >  1:46435/1‑148 (MQ=255)
                            aGGCTGTAACGACTTTTCCATCGGCATCGAGCTGGAAGGCACCGACGACCTGCCCTATACCGATGCCCAGTATGCGGTGCTGGAGCAACTCACCCGGCACATCCGCAGTGCCTGGCCAGCCATCGACCTGAACCGCATTCAGGGCcaca                                                                                                                     <  2:323335/149‑1 (MQ=255)
                             ggCTGTAACGACTTTTCCATCGGCATCGAGCTGGAAGGCACCGACGACCTGCCCTATACCGATGCCCAGTATGCGGTGCTGGAGCAACTCACCCGGCACATCCGCAGTGCCTGGCCAGCCATCGACCTGAACCGCATTCAGGGCCACAg                                                                                                                    <  3:302701/149‑1 (MQ=255)
                                    tcagctcttCCAACGGCATCGAGCTGGAAGGCACCGACGACCTGCCCTATACCGATGCCCAGTATGCGGTGCTGGAGCAACTCACCCGGCACATCCGCAGTGCCTGGCCAGCCATCGACCTGAACCGCATTCAGGGCCACAGCGACAtt                                                                                                             >  8:164938/8‑149 (MQ=255)
                                    aCGACTTTTCCATCGGCATCGAGCTGGAAGGCACCGACGACCTGCCCTATACCGATGCCCAGTATGCGGTGCTGGAGCAACTCACCCGGCACATCCGCAGTGCCTGGCCAGCCATCGACCTGAACCGCATTCAGGGCCACAGCGACAtt                                                                                                             >  1:117516/1‑149 (MQ=255)
                                    aCGACTTTTCCATCGGCATCGAGCTGGAAGGCACCGACGACCTGCCCTATACCGATGCCCAGTATGCGGTGCTGGAGCAACTCACCCGGCACATCCGCAGTGCCTGGCCAGCCATCGACCTGAACCGCATTCAGGGCCACAGCGACAtt                                                                                                             <  3:57609/149‑1 (MQ=255)
                                     cGACTTTTCCATCGGCATCGAGCTGGAAGGCACCGACGACCTGCCCTATACCGATGCCCAGTATGCGGTGCTGGAGCAACTCACCCGGCACATCCGCAGTGCCTGGCCAGCCATCGACCTGAACCGCATTCAGGGCCACAGCGACATTg                                                                                                            >  1:36062/1‑149 (MQ=255)
                                      gACTTTTCCATCGGCATCGAGCTGGAAGGCACCGACGACCTGCCCTATACCGATGCCCAGTATGCGGTGCTGGAGCAACTCACCCGGCACATCCGCAGTGCCTGGCCAACCATCGACCTGAACCGCATTCAGGGCCACAGCGACATTGc                                                                                                           >  5:111931/1‑149 (MQ=255)
                                         ttttCCATCGGCATCGAGCTGGAAGGCACCGACGACCTGCCCTATACCGATGCCCAGTATGCGGTGCTGGAGCAACTCACCCGGCACATCCGCAGTGCCTGGCCAGCCATCGACCTGAACCGCATTc                                                                                                                              <  1:266880/127‑1 (MQ=255)
                                         ttttCCATCGGCATCGAGCTGGAAGGCACCGACGACCTGCCCTATACCGATGCCCAGTATGCGGTGCTGGAGCAACTCACCCGGCACATCCGCAGTGCCTGGCCAGCCATCGACCTGAACCGCATTc                                                                                                                              >  2:266880/1‑127 (MQ=255)
                                                  ggCATCGAGCTGGAAGGCACCGACGACCTGCCCTATACCGATGCCCAGTATGCGGTGCTGGAGCAACTCACCCGGCACATCCGCAGTGCCTGGCCAGCCATCGACCTGAACCGCATTCAGGGCCACAGCGACATTGCCCCGCAGCGCaa                                                                                               <  3:315551/149‑1 (MQ=255)
                                                   gCATCGAGATGGAAGGCACCGACGACCTGCCCTATACCGATGCCCAGTATGCGGAGCTGGAGCAACTCACCCGGCACATCCGCAGTGCCTGGCCAGACATCCACCTGAACCGCATTCAGGGCCACAGCGACATTGCCCACCAgcgc                                                                                                 >  4:129193/1‑146 (MQ=255)
                                                                aGGCACCGACGACCTGCCCTATACCGATGCCCAGTATGCGGTGCTGGAGCAACTCACCCGGCACATCCGCAGTGCCTGGCCAGCCATCGa                                                                                                                                            >  8:244703/1‑90 (MQ=255)
                                                                aGGCACCGACGACCTGCCCTATACCGATGCCCAGTATGCGGTGCTGGAGCAACTCACCCGGCACATCCGCAGTGCCTGGCCAGCCATCGa                                                                                                                                            <  7:244703/90‑1 (MQ=255)
                                                                   cACCGACGACCTGCCCTATACCGATGCCCAGTATGCGGTGCTGGAGCAACTCACCCGGCACATCCGGAGTGCCTGGCCAGCCATCGACCTGAACCGCATTCAGGGCCACAGCGACATTGCCCCGCAGCGCAAGACCGACcccggcccgg                                                                              <  2:147277/149‑1 (MQ=255)
                                                                      cgacgaCCTGCCCTATACCGATGCCCAGTATGCGGTGCTGGAGCAACTCACACGGCACATCCGCAGTGCCTGGCCAGCCATCGACCTGAACCGCATTCAGGGCCACAGCGACATTGCCCCGCAGCGCAAGACCGACCCCGGCCCGGCCt                                                                           >  6:131350/1‑149 (MQ=255)
                                                                       gacgacCTGCCCTATACCGATGCCCAGTATGCGGTGCTGGAGCAACTCACCCGGCACATCCGCAGTGCCTGGCCAGCCATCGACCTGAACCGCATTCAGGGCCACAGCGACATTGCCCCGCAGCGCAAGACCGACCCCGGCCCGGCCtt                                                                          <  2:117516/149‑1 (MQ=255)
                                                                        acgacCTGCCCTATACCGATGCCCAGTATGCGGTGCTGGAGCAACTCACCCGGCACATCCGCAGTGCCTGGCCAGCCATCGACCTGAACCGCATTCAGGGCCACAGCGACATTGCCCCGCAGCGCAAGACCGACCCCGGCCCGGCCTTc                                                                         <  2:36062/149‑1 (MQ=255)
                                                                           acCTGCCCTATACCGATGCCCAGTATGCGGTGCTGGAGCAACTCACCCGGCACATCCGCAGTGCCTGGCCAACCATCGACCTGAACCGCATTCAGGGCCACAGCGACATTGCCCCGCAGCGCAAGACCGACCCCGGCCCGGCCTTCGAc                                                                      <  6:111931/149‑1 (MQ=255)
                                                                            ccTGCCCTATACCGATGCCCAGTATGCGGTGCTTGAGCAACTCACGAGGCACATCCGGAGTGCCTGGCCAGCCATCGTCCTGAACCGCATTAAGGGCCACAGCGACATTTCGCCGCAGCGCAATACCGACCCCGGCCCGGCCTTCGACt                                                                     <  4:290730/149‑1 (MQ=255)
                                                                             cTGCCCTATACCGATGCCCAGTATGCGGTGCTGGAGCAACTCACCCGGCACATCCGCAGTGCCTGGCCAGCCATCGACCTGAACCGCATTCAGGGCCACAGCGACATTGCCCCGCAGCGCAAGACCGACCCCGGCCCGGCCTTCGACTg                                                                    <  2:46435/149‑1 (MQ=255)
                                                                                  cTATACCGATGCCCAGTATGCGGTGCTGGAGCAACTCACCCGGCACATCCGCAGTGCCTGGCCAGCCATCGACCTGAACCGCATTCAGGGCCACAGCGACATTGCCCCGCAGCGCAAGACCGACCCCGGCCCGGCCTTCGACTGGTcgc                                                               >  7:144058/1‑149 (MQ=255)
                                                                                    ataCCGATGCCCAGTATGCGGTGCTGGAGCAACTCACCCGGCACATCCGCAGTGCCTGGCCAGCCATCGACCTGAACCGCATTCAGGGCCACAGCGACATTGCCCCGCAGCGCAAGACCGACCCCGGCCCGGCCTTCGACTGGTCgcgc                                                             >  1:168312/1‑149 (MQ=255)
                                                                                           tGCCCAGTATGCGGTGCTGGAGCAACTCACCCGGCACATCCGCAGTGCCTGGCCAGCCATCGACCTGAACCGCATTCAGGGCCACAGCGACATTGCCCCGCAGCGCAAGACCGACCCCGGCCCGGCCTTCGACTGGTcgc                                                               >  6:107674/1‑140 (MQ=255)
                                                                                           tGCCCAGTATGCGGTGCTGGAGCAACTCACCCGGCACATCCGCAGTGCCTGGCCAGCCATCGACCTGAACCGCATTCAGGGCCACAGCGACATTGCCCCGCAGCGCAAGACCGACCCCGGCCCGGCCTTCGACTGGTcgc                                                               <  5:107674/140‑1 (MQ=255)
                                                                                            gCCCAGTATGCGGTGCTGGAGCAACTCACCCGGCACATCCGCAGTGCCTGGCCAGCCATCGACCTGAACCGCATTCAGGGCCACAGCGACATTGCCCCGCAGCGCAAGACCGACCCCGGCCCGGCCTTCGACTGGTCGCGCTAc                                                          >  3:43726/1‑144 (MQ=255)
                                                                                            gCCCAGTATGCGGTGCTGGAGCAACTCACCCGGCACATCCGCAGTGCCTGGCCAGCCATCGACCTGAACCGCATTCAGGGCCACAGCGACATTGCCCCGCAGCGCAAGACCGACCCCGGCCCGGCCTTCGACTGGTCGCGCTAc                                                          <  4:43726/144‑1 (MQ=255)
                                                                                                    tgctgtgctgGAGCAACTCACCCGGCACATCCGCAGTGCCTGGCTAGCCATCGACCTGAACCGCATTCAGGGCCACAGCGACATTGCCCCGCAGCGCAAGACCGACCCCGGCCCGGCCTTCGACTGGTCGCGCTACCGCCAGGCGCTg                                              <  7:313509/148‑1 (MQ=255)
                                                                                                    tGCGGTGCTGGAGCAACTCACCCGGCACATCCGCAGTGCCTGGCCAGCCATCGACCTGAACCGCATTCAGGGCCACAGCGACATTGCCCCGCAGCGCAAGACCGACCCCGGCCCGGCCTTCGACTGGTCGCGCTACCGCCAGGCGCTg                                              <  8:144058/148‑1 (MQ=255)
                                                                                                       ggTGCTGGAGCAACTCACCCGGCACATCCGCAGTGCCTGGCCAACCATCGACCTGAACCGCATTCAGGGCCACAGCGACATTGCCCCGCAGCGCAAGACCGACCCCGGCCCGGCCTTCGACTGGTCGCGCTACCg                                                        >  3:118367/1‑135 (MQ=255)
                                                                                                       ggTGCTGGAGCAACTCACCCGGCACATCCGCAGTGCCTGGCCAACCATCGACCTGAACCGCATTCAGGGCCACAGCGACATTGCCCCGCAGCGCAAGACCGACCCCGGCCCGGCCTTCGACTGGTCGCGCTACCg                                                        <  4:118367/135‑1 (MQ=255)
                                                                                                            tGGAGAAACTCACCCGGCACATCCGCAGTGACTGGCCAGACATCGACCTGAACAGCATTAAGGGCCACAGCGACATTGACCAGAAGCGCAAGACAGACCCAGGACCGGCCTTCGACTGGTCGCGCTAACGCCAGGCGCTGCAGCACAAc                                     >  4:285811/1‑149 (MQ=255)
                                                                                                             ggAGCAACTCACCCGGCACATCCGCAGTGCCTGGCCAGCCATCGACCTGAACCGCATTCAGGGCCACAGCGACATTGCCCCGCAGCGCAAGACCGACCCCGGCCCGGCCTTCGACTGGTCGCGCTACCGCCAGGCGCTGCAGCACAAc                                     >  3:295171/1‑148 (MQ=255)
                                                                                                              gAGCAACTCACCCGGCACATCCGCAGTGCCTGGCCAGCCATCGACCTGAACCGCATTCAGGGCCACAGCGACATTGCCCCGCAGCGCAAGACCGACCCCGGCCCGGCCTTCGACTGGTCGCGCTACCGCCAGGCGCTGCAGCACAAc                                     >  7:97894/1‑147 (MQ=255)
                                                                                                                   aCTCACCCGGCACATCCGCAGTGCCTGGCCAGCCATCGACCTGAACCGCATTCAGGGCCACAGCGACATTGCCCCGCAGCGCAAGACCGACCCCGGCCCGGCCTTCGACTGGTCGCGCTACCGCCAgg                                                   <  3:138271/128‑1 (MQ=255)
                                                                                                                   aCTCACCCGGCACATCCGCAGTGCCTGGCCAGCCATCGACCTGAACCGCATTCAGGGCCACAGCGACATTGCCCCGCAGCGCAAGACCGACCCCGGCCCGGCCTTCGACTGGTCGCGCTACCGCCAgg                                                   >  4:138271/1‑128 (MQ=255)
                                                                                                                   aCTCACCCGGCACATCCGCAGTGCCTGGCCAGCCATCGACCTGAACCGCATTCAGGGCCACAGCGACATTGCCCCGCAGCGCAAGACCGACCCCGGCCCGGCCTTCGACTGGTCGCGCTACCGCCAGGCGCTGCAGCACAACGAGGAc                               >  3:300188/1‑148 (MQ=255)
                                                                                                                   aCTCACCCGGCACATCCGCAGTGCCTGGCCAGCCATCGACCTGAACCGCATTCAGGGCCACAGCGACATTGCCCCGCAGCGCAAGACCGACCCCGGCCCGGCCTTCGACTGGTCGCGCTACCGCCAGGCGCTGCAGCACAACGAGGAc                               >  6:178006/1‑148 (MQ=255)
                                                                                                                     tCACCCGGCACATCCGCAGTGCCTGGCCAGCCATCGACCTGAACCGCATTCAGGGCCACAGCGACATTGCCCCGCAGCGCAAGACCGACCCCGGCCCGGCCTTCGACTGGTCGCGCTACCGCCAGGCGCTGCAGCACAACGAGGACAAg                            <  3:129193/149‑1 (MQ=255)
                                                                                                                       aCCCGGCACATCCGCAGTGCCTGGCCAGCCNTCGACCTGAACCGCATTCAGGGCCACAGCGACATTGCCCCGCAGCGCAAGACCGACCCCGGCCCGGCCTTCGACTGGTCGCGCTACCGCCAGGCGCTGCAGCACAACGAGGACAAgg                           >  1:29871/1‑148 (MQ=255)
                                                                                                                         ccGGCACATCCGCAGTGCCTGGCCAGCCATCGACCTGAACCGCATTCAGGGCCACAGCGACATTGCCCCGCAGCGCAAGACCGACCCCGGCCCGGCCTTCGACTGGTCGCGCTACCGCCAGGCGCTGCAGCACAACGAGGACAAGGCAt                        <  2:168312/149‑1 (MQ=255)
                                                                                                                         ccGGCACATCCGCAGTGCCTGGCCAGCCATCGACCTGAACCGCATTCAGGGCCACAGCGACATTGCCCCGCAGCGCAAGACCGACCCCGGCCCGGCCTTCGACTGGTCGCGCTACCGCCAGGCGCTGCAGCACAACGAGGACAAGGCAt                        <  7:164938/149‑1 (MQ=255)
                                                                                                                           ggCACATCCGCAGTGCCTGGCCAGCCATCGACCTGAACCGCATTCAGGGCCACAGCGACATTGCCCCGCAGCGCAAGACCGACCCCGGCCCGGCCTTCGACTGGTCGCGCTACCGCCAGGCGCTGCAGCACAACGAGGACAAGGCATGa                      <  2:3799/149‑1 (MQ=255)
                                                                                                                                 tCCGCAGTGCCTGGCCAGCCATCGACCTGAACCGCATTCAGGGCCACAGCGACATTGCCCCGCAGCGCAAGACCGACCCCGGCCCGGCCTTCGACTGGTCGCGCTACCGCCAGGCGCTGCAGCACAACGAGGACAAGGCATGAGTTTTc                >  7:331772/1‑149 (MQ=255)
                                                                                                                                    gCAGTGCCTGGCCAGCCATCGACCTGAACCGCATTCAGGGCCACAGCGACATTGCCCCGCAGCGCAAGACCGACCCCGGCCCGGCCTTCGACTGGTCGCGCTACCGCCAGGCGCTGCAGCACAACGAGGACAAGGCATGAGTTTTCTgg             >  1:218999/1‑149 (MQ=255)
                                                                                                                                    gCAGTGCCTGGCCAGCCATCGACCTGAACCGCATTCAGGGCCACAGCGACATTGCCCCGCAGCGCAAGACCGACCCCGGCCCGGCCTTCGACTGGTCGCGCTACCGCCAGGCGCTGCAGCACAACGAGGACAAGGCATGAGTTTTCTgg             >  1:240594/1‑149 (MQ=255)
                                                                                                                                               ccagccaTCGACCTGAACCGCATTCAGGGCCACAGCGACATTGCCCCGCAGCGCAAGACCGACCCCGGCCCGGCCTTCGACTGGTCGCGCTACCGCCAGGCGCTGCAGCACAACGAGGACAAGGCATGAGTTTTCTGGTGTTGTTGCTg  >  4:272269/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CGGGTGTGTCGAGCTTTGACGGGCGCGAAGGCTGTAACGACTTTTCCATCGGCATCGAGCTGGAAGGCACCGACGACCTGCCCTATACCGATGCCCAGTATGCGGTGCTGGAGCAACTCACCCGGCACATCCGCAGTGCCTGGCCAGCCATCGACCTGAACCGCATTCAGGGCCACAGCGACATTGCCCCGCAGCGCAAGACCGACCCCGGCCCGGCCTTCGACTGGTCGCGCTACCGCCAGGCGCTGCAGCACAACGAGGACAAGGCATGAGTTTTCTGGTGTTGTTGCTG  >  NC_002947/906291‑906582

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