Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 648,692 C→T 100% R69H (CGC→CAC)  mrdB ← cell wall shape‑determining protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL606648,6920C→T100.0% 405.4 143R69H (CGC→CAC) mrdBcell wall shape‑determining protein
Reads supporting (aligned to +/- strand):  new base (98/45):  ref base (0/0):  total (98/45)

Only 100 of 143 total aligned reads displayed.

TATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGAT  >  REL606/648660‑648723
                                |                               
tataGAGATAGGGGGCCCAGCCTTCATAAACGTGt                               <  1:17323571/35‑1 (MQ=255)
tataGAGATAGGGGGCCCAGCCTTCATAAACGTGt                               >  1:11839953/1‑35 (MQ=255)
tataGAGATAGGGGGCCCAGCCTTCATAAACGTGt                               <  1:11411610/35‑1 (MQ=255)
tataGAGATAGGGGGCCCAGCCTTCATAAACGTGt                               <  1:5443889/35‑1 (MQ=255)
tataGAGATAGGGGGCCCAGCCTTCATAAACGTGt                               >  1:387541/1‑35 (MQ=255)
 ataGAGATAGGGGGCCCAGCCTTCATAAACGTGTg                              >  1:6961702/1‑35 (MQ=255)
 ataGAGATAGGGGGCCCAGCCTTCATAAACGTGTg                              <  1:10295550/35‑1 (MQ=255)
 ataGAGATAGGGGGCCCAGCCTTCATAAACGTGTg                              >  1:4246071/1‑35 (MQ=255)
 ataGAGATAGGGGGCCCAGCCTTCATAAACGTGTg                              >  1:8523110/1‑35 (MQ=255)
  taGAGATAGGGGGCCCAGCCTTCATAAACGTGTgg                             >  1:17471525/1‑35 (MQ=255)
  taGAGATAGGGGGCCCAGCCTTCATAAACGTGTgg                             >  1:5412007/1‑35 (MQ=255)
  taGAGATAGGGGGCCCAGCCTTCATAAACGTGTgg                             >  1:10101998/1‑35 (MQ=255)
    gagaTAGGGGGCCCAGCCTTCATAAACGTGTggag                           >  1:5156557/1‑35 (MQ=255)
    gagaTAGGGGGCCCAGCCTTCATAAACGTGTggag                           >  1:10815662/1‑35 (MQ=255)
    gagaTAGGGGGCCCAGCCTTCATAAACGTGTggag                           >  1:17091373/1‑35 (MQ=255)
     agaTAGGGGGCCCAGCCTTCATAAACGTGTggagg                          <  1:9264225/35‑1 (MQ=255)
     agaTAGGGGGCCCAGCCTTCATAAACGTGTggagg                          >  1:5310412/1‑35 (MQ=255)
     agaTAGGGGGCCCAGCCTTCATAAACGTGTggagg                          >  1:14779614/1‑35 (MQ=255)
      gaTAGGGGGCCCAGCCTTCATAAACGTGTggagga                         >  1:1399648/1‑35 (MQ=255)
      gaTAGGGGGCCCAGCCTTCATAAACGTGTggagga                         >  1:4694044/1‑35 (MQ=255)
      gaTAGGGGGCCCAGCCTTCATAAACGTGTggagga                         <  1:10360297/35‑1 (MQ=255)
       aTAGGGGGCCCAGCCTTCATAAACGTGTGGAGGaa                        >  1:8474427/1‑35 (MQ=255)
       aTAGGGGGCCCAGCCTTCATAAACGTGTGGAGGaa                        <  1:685809/35‑1 (MQ=255)
       aTAGGGGGCCCAGCCTTCATAAACGTGTGGAGGaa                        >  1:10625865/1‑35 (MQ=255)
       aTAGGGGGCCCAGCCTTCATAAACGTGTGGAGGaa                        <  1:11412834/35‑1 (MQ=255)
        tAGGGGGCCCAGCCTTCATAAACGTGTGGAGTAAt                       >  1:15612182/1‑35 (MQ=38)
        tAGGGGGCCCAGCCTTCATAAACGTGTGGAGTAAt                       >  1:20614053/1‑35 (MQ=38)
        tAGGGGGCCCAGCCTTCATAAACGTGTGGAGGAAt                       >  1:14290039/1‑35 (MQ=255)
        tAGGGGGCCCAGCCTTCATAAACGTGTGGAGGAAt                       >  1:8505700/1‑35 (MQ=255)
        tAGGGGGCCCAGCCTTCATAAACGTGTGGAGGAAt                       >  1:18534264/1‑35 (MQ=255)
        tAGGGGGCCCAGCCTTCATAAACGTGTGGAGGAAt                       >  1:2820839/1‑35 (MQ=255)
        tAGGGGGCCCAGCCTTCATAAACGTGTGGAGGAAt                       <  1:2943721/35‑1 (MQ=255)
        tAGGGGGCCCAGACTTCATAAACGTGTGGAGGAAt                       <  1:7674063/35‑1 (MQ=38)
         aGGGGGCCCAGCCTTCATAAACGTGTGGAGGAAtt                      <  1:17706035/35‑1 (MQ=255)
         aGGGGGCCCAGCCTTCATAAACGTGTGGAGGAAtt                      <  1:12694529/35‑1 (MQ=255)
         aGGGGGCCCAGCATTCATAAACGTGTGGAGGAAtt                      <  1:2121905/35‑1 (MQ=38)
          gggggCCCAGCCTTCATAAACGTGTTGAGGAAttt                     >  1:1806725/1‑35 (MQ=38)
          gggggCCCAGCCTTCATAAACGTGTGGAGGAAttt                     >  1:17638853/1‑35 (MQ=255)
          gggggCCCAGCCTTCATAAACGTGTGGAGGAAttt                     >  1:20076615/1‑35 (MQ=255)
          gggggCCCAGCCTTCATAAACGTGTGGAGGAAttt                     <  1:11447387/35‑1 (MQ=255)
          gggggCCCAGCCTTCATAAACGTGTGGAGGAAttt                     <  1:6709415/35‑1 (MQ=255)
          gggggCCCAGCCTTCATAAACGTGTGGAGGAAttt                     <  1:14548818/35‑1 (MQ=255)
          gggggCCCAGCCTTCATAAACGTGTGGAGGAAttt                     >  1:11919805/1‑35 (MQ=255)
           ggggCCCAGCCTTCATAAACGTGTGGAGGAATTTg                    <  1:11648427/35‑1 (MQ=255)
           ggggCCCAGCCTTCATAAACGTGTGGAGGAATTTg                    <  1:14190023/35‑1 (MQ=255)
           ggggCCCAGCCTTCATAAACGTGTGGAGGAATTTg                    >  1:12990226/1‑35 (MQ=255)
            gggCCCAGCCTTCATAAACGTGTGGAGGAATTTgc                   >  1:563216/1‑35 (MQ=255)
             ggCCCAGCCTTCATAAACGTGTGGAGGAATTTgcg                  >  1:12501271/1‑35 (MQ=255)
             ggCCCAGCCTTCATAAACGTGTGGAGGAATTTgcg                  >  1:20525927/1‑35 (MQ=255)
             ggCCCAGCCTTCATAAACGTGTGGAGGAATTTgcg                  <  1:13360987/35‑1 (MQ=255)
             ggCCCAGCCTTCATAAACGTGTGGAGGAATTTgcg                  <  1:16636017/35‑1 (MQ=255)
             ggCCCAGCCTTAATAAACGTGTGGAGGAATTTgcg                  <  1:18043791/35‑1 (MQ=37)
              gCCCAGCCTTCATAAACGTGTGGAGGAATTTgcgc                 <  1:5074539/35‑1 (MQ=255)
              gCCCAGCCTTCATAAACGTGTGGAGGAATTTgcgc                 <  1:3414544/35‑1 (MQ=255)
              gCCCAGCCTTCATAAACGTGTGGAGGAATTTgcgc                 >  1:424359/1‑35 (MQ=255)
              gCCCAGCCTTCATAAACGTGTGGAGGAATTTgcgc                 >  1:15518405/1‑35 (MQ=255)
              gCCCAGCCTTCATAAACGTGTGGAGGAATTTgcgc                 <  1:19765092/35‑1 (MQ=255)
              gCCCAGCCTTCATAAACGTGTGGAGGAATTTgcgc                 >  1:19977080/1‑35 (MQ=255)
               cccAGCCTTCATAAACGTGTGGAGGAATTTGCGcc                >  1:12321731/1‑35 (MQ=255)
               cccAGCCTTCATAAACGTGTGGAGGAATTTGCGcc                >  1:9410378/1‑35 (MQ=255)
                ccAGCCTTCATAAACGTGTGGAGGAATTTGCGCCa               >  1:18515053/1‑35 (MQ=255)
                ccAGCCTTCATAAACGTGTGGAGGAATTTGCGCCa               <  1:14056255/35‑1 (MQ=255)
                ccAGCCTTCATAAACGTGTGGAGGAATTTGCGCCa               <  1:2979121/35‑1 (MQ=255)
                 cAGCCTTCATAAACGTGTGGAGGAATTTGCGCCAt              <  1:4110071/35‑1 (MQ=255)
                 cAGCCTTCATAAACGTGTGGAGGAATTTGCGCCAt              <  1:5269807/35‑1 (MQ=255)
                 cAGCCTTCATAAACGTGTGGAGGAATTTGCGCCAt              <  1:7053715/35‑1 (MQ=255)
                 cAGCCTTCATAAACGTGTGGAGGAATTTGCGCCAt              >  1:10075517/1‑35 (MQ=255)
                 cAGCCTTCATAAACGTGTGGAGGAATTTGCGCCAt              >  1:18812846/1‑35 (MQ=255)
                  aGCCTTCATAAACGTGTGGAGGAATTTGCGCCATc             >  1:7811866/1‑35 (MQ=255)
                  aGCCTTCATAAACGTGTGGAGGAATTTGCGCCATc             <  1:14586106/35‑1 (MQ=255)
                  aGCCTTCATAAACGTGTGGAGGAATTTGCGCCATc             >  1:9331181/1‑35 (MQ=255)
                  aGCCTTCATAAACGTGTGGAGGAATTTGCGCCATc             <  1:4790383/35‑1 (MQ=255)
                  aGCCTTCATAAACGTGTGGAGGAATTTGCGCCATc             >  1:17273113/1‑35 (MQ=255)
                   gCCTTCATAAACGTGTGGAGGAATTTGCGCCATca            >  1:10838556/1‑35 (MQ=255)
                   gCCTTCATAAACGTGTGGAGGAATTTGCGCCATca            >  1:15181818/1‑35 (MQ=255)
                    ccTTCATAAACGTGTGGAGTAATTTGCGCCATcac           >  1:6970546/1‑35 (MQ=37)
                    ccTTCATAAACGTGTGGAGGAATTTGCGCCATcac           <  1:19100836/35‑1 (MQ=255)
                    ccTTCATAAACGTGTGGAGGAATTTGCGCCATcac           >  1:9038737/1‑35 (MQ=255)
                    ccTTCATAAACGTGTGGAGGAATTTGCGCCATcac           >  1:11761606/1‑35 (MQ=255)
                    ccTTCATAAACGTGTGGAGGAATTTGCGCCATcac           <  1:12615890/35‑1 (MQ=255)
                     cTTCATAAACGTGTGGAGGAATTTGCGCCATcacc          >  1:13801042/1‑35 (MQ=255)
                      ttCATAAACGTGTGGAGGAATTTGCGCCATcacca         <  1:20565624/35‑1 (MQ=255)
                       tCATAAACGTGTGGAGGAATTTGCGCCATcaccac        >  1:17792669/1‑35 (MQ=255)
                       tCATAAACGTGTGGAGGAATTTGCGCCATcaccac        >  1:13787576/1‑35 (MQ=255)
                        cATAAACGTGTGGAGGAATTTGCGCCATCaccacc       >  1:19684716/1‑35 (MQ=255)
                        cATAAACGTGTGGAGGAATTTGCGCCATCaccacc       >  1:3063676/1‑35 (MQ=255)
                         aTAAACGTGTGGAGGAATTTGCGCCATCAccacca      >  1:9557435/1‑35 (MQ=255)
                         aTAAACGTGTGGAGGAATTTGCGCCATCAccacca      >  1:17668839/1‑35 (MQ=255)
                          tAAACGTGTGGAGGAATTTGCGCCATCACCACCat     >  1:7009512/1‑35 (MQ=255)
                           aaaCGTGTGGAGGAATTTGCGCCATCACCACCatg    >  1:17924219/1‑35 (MQ=255)
                           aaaCGTGTGGAGGAATTTGCGCCATCACCACCatg    >  1:9041512/1‑35 (MQ=255)
                           aaaCGTGTGGAGGAATTTGCGCCATCACCACCatg    >  1:18372072/1‑35 (MQ=255)
                           aaaCGTGTGGAGGAATTTGCGCCATCACCACCata    >  1:4210560/1‑34 (MQ=255)
                            aaCGTGTGGAGGAATTTGCGCCATCACCACCatga   >  1:14397892/1‑35 (MQ=255)
                            aaCGTGTGGAGGAATTTGCGCCATCACCACCatga   >  1:13417338/1‑35 (MQ=255)
                            aaCGTGTGGAGGAATTTGCGCCATCACCACCatga   >  1:1905586/1‑35 (MQ=255)
                             aCGTGTGGAGGAATTTGCGCCATCACCACCatgat  >  1:13831423/1‑35 (MQ=255)
                             aCGTGTGGAGGAATTTGCGCCATCACCACCatgat  >  1:17123615/1‑35 (MQ=255)
                             aCGTGTGGAGGAATTTGCGCCATCACCACCatgat  >  1:17941023/1‑35 (MQ=255)
                             aCGTGTGGAGGAATTTGCGCCATCACCACAatgat  <  1:7066709/35‑1 (MQ=37)
                                |                               
TATAGAGATAGGGGGCCCAGCCTTCATAAACGCGTGGAGGAATTTGCGCCATCACCACCATGAT  >  REL606/648660‑648723

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.