Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 645,855 A→G L231S (TTA→TCA)  dacA ← D‑alanyl‑D‑alanine carboxypeptidase (penicillin‑binding protein 5)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606645,8550AG100.0% 260.1 / NA 61L231S (TTA→TCA) dacAD‑alanyl‑D‑alanine carboxypeptidase (penicillin‑binding protein 5)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/34);  total (27/34)

AGCAACAAGGTTGTAACCTGCTTTGTCAGTGTGTCCGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATAACAGGCCGTTACGGTTCAGCTGGCGAATACCGTTAAACGTAAATTCTTTTTCTTTATAGATCGAGTATTCA  >  REL606/645782‑645926
                                                                         |                                                                       
aGCAACAAGGTTGTAACCTGCTTTGTCAGTGTGTCCGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGAc                                                                       <  1:1654214/76‑1 (MQ=255)
aGCAACAAGGTTGTAACCTGCTTTGTCAGTGTGTCCGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGAc                                                                       <  2:1399159/76‑1 (MQ=255)
 gCAACAAGGTTGTAACCTGCTTTGTCAGTGTGTCCGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGACa                                                                      <  1:76923/76‑1 (MQ=255)
   aacaaGGTTGTAACCTGCTTTGTCAGTGTGTCCGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGACAgg                                                                    >  1:47691/1‑76 (MQ=255)
   aacaaGGTTGTAACCTGCTTTGTCAGTGTGTCCGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGACAgg                                                                    >  1:1325132/1‑76 (MQ=255)
   aacaaGGTTGTAACCTGCTTTGTCAGTGTGTCCGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGACAgg                                                                    <  1:1774375/76‑1 (MQ=255)
     caaGGTTGTAACCTGCTTTGTCAGTGTGTCCGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGACAGGcc                                                                  <  1:402655/76‑1 (MQ=255)
      aaGGTTGTAACCTGCTTTGTCAGTGTGTCCGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCg                                                                 >  2:813230/1‑76 (MQ=255)
      aaGGTTGTAACCTGCTTTGTCAGTGTGTCCGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCg                                                                 >  2:325242/1‑76 (MQ=255)
       aGGTTGTAACCTGCTTTGTCAGTGTGTCCGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCGt                                                                >  2:1310223/1‑76 (MQ=255)
       aGGGTGTAACCTGCTTTGTCAGTGTGTCCGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCGt                                                                >  2:1310239/1‑76 (MQ=255)
           tGTAACCTGCTTTGTCAGTGTGTCCGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACg                                                            >  1:1103505/1‑76 (MQ=255)
           tGTAACCTGCTTTGTCAGTGTGTCCGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACg                                                            >  1:65928/1‑76 (MQ=255)
             tAACCTGCTTTGTCAGTGTGTCCGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGt                                                          >  2:1267527/1‑76 (MQ=255)
             tAACCTGCTTTGTCAGTGTGTCCGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGt                                                          >  1:772068/1‑76 (MQ=255)
               cgatGCTTTGTCAGTGTGTCCGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTc                                                        >  2:851893/4‑76 (MQ=255)
                  tGCTTTGTCAGTGTGTCCGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGc                                                     >  2:172958/1‑76 (MQ=255)
                   gCTTTGTCAGTGTGTCCGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGCt                                                    >  2:1670277/1‑76 (MQ=255)
                       tGTCAGTGTGTCCGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCg                                                <  1:291596/76‑1 (MQ=255)
                       tGTCAGTGTGTCCGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCg                                                >  1:1605255/1‑76 (MQ=255)
                          cAGTGTGTCCGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCGAAt                                             <  2:85137/76‑1 (MQ=255)
                          cAGTGTGTCCGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCGAAt                                             >  1:684367/1‑76 (MQ=255)
                           aGTGTGTCCGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCGAATa                                            <  2:1342106/76‑1 (MQ=255)
                           aGTGTGTCCGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCGAATa                                            >  2:1105222/1‑76 (MQ=255)
                           aGTGGGGCCGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCGAATa                                            <  2:950829/76‑1 (MQ=255)
                              gtgtCCGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCGAATACCg                                         <  2:1932140/76‑1 (MQ=255)
                                gtCCGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCGAATACCGtt                                       >  1:820919/1‑76 (MQ=255)
                                gtCCGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCGAATACCGtt                                       >  2:1433753/1‑76 (MQ=255)
                                gtCCGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCGAATACCGtt                                       <  1:1469330/76‑1 (MQ=255)
                                 tCCGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCGAATACCGTTa                                      <  2:1299482/76‑1 (MQ=255)
                                 tCCGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCGAATACCGTTa                                      >  1:396443/1‑76 (MQ=255)
                                 tCCGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCGAATACCGTTa                                      <  2:1324608/76‑1 (MQ=255)
                                 tCCGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCGAATACCGTTa                                      <  2:1920761/76‑1 (MQ=255)
                                 tCCGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCGAATACCGTTa                                      <  2:1700805/76‑1 (MQ=255)
                                 tCCGGTTTTGATGCCGTCGACATGCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCGAATACCGTTa                                      <  2:1299466/76‑1 (MQ=255)
                                  ccGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCGAATACCGTTaa                                     <  2:618432/76‑1 (MQ=255)
                                  ccGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCGAATACCGTTaa                                     <  2:517767/76‑1 (MQ=255)
                                  ccGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCGAATACCGTTaa                                     <  2:417879/76‑1 (MQ=255)
                                  ccGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCGAATACCGTTaa                                     <  2:1392701/76‑1 (MQ=255)
                                   cGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCGAATACCGTTaaa                                    <  1:529526/76‑1 (MQ=255)
                                     gTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCGAATACCGTTAAACg                                  >  1:1213453/1‑76 (MQ=255)
                                      ttttGATGCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCGAATACCGTTAAACGt                                 >  2:649071/1‑76 (MQ=255)
                                         tGATGCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCGAATACCGTTAAACGTaaa                              >  1:1519312/1‑76 (MQ=255)
                                         tGATGCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCGAATACCGTTAAACGTaaa                              <  1:294045/76‑1 (MQ=255)
                                            tGCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCGAATACCGTTAAACGTAAAttc                           <  2:47691/76‑1 (MQ=255)
                                            tGCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCGAATACCGTTAAACGTAAAttc                           <  2:1325132/76‑1 (MQ=255)
                                             gCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCGAATACCGTTAAACGTAAAttct                          <  1:1133213/76‑1 (MQ=255)
                                             gCCGTCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCGAATACCGTTAAACGTAAAttct                          <  2:820919/76‑1 (MQ=255)
                                              ccGTCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCGAATACCGTTAAACGTAAAttctt                         >  2:757031/1‑76 (MQ=255)
                                                 tCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCGAATACCGTTAAACGTAAAttcttttt                      <  2:32007/76‑1 (MQ=255)
                                                 tCGACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCGAATACCGTTAAACGTAAAttcttttt                      <  2:785000/76‑1 (MQ=255)
                                                   gACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCGAATACCGTTAAATGTAAAttcttttttt                    >  1:110153/1‑74 (MQ=255)
                                                   gACATTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCGAATACCGTTAAACGCAAAttctttttct                    >  2:133805/1‑76 (MQ=255)
                                                      aTTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCGAATACCGTTAAACGTAAATTCTTTTTCTTta                 >  1:1377841/1‑76 (MQ=255)
                                                      aTTCAGGCTGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCGAATACCGTTAAACGTAAATTCTTTTTCTTta                 <  2:794000/76‑1 (MQ=255)
                                                           ggCTGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCGAATACCGTTAAACGTAAATTCTTTTTCTTTATAGAt            <  2:1666136/76‑1 (MQ=255)
                                                              tGTTATCCCATGACAGGCCGTTACGGTTCAGCTGGCGAATACCGTTAAACGTAAATTCTTTTTCTTTATAGATCGa         >  1:364548/1‑76 (MQ=255)
                                                                ttATCCCATGACAGGCCGTTACGGTTCAGCTGGCGAATACCGTTAAACGTAAATTCTTTTTCTTTATAGATCGAGt       <  1:264427/76‑1 (MQ=255)
                                                                  aTCCCATGACAGGCCGTTACGGTTCAGCTGGCGAATACCGTTAAACGTAAATTCTTTTTCTTTATAGATCGAGTat     <  2:1595787/76‑1 (MQ=255)
                                                                   tCCCATGACAGGCCGTTACGGTTCAGCTGGCGAATACCGTTAAACGTAAATTCTTTTTCTTTATAGATCGAGTatt    <  2:1951950/76‑1 (MQ=255)
                                                                     ccATGACAGGCCGTTACGGTTCAGCTGGCGAATACCGTTAAACGTAAATTCTTTTTCTTTATAGATCGAGTattca  >  2:664831/1‑76 (MQ=255)
                                                                     ccATGACAGGCCGTTACGGTTCAGCTGGCGAATACCGTTAAACGTAAATTCTTTTTCTTTATAGATCGAGTattca  <  1:1433753/76‑1 (MQ=255)
                                                                         |                                                                       
AGCAACAAGGTTGTAACCTGCTTTGTCAGTGTGTCCGGTTTTGATGCCGTCGACATTCAGGCTGTTATCCCATAACAGGCCGTTACGGTTCAGCTGGCGAATACCGTTAAACGTAAATTCTTTTTCTTTATAGATCGAGTATTCA  >  REL606/645782‑645926

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