New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	= 582186	107 (1.210)	6 (0.070)	5/240	NT	5.3%	coding (49/1923 nt)	selB	selenocysteyl‑tRNA‑specific translation elongation factor
	?	NC_002947	= 582220	108 (1.250)					coding (83/1923 nt)	selB	selenocysteyl‑tRNA‑specific translation elongation factor

TGGGCCGCATTGACGACGATGCCCTGTGGCTGGACCTGCGCCAGCTGGATGACGAAGCCCAGTGGCTGGCCCAGTTGCCAGCCCTGCAGCTGGGGCCGGTGCAGTGATCGTTGGTACCGCTGGCCACATCGACCACGGCAAGACCGCACTGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/582035‑582186 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgcCGCTGGTCGCCGGCCTGGCCGGTCAGCGCCAGCAGCAGT  <  NC_002947/582220‑582182                                                                                                                                                                                                 TGGGCCGCATTGACGACGATGCCCTGTGGCTGGACCTGCGCCAGCTGGATGACGAAGCCCAGTGGCTGGCCCAGTTGCCAGCCCTGCAGCTGGGGCCGGTGCAGTGATCGTTGGTACCGCTGGCCACATCGACCACGGCAAGACCGCACT                                           >  1:1663485/1‑150                                         CCAGCTGGATGACGAAGCCCAGTGGCTGGCCCAGTTGCCAGCCCTGCAGCTGGGGGCGGTGCAGTGATCGTTGGTACCGCTGGCCACATCGACCACGGCAAGACCGCACTGCCGCTGGTCGCCGGCCTGGCCGGTCAGCGCCAGCAGCAG   <  2:1219428/150‑1                                          CAGCTGGATGACGAAGCCCAGTGGCTGGCCCAGTTGCCAGCCCTGCAGCTGGGGCCGGTGCAGTGATCGTTGGTACCGCTGGCCACATCGACCACGGCAAGACCGCACTGCCGCTGGTCGCCGGCCTGGCCGGTCAGCGCCAGCAGCAGT  <  1:1915154/150‑1                                                                                AGCCCTGCAGCTGGGGCCGGTGCAGTGATCGTTGGTACCGCTGGCCACATCGACCACGGCAAGACCGCACTGCCGCTGGTCGCCGGCCTGGCCGGTCAGCGCCAGCAG      <  1:1000994/108‑1                                                                                AGCCCTGCAGCTGGGGCCGGTGCAGTGATCGTTGGTACCGCTGGCCACATCGACCACGGCAAGACCGCACTGCCGCTGGTCGCCGGCCTGGCCGGTCAGCGCCAGCAG      >  2:1000994/1‑108                                                                                                                TGGTACCGCTGGCCACATCGACCACGGCAAGACCGCACTGCCGCTGGTCGCCGGCCTGGCCGGTCAGCGCCAGCAGCAGT  >  1:1334320/1‑80                                                                                                                TGGTACCGCTGGCCACATCGACCACGGCAAGACCGCACTGCCGCTGGTCGCCGGCCTGGCCGGTCAGCGCCAGCAGCAGT  <  2:1334320/80‑1                                                                                                                                                                                                 TGGGCCGCATTGACGACGATGCCCTGTGGCTGGACCTGCGCCAGCTGGATGACGAAGCCCAGTGGCTGGCCCAGTTGCCAGCCCTGCAGCTGGGGCCGGTGCAGTGATCGTTGGTACCGCTGGCCACATCGACCACGGCAAGACCGCACTGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/582035‑582186 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgcCGCTGGTCGCCGGCCTGGCCGGTCAGCGCCAGCAGCAGT  <  NC_002947/582220‑582182

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