New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	AP009048	4205005 =	48 (0.850)	3 (0.060)	3/96	NT	6.1%	coding (205/1290 nt)	rsmB	16S rRNA m5C967 methyltransferase, S‑adenosyl‑L‑methionine‑dependent
	?	AP009048	4205036 =	47 (0.890)					coding (174/1290 nt)	rsmB	16S rRNA m5C967 methyltransferase, S‑adenosyl‑L‑methionine‑dependent

CAAGAGTTGTGCTTTGGCGTACTGCGTACGCTTTCGCAGTTAGACTGGCTGATTAATAAGTTAATGGCCCGTC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AP009048/4205077‑4205005 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gtcTAACTGCGAAAGCGTACGCAGTACGCCAAAGCACAACTCTTGAAGAAGTGCTTTGTCTTTATCGGA  >  AP009048/4205036‑4205101                                                                                                                                             CAAGAGTTGTGCTTTGGCGTACTGCGTACGCTTTCGCAGTTAGACTGGCTGATTAATAAGTTAATGGCCCG                                                                      >  1:3697272/1‑71                                TTTCGCAGTTAGACTGGCTGATTAATAAGTTAATGGCCCGTC                                                                    <  1:1512099/42‑1                                             CTGGCTGATTAATAAGTTAATGGCCCGTCTAACTGCGAAAGCGTACGCAGTACGCCA                                        <  1:3178399/57‑1                                               GGCTGATTAATAAGTTAATGGCCCGTCTAACTGCGAAAGCGTACGCAGTACGCCAAAGCACAACTCTTG                          >  1:44475/1‑69                                                           AGTTAATGGCCCGTCTAACTGCGAAAGCGTACGCAGTACGCCAAAGCACAACTCTTGAAGAAGTGCTTTGT            <  1:550675/71‑1                                                                       GTCTAACTGCGAAAGCGTACGCAGTACGCCAAAGCACAACTCTTG                          <  1:142646/45‑1                                                                       GTCTAACTGCGAAAGCGTACGCAGTACGCAAAAGCACAACTCTTGAAGAAGTGCTTTGTCTTTATCGGA  <  1:2156930/69‑1                                                                        TCTAACTGCGAAAGCGTACGCAGTACGCCAAAGCACAACTCTTGAAGAAGTGC                 <  1:3254070/53‑1                                                                        TCTAACTGCGAAAGCGTACGCAGTACGCCAAAGCACAACTCTTGAAGAAGT                   >  1:4720780/1‑51                                                                        TCTAACTGCGAAAGCGTACGCAGTACGCCAAAGCACAACTC                             >  1:4425253/1‑41                                                                         CTAACTGCGAAAGCGTACGCAGTACGCCAAAGCACAACTCTTGAAGAAG                    >  1:4480873/1‑49                                                                                                                                             CAAGAGTTGTGCTTTGGCGTACTGCGTACGCTTTCGCAGTTAGACTGGCTGATTAATAAGTTAATGGCCCGTC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  AP009048/4205077‑4205005 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gtcTAACTGCGAAAGCGTACGCAGTACGCCAAAGCACAACTCTTGAAGAAGTGCTTTGTCTTTATCGGA  >  AP009048/4205036‑4205101

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