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,3910TG96.4% 456.5 / 2.2 110I471L (ATC→CTC) mrdAtranspeptidase involved in peptidoglycan synthesis (penicillin‑binding protein 2)
Reads supporting (aligned to +/- strand):  ref base T (4/0);  new base G (57/49);  total (61/49)
Fisher's exact test for biased strand distribution p-value = 1.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.32e-01

Only 100 of 125 total aligned reads displayed.

GCACCTGTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGATCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCGGAATGGTGTC  >  REL606/649319‑649463
                                                                        |                                                                        
gccccTTTTTGCTGTCTTCGGCGGTGCTCATCAGCAAACGCGGAACCTTCACGATACCGTCATTAATCAGGagcat                                                                       <  1:2576069/73‑1 (MQ=255)
gCACCTGTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGagcat                                                                       >  2:3170885/1‑76 (MQ=255)
gCACCTGTTTGCCGTCTGCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGagcat                                                                       <  1:409239/76‑1 (MQ=255)
 cACCTGTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGagcatc                                                                      <  2:4768211/76‑1 (MQ=255)
  cccTGTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAgcatca                                                                     <  1:3387175/75‑1 (MQ=255)
  aCCTGTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAgcatca                                                                     <  2:2565857/76‑1 (MQ=255)
  aCCTGTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAgcatca                                                                     <  1:2868274/76‑1 (MQ=255)
  aCCTGTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAgcatca                                                                     <  1:2473519/76‑1 (MQ=255)
  aCCTGTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAgcatca                                                                     >  2:3973527/1‑76 (MQ=255)
    cTGTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGt                                                                   <  1:100421/76‑1 (MQ=255)
    cTGTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGt                                                                   >  1:916705/1‑76 (MQ=255)
      gTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGc                                                                 >  1:3235590/1‑76 (MQ=255)
      gTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGc                                                                 >  2:4364119/1‑76 (MQ=255)
      gTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGc                                                                 >  1:4593511/1‑76 (MQ=255)
       tttGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGcc                                                                >  1:291389/1‑76 (MQ=255)
         tGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCtt                                                              >  2:2933462/1‑76 (MQ=255)
         tGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCtt                                                              <  2:4466152/76‑1 (MQ=255)
          gCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTa                                                             <  2:3744283/76‑1 (MQ=255)
           ccGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTGACGATACCGTCATTAATCAGGAGCATCAGTGCCTTAc                                                            <  2:322391/76‑1 (MQ=255)
             gTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTc                                                          >  1:4293199/1‑76 (MQ=255)
              tCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCa                                                         <  1:3275600/76‑1 (MQ=255)
               cTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGATCATCAGTGCCTTACTCAt                                                        >  1:1295516/1‑76 (MQ=255)
               cTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCAt                                                        >  2:2448352/1‑76 (MQ=255)
               cTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCAt                                                        >  2:468176/1‑76 (MQ=255)
                ttCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATc                                                       >  1:1254456/1‑76 (MQ=255)
                ttCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATc                                                       >  2:361252/1‑76 (MQ=255)
                 tCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCt                                                      <  2:85781/76‑1 (MQ=255)
                 tCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCt                                                      <  1:3711275/76‑1 (MQ=255)
                   ggcggTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTgg                                                    <  1:3681416/76‑1 (MQ=255)
                   ggcggTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTgg                                                    >  2:3738864/1‑76 (MQ=255)
                    gcggTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGa                                                   >  2:4625243/1‑76 (MQ=255)
                     cggTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGAt                                                  >  2:4512850/1‑76 (MQ=255)
                      ggTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGAtt                                                 <  2:2538870/76‑1 (MQ=255)
                        tGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTgg                                               >  2:2753300/1‑76 (MQ=255)
                         gCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTggg                                              >  1:1008362/1‑76 (MQ=255)
                           tcatcaGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGt                                            <  1:637741/76‑1 (MQ=255)
                            catcaGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTc                                           <  2:4692374/76‑1 (MQ=255)
                            catcaGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTc                                           >  2:4112189/1‑76 (MQ=255)
                            catcaGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTc                                           <  2:13048/76‑1 (MQ=255)
                            catcaGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTc                                           >  2:1754770/1‑76 (MQ=255)
                            catcaGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTc                                           <  2:2592361/76‑1 (MQ=255)
                             atcaGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCg                                          >  2:1677070/1‑76 (MQ=255)
                              acaGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGc                                         <  2:376357/75‑1 (MQ=255)
                               caGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGATCATCAGTGCCTTACTCATCTGGATTGGGGTCGCt                                        >  1:704595/1‑76 (MQ=255)
                               caGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCt                                        <  2:3909470/76‑1 (MQ=255)
                                aGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTg                                       <  1:3418279/76‑1 (MQ=255)
                                  cAAATGCGGAACCTTCACGATACCGTCATTAATCAGGATCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTc                                     >  1:621056/1‑76 (MQ=255)
                                  cAAATGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTc                                     >  2:1698262/1‑76 (MQ=255)
                                   aaaTGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTcc                                    <  2:1173111/76‑1 (MQ=255)
                                   aaaTGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTcc                                    <  2:752225/76‑1 (MQ=255)
                                   aaaTGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTcc                                    >  2:644673/1‑76 (MQ=255)
                                    aaTGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCa                                   >  2:2208282/1‑76 (MQ=255)
                                    aaTGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCa                                   <  2:4009558/76‑1 (MQ=255)
                                      tGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGt                                 >  1:3204954/1‑76 (MQ=255)
                                      tGCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGt                                 <  1:1570378/76‑1 (MQ=255)
                                       gCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTa                                <  1:5013671/76‑1 (MQ=255)
                                       gCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTa                                >  1:1476532/1‑76 (MQ=255)
                                       gCGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTa                                <  2:2478024/76‑1 (MQ=255)
                                        cGGAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTaa                               <  1:1534044/76‑1 (MQ=255)
                                        aggAACCTCCACGATACCGTCATTAATCAGAAGCATCAGTGCCTTACTCATCTGGATTGGGGTCTCTGTCCAGTaa                               <  2:1642918/75‑1 (MQ=255)
                                         ggAACCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAAc                              >  1:1697264/1‑76 (MQ=255)
                                           aaCCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAAccc                            >  1:1055268/1‑76 (MQ=255)
                                            aCCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCt                           <  1:1743696/76‑1 (MQ=255)
                                            aCCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCt                           <  1:4611816/76‑1 (MQ=255)
                                            aCCTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCt                           <  1:4618532/76‑1 (MQ=255)
                                             ccTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTg                          <  1:4901635/76‑1 (MQ=255)
                                             ccTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTg                          <  2:1883742/76‑1 (MQ=255)
                                              cTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGa                         >  1:1322486/1‑76 (MQ=255)
                                              cTTCACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGa                         >  1:4012166/1‑76 (MQ=255)
                                                 cACGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTCACCCTGACCg                      >  2:2882684/1‑76 (MQ=255)
                                                  aCGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGa                     >  1:4808882/1‑76 (MQ=255)
                                                  aCGATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGa                     >  1:3703436/1‑76 (MQ=255)
                                                    gATACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATa                   <  1:1195731/76‑1 (MQ=255)
                                                     aTACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATAc                  >  2:246312/1‑76 (MQ=255)
                                                     aTACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATAc                  >  2:3461755/1‑76 (MQ=255)
                                                     aTACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATAc                  >  1:4005139/1‑76 (MQ=255)
                                                      tACCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATAcc                 >  1:1001117/1‑76 (MQ=255)
                                                       aCCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCa                >  1:1497489/1‑76 (MQ=255)
                                                       aCCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCa                >  1:4327725/1‑76 (MQ=255)
                                                       aCCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCa                <  1:3868422/76‑1 (MQ=255)
                                                       aCCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCa                <  1:3011245/76‑1 (MQ=255)
                                                       aCCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCa                <  1:2358952/76‑1 (MQ=255)
                                                       aCCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCa                <  1:4773787/76‑1 (MQ=255)
                                                       aCCGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCa                <  2:1867920/76‑1 (MQ=255)
                                                       aCCGTCATTAATCAGGAGCAGCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCa                <  1:4401289/76‑1 (MQ=255)
                                                        ccGTCATTAATCAGGATCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCaa               >  1:836963/1‑76 (MQ=255)
                                                         cGTCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAAc              >  1:3995252/1‑76 (MQ=255)
                                                           tCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCg            <  2:2972204/76‑1 (MQ=255)
                                                           tCATTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCg            >  1:2445147/1‑76 (MQ=255)
                                                             aTTAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCGGa          >  1:2179974/1‑76 (MQ=255)
                                                              ttAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCGGaa         >  2:1501511/1‑76 (MQ=255)
                                                              ttAATCAGGAGCATCAGTGCCTTACTCATCCGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCGGaa         >  2:797274/1‑76 (MQ=255)
                                                               tAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCGGAAt        <  1:3078600/76‑1 (MQ=255)
                                                               tAATCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCGGAAt        >  2:4748754/1‑76 (MQ=255)
                                                               tAAGCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCGGAAt        <  2:1604433/76‑1 (MQ=255)
                                                                  tCAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCGGAATGgt     >  1:4671906/1‑76 (MQ=255)
                                                                   cAGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCGGAATGgtg    >  2:3095365/1‑76 (MQ=255)
                                                                    aGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCGGAATGgtgt   <  2:1103349/76‑1 (MQ=255)
                                                                    aGGAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCGGAATGgtgt   <  2:3306768/76‑1 (MQ=255)
                                                                     ggAGCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCGGAATGGTGTc  <  1:2240570/76‑1 (MQ=255)
                                                                        |                                                                        
GCACCTGTTTGCCGTCTTCGGCGGTGCTCATCAGCAAATGCGGAACCTTCACGATACCGTCATTAATCAGGATCATCAGTGCCTTACTCATCTGGATTGGGGTCGCTGTCCAGTAACCCTGACCGATACCAACCGGAATGGTGTC  >  REL606/649319‑649463

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: