Predicted mutation
evidence seq id position mutation freq annotation gene description
MC JC REL606 547,700 Δ8,224 bp 100% IS1‑mediated [nmpC][ECB_00513] [nmpC], ybcR, ybcS, ybcT, ybcU, ECB_00510, nohB, ECB_00512, [ECB_00513]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ REL606 546952–547700 555923–555877 8178–8972 24 [20] [1] 87 [insB‑6]–[ECB_00513] [insB‑6], insA‑6, nmpC, ybcR, ybcS, ybcT, ybcU, ECB_00510, nohB, ECB_00512, [ECB_00513]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 555924 =1 (0.020)59 (1.080) 34/62 NT 98.3% coding (1209/2346 nt) ECB_00513 conserved hypothetical protein
?REL606 = 2035093 NA (NA)noncoding (1/768 nt) IS1 repeat region

TGGAGAAACCGAAACCGGTAATGACTCCAACTTATTGATAGTGTTTTA  >  REL606/546916‑546963
                                   |            
tGGAGAAACCGAAACCGGTAATGACTCCAACTTAtt              >  1:2621145/1‑36 (MQ=255)
tGGAGAAACCGAAACCGGTAATGACTCCAACTTAtt              >  1:4436568/1‑36 (MQ=255)
tGGAGAAACCGAAACCGGTAATGACTCCAACTTAtt              <  1:4067630/36‑1 (MQ=255)
 ggAGAAACCGAAACCGGTAATGACTCCAACTTATTg             >  1:1174315/1‑36 (MQ=255)
 ggAGAAACCGAAACCGGTAATGACTCCAACTTATTg             >  1:7678888/1‑36 (MQ=255)
 ggAGAAACCGAAACCGGTAATGACTCCAACTTATTg             <  1:632602/36‑1 (MQ=255)
 ggAGAAACCGAAACCGGGAATGACTCCAACTTATTg             >  1:7759597/1‑36 (MQ=255)
  gagaAACCGAAACNGGTAATNACTCCAACTTATTNa            >  1:613384/1‑36 (MQ=21)
    gaaaccgaaaccGGTAATGACTCCAACTTATTGATa          >  1:5112371/1‑36 (MQ=255)
    gaaaccgaaaccGGTAATGACTCCAACTTATTGATa          >  1:3477781/1‑36 (MQ=255)
    gaaaccgaaaccGGTAATGACTCCAACTTATTGATa          >  1:1939140/1‑36 (MQ=255)
    gaaaccgaaaccGGTAATGACTCCAACTTATTGATa          >  1:1009177/1‑36 (MQ=255)
    gaaaccgaaaccGGGAATGACTCCAACTTATTGATa          >  1:5657225/1‑36 (MQ=255)
     aaaccgaaaccgGTAATGACTCCAACTTATTGATAg         >  1:4320809/1‑36 (MQ=255)
      aaccgaaaccgGTAATGACTCCGCCTTATTGATAgt        <  1:3633769/36‑1 (MQ=18)
      aaccgaaaccgGTAATGACTCCAACTTATTGATAgt        <  1:3799278/36‑1 (MQ=255)
       accgaaaccgGTAATGACTCCAACTTATTGATAgtg       >  1:4085908/1‑36 (MQ=255)
       accgaaaccgGTAATGACTCCAACTTATTGATAgtg       >  1:1829683/1‑36 (MQ=255)
        ccgaaaccgGTAATGACTCCAACTTTTTGATAgtgt      >  1:4463178/1‑36 (MQ=255)
         cgaaaccgGTAATGACTCCAACTTATTGATAGTGtt     >  1:6350492/1‑36 (MQ=255)
         cgaaaccgGTAATGACTCCAACTTATTGATAGTGtt     >  1:8126918/1‑36 (MQ=255)
            aaccgGTAATGACTCCAACTTATTGATAGTGTTTTa  >  1:3299216/1‑36 (MQ=35)
            aaccgGTAATGACTCCAACTTATTGATAGTGTTTTa  >  1:7852363/1‑36 (MQ=35)
            aaccgGTAATGACTCCAACTTATTGATAGGGTTTTa  >  1:7831177/1‑36 (MQ=255)
                                   |            
TGGAGAAACCGAAACCGGTAATGACTCCAACTTATTGATAGTGTTTTA  >  REL606/546916‑546963

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 26 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG