Predicted mutation
evidence seq id position mutation annotation gene description
MC JC REL606 547,700 Δ8,125 bp 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 546959–547698 555824–555805 8108–8866 34 [30] [0] 78 [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 = 547699NA (NA)78 (0.980) 52/100 0.2 100% noncoding (1/768 nt) IS1 repeat region
?REL606 555825 = 0 (0.000)coding (1110/2346 nt) ECB_00513 conserved hypothetical protein

GTATAAATATGTCTTACAGCCACTGCCTGATGTGTGGATACCGTTTAACGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATAAGCAAATAA  >  REL606/555777‑555874
                                                |                                                 
cctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGa                                                 <  1:240348‑M2/3‑1 (MQ=255)
cctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGa                                                 <  1:3832361‑M2/3‑1 (MQ=255)
cctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGa                                                 >  1:7273152‑M2/49‑51 (MQ=255)
cctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGa                                                 <  1:7206043‑M2/3‑1 (MQ=255)
cctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGa                                                 >  1:5652477‑M2/49‑51 (MQ=255)
cctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGa                                                 <  1:5820518‑M2/3‑1 (MQ=255)
cctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGa                                                 >  1:6032036‑M2/49‑51 (MQ=255)
 ctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGa                                                 <  1:5339433‑M2/3‑1 (MQ=255)
  tcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGAtt                                               <  1:74676‑M2/5‑1 (MQ=255)
   caaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGATTc                                              >  1:1621587‑M2/46‑51 (MQ=255)
   caaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGATTc                                              >  1:6294574‑M2/46‑51 (MQ=255)
   caaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGATTc                                              <  1:1736209‑M2/6‑1 (MQ=255)
    aaaaacaccatcatacactaaatcagtaagttggcagcatcaccCGATTCa                                             <  1:2591578‑M2/7‑1 (MQ=255)
       aacaccatcatacactaaatcagtaagttggcagcatcaccCGATTCACTg                                          <  1:2593732‑M2/10‑1 (MQ=255)
       aacaccatcatacactaaatcagtaagttggcagcatcaccCGATTCACTg                                          <  1:262037‑M2/10‑1 (MQ=255)
        acaccatcatacactaaatcagtaagttggcagcatcaccCGATTCACTgg                                         <  1:407544‑M2/11‑1 (MQ=255)
           ccatcatacactaaatcagtaagttggcagcatcaccCGATTCACTGGata                                      >  1:2315380‑M2/38‑51 (MQ=255)
           ccatcatacactaaatcagtaagttggcagcatcaccCGATTCACTGGata                                      >  1:246036‑M2/38‑51 (MQ=255)
           ccatcatacactaaatcagtaagttggcagcatcaccCGATTCACTGGata                                      <  1:4264297‑M2/14‑1 (MQ=255)
           ccatcatacactaaatcagtaagttggcagcatcaccCGATTCACTGGata                                      <  1:5116479‑M2/14‑1 (MQ=255)
            catcatacactaaatcagtaagttggcagcatcaccCGATTCACTGGatat                                     <  1:1310880‑M2/15‑1 (MQ=255)
            catcatacactaaatcagtaagttggcagcatcaccCGATTCACTGGatat                                     >  1:5107158‑M2/37‑51 (MQ=255)
            catcatacactaaatcagtaagttggcagcatcaccCGATTCACTGGatat                                     <  1:3978640‑M2/15‑1 (MQ=255)
            catcatacactaaatcagtaagttggcagcatcaccCGATTCACTGGatat                                     <  1:1598798‑M2/15‑1 (MQ=255)
             atcatacactaaatcagtaagttggcagcatcaccCGATTCACTGGATAtg                                    <  1:1947899‑M2/16‑1 (MQ=255)
             atcatacactaaatcagtaagttggcagcatcaccCGATTCACTGGATAtg                                    <  1:3233192‑M2/16‑1 (MQ=255)
              tcatacactaaatcagtaagttggcagcatcaccCGATTCACTGGATAtgc                                   <  1:7073509‑M2/17‑1 (MQ=255)
                atacactaaatcagtaagttggcagcatcaccCGATTCACTGGATAtgctt                                 <  1:5479435‑M2/19‑1 (MQ=255)
                atacactaaatcagtaagttggcagcatcaccCGATTCACTGGATAtgctt                                 >  1:5127369‑M2/33‑51 (MQ=255)
                 tacactaaatcagtaagttggcagcatcaccCGATTCACTGGATAtgcttg                                >  1:3449045‑M2/32‑51 (MQ=255)
                 tacactaaatcagtaagttggcagcatcaccCGATTCACTGGATAtgcttg                                <  1:5859670‑M2/20‑1 (MQ=255)
                  acactaaatcagtaagttggcagcatcaccCGATTCACTGGATAtgcttgc                               >  1:6688990‑M2/31‑51 (MQ=255)
                   cactaaatcagtaagttggcagcatcaccCGATTCACTGGATAtgcttgct                              >  1:4130091‑M2/30‑51 (MQ=255)
                    actaaatcagtaagttggcagcatcaccCGATTCACTGGATATGCTTGCTg                             <  1:2481877‑M2/23‑1 (MQ=255)
                    actaaatcagtaagttggcagcatcaccCGATTCACTGGATATGCTTGCTg                             >  1:3429004‑M2/29‑51 (MQ=255)
                    actaaatcagtaagttggcagcatcaccCGATTCACTGGATATGCTTGCTg                             >  1:376466‑M2/29‑51 (MQ=255)
                     ctaaatcagtaagttggcagcatcaccCGATTCACTGGATATGCTTGCTgg                            >  1:4468490‑M2/28‑51 (MQ=255)
                      taaatcagtaagttggcagcatcaccCGATTCACTGGATATGCTTGCTGGc                           <  1:2233541‑M2/25‑1 (MQ=255)
                      taaatcagtaagttggcagcatcaccCGATTCACTGGATATGCTTGCTGGc                           >  1:1820596‑M2/27‑51 (MQ=255)
                       aaatcagtaagttggcagcatcaccCGATTCACTGGATATGCTTGCTGGCt                          >  1:6790648‑M2/26‑51 (MQ=255)
                       aaatcagtaagttggcagcatcaccCGATTCACTGGATATGCTTGCTGGCt                          <  1:1665137‑M2/26‑1 (MQ=255)
                        aatcagtaagttggcagcatcaccCGATTCACTGGATATGCTTGCTGGCtt                         <  1:3942314‑M2/27‑1 (MQ=255)
                         atcagtaagttggcagcatcaccCGATTCACTGGATATGCTTGCTGGCttt                        <  1:5020752‑M2/28‑1 (MQ=255)
                         atcagtaagttggcagcatcaccCGATTCACTGGATATGCTTGCTGGCttt                        <  1:4685838‑M2/28‑1 (MQ=255)
                         atcagtaagttggcagcatcaccCGATTCACTGGATATGCTTGCTGGCttt                        >  1:2136808‑M2/24‑51 (MQ=255)
                          tcagtaagttggcagcatcaccCGATTCACTGGATATGCTTGCTGGCtttt                       <  1:2964389‑M2/29‑1 (MQ=255)
                           cagtaagttggcagcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTc                      <  1:5203816‑M2/30‑1 (MQ=255)
                           cagtaagttggcagcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTc                      >  1:133211‑M2/22‑51 (MQ=255)
                           cagtaagttggcagcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTc                      <  1:4051276‑M2/30‑1 (MQ=255)
                            agtaagttggcagcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTCg                     <  1:1358649‑M2/31‑1 (MQ=255)
                             gtaagttggcagcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGc                    >  1:738546‑M2/20‑51 (MQ=255)
                              taagttggcagcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGcc                   <  1:2192818‑M2/33‑1 (MQ=255)
                                agttggcagcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTg                 <  1:3117868‑M2/35‑1 (MQ=255)
                                 gttggcagcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTgg                <  1:6945203‑M2/36‑1 (MQ=255)
                                 gttggcagcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTgg                >  1:3509874‑M2/16‑51 (MQ=255)
                                  ttggcagcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGt               <  1:7134746‑M2/37‑1 (MQ=255)
                                   tggcagcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGtt              <  1:904867‑M2/38‑1 (MQ=255)
                                    ggcagcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTta             >  1:1365855‑M2/13‑51 (MQ=255)
                                    ggcagcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTta             >  1:1620675‑M2/13‑51 (MQ=255)
                                    ggcagcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTta             <  1:6404339‑M2/39‑1 (MQ=255)
                                     gcagcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTtat            <  1:1347727‑M2/40‑1 (MQ=255)
                                       agcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATaa          <  1:4195708‑M2/42‑1 (MQ=255)
                                       agcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATaa          <  1:7297584‑M2/42‑1 (MQ=255)
                                       agcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATaa          <  1:472533‑M2/42‑1 (MQ=255)
                                       agcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATaa          <  1:3997086‑M2/42‑1 (MQ=255)
                                        gcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTgg                <  1:5791901‑M2/36‑1 (MQ=255)
                                        gcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATAAg         >  1:3849039‑M2/9‑51 (MQ=255)
                                        gcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATAAg         <  1:7028763‑M2/43‑1 (MQ=255)
                                        gcatcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATAAg         <  1:435647‑M2/43‑1 (MQ=255)
                                         catcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATAAGc        >  1:809779‑M2/8‑51 (MQ=255)
                                          atcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATAAGCa       <  1:4356703‑M2/45‑1 (MQ=255)
                                          atcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATAAGCa       >  1:3778161‑M2/7‑51 (MQ=255)
                                          atcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATAAGCa       >  1:5825882‑M2/7‑51 (MQ=255)
                                           tcaccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATAAGCaa      <  1:3713872‑M2/46‑1 (MQ=255)
                                             accCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATAAGCAAAt    <  1:6951094‑M2/48‑1 (MQ=255)
                                              ccCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATAAGCAAATa   <  1:914973‑M2/49‑1 (MQ=255)
                                               cCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATAAGCAAATaa  >  1:6210412‑M2/2‑51 (MQ=255)
                                               cCGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATAAGCAAATaa  <  1:2743902‑M2/50‑1 (MQ=255)
                                                |                                                 
GTATAAATATGTCTTACAGCCACTGCCTGATGTGTGGATACCGTTTAACGATTCACTGGATATGCTTGCTGGCTTTTCGCCTGGTTATAAGCAAATAA  >  REL606/555777‑555874

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