Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,326,566 A→T D67E (GAT→GAA mreD ← cell wall structural complex MreBCD transmembrane component MreD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,326,5660AT98.7% 619.2 / ‑0.0 152D67E (GAT→GAAmreDcell wall structural complex MreBCD transmembrane component MreD
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base T (101/49);  total (102/50)
Fisher's exact test for biased strand distribution p-value = 5.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

Only 100 of 153 total aligned reads displayed.

ATTTCAGCGCCACCAGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGATCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAAC  >  REL606/3326493‑3326638
                                                                         |                                                                        
aTTTCAGCGCCACCAGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTc                                                                        >  2:4384868/1‑76 (MQ=255)
 tttCAGCGCCACCAGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTcc                                                                       >  2:4221630/1‑76 (MQ=255)
  ttCAGCGCCACCAGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCa                                                                      >  1:3538800/1‑76 (MQ=255)
   tcggcgcCACCAGGAACGCAATGAGGCACATCGCAAATCCGGGTCCCCAAGGCGTGGAGCCGCTGATCAGTTCCAg                                                                     <  2:2797052/73‑1 (MQ=255)
   tCAGCGCCACCAGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAg                                                                     >  1:3864669/1‑76 (MQ=255)
   tCAGCGCCACCAGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAg                                                                     >  2:6450068/1‑76 (MQ=255)
    cAGCGCCACCAGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGt                                                                    >  2:2343803/1‑76 (MQ=255)
     aGCGCCACCAGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTa                                                                   >  1:5056147/1‑76 (MQ=255)
     aGCGCCACCAGGAAAGAAAGAATGCTCATCGCCAATGCGCGTACGCCAAGCGGGGAGCCGCTGATCAGTTCCAGTa                                                                   <  2:5441482/76‑1 (MQ=255)
       ccccaccaGGTAAGCAATGATGCTCACCGCCAAGCCGCGCACGCCAAGCGGCGAGCCGCTGATCAGTTCCAGTATg                                                                 <  1:3214774/74‑1 (MQ=255)
        gcCACCAGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATgg                                                                >  1:4177636/1‑76 (MQ=255)
        cccacaaGGTAAGCAAGGATGCCCACCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATgg                                                                <  1:1578188/75‑1 (MQ=255)
         ccaccaGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGc                                                               >  2:4425415/1‑76 (MQ=255)
          cacgaGGAAAGAAATGATCCTCACCGCAAATACGGGTACCCAAGGCGGGGAGCCGCTGATCAGTTCCAGTATGGCa                                                              <  2:4512353/76‑1 (MQ=255)
          caccaGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCa                                                              >  2:4707386/1‑76 (MQ=255)
          caccaGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCa                                                              >  1:5673387/1‑76 (MQ=255)
          caccaGGAAAGCAATGTTGCCCATCGCCAATCCGCGTACCCAAAGCGGGGAGCCGCTGATCAGTTCCAGTATGGCa                                                              <  2:857731/76‑1 (MQ=255)
           accaGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCAc                                                             >  2:2434254/1‑76 (MQ=255)
           accaGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCAc                                                             >  1:6200149/1‑76 (MQ=255)
            ccaGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCAcc                                                            >  2:6941105/1‑76 (MQ=255)
              aGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCa                                                          >  1:5013436/1‑76 (MQ=255)
                gTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATc                                                        >  1:1881698/1‑76 (MQ=255)
                 tAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGAACCCATca                                                       >  2:2185357/1‑76 (MQ=255)
                 tAAGAAATGAGCCGCACCCCCAATGCGCGTACGCCAGGCGGCGAGCCGCTGATCAGTTCCAGTATGGCACCCATca                                                       <  1:6899585/76‑1 (MQ=255)
                 gaaGAAATGATCCTCATCGCCAATCCCCGTACGCCAAGCGGGGAGCCGCTGATCAGTTCCAGTATGGCACCCATca                                                       <  2:6904446/75‑1 (MQ=255)
                     cAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACaaa                                                   >  1:4543525/1‑76 (MQ=255)
                      atggatgCTCAGCGCCAATCCGCGTCCGCCAGGCGTGGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACaaaa                                                  <  2:6255459/73‑1 (MQ=255)
                      aaTGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACaaaa                                                  >  2:1815449/1‑76 (MQ=255)
                       atgatgCTCATCGCCAATACGCGTACGCCAAGCGTCGGGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAAc                                                 >  2:6530002/1‑76 (MQ=255)
                       atgatgCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAAc                                                 >  1:6036042/1‑76 (MQ=255)
                        tgatgCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACGCATCACAAAAcc                                                >  1:1745030/1‑76 (MQ=255)
                        tgatgCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAAcc                                                >  1:3867841/1‑76 (MQ=255)
                            gCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCtgtg                                            >  2:2340783/1‑76 (MQ=255)
                             cTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGc                                           >  1:4928700/1‑76 (MQ=255)
                              tCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGcc                                          >  2:5295158/1‑76 (MQ=255)
                                aTCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCca                                        >  2:1726224/1‑76 (MQ=255)
                                aTCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCca                                        >  1:4128792/1‑76 (MQ=255)
                                 tCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCcac                                       >  2:3590340/1‑76 (MQ=255)
                                 tCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCcac                                       >  2:623419/1‑76 (MQ=255)
                                  cGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCcaca                                      >  1:2527186/1‑76 (MQ=255)
                                  cGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCcaca                                      >  2:2999454/1‑76 (MQ=255)
                                   ccaaaTCCGCGTACCCAAGGGGGCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACAt                                     <  2:4250300/73‑1 (MQ=255)
                                    caaaTGCGCGTCCCCAAAGCGGCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACAtt                                    <  2:3956066/74‑1 (MQ=255)
                                     cAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACAttt                                   >  2:6559136/1‑76 (MQ=255)
                                     cAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACAttt                                   >  2:2055312/1‑76 (MQ=255)
                                     cAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACAttt                                   >  2:2424531/1‑76 (MQ=255)
                                     cAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGATCCAGTATGGCACCCATCACAAAACCTGTGCCCACAttt                                   >  1:4862987/1‑76 (MQ=255)
                                      aaTCCGCGGACGCCAGGGGTGGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTa                                  <  1:3372694/76‑1 (MQ=255)
                                      aaTCCGCGCACGCAAGGGGTGGAGCCGCTGATCAGATCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTa                                  <  1:3540784/76‑1 (MQ=255)
                                      aaTACGCGTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTa                                  >  2:578616/1‑76 (MQ=255)
                                          cgcgTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTAcgcg                              >  2:6942559/1‑76 (MQ=255)
                                          cgcgTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTAcgcg                              >  2:2995297/1‑76 (MQ=255)
                                           gcgTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGa                             >  2:1237633/1‑76 (MQ=255)
                                           gcgTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGCATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGa                             >  1:742285/1‑76 (MQ=255)
                                             gTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATg                           >  2:55982/1‑76 (MQ=255)
                                             gTACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATg                           >  2:102439/1‑76 (MQ=255)
                                              tACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGa                          >  2:4606898/1‑76 (MQ=255)
                                              tACGCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGa                          >  1:6595556/1‑76 (MQ=255)
                                                 gCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGc                       >  1:1008598/1‑76 (MQ=255)
                                                 gCCAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGc                       >  1:6105901/1‑76 (MQ=255)
                                                  ccAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCa                      >  2:4991549/1‑76 (MQ=255)
                                                   cAAGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCaa                     >  1:3495850/1‑76 (MQ=255)
                                                    aaGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAg                    >  2:2661428/1‑76 (MQ=255)
                                                    aaGCGTCGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAg                    >  1:237431/1‑76 (MQ=255)
                                                      gCGTGGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGc                  <  2:4658567/76‑1 (MQ=255)
                                                        ggggAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCa                <  2:6743943/73‑1 (MQ=255)
                                                        ggggAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCa                <  2:3229007/73‑1 (MQ=255)
                                                        gcggAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCa                <  1:1213936/73‑1 (MQ=255)
                                                         gggAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAg               <  2:318313/74‑1 (MQ=255)
                                                          ggAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAgg              <  1:2138038/75‑1 (MQ=255)
                                                          cGAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAgg              >  2:5302780/1‑76 (MQ=255)
                                                           gAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGa             >  1:1308729/1‑76 (MQ=255)
                                                           gAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGa             >  1:2237889/1‑76 (MQ=255)
                                                           gAGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGa             >  2:6890122/1‑76 (MQ=255)
                                                            aGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGAt            <  1:5154273/76‑1 (MQ=255)
                                                            aGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGAt            <  1:4654358/76‑1 (MQ=255)
                                                            aGCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGAt            <  1:2340908/76‑1 (MQ=255)
                                                             gCCGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATc           <  1:1698556/76‑1 (MQ=255)
                                                              ccGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATcc          <  1:6450068/76‑1 (MQ=255)
                                                              ccGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATcc          >  1:6817648/1‑76 (MQ=255)
                                                               cGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCa         >  2:5731917/1‑76 (MQ=255)
                                                               cGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCa         <  2:5089947/76‑1 (MQ=255)
                                                               cGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCa         <  2:158732/76‑1 (MQ=255)
                                                               cGCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCa         >  2:3965859/1‑76 (MQ=255)
                                                                gCTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCaa        >  2:1390448/1‑76 (MQ=255)
                                                                 cTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAAt       >  2:5981403/1‑76 (MQ=255)
                                                                 cTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAAt       >  2:4635371/1‑76 (MQ=255)
                                                                 cTGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAAt       >  1:3744409/1‑76 (MQ=255)
                                                                  tGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATa      >  2:4788580/1‑76 (MQ=255)
                                                                  tGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATa      >  1:1573432/1‑76 (MQ=255)
                                                                  tGATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATa      <  1:777985/76‑1 (MQ=255)
                                                                   gATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATac     >  1:624395/1‑76 (MQ=255)
                                                                   gATCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGACCCAATac     >  1:6895011/1‑76 (MQ=255)
                                                                    aTCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATaca    <  1:3739052/76‑1 (MQ=255)
                                                                    aTCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATaca    <  1:573065/76‑1 (MQ=255)
                                                                    aTCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATaca    <  1:1772623/76‑1 (MQ=255)
                                                                    aTCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATaca    <  1:6886139/76‑1 (MQ=255)
                                                                     tCAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATacaa   <  1:192531/76‑1 (MQ=255)
                                                                      cAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATacaac  <  2:1035005/76‑1 (MQ=255)
                                                                      cAGTTCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATacaac  <  1:1059633/76‑1 (MQ=255)
                                                                         |                                                                        
ATTTCAGCGCCACCAGGTAAGCAATGATGCTCATCGCCAATACGCGTACGCCAAGCGTCGAGCCGCTGATCAGATCCAGTATGGCACCCATCACAAAACCTGTGCCCACATTTACGCGATGAGGCAAGGCCAGGATCCAATACAAC  >  REL606/3326493‑3326638

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