New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3091309 =	11 (0.550)	4 (0.220)	3/118	1.8	25.1%	coding (176/732 nt)	rsmE	16S rRNA m(3)U1498 methyltransferase, SAM‑dependent
	?	NC_000913	3091468 =	14 (0.790)					coding (335/732 nt)	rsmE	16S rRNA m(3)U1498 methyltransferase, SAM‑dependent

GGGCAGGCGTTGCAATTGTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGC  >  NC_000913/3091242‑3091368                                                                    |                                                            atcttctgccactgctgaagcttcttgttcagacgttcactatccagtttaacgccgcagcgctcaga                                                             >  2:560557‑M1/68‑68 (MQ=255)  ggCAGGCGTTGCAATTGTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCaaa                                                            >  1:177895/1‑68 (MQ=255)        cGTTGCAATTGTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGc                                                      >  1:417601/1‑68 (MQ=255)                 tgattGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCTAAAAAAGCGTGGAAGTg                                             >  2:679672/4‑68 (MQ=255)                 tGTTTGACGGTAGCGACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTg                                             <  2:463659/68‑1 (MQ=255)                        cGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGc                                      >  1:513479/1‑68 (MQ=255)                               aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAgg                               >  2:353415/1‑68 (MQ=255)                                 ccAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGcc                             >  1:644575/1‑68 (MQ=255)                                 ccAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGcc                             >  2:431553/1‑68 (MQ=255)                                 caaGGTCTTTGACGCCGAAATTACAAGCGACAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGcc                             >  2:11903/3‑68 (MQ=255)                                                    acgccgcagcgctcagAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAAtc          <  1:359761‑M1/52‑1 (MQ=255)                                                    acgccgcagcgctcagAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAAtc          <  1:573692‑M1/52‑1 (MQ=255)                                                      gccgcagcgctcagAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATc                 <  1:197714‑M1/45‑1 (MQ=255)                                                      gccgcagcgctcagAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATc                 >  2:197714‑M1/15‑59 (MQ=255)                                                           gctccAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTg   <  2:17375/65‑1 (MQ=255)                                                            cgcCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGc  >  1:411815/1‑68 (MQ=255)                                                                    |                                                            GGGCAGGCGTTGCAATTGTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGATCGCGAATCTCCGCTGC  >  NC_000913/3091242‑3091368

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