New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	4398991 =	42 (2.120)	3 (0.180)	3/136	1.5	12.1%	intergenic (+79/‑101)	rnr/rlmB	ribonuclease R/23S rRNA (guanosine(2251)‑2'‑O)‑methyltransferase RlmB
	?	NZ_CP009273	4399052 =	9 (0.550)					intergenic (+140/‑40)	rnr/rlmB	ribonuclease R/23S rRNA (guanosine(2251)‑2'‑O)‑methyltransferase RlmB
Rejected: Frequency below/above cutoff threshold.

ATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCT  >  NZ_CP009273/4398968‑4399135                                                                                     |                                                                                    acgttttgcctgatgcactacgtttatccggcctacatgatctctgcaatatattgaatttgcgtgattatgtaggccggataaACGTAg                                                                                <  1:276121‑M2/6‑1 (MQ=255)                  ctacgtTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGt                                                               <  1:276121‑M1/84‑1 (MQ=255)                        TTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAaacaa                                                         >  1:225951‑M1/1‑90 (MQ=255)                          ATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACg                                                       >  2:55782‑M1/1‑90 (MQ=255)                          ATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACg                                                       >  2:287662‑M1/1‑90 (MQ=255)                           TCCGGCCTACATAATCACGCAA‑TTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAg                                                     >  1:142166‑M1/1‑90 (MQ=255)                              ggcctacatgatctctgcaatatattgaatttgcgtgattatgtaggccggataaACGTAGTGCATCAGGCAAAACGTAAACAACGAGTa                                                   >  1:25783‑M2/56‑90 (MQ=255)                                           aCGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGa                                      >  1:555280/1‑90 (MQ=255)                                             tgcaatatattgaatttgcgtgattatgtaggccggataaACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGaaa                                    <  2:379417‑M2/50‑1 (MQ=255)                                             gCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGaaa                                    <  2:478454/90‑1 (MQ=255)                                             gCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGaaa                                    >  2:460774/1‑90 (MQ=255)                                             gCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGaaa                                    >  2:19230/1‑90 (MQ=255)                                                      atatatTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACg                           <  2:341260/90‑1 (MQ=255)                                                          atTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCAt                       >  1:335981/1‑90 (MQ=255)                                                           ttGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATc                      >  2:224623/1‑90 (MQ=255)                                                                gagaTCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGc                 <  1:24248/90‑1 (MQ=255)                                                                gagaTCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGc                 <  1:227977/90‑1 (MQ=255)                                                                            gccggataaACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCct     >  2:449814‑M2/10‑90 (MQ=255)                                                                            gccggataaACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCct     >  1:574313‑M2/10‑90 (MQ=255)                                                                            gccggataaACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCct     >  2:370879‑M2/10‑90 (MQ=255)                                                                              cggataaACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCctgc   <  1:449814‑M2/83‑1 (MQ=255)                                                                               ggataaACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGc        <  2:156346‑M2/78‑1 (MQ=255)                                                                               ggataaACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGc        >  1:156346‑M2/7‑84 (MQ=255)                                                                               ggataaACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCctgct  >  1:536563‑M2/7‑90 (MQ=255)                                                                                     |                                                                                    ATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACATAATCACGCAAATTCAATATATTGCAGAGATCATGTAGGCCGGATAAACGTAGTGCATCAGGCAAAACGTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCT  >  NZ_CP009273/4398968‑4399135

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