Predicted mutation
evidence seq id position mutation freq annotation gene description
MC JC REL606 547,700 Δ2,651 bp 100% IS1‑mediated [nmpC][ybcU] [nmpC], ybcR, ybcS, ybcT, [ybcU]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ REL606 546953–547700 550350 2651–3398 55 [48] [7] 94 [insB‑6]–[ybcU] [insB‑6], insA‑6, nmpC, ybcR, ybcS, ybcT, [ybcU]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 547699NA (NA)59 (0.570) 34/48 NT 93.7% noncoding (1/768 nt) IS1 repeat region
?REL606 550351 = 4 (0.040)coding (12/294 nt) ybcU predicted lipoprotein

GAGAAACCGAAACCGGTAATGACTCCAACTTATTGATAGTGTTTTA  >  REL606/546918‑546963
                                  |           
gagaAACCGAAACCGGTAATGACTCCAACTTATTg             <  1:8614514/35‑1 (MQ=255)
gagaAACCGAAACCGGTAATGACTCCAACTTATTg             >  1:7308418/1‑35 (MQ=255)
gagaAACCGAAACCGGTAATGACTCCAACTTATTg             <  1:16229757/35‑1 (MQ=255)
gagaAACCGAAACCGGTAATGACTCCAACTTATTg             <  1:12780664/35‑1 (MQ=255)
gagaAACCGAAACCGGTAATGACTCCAACTTATTg             <  1:19135881/35‑1 (MQ=255)
gagaAACCGAAACCGGTAATGACTCCAACTTATTg             >  1:3614739/1‑35 (MQ=255)
gagaAACCGAAACCGGTAATGACTCCAACTTATTg             >  1:14368138/1‑35 (MQ=255)
 agaAACCGAAACCGGTAATTACTCCAACTTATTGa            >  1:18839010/1‑35 (MQ=255)
 agaAACCGAAACCGGTAATGACTCCAACTTATTGa            >  1:8403788/1‑35 (MQ=255)
 agaAACCGAAACCGGTAATGACTCCAACTTATTGa            <  1:16903870/35‑1 (MQ=255)
 agaAACCGAAACCGGTAATGACTCCAACTTATTGa            >  1:5866492/1‑35 (MQ=255)
 agaAACCGAAACCGGTAATGACTCCAACTTATTGa            >  1:1186195/1‑35 (MQ=255)
 agaAACCGAAACCGGTAATGACTCCAACTTATTGa            <  1:5536993/35‑1 (MQ=255)
 agaAACCGAAACCGGTAATGACTCCAACTTATTGa            <  1:12777638/35‑1 (MQ=255)
  gaaaccgaaaccTGTAATGACTCCAACTTATTGAt           >  1:16850372/1‑35 (MQ=255)
  gaaaccgaaaccGGTAATGACTCCAACTTATTGAt           <  1:20598048/35‑1 (MQ=255)
   aaaccgaaaccgGTAATGACTCCAACTTATTGATa          >  1:3784976/1‑35 (MQ=255)
   aaaccgaaaccgGTAATGACTCCAACTTATTGATa          <  1:5164714/35‑1 (MQ=255)
   aaaccgaaaccgGTAATGACTCCAACTTATTGATa          >  1:18035627/1‑35 (MQ=255)
   aaaccgaaaccgGTAATGACTCCAACTTATTGATa          >  1:7141568/1‑35 (MQ=255)
   aaaccgaaaccgGTAATGACTCCAACTTATTGATa          >  1:7278149/1‑35 (MQ=255)
   aaaccgaaaccgGTAATGACTCCAACTTATTGATa          <  1:9907751/35‑1 (MQ=255)
   aaaccgaaaccgGTAATGACTCCAACTTATTGATa          >  1:16422188/1‑35 (MQ=255)
   aaaccgaaaccgGTAATGACTCCAACTTATTGATa          <  1:10646496/35‑1 (MQ=255)
    aaccgaaaccgGTAATGACTCCAACTTATTGATAg         >  1:13170113/1‑35 (MQ=255)
    aaccgaaaccgGTAATGACTCCAACTTATTGATAg         <  1:15033653/35‑1 (MQ=255)
    aaccgaaaccgGTAATGACTCCAACTTATTGATAg         >  1:19320979/1‑35 (MQ=255)
     accgaaaccgGTAATGACTCCAACTTATTGATAgt        >  1:12144583/1‑35 (MQ=255)
     accgaaaccgGTAATGACTCCAACTTATTGATAgt        >  1:10351298/1‑35 (MQ=255)
     accgaaaccgGTAATGACTCCAACTTATTGATAgt        >  1:19487813/1‑35 (MQ=255)
      ccgaaaccgGTAATGACTCCAACTTATTGATAgtg       >  1:6203846/1‑35 (MQ=255)
      ccgaaaccgGTAATGACTCCAACTTATTGATAgtg       <  1:2444030/35‑1 (MQ=255)
      ccgaaaccgGTAATGACTCCAACTTATTGATAgtg       >  1:13215653/1‑35 (MQ=255)
      ccgaaaccgGTAATGACTCCAACTTATTGATAgtg       <  1:11627884/35‑1 (MQ=255)
       tctaaccgGTAATGACTCCAACTTATTGATAgtgt      <  1:841545/32‑1 (MQ=255)
       cgaaaccgGTAATGACTCCAACTTATTGATAgtgt      <  1:5972027/35‑1 (MQ=255)
       cgaaaccgGTAATGACTCCAACTTATTGATAgtgt      >  1:17112911/1‑35 (MQ=255)
       cgaaaccgGTAATGACTCCAACTTATTGATAgcgt      >  1:8796562/1‑35 (MQ=255)
       cgaaaccgGTAATGACTCCAAATTATTGATAgtgt      <  1:10242962/35‑1 (MQ=255)
        gaaaccgGTAATGACTCCAACTTATTGATAGTGtt     <  1:8485148/35‑1 (MQ=40)
        gaaaccgGTAATGACTCCAACTTATTGATAGTGtt     <  1:10175165/35‑1 (MQ=40)
        gaaaccgGTAATGACTCCAACTTATTGATAGTGtt     <  1:2799318/35‑1 (MQ=40)
        gaaaccgGTAATGACTCCAACTTATTGATAGTGtt     <  1:18946105/35‑1 (MQ=40)
         aaaccgGTAATGACTCCAACTTATTGATAGTGttt    <  1:5813792/35‑1 (MQ=38)
         aaaccgGTAATGACTCCAACTTATTGATAGTGttt    <  1:11462927/35‑1 (MQ=38)
         aaaccgGTAATGACTCCAACTTATTGATAGTGttt    <  1:19879364/35‑1 (MQ=38)
         aaaccgGTAATGACTCCAACTTATTGATAGTGttt    >  1:14968354/1‑35 (MQ=38)
          aaccgGTAATGACTCCAACTTATTGATAGTGtttt   <  1:5911184/35‑1 (MQ=35)
          aaccgGTAATGACTCCAACTTATTGATAGTGtttt   <  1:19751573/35‑1 (MQ=35)
          aaccgGTAATGACTCAAACTTATTGATAGTGtttt   <  1:16650162/35‑1 (MQ=255)
           accgGTAATGACTCCAACTTATTGATAGTGTTTTa  >  1:12765254/1‑35 (MQ=32)
           accgGTAATGACTCCAACTTATTGATAGTGTTTTa  <  1:13415412/35‑1 (MQ=32)
           accgGTAATGACTCCAACTTATTGATAGTGTTTTa  >  1:18627561/1‑35 (MQ=32)
           accgGTAATGACTCCAACTTATTGATAGTGTTTTa  >  1:15280572/1‑35 (MQ=32)
           accgGTAATGACTCCAACTTATTGATAGTGGTTTa  >  1:11444538/1‑35 (MQ=12)
                                  |           
GAGAAACCGAAACCGGTAATGACTCCAACTTATTGATAGTGTTTTA  >  REL606/546918‑546963

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 33 ≤ 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.