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

TGTAATAAGCATTGCCAGAGCGGCAGAAAACAACATTTTTTTCATCTTATTATCCTGCATTGTTAAAA  >  REL606/550318‑550385
                                 |                                  
tGTAATAAGCATTGCCAGAGCGGCAGAAAacaaca                                   >  1:10716020/1‑35 (MQ=255)
catacactaaatcagtaagttggcagcatcaccCa                                   <  1:16639672‑M2/2‑1 (MQ=255)
catacactaaatcagtaagttggcagcatcaccCa                                   <  1:19090544‑M2/2‑1 (MQ=255)
catacactaaatcagtaagttggcagcatcaccCa                                   <  1:670116‑M2/2‑1 (MQ=255)
catacactaaatcagtaagttggcagcatcaccCa                                   <  1:20218424‑M2/2‑1 (MQ=255)
catacactaaatcagtaagttggcagcatcaccCa                                   <  1:14389154‑M2/2‑1 (MQ=255)
catacactaaatcagtaagttggcagcatcaccCa                                   <  1:561119‑M2/2‑1 (MQ=255)
catacactaaatcagtaagttggcagcatcaccCa                                   >  1:2683632‑M2/34‑35 (MQ=255)
catacactaaatcagtaagttggcagcatcaccCa                                   <  1:11746488‑M2/2‑1 (MQ=255)
catacactaaatcagtaagttggcagcatcaccCa                                   >  1:3034328‑M2/34‑35 (MQ=255)
catacactaaatcagtaagttggcagcatcaccCa                                   >  1:3799628‑M2/34‑35 (MQ=255)
catacactaaatcagtaagttggcagcatcaccCa                                   >  1:10609978‑M2/34‑35 (MQ=255)
catacactaaatcagtaagttggcagcatcaccCa                                   <  1:4511867‑M2/2‑1 (MQ=255)
catacactaaatcagtaagttggcagcatcaccCa                                   <  1:14457173‑M2/2‑1 (MQ=255)
catacactaaatcagtaagttggcagcatcacaCa                                   <  1:9220841‑M2/2‑1 (MQ=255)
 atacactaaatcagtaagttggcagcatcaccCAt                                  >  1:6999817‑M2/33‑35 (MQ=255)
 atacactaaatcagtaagttggcagcatcaccCAt                                  >  1:19336272‑M2/33‑35 (MQ=255)
  tacactaaatcagtaagttggcagcatcaccCatg                                 >  1:7082398‑M2/32‑34 (MQ=255)
  tacactaaatcagtaagttggcagcatcaccCAtt                                 <  1:6522136‑M2/4‑1 (MQ=255)
  tacactaaatcagtaagttggcagcatcaccCAtt                                 <  1:7148470‑M2/4‑1 (MQ=255)
  tacactaaatcagtaagttggcagcatcaccCAtt                                 >  1:4200401‑M2/32‑35 (MQ=255)
   acactaaatcagtaagttggcagcatcaccCAttt                                <  1:477804‑M2/5‑1 (MQ=255)
   acactaaatcagtaagttggcagcatcaccCAttt                                <  1:10049345‑M2/5‑1 (MQ=255)
    cactaaatcagtaagttggcagcatcaccCAtttt                               >  1:2350079‑M2/30‑35 (MQ=255)
    cactaaatcagtaagttggcagcatcaccCAtttt                               <  1:17801198‑M2/6‑1 (MQ=255)
    cactaaatcagtaagttggcagcatcaccCAtttt                               <  1:9063694‑M2/6‑1 (MQ=255)
      ctaaatcagtaagttggcagcatcaccCAtttttt                             >  1:6648175‑M2/28‑35 (MQ=255)
      ctaaatcagtaagttggcagcatcaccCAtttttt                             <  1:11218900‑M2/8‑1 (MQ=255)
       taaatcagtaagttggcagcatcagcCAttttttt                            <  1:15599576‑M2/9‑1 (MQ=255)
       taaatcagtaagttggcagcatcaccCAttttttt                            <  1:18708835‑M2/9‑1 (MQ=255)
       taaatcagtaagttggcagcatcaccCAttttttt                            <  1:5470538‑M2/9‑1 (MQ=255)
        gCATTGCCAGAGCGGCAGAAAACAACATTTTTTTc                           <  1:19901404/35‑1 (MQ=255)
        aaatcagtaagttggcagcatcaccCATTTTTTTc                           >  1:20396147‑M2/26‑35 (MQ=255)
        aaatcagtaagttggcagcatcaccCATTTTTTTc                           >  1:12602237‑M2/26‑35 (MQ=255)
         cATTGCCAGAGCGGCAGAAAACAACATTTTTTTCa                          <  1:20023630/35‑1 (MQ=255)
         aatcagtaagttggcagcatcaccCATTTTTTTCa                          >  1:12250178‑M2/25‑35 (MQ=255)
         aatcagtaagttggcagcatcaccAATTTTTTTCa                          <  1:12836556‑M2/11‑1 (MQ=255)
          atcagtaagttggcagcatcaccCATTTTTTTCAt                         <  1:1432407‑M2/12‑1 (MQ=255)
          atcagtaagttggcagcatcaccCATTTTTTTCAt                         <  1:1104505‑M2/12‑1 (MQ=255)
          atcagtaagttggcagcatcaccCATTTTTTTCAt                         >  1:18133760‑M2/24‑35 (MQ=255)
           tcagtaagttggcagcatcaccCATTTTTTTCATc                        >  1:6910058‑M2/23‑35 (MQ=255)
           tcagtaagttggcagcatcaccCATTTTTTTCATc                        <  1:4882655‑M2/13‑1 (MQ=255)
           tcagtaagttggcagcatcaccCATTTTTTTCATc                        >  1:18967033‑M2/23‑35 (MQ=255)
           tcagtaagttggcagcatcaccCATTTTTTTCATc                        >  1:4044685‑M2/23‑35 (MQ=255)
           tcagtaagttggcagcatcaccCATTTTTTTCATc                        >  1:10598480‑M2/23‑35 (MQ=255)
            cagtaagttggcagcatcaccCATTTTTTTCATCt                       <  1:19044861‑M2/14‑1 (MQ=255)
            cagtaagttggcagcatcaccCATTTTTTTCATCt                       <  1:653772‑M2/14‑1 (MQ=255)
            cagtaagttggcagcataaccCATTTTTTTCATCt                       <  1:6249004‑M2/14‑1 (MQ=255)
             gCCAGAGCGGCAGAAAACAACATTTTTTTCATCtt                      <  1:3196388/35‑1 (MQ=255)
             agtaagtttgcagcatcaccCATTTTTTTCATCtt                      >  1:1699091‑M2/21‑35 (MQ=255)
             agtaagttggcagcatcaccCATTTTTTTCATCtt                      <  1:12901107‑M2/15‑1 (MQ=255)
               taagttggcagcatcaccCATTTTTTTCATCttat                    >  1:13470782‑M2/19‑35 (MQ=255)
               taagttggcagcataaccCATTTTTTTCATCttat                    <  1:8056764‑M2/17‑1 (MQ=255)
                aagttggcagcatcaccCATTTTTTTCATCttatt                   >  1:2861219‑M2/18‑35 (MQ=255)
                aagttggcagcatcaccCATTTTTTTCATCttatt                   <  1:3673967‑M2/18‑1 (MQ=255)
                aagttggcagcatcaccCATTTTTTTCATCttatt                   >  1:1398461‑M2/18‑35 (MQ=255)
                 agttggcagcatcaccCATTTTTTTCATCttatta                  >  1:2798034‑M2/17‑35 (MQ=255)
                 agttggcagcatcaccCATTTTTTTCATCttatta                  >  1:6731446‑M2/17‑35 (MQ=255)
                 agttggcagcatcaccCATTTTTTTCATCttatta                  >  1:6748301‑M2/17‑35 (MQ=255)
                  gttggcagcatcaccCATTTTTTTCATCTtattat                 <  1:6050896‑M2/20‑1 (MQ=255)
                  gttggcagcatcaccCATTTTTTTCATCTtattat                 <  1:3845981‑M2/20‑1 (MQ=255)
                  gttggcagcatcaccCATTTTTTTCATCTtattat                 >  1:14832810‑M2/16‑35 (MQ=255)
                    tggcagcatcaccCATTTTTTTCATCTTATTATcc               >  1:2571021‑M2/14‑35 (MQ=255)
                    tggcagcatcaccCATTTTTTTCATCTTATTATcc               <  1:18849183‑M2/22‑1 (MQ=255)
                     ggcagcatcaccCATTTTTTTCATCTTATTATCCt              >  1:14876006‑M2/13‑35 (MQ=255)
                     ggcagcatcaccCATTTTTTTCATCTTATTATCCt              >  1:4022064‑M2/13‑35 (MQ=255)
                     ggcagcatcaccCATTTTTTTCATCTTATTATCCt              <  1:803152‑M2/23‑1 (MQ=255)
                     ggcagcatcaccCATTTTTTTCATCTTATTATCCt              <  1:20201317‑M2/23‑1 (MQ=255)
                      gcagcatcaccCATTTTTTTCATCTTATTATCCTg             <  1:3733042‑M2/24‑1 (MQ=255)
                      gcagcatcaccCATTTTTTTCATCTTATTATCCTg             >  1:14404573‑M2/12‑35 (MQ=255)
                      gcagcatcaccCATTTTTTTCATCTTATTATCCTg             >  1:15302191‑M2/12‑35 (MQ=255)
                      gcagcatcaccCATTTTTTTCATCTTATTATCCTg             <  1:18509552‑M2/24‑1 (MQ=255)
                       cagcatcaccCATTTTTTTCATCTTATTATCCTGc            <  1:7342326‑M2/25‑1 (MQ=255)
                       cagcatcaccCATTTTTTTCATCTTATTATCCTGc            >  1:8293703‑M2/11‑35 (MQ=255)
                       cagcatcaccAATTTTTTTCATCTTATTATCCTGc            <  1:741897‑M2/25‑1 (MQ=255)
                        agcatcaccCATTTTTTTCATCTTATTATCCTGCa           >  1:5623083‑M2/10‑35 (MQ=255)
                        agcatcaccCATTTTTTTCATCTTATTATCCTGCa           >  1:5563768‑M2/10‑35 (MQ=255)
                          catcaccCATTTTTTTCATCTTATTATCCTGCAtt         >  1:1258776‑M2/8‑35 (MQ=255)
                          catcaccCATTTTTTTCATCTTATTATCCTGCAtt         <  1:13213322‑M2/28‑1 (MQ=255)
                          catcaccCATTTTTTTCATCTTATTATCCTGCAtt         <  1:7238959‑M2/28‑1 (MQ=255)
                          catcaccCATTTTTTTCATCTTATTATCCTGCAtt         <  1:15149477‑M2/28‑1 (MQ=255)
                          catcaccCATTTTTTTCATCTTATTATCCTGCAtt         <  1:3420287‑M2/28‑1 (MQ=255)
                          catcaccCATTTTTTTCATCTTATTATCCTGCAtt         <  1:4314625‑M2/28‑1 (MQ=255)
                           cataccCATTTTTTTCATCTTATTATCCTGCATTg        <  1:465853‑M2/29‑1 (MQ=255)
                           atcaccCATTTTTTTCATCTTATTATCCTGCATTg        <  1:2245245‑M2/29‑1 (MQ=255)
                           atcaccCATTTTTTTCATCTTATTATCCTGCATTg        <  1:12131796‑M2/29‑1 (MQ=255)
                           aaaCAACATTTTTTTCATCTTATTATCCTGCATTg        >  1:10803757/1‑35 (MQ=255)
                            tcacgCATTTTTTTCATCTTATTATCCTGCATTGt       >  1:14960340‑M2/6‑35 (MQ=255)
                            aacaacATTTTTTTCATCTTATTATCCTGCATTGt       >  1:8036022/1‑35 (MQ=255)
                             caccCATTTTTTTCATCTTATTATCCTGCATTGtt      <  1:5612833‑M2/31‑1 (MQ=255)
                               ccCATTTTTTTCATCTTATTATCCTGCATTGTTaa    <  1:2873911‑M2/33‑1 (MQ=255)
                               ccCATTTTTTTCATCTTATTATCCTGCATTGTTaa    <  1:10867519‑M2/33‑1 (MQ=255)
                                 caTTTTTTTCATCTTATTATCCTGCATTGTTaaaa  >  1:19395556/1‑35 (MQ=255)
                                 caTTTTTTTCATCTTATTATCCTGCATTGTTaaaa  <  1:9921012/35‑1 (MQ=255)
                                 |                                  
TGTAATAAGCATTGCCAGAGCGGCAGAAAACAACATTTTTTTCATCTTATTATCCTGCATTGTTAAAA  >  REL606/550318‑550385

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