Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 649,391 T→G I471L (ATC→CTC)  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,3910TG98.8% 680.5 / ‑2.8 165I471L (ATC→CTC) mrdAtranspeptidase involved in peptidoglycan synthesis (penicillin‑binding protein 2)
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base G (89/74);  total (91/74)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

Only 100 of 179 total aligned reads displayed.

GCACCTGTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGATCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCGGAATGGTGT  >  REL606/649319‑649462
                                                                        |                                                                       
gCACCTGTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGagcat                                                                      <  1:679626/76‑1 (MQ=255)
gCACCTGTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGagcat                                                                      <  1:1908704/76‑1 (MQ=255)
gCACCTGTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGagcat                                                                      <  1:1826439/76‑1 (MQ=255)
 cACCTGTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGagcatc                                                                     <  1:3228338/76‑1 (MQ=255)
  aCCTGTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAgcatca                                                                    <  2:800871/76‑1 (MQ=255)
  aCCTGTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAgcatca                                                                    <  2:4294055/76‑1 (MQ=255)
   ccTGTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAg                                                                   <  2:4110522/76‑1 (MQ=255)
   ccTGTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAg                                                                   >  1:2220690/1‑76 (MQ=255)
    cTGTTTGCCGTCTTCGGCGGTGCTCATCAGCACATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGt                                                                  >  2:490457/1‑76 (MQ=255)
    cTGTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGt                                                                  >  2:5895936/1‑76 (MQ=255)
     tGTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTg                                                                 <  1:882783/76‑1 (MQ=255)
      gTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGc                                                                >  1:38961/1‑76 (MQ=255)
      gTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGc                                                                >  1:422138/1‑76 (MQ=255)
       tttGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGcc                                                               >  2:3096382/1‑76 (MQ=255)
       tttGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGcc                                                               <  1:2860634/76‑1 (MQ=255)
        ttGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCt                                                              >  2:1493179/1‑76 (MQ=255)
        ttGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCt                                                              <  1:62726/76‑1 (MQ=255)
         tGCCGTCTTCGGCGGTGCTCGTCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCtt                                                             <  2:88842/76‑1 (MQ=255)
         tGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGGACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCtt                                                             >  2:2845586/1‑76 (MQ=255)
         tGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCtt                                                             >  1:2891153/1‑76 (MQ=255)
         tGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCtt                                                             >  1:2785685/1‑76 (MQ=255)
         tGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCtt                                                             >  2:5440246/1‑76 (MQ=255)
         tGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCtt                                                             >  1:3642494/1‑76 (MQ=255)
         tGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCtt                                                             <  2:4655024/76‑1 (MQ=255)
          gCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTa                                                            >  2:260277/1‑76 (MQ=255)
          gCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTa                                                            <  2:1973419/76‑1 (MQ=255)
          gCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTa                                                            >  2:2428245/1‑76 (MQ=255)
          gCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTa                                                            >  2:4730466/1‑76 (MQ=255)
           ccGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTAc                                                           >  2:4333971/1‑76 (MQ=255)
           ccGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTAc                                                           >  2:6055832/1‑76 (MQ=255)
             gTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTc                                                         >  2:4542987/1‑76 (MQ=255)
              tCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCa                                                        <  1:965833/76‑1 (MQ=255)
               cTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCAt                                                       <  2:4344805/76‑1 (MQ=255)
               cTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCAt                                                       >  1:6154128/1‑76 (MQ=255)
               cTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCAt                                                       >  1:132128/1‑76 (MQ=255)
               cTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCAt                                                       <  2:3970852/76‑1 (MQ=255)
                 tCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCt                                                     <  1:2651859/76‑1 (MQ=255)
                 tCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCt                                                     >  2:2547573/1‑76 (MQ=255)
                   ggcggTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTgg                                                   <  1:2519991/76‑1 (MQ=255)
                   ggcggTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTgg                                                   <  1:5635507/76‑1 (MQ=255)
                    gcggTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGa                                                  <  1:3935861/76‑1 (MQ=255)
                    gcggTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGa                                                  <  2:236309/76‑1 (MQ=255)
                      ggTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGAtt                                                <  2:6300084/76‑1 (MQ=255)
                      ggTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGAtt                                                <  2:2740982/76‑1 (MQ=255)
                       gTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTg                                               <  2:1839341/76‑1 (MQ=255)
                         gCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTggg                                             >  1:367159/1‑76 (MQ=255)
                           tcatcaGCAAATGCGGAACCTTCACGATACCTTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGt                                           <  2:2843205/76‑1 (MQ=255)
                           tcatcaGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGt                                           >  2:1752046/1‑76 (MQ=255)
                            catcaGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTc                                          <  2:4993864/76‑1 (MQ=255)
                             atcaGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGATCATCAGTGCCTTACTCATCTGGATTGGGGTCg                                         >  1:4636451/1‑76 (MQ=255)
                             atcaGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCg                                         <  1:1990982/76‑1 (MQ=255)
                              tcaGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGc                                        >  2:3008896/1‑76 (MQ=255)
                                aGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTg                                      >  2:4732716/1‑76 (MQ=255)
                                aGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTg                                      >  1:5284151/1‑76 (MQ=255)
                                  cAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTc                                    >  1:5173101/1‑76 (MQ=255)
                                   aaaTGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTcc                                   <  2:4040124/76‑1 (MQ=255)
                                   aaaTGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTcc                                   >  1:4687899/1‑76 (MQ=255)
                                    aaTGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCa                                  >  2:3473928/1‑76 (MQ=255)
                                    aaTGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCa                                  >  1:3818019/1‑76 (MQ=255)
                                      tGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGt                                <  1:104025/76‑1 (MQ=255)
                                      tGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGt                                >  1:3973100/1‑76 (MQ=255)
                                       gCGGAACCTTCACGATACCGTCATTAATCAGGATCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTa                               >  1:5927281/1‑76 (MQ=255)
                                       gCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGGCGCTGTCCAGTa                               <  2:6227083/76‑1 (MQ=255)
                                        cGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTaa                              <  1:351688/76‑1 (MQ=255)
                                        cGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTaa                              <  1:4522045/76‑1 (MQ=255)
                                        cGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTaa                              <  1:5235957/76‑1 (MQ=255)
                                         ggAGCCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTacc                             <  2:812321/76‑3 (MQ=255)
                                         ggAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAAc                             >  2:1291376/1‑76 (MQ=255)
                                         ggAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAAc                             >  2:3190215/1‑76 (MQ=255)
                                          gAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAAcc                            >  2:5568975/1‑76 (MQ=255)
                                           aaCCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAAccc                           >  1:1010791/1‑76 (MQ=255)
                                            aCCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCt                          <  1:490457/76‑1 (MQ=255)
                                            aCCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCt                          <  1:4847194/76‑1 (MQ=255)
                                            aCCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCt                          <  1:1421557/76‑1 (MQ=255)
                                             ccTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTg                         <  2:5875238/76‑1 (MQ=255)
                                              cTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGa                        <  2:6007341/76‑1 (MQ=255)
                                               ttCACGTTACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGAc                       <  2:235778/76‑1 (MQ=255)
                                               ttCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGAc                       >  1:429164/1‑76 (MQ=255)
                                                 cACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCg                     <  1:5531349/76‑1 (MQ=255)
                                                   cGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGAt                   <  2:4875658/76‑1 (MQ=255)
                                                    gATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATa                  <  1:5899678/76‑1 (MQ=255)
                                                     aTACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATAc                 <  2:6191829/76‑1 (MQ=255)
                                                     aTACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATAc                 >  1:2730286/1‑76 (MQ=255)
                                                     aTACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATAc                 <  2:1051778/76‑1 (MQ=255)
                                                      tACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATAcc                >  1:5068276/1‑76 (MQ=255)
                                                      tACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATAcc                <  2:1324493/76‑1 (MQ=255)
                                                         cGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAAc             <  2:5350609/76‑1 (MQ=255)
                                                          gTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAAcc            >  1:2782629/1‑76 (MQ=255)
                                                           tCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCg           >  2:5199770/1‑76 (MQ=255)
                                                           tCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCg           >  1:1337009/1‑76 (MQ=255)
                                                              ttAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCGGaa        >  1:3307613/1‑76 (MQ=255)
                                                              ttAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCGGaa        <  2:1623076/76‑1 (MQ=255)
                                                               tAATCAGGAGCATCAGTGCCTTCCTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCGGAAt       <  1:1735739/76‑1 (MQ=255)
                                                               tAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCGGAAt       >  1:1221265/1‑76 (MQ=255)
                                                                aaTCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCGGAATg      >  2:2497944/1‑76 (MQ=255)
                                                                 aTCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCGGAATgg     >  2:3738122/1‑76 (MQ=255)
                                                                 aTCAGGAGCATCAGTCCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCGGAATgg     <  1:4320416/76‑1 (MQ=255)
                                                                   cAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCGGAATGgtg   >  1:5935834/1‑76 (MQ=255)
                                                                    aGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCGGAATGgtgt  <  2:5243327/76‑1 (MQ=255)
                                                                    aGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCGGAATGgtgt  >  2:3780172/1‑76 (MQ=255)
                                                                        |                                                                       
GCACCTGTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGATCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCGGAATGGTGT  >  REL606/649319‑649462

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