Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 649,574 T→C K410E (AAA→GAA)  mrdA ← transpeptidase involved in peptidoglycan synthesis (penicillin‑binding protein 2)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606649,5740TC99.2% 490.0 / ‑4.2 118K410E (AAA→GAA) mrdAtranspeptidase involved in peptidoglycan synthesis (penicillin‑binding protein 2)
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (64/53);  total (65/53)
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.37e-01

Only 100 of 125 total aligned reads displayed.

CATTCGCGGGTAGGCATGTTGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTC  >  REL606/649501‑649646
                                                                         |                                                                        
cATTCGCGGGTAGGCATGTTGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGGAATGACCATAACCGAATTCAc                                                                        >  1:4824211/1‑76 (MQ=255)
 aTTCGCGGGTAGGCATGTTGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTCAcc                                                                       <  2:217680/76‑1 (MQ=255)
  ttCGCGGGTAGGCATGTTGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTCAccc                                                                      <  2:2699481/76‑1 (MQ=255)
      agggTAGGGTTTTTGCAGGAGCGGTCTTCCGCCAGGGCGATACCGGTGTAATGACCATACCCGAATTCACgcgtcc                                                                  <  2:3250616/75‑1 (MQ=255)
          tAGGCATGTTGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTCACCTATCCATTc                                                              >  2:4418596/1‑76 (MQ=255)
          tAGGCATGTTGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTc                                                              <  1:4242462/76‑1 (MQ=255)
           aGGCATGTTGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCg                                                             <  2:4617214/76‑1 (MQ=255)
           aGGCATGTTGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCg                                                             <  2:2064472/76‑1 (MQ=255)
           aGGCATGTTGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCg                                                             >  2:3998150/1‑76 (MQ=255)
           aGGCATGTTGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCg                                                             >  1:1923748/1‑76 (MQ=255)
            ggCATGTTGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCgg                                                            <  1:2664685/76‑1 (MQ=255)
             gCATGTTGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGga                                                           >  2:3141257/1‑76 (MQ=255)
              cATGTTGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGgag                                                          >  1:4044153/1‑76 (MQ=255)
                tGTTGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGagag                                                        >  1:2364619/1‑76 (MQ=255)
                tGTTGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGagag                                                        <  2:608846/76‑1 (MQ=255)
                tGTTGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGagag                                                        >  2:1185842/1‑76 (MQ=255)
                  ttGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGc                                                      >  1:5031825/1‑76 (MQ=255)
                     ccGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGcgat                                                   >  1:5108168/1‑76 (MQ=255)
                     ccGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGcgat                                                   >  2:4000550/1‑76 (MQ=255)
                       ggAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGcgatcg                                                 >  2:4282594/1‑76 (MQ=255)
                        gAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGcgatcga                                                >  1:2720860/1‑76 (MQ=255)
                        gAACGTTCTTCCGCCAGGTCGATACCAGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGcgatcga                                                >  1:4486901/1‑76 (MQ=255)
                             ttcttcCGCCAGGTGGATACCGGTGTAATGACCATAACGGAATTCACTCATCCATTCGGAGAGGCGATGGATtccc                                           >  2:2083794/1‑76 (MQ=255)
                             ttcttcCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATctcc                                           >  2:4737280/1‑76 (MQ=255)
                             ttcttcCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATcccc                                           >  2:3347077/1‑76 (MQ=255)
                             ttcttcCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATcccc                                           <  2:958435/76‑1 (MQ=255)
                              tcttcCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCa                                          >  1:3481713/1‑76 (MQ=255)
                              tcttcCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCa                                          >  1:1716111/1‑76 (MQ=255)
                              tcttcCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCa                                          >  1:1889727/1‑76 (MQ=255)
                               cttcCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCat                                         <  2:2918684/76‑1 (MQ=255)
                               cttcCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCat                                         >  1:2117385/1‑76 (MQ=255)
                                ttcCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCata                                        <  2:4846873/76‑1 (MQ=255)
                                 tcCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCatat                                       >  1:313002/1‑76 (MQ=255)
                                 tcCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCatat                                       >  1:3675091/1‑76 (MQ=255)
                                 tcCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCatat                                       <  1:3814497/76‑1 (MQ=255)
                                 tcCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCatat                                       <  2:3174524/76‑1 (MQ=255)
                                  ccGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATc                                      <  2:4295540/76‑1 (MQ=255)
                                  ccGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATc                                      >  1:4668282/1‑76 (MQ=255)
                                  ccGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATc                                      <  1:4031398/76‑1 (MQ=255)
                                  ccGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATc                                      >  2:4950770/1‑76 (MQ=255)
                                   cGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCg                                     <  1:3629196/76‑1 (MQ=255)
                                      cAGGTCGATACCGGTGTAATGACCATAACCGAATTCCCCCATCCAATCGGAGAGGCGATCGATCCCCATATCGTAg                                  <  2:3764204/76‑1 (MQ=255)
                                      cAGGTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAg                                  >  2:1546598/1‑76 (MQ=255)
                                      cAGGTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAg                                  >  1:2064204/1‑76 (MQ=255)
                                       aGGTCGATACCGGTGTAATGACCATAACCGAATTCCCCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAgg                                 <  1:1033223/76‑1 (MQ=255)
                                         gTCGATACCGGTGTAATGCCCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGcc                               <  1:603555/76‑1 (MQ=255)
                                         gTCGATACCGGTGTAATGACCATAACCGAATTCGCCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGcc                               <  1:5104492/76‑1 (MQ=255)
                                         gTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGcc                               <  1:2022046/76‑1 (MQ=255)
                                         gTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGcc                               >  1:3748833/1‑76 (MQ=255)
                                         gTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGcc                               >  1:3746530/1‑76 (MQ=255)
                                         gTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGcc                               <  1:3426532/76‑1 (MQ=255)
                                         gTCGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGcc                               >  2:1135395/1‑76 (MQ=255)
                                           cGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCAc                             >  1:2750043/1‑76 (MQ=255)
                                           cGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCAc                             >  2:4196422/1‑76 (MQ=255)
                                           cGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCAc                             >  1:3812896/1‑76 (MQ=255)
                                           cGATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCAc                             >  1:1603518/1‑76 (MQ=255)
                                            gATACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCAcc                            >  1:2120258/1‑76 (MQ=255)
                                              tACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTg                          <  2:4421504/76‑1 (MQ=255)
                                              tACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTg                          <  1:340793/76‑1 (MQ=255)
                                              tACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTg                          <  1:3930135/76‑1 (MQ=255)
                                              tACCGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTg                          <  1:4460197/76‑1 (MQ=255)
                                                 cGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATa                       >  1:474794/1‑76 (MQ=255)
                                                 cGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATa                       <  2:4925391/76‑1 (MQ=255)
                                                 cGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATa                       <  1:1871376/76‑1 (MQ=255)
                                                 cGGTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATa                       >  2:988881/1‑76 (MQ=255)
                                                  ggTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATag                      <  1:3170555/76‑1 (MQ=255)
                                                  ggTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATag                      <  2:3365507/76‑1 (MQ=255)
                                                  ggTGTAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATag                      <  1:1233233/76‑1 (MQ=255)
                                                   gtgtAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATaga                     <  2:792516/76‑1 (MQ=255)
                                                    tgtAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATagaa                    >  2:3871185/1‑76 (MQ=255)
                                                     gtAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATagaag                   >  1:4463618/1‑76 (MQ=255)
                                                     gtAATGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATagaag                   <  1:1546598/76‑1 (MQ=255)
                                                       aaTGACCATAACCGATTTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCTCCTGATAgaagaa                 <  2:2999250/76‑1 (MQ=255)
                                                       aaTGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATAgaagaa                 <  1:5133385/76‑1 (MQ=255)
                                                       aaTGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATAgaagaa                 <  1:137718/76‑1 (MQ=255)
                                                       aaTGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATAgaagaa                 <  2:1761439/76‑1 (MQ=255)
                                                       aaTGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATAgaagaa                 >  2:1634980/1‑76 (MQ=255)
                                                       aaTGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATAgaagaa                 >  2:1434715/1‑76 (MQ=255)
                                                        aTGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATAGaagaag                >  2:2353873/1‑76 (MQ=255)
                                                        aTGACAATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATAGaagaag                >  2:962772/1‑76 (MQ=255)
                                                         tGACCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAgg               <  2:3865420/76‑1 (MQ=255)
                                                           aCCATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAGGTa             >  1:1591439/1‑76 (MQ=255)
                                                            ccATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAGGTAt            >  1:1048627/1‑76 (MQ=255)
                                                            ccATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAGGTAt            <  1:2802706/76‑1 (MQ=255)
                                                            ccATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAGGTAt            >  1:353433/1‑76 (MQ=255)
                                                             cATAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATc           >  1:1711497/1‑76 (MQ=255)
                                                              aTAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATcc          >  2:212302/1‑76 (MQ=255)
                                                              aTAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATcc          >  1:1526416/1‑76 (MQ=255)
                                                               tAACCGAATTCCCCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCg         <  2:1896160/76‑1 (MQ=255)
                                                               tAACCGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCg         >  1:436896/1‑76 (MQ=255)
                                                                  ccGAATTCCCCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAg      <  2:4424918/76‑1 (MQ=255)
                                                                  ccGAATTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAg      >  2:2898904/1‑76 (MQ=255)
                                                                    gAATTCCCCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGAt    <  2:3748833/76‑1 (MQ=255)
                                                                    gAATTCACCCATCCATGGGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGAt    <  2:2187493/76‑1 (MQ=255)
                                                                     aaTTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGAtt   >  2:4830629/1‑76 (MQ=255)
                                                                     aaTTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGAtt   >  1:1969153/1‑76 (MQ=255)
                                                                     aaTTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGAtt   <  2:5166638/76‑1 (MQ=255)
                                                                      aTTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGAttc  <  1:922023/76‑1 (MQ=255)
                                                                      aTTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGAttc  <  1:2403258/76‑1 (MQ=255)
                                                                      aTTCACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGAttc  >  2:2334579/1‑76 (MQ=255)
                                                                         |                                                                        
CATTCGCGGGTAGGCATGTTGCCGGAACGTTCTTCCGCCAGGTCGATACCGGTGTAATGACCATAACCGAATTTACCCATCCATTCGGAGAGGCGATCGATCCCCATATCGTAGGCCACCTGATAGAAGAAGGTATCCGCAGATTC  >  REL606/649501‑649646

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: