Predicted mutation
evidence seq id position mutation annotation gene description
MC JC REL606 1,974,040 Δ8,273 bp IS1‑mediated [yedU][yeeI] [yedU], yedV, yedW, yedX, yedY, yedZ, yodA, yodB, serU, [yeeI]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ REL606 1973304–1974039 1982312 8274–9009 34 [32] [1] 74 insA‑13–[yeeI] insA‑13,insB‑13,yedU,yedV,yedW,yedX,yedY,yedZ,yodA,yodB,serU,[yeeI]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 1974039NA (NA)73 (0.920) 49/100 0.4 98.7% noncoding (768/768 nt) IS1 repeat region
?REL606 1982313 = 1 (0.010)coding (664/798 nt) yeeI hypothetical protein

TGCCAGTGATCCTGCTGAATGTTTTGCCGTACTTTCTGAATATTTCTTTAGCGCCCCAGAACTTTTTGCTCCTCGTTTCCCTTCATTGTGGCAACGTTT  >  REL606/1982264‑1982362
                                                 |                                                 
cgggcattatctgaacataaaacactatcaataagttggagtcattaccAg                                                  <  1:6196500‑M2/2‑1 (MQ=255)
 gggcattatctgaacataaaacactatcaataagttggagtcattaccAgc                                                 >  1:182212‑M2/49‑51 (MQ=255)
 gggcattatctgaacataaaacactatcaataagttggagtcattaccAgc                                                 <  1:1188058‑M2/3‑1 (MQ=255)
  ggcattatctgaacataaaacactatcaataagttggagtcattaccAgcg                                                >  1:233420‑M2/48‑51 (MQ=255)
  ggcattatctgaacataaaacactatcaataagttggagtcattaccAgcg                                                <  1:5028880‑M2/4‑1 (MQ=255)
  ccAGTGATCCTGCTGAATGTTTTGCCGTACTTTCTGAATATTTCTTTAgcg                                                >  1:1189130/1‑51 (MQ=255)
   gcattatctgaacataaaacactatcaataagttggagtcattaccAgcgc                                               >  1:2751678‑M2/47‑51 (MQ=255)
    cattatctgaacataaaacactatcaataagttggagtcattaccAGCGcc                                              >  1:3329234‑M2/46‑51 (MQ=255)
    cattatctgaacataaaacactatcaataagttggagtcattaccAGCGcc                                              >  1:1008426‑M2/46‑51 (MQ=255)
     attatctgaacataaaacactatcaataagttggagtcattaccAGCGccc                                             >  1:4682073‑M2/45‑51 (MQ=255)
     attatctgaacataaaacactatcaataagttggagtcattaccAGCGccc                                             >  1:2804356‑M2/45‑51 (MQ=255)
     attatctgaacataaaacactatcaataagttggagtcattaccAGCGccc                                             <  1:1515630‑M2/7‑1 (MQ=255)
       tatctgaacataaaacactatcaataagttggagtcattaccAGCGCCCCa                                           >  1:2798851‑M2/43‑51 (MQ=255)
        atctgaacataaaacactatcaataagttggagtcattaccAGCGCCCCAg                                          <  1:1452184‑M2/10‑1 (MQ=255)
         tctgaacataaaacactatcaataagttggagtcattaccAGCGCCCCAGa                                         >  1:6217084‑M2/41‑51 (MQ=255)
         tctgaacataaaacactatcaataagttggagtcattaccAGCGCCCCAGa                                         >  1:222290‑M2/41‑51 (MQ=255)
         tctgaacataaaacactatcaataagttggagtcattaccAGCGCCCCAGa                                         >  1:808947‑M2/41‑51 (MQ=255)
          ctgaacataaaacactatcaataagttggagtcattaccAGCGCCCCAGaa                                        <  1:5723445‑M2/12‑1 (MQ=255)
          ctgaacataaaacactatcaataagttggagtcattaccAGCGCCCCAGaa                                        <  1:3245647‑M2/12‑1 (MQ=255)
          ctgaacataaaacactatcaataagttggagtcattaccAGCGCCCCAGaa                                        <  1:5767743‑M2/12‑1 (MQ=255)
          ctgaacataaaacactatcaataagttggagtcattaccAGCGCCCCAGaa                                        >  1:332875‑M2/40‑51 (MQ=255)
           tgaacataaaacactatcaataagttggagtcattaccAGCGCCCCAGAAc                                       >  1:3202435‑M2/39‑51 (MQ=255)
            gaacataaaacactatcaataagttggagtcattaccAGCGCCCCAGAACt                                      >  1:5903996‑M2/38‑51 (MQ=255)
             aacataaaacactatcaataagttggagtcattaccAGCGCCCCAGAACtt                                     >  1:4822861‑M2/37‑51 (MQ=255)
              acataaaacactatcaataagttggagtcattaccAGCGCCCCAGAACttt                                    >  1:4748143‑M2/36‑51 (MQ=255)
              acataaaacactatcaataagttggagtcattaccAGCGCCCCAGAACttt                                    >  1:2374903‑M2/36‑51 (MQ=255)
               cataaaacactatcaataagttggagtcattaccAGCGCCCCAGAACtttt                                   <  1:1147173‑M2/17‑1 (MQ=255)
               cataaaacactatcaataagttggagtcattaccAGCGCCCCAGAACtttt                                   >  1:6019114‑M2/35‑51 (MQ=255)
               cataaaacactatcaataagttggagtcattaccAGCGCCCCAGAACtttt                                   <  1:2172790‑M2/17‑1 (MQ=255)
                ataaaacactatcaataagttggagtcattaccAGCGCCCCAGAACttttt                                  >  1:843407‑M2/34‑51 (MQ=255)
                ataaaacactatcaataagttggagtcattaccAGCGCCCCAGAACttttt                                  >  1:83351‑M2/34‑51 (MQ=255)
                ataaaacactatcaataagttggagtcattaccAGCGCCCCAGAACttttt                                  <  1:1250196‑M2/18‑1 (MQ=255)
                  aaaacactatcaataagttggagtcattaccAGCGCCCCAGAACTTTTTGc                                >  1:4068990‑M2/32‑51 (MQ=255)
                        ctatcaataagttggagtcattaccAGCGCCCCAGAACTTTTTGCTCCTCg                          >  1:1882364‑M2/26‑51 (MQ=255)
                         tntcaataagttggagtcattaccAGCGCCCCAGAACTTTTTGCTCCTCGt                         >  1:4817775‑M2/25‑51 (MQ=255)
                         tatcaataagttggagtcattaccAGCGCCCCAGAACTTTTTGCTCCTCGt                         >  1:3149665‑M2/25‑51 (MQ=255)
                          atcaataagttggagtcattaccAGCGCCCCAGAACTTTTTGCTCCTCGtt                        <  1:3773929‑M2/28‑1 (MQ=255)
                            caataagttggagtcattaccAGCGCCCCAGAACTTTTTGCTCCTCGTTTc                      >  1:5771546‑M2/22‑51 (MQ=255)
                            caataagttggagtcattaccAGCGCCCCAGAACTTTTTGCTCCTCGTTTc                      >  1:2423343‑M2/22‑51 (MQ=255)
                             aataagttggagtcattaccAGCGCCCCAGAACTTTTTGCTCCTCGTTTcc                     >  1:3963911‑M2/21‑51 (MQ=255)
                             aataagttggagtcattaccAGCGCCCCAGAACTTTTTGCTCCTCGTTTcc                     >  1:6566975‑M2/21‑51 (MQ=255)
                              ataagttggagtcattaccAGCGCCCCAGAACTTTTTGCTCCTCGTTTccc                    >  1:6148846‑M2/20‑51 (MQ=255)
                               taagttggagtcattaccAGCGCCCCAGAACTTTTTGCTCCTCGTTTCCCt                   >  1:5946799‑M2/19‑51 (MQ=255)
                               taagttggagtcattaccAGCGCCCCAGAACTTTTTGCTCCTCGTTTCCCt                   >  1:1116612‑M2/19‑51 (MQ=255)
                                aagttggagtcattaccAGCGCCCCAGAACTTTTTGCTCCTCGTTTCCCtt                  <  1:4651380‑M2/34‑1 (MQ=255)
                                aagttggagtcattaccAGCGCCCCAGAACTTTTTGCTCCTCGTTTCCCtt                  >  1:3389779‑M2/18‑51 (MQ=255)
                                aagttggagtcattaccAGCGCCCCAGAACTTTTTGCTCCTCGTTTCCCtt                  <  1:58994‑M2/34‑1 (MQ=255)
                                 agttggagtcattaccAGCGCCCCAGAACTTTTTGCTCCTCGTTTCCCTTc                 >  1:2365318‑M2/17‑51 (MQ=255)
                                 agttggagtcattaccAGCGCCCCAGAACTTTTTGCTCCTCGTTTCCCTTc                 >  1:6222643‑M2/17‑51 (MQ=255)
                                   ttggagtcattaccAGCGCCCCAGAACTTTTTGCTCCTCGTTTCCCTTCAt               >  1:5848261‑M2/15‑51 (MQ=255)
                                   ttggagtcattaccAGCGCCCCAGAACTTTTTGCTCCTCGTTTCCCTTCAt               >  1:6388930‑M2/15‑51 (MQ=255)
                                    tggagtcattaccAGCGCCCCAGAACTTTTTGCTCCTCGTTTCCCTTCAtt              <  1:2822013‑M2/38‑1 (MQ=255)
                                    tggagtcattaccAGCGCCCCAGAACTTTTTGCTCCTCGTTTCCCTTCAtt              >  1:5214777‑M2/14‑51 (MQ=255)
                                     ggagtcattaccAGCGCCCCAGAACTTTTTGCTCCTCGTTTCCCTTCATtg             >  1:6179273‑M2/13‑51 (MQ=255)
                                     ggagtcattaccAGCGCCCCAGAACTTTTTGCTCCTCGTTTCCCTTCATtg             <  1:746082‑M2/39‑1 (MQ=255)
                                       agtcattaccAGCGCCCCAGAACTTTTTGCTCCTCGTTTCCCTTCATtgtg           >  1:3433585‑M2/11‑51 (MQ=255)
                                       agtcattaccAGCGCCCCAGAACTTTTTGCTCCTCGTTTCCCTTCATtgtg           >  1:6781420‑M2/11‑51 (MQ=255)
                                         tcattaccAGCGCCCCAGAACTTTTTGCTCCTCGTTTCCCTTCATTGTGGc         <  1:6895830‑M2/43‑1 (MQ=255)
                                         tcattaccAGCGCCCCAGAACTTTTTGCTCCTCGTTTCCCTTCATTGTGGc         >  1:719945‑M2/9‑51 (MQ=255)
                                          cattaccAGCGCCCCAGAACTTTTTGCTCCTCGTTTCCCTTCATTGTGGCa        >  1:5824365‑M2/8‑51 (MQ=255)
                                          cattaccAGCGCCCCAGAACTTTTTGCTCCTCGTTTCCCTTCATTGTGGCa        >  1:3440301‑M2/8‑51 (MQ=255)
                                           attaccAGCGCCCCAGAACTTTTTGCTCCTCGTTTCCCTTCATTGTGGCaa       <  1:5926969‑M2/45‑1 (MQ=255)
                                           attaccAGCGCCCCAGAACTTTTTGCTCCTCGTTTCCCTTCATTGTGGCaa       <  1:1680598‑M2/45‑1 (MQ=255)
                                           attaccAGCGCCCCAGAACTTTTTGCTCCTCGTTTCCCTTCATTGTGGCaa       >  1:1924762‑M2/7‑51 (MQ=255)
                                            ttaccAGCGCCCCAGAACTTTTTGCTCCTCGTTTCCCTTCATTGTGGCAAc      >  1:3736944‑M2/6‑51 (MQ=255)
                                            ttaccAGCGCCCCAGAACTTTTTGCTCCTCGTTTCCCTTCATTGTGGCAAc      <  1:365447‑M2/46‑1 (MQ=255)
                                            ttaccAGCGCCCCAGAACTTTTTGCTCCTCGTTTCCCTTCATTGTGGCAAc      <  1:5024645‑M2/46‑1 (MQ=255)
                                            ttaccAGCGCCCCAGAACTTTTTGCTCCTCGTTTCCCTTCATTGTGGCAAc      >  1:7117172‑M2/6‑51 (MQ=255)
                                             taccAGCGCCCCAGAACTTTTTGCTCCTCGTTTCCCTTCATTGTGGCAACg     >  1:3819920‑M2/5‑51 (MQ=255)
                                             taccAGCGCACCAGAACTTTTTGCTCCTCGTTTCCCTTCATTGTGGCAACg     >  1:1011860‑M2/5‑51 (MQ=255)
                                                cAGCGCCCCAGAACTTTTTGCTCCTCGTTTCCCTTCATTGTGGCAACGttt  <  1:4709880‑M2/50‑1 (MQ=255)
                                                cAGCGCCCCAGAACTTTTTGCTCCTCGTTTCCCTTCATTGTGGCAACGttt  <  1:6097765‑M2/50‑1 (MQ=255)
                                                cAGCGCCCCAGAACTTTTTGCTCCTCGTTTCCCTTCATTGTGGCAACGttt  >  1:2060726‑M2/2‑51 (MQ=255)
                                                cAGCGCCCCAGAACTTTTTGCTCCTCGTTTCCCTTCATTGTGGCAACGttt  <  1:5733089‑M2/50‑1 (MQ=255)
                                                 |                                                 
TGCCAGTGATCCTGCTGAATGTTTTGCCGTACTTTCTGAATATTTCTTTAGCGCCCCAGAACTTTTTGCTCCTCGTTTCCCTTCATTGTGGCAACGTTT  >  REL606/1982264‑1982362

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