New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AP009048 = 7126587 (1.070)4 (0.070)
+AGGGTCCCTAAGTC		 3/70 NT 8.3% coding (879/879 nt) araC DNA‑binding transcriptional dual regulator
?AP009048 = 3464920 37 (0.450)intergenic (‑296/+5) murB/rrfB UDP‑N‑acetylenolpyruvoylglucosamine reductase, FAD‑binding/5S ribosomal RNA

GAATGATGTAGCCGTCAAGTTGTCATAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AP009048/71238‑71265
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TAAAACGAAAGGCTCAGTCGAAAGACTGGGCC  <  AP009048/3464920‑3464889
                            ||||||||||||||                                
GAATGATGTAGCCGTCAAGTTGTCATAAAGGGTCCCTAAGTCTAAAACGAAAGGCTCAGTCGAAAGACTGG     >  1:125376/1‑71
   TGATGTAGCCGTCAAGTTGTCATAAAGGGTCCCTAAGTCTAAAACGAAAGGCTCAGTCGAAAGACTGGGCC  >  1:625004/1‑71
   TGATGTAGCCGTCAAGTTGTCATAAAGGGTCCCTAAGTCTAAAACGAAAGGCTCAGTCGAAAGACTGGGCC  >  1:7450081/1‑71
      TGTAGCCGTCAAGTTGTCATAAAGGGTCCCTAAGTCTAAAACGAAAGGCTCAGTCGAAAGACTGGGC   >  1:7670336/1‑67
                              GGTCCCTAAGTCTAAAACGAAAGGCTCAGTCGAAAGACTGGGC   >  1:2670660/1‑43
                              GGTCCCTAAGTCTAAAACGAAAGGCTCAGTCGAAAGACTGGGC   >  1:3954350/1‑43
                              GGTCCCTAAGTCTAAAACGAAAGGCTCAGTCGAAAGACTGGGC   >  1:6866585/1‑43
                            ||||||||||||||                                
GAATGATGTAGCCGTCAAGTTGTCATAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AP009048/71238‑71265
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TAAAACGAAAGGCTCAGTCGAAAGACTGGGCC  <  AP009048/3464920‑3464889

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