Predicted mutation
evidence seq id position mutation annotation gene description
MC JC REL606 3,894,997 Δ6,934 bp IS150‑mediated rbsD[yieO] rbsD, rbsA, rbsC, rbsB, rbsK, rbsR, [yieO]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ REL606 3893611–3894994 3901930 6937–8320 21 [20] [0] 70 [insJ‑5]–[yieO] [insJ‑5],insK‑5,rbsD,rbsA,rbsC,rbsB,rbsK,rbsR,[yieO]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 = 3652533NA (NA)67 (1.160) 51/146 0.1 100% noncoding (1443/1443 nt) IS150 repeat region
?REL606 3901931 = 0 (0.000)coding (489/1428 nt) yieO predicted multidrug or homocysteine efflux system

CTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGTCTCAAGC  >  REL606/3652460‑3652608
                                                                         |                                                                           
cTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACaaa                                                                           >  1:771142/1‑76 (MQ=31)
cTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACaaa                                                                           >  1:1238077/1‑76 (MQ=31)
cTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACaaa                                                                           <  1:1818970/76‑1 (MQ=31)
 tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAAc                                                                          <  1:1082324/76‑1 (MQ=32)
 tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAAc                                                                          >  2:624739/1‑76 (MQ=32)
 tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAAc                                                                          <  1:623212/76‑1 (MQ=32)
 tGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAAc                                                                          <  2:379408/76‑1 (MQ=32)
  gACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACt                                                                         <  2:1063321/76‑1 (MQ=34)
  gACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACt                                                                         <  1:5929/76‑1 (MQ=34)
  gACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACt                                                                         <  1:5907/76‑1 (MQ=34)
  gACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACt                                                                         >  1:1558037/1‑76 (MQ=34)
   aCTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACtt                                                                        >  1:825582/1‑76 (MQ=35)
      ccAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACtttaa                                                                     >  2:1400991/1‑74 (MQ=38)
      ccAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGa                                                                     >  2:423985/1‑76 (MQ=40)
      ccAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGa                                                                     <  2:1369145/76‑1 (MQ=40)
       cAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGAt                                                                    >  2:925254/1‑76 (MQ=255)
         aTTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGAttt                                                                  >  1:565534/1‑76 (MQ=255)
          ttGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTta                                                                 <  2:615509/76‑1 (MQ=255)
           tGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTtat                                                                >  1:870559/1‑76 (MQ=255)
              atatCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGt                                                             >  1:824633/1‑76 (MQ=255)
               tatCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTc                                                            <  2:516971/76‑1 (MQ=255)
                atCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCa                                                           <  1:598737/76‑1 (MQ=255)
                atCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCa                                                           >  2:1194940/1‑76 (MQ=255)
                atCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCa                                                           <  2:1317928/76‑1 (MQ=255)
                   ggAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGt                                                        >  1:1004363/1‑76 (MQ=255)
                    gAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTg                                                       >  2:1424362/1‑76 (MQ=255)
                      aTCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTggg                                                     <  1:551625/76‑1 (MQ=255)
                       tCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTgggg                                                    >  2:1455698/1‑76 (MQ=255)
                       tCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTgggg                                                    >  2:364299/1‑76 (MQ=255)
                        cacaCCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGc                                                   >  2:40970/4‑76 (MQ=255)
                         aGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCt                                                  >  1:1687498/1‑76 (MQ=255)
                             cTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTtc                                              >  1:453329/1‑76 (MQ=255)
                              tataTGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTtct                                             <  2:1302836/76‑1 (MQ=255)
                               atatGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTtctc                                            >  1:714031/1‑76 (MQ=255)
                               atatGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTtctc                                            >  2:232421/1‑76 (MQ=255)
                               atatGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTtctc                                            >  2:1643678/1‑76 (MQ=255)
                               atatGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTtctc                                            >  1:452981/1‑76 (MQ=255)
                               atatGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTtctc                                            >  1:1835805/1‑76 (MQ=255)
                                tatGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTctct                                           <  2:437959/76‑1 (MQ=255)
                                 atGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCtctg                                          >  2:459138/1‑76 (MQ=255)
                                 atGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCtctg                                          >  1:1152782/1‑76 (MQ=255)
                                  tGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCtctgt                                         >  1:916576/1‑76 (MQ=255)
                                        gtgtTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCt                                   <  1:1291259/76‑1 (MQ=255)
                                          gtTTAACTGTCCAACTTTTTGNGGTCGGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCttt                                 <  2:296648/76‑1 (MQ=255)
                                          gtTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCttt                                 >  1:926927/1‑76 (MQ=255)
                                           tttAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTc                                >  1:134771/1‑76 (MQ=255)
                                           tttAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTc                                >  1:720792/1‑76 (MQ=255)
                                           tttAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTc                                >  1:696598/1‑76 (MQ=255)
                                            ttAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCg                               >  1:1765765/1‑76 (MQ=255)
                                             tAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCgg                              >  1:1283526/1‑76 (MQ=255)
                                             tAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCgg                              <  1:1554267/76‑1 (MQ=255)
                                             tAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCgg                              >  2:1524773/1‑76 (MQ=255)
                                              aaCTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGt                             >  1:1640534/1‑76 (MQ=255)
                                              aaCTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGt                             >  2:1734932/1‑76 (MQ=255)
                                              aaCTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGt                             >  1:272454/1‑76 (MQ=255)
                                              aaCTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGt                             >  1:770956/1‑76 (MQ=255)
                                                cTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGa                           >  1:838977/1‑76 (MQ=255)
                                                  gTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAAt                         >  2:997471/1‑76 (MQ=255)
                                                  gTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAAt                         <  1:833922/76‑1 (MQ=255)
                                                  gTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAAt                         >  1:972577/1‑76 (MQ=255)
                                                  gTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAAt                         >  2:1423505/1‑76 (MQ=255)
                                                       aCTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTg                    <  2:1261363/76‑1 (MQ=255)
                                                        cTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTga                   <  2:904705/76‑1 (MQ=255)
                                                          ttttGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTgaga                 >  1:120092/1‑76 (MQ=255)
                                                           tttGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGAc                >  1:49193/1‑76 (MQ=255)
                                                           tttGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGAc                >  2:870212/1‑76 (MQ=255)
                                                            ttGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACa               >  1:1118090/1‑76 (MQ=255)
                                                             tGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACaa              <  2:287496/76‑1 (MQ=255)
                                                             tGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACaa              >  1:1430523/1‑76 (MQ=255)
                                                                 gTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAg          >  1:1366600/1‑76 (MQ=255)
                                                                   cAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGtc        >  2:867335/1‑76 (MQ=255)
                                                                     gTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGtctc      <  2:1130837/76‑1 (MQ=255)
                                                                     gTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGtctc      <  2:1761236/76‑1 (MQ=255)
                                                                       aCAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGTCTCaa    >  2:579208/1‑76 (MQ=255)
                                                                       aCAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGTCTCaa    <  2:628399/76‑1 (MQ=255)
                                                                         aaaCTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGTCTCAAGc  <  1:1262045/76‑1 (MQ=255)
                                                                         aaaCTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGTCTCAAGc  >  2:625382/1‑76 (MQ=255)
                                                                         aaaCTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGTCTCAAGc  <  1:1124374/76‑1 (MQ=255)
                                                                         aaaCTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGTCTCAAGc  >  2:1069829/1‑76 (MQ=255)
                                                                         aaaCTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGTCTCAAGc  >  2:158824/1‑76 (MQ=255)
                                                                         aaaCTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGTCTCAAGc  >  2:1687087/1‑76 (MQ=255)
                                                                         |                                                                           
CTGACTCCAATTGAATATCGGAATCAGACCTATATGCCTCGTGTTTAACTGTCCAACTTTTTGGGGTCAGTACAAACTTTGATTTATAGTCAGGTGGGGCTTTTCTCTGTCTGCCTTTCGGTGAATACCTGAGACAAACAGTCTCAAGC  >  REL606/3652460‑3652608

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