New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 232533730 (0.740)3 (0.080) 3/154 2.8 9.5% pseudogene (2455/4515 nt) BW25113_RS11635 alpha‑2‑macroglobulin family protein
?NZ_CP009273 = 2325338 30 (0.820)pseudogene (2454/4515 nt) BW25113_RS11635 alpha‑2‑macroglobulin family protein
Rejected: Frequency below/above cutoff threshold.

TGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NZ_CP009273/2325245‑2325337
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gttgcgccTGGCGCAACTAACCGCAGCGAGCGGCGAGTAAAAATGCTTGAACGTCCACGGCGTGAAGAGGTGGATACCGCGGCATGGATGCCGT  <  NZ_CP009273/2325338‑2325253
                                                                                                                                                                                   
TGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGC                                                                                           >  2:874806/1‑90
  TGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGC                                                                                         <  1:482831/90‑1
  TGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGC                                                                                         <  2:916182/90‑1
                                     TCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCTGGCGCAACTAACCGCAGCGAGCGGCGAGTAAAA                                                      >  2:668819/1‑90
                                             TGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCCTGGCGCAACTAACCGCAGCGAGCGGCGAGTAAAAATGCTTGA                                              >  1:502291/1‑90
                                                                                  TTAGTTGCGCCTGGCGCAACTAACCGCAGCGAGCGGCGAGTAAAAATGCTTGAACGTCCACGGCGTGAAGAGGTGGATACCGCGGCATGG         >  1:4248/1‑90
                                                                                         CGCCTGGCGCAACTAACCGCAGCGAGCGGCGAGTAAAAATGCTTGAACGTCCACGGCGTGAAGAGGTGGATACCGCGGCATGGATGCCGT  <  1:440354/90‑1
                                                                                                                                                                                   
TGTGAGTGACGGCATCCATGCCGCGGTATCCACCTCTTCACGCCGTGGACGTTCAAGCATTTTTACTCGCCGCTCGCTGCGGTTAGTTGCGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NZ_CP009273/2325245‑2325337
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gttgcgccTGGCGCAACTAACCGCAGCGAGCGGCGAGTAAAAATGCTTGAACGTCCACGGCGTGAAGAGGTGGATACCGCGGCATGGATGCCGT  <  NZ_CP009273/2325338‑2325253

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