New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1482256 =0 (0.000)3 (0.150)
+GTATGGCTGCGAT
3/70 2.1 28.3% coding (605/606 nt) azoR FMN dependent NADH:quinone oxidoreductase
?NC_000913 = 3693451 21 (0.770)coding (1786/2619 nt) bcsA cellulose synthase catalytic subunit

TGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTT  >  NC_000913/1482208‑1482304
                                                |                                                
cctatcggcagaccgacaagataccacgccagcacatcgcagccatact                                                  <  1:478480‑M2/1‑1 (MQ=255)
   atcggcagaccgacaagataccacgccagcacatcgcagccatacTATg                                               >  1:1898400‑M2/46‑49 (MQ=255)
      ggcagaccgacaagataccacgccagcacatcgcagccatacTATg                                               >  1:1882518‑M2/43‑46 (MQ=255)
        cagaccgacaagataccacgccagcacatcgcagccatacTATGCAGa                                           >  1:1371263‑M2/41‑48 (MQ=255)
           accgacaagataccacgccagcacatcgcagccatacTATGCAGAAAc                                        >  1:1280449‑M2/38‑48 (MQ=255)
                caagataccacgccagcacatcgcagccatacTATGCAGAAACAATGCt                                  >  1:1152491‑M2/33‑49 (MQ=255)
                 aagataccacgccagcacatcgcagccatacTATGCAGAAACAATGCTg                                 >  1:192147‑M2/32‑49 (MQ=255)
                  agataccacgccagcacatcgcagccatacTATGCAGAAACAATGCTGt                                <  1:1359649‑M2/19‑1 (MQ=255)
                    ataccacgccagcacatcgcagccatacTATGCAGAAACAATGCTGTc                               >  1:163659‑M2/29‑48 (MQ=255)
                    ataccacgccagcacatcgcagccatacTATGCAGAAACAATGCTGTCg                              <  1:1305018‑M2/21‑1 (MQ=255)
                    ataccacgccagcacatcgcagccatacTATGCAGAAACAATGCTGTCg                              <  1:1207274‑M2/21‑1 (MQ=255)
                      accacgccagcacatcgcagccatacTATGCAGAAACAATGCTGTCGAt                            >  1:1192822‑M2/27‑49 (MQ=255)
                      accacgccagcacatcgcagccatacTATGCAGAAACAATGCTGTCGAt                            >  1:424227‑M2/27‑49 (MQ=255)
                           gccagcacatcgcagccatacTATGCAGAAACAATGCTGTCGATGgctg                       >  1:1598266‑M2/22‑49 (MQ=255)
                           gccagcacatcgcagccatacTATGCAGAAACAATGCTGTCGATGgctg                       >  1:29936‑M2/22‑49 (MQ=255)
                               gcacatcgcagccatacTATGCAGAAACAATGCTGTCGATGGCTGCttt                   >  1:958181‑M1/18‑49 (MQ=255)
                                cacatcgcagccatacTATGCAGAAACAATGCTGTCGATGGCTGCtttt                  >  1:2029381‑M1/17‑49 (MQ=255)
                                 acatcgcagccatacTATGCAGAAACAATGCTGTCGATGGCTGCTTTTg                 <  1:1524416‑M1/34‑1 (MQ=255)
                                      gcagccatacTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGt              >  1:1757041‑M2/11‑47 (MQ=255)
                                       cagccatacTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAg           <  1:1173482‑M2/40‑1 (MQ=255)
                                        agccatacTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGa          >  1:727759‑M2/9‑49 (MQ=255)
                                        agccatacTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGa          <  1:1074811‑M2/41‑1 (MQ=255)
                                         gccatacTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGAc         >  1:241531‑M2/8‑49 (MQ=255)
                                            atacTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACtgt      <  1:109502‑M2/45‑1 (MQ=255)
                                            atacTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACtgt      <  1:126668‑M2/45‑1 (MQ=255)
                                             tacTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACtgtg     <  1:1961050‑M2/46‑1 (MQ=255)
                                              acTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGc    >  1:994448‑M2/3‑49 (MQ=255)
                                                tATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCtt  >  1:1226965/1‑49 (MQ=255)
                                                |                                                
TGAGATGTGGACTGAATATCTACAGTCCACATCAAGACCGTGTCCGGTTATGCAGAAACAATGCTGTCGATGGCTGCTTTTGCGTCAGACTGTGCTT  >  NC_000913/1482208‑1482304

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