New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	4399442 =	23 (0.850)	3 (0.130)	3/150	2.1	12.4%	coding (351/732 nt)	rlmB	23S rRNA (guanosine(2251)‑2'‑O)‑methyltransferase RlmB
	?	NZ_CP009273	4399525 =	22 (0.920)					coding (434/732 nt)	rlmB	23S rRNA (guanosine(2251)‑2'‑O)‑methyltransferase RlmB
Rejected: Frequency below/above cutoff threshold.

GATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCG  >  NZ_CP009273/4399355‑4399524                                                                                        |                                                                                   gATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGAcgc                                                                                  >  2:209367/1‑90 (MQ=255)    cgacgatccagatattctcttcctgcaacatacgcatggtgcgcgcaaggttagtcacccgaatcagcggaacgctttctgccgcgccgc                                                                               >  2:209523‑M1/85‑90 (MQ=255)      cttcgctCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAg                                                                             <  2:616341/90‑1 (MQ=255)        tcgctCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGc                                                                           <  2:706578/90‑1 (MQ=255)        tcgctCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGc                                                                           <  1:283041/90‑1 (MQ=255)               aTCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATg                                                                    >  1:341781/1‑90 (MQ=255)                         ccTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATtgt                                                          <  2:428553/90‑1 (MQ=255)                            gctgATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGcc                                                       <  1:721959/90‑1 (MQ=255)                              tgATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGa                                                     >  1:746917/1‑90 (MQ=255)                               gATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGaa                                                    <  2:252076/90‑1 (MQ=255)                                   cTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGAt                                                <  2:503112/90‑1 (MQ=255)                                        cGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTc                                           <  2:642440/90‑1 (MQ=255)                                         ggCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTcc                                          >  2:135850/1‑90 (MQ=255)                                          gCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCg                                         >  1:713359/1‑90 (MQ=255)                                             tAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGcac                                      <  1:716765/90‑1 (MQ=255)                                             tAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGcac                                      <  2:341781/90‑1 (MQ=255)                                                    tCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCaa                               >  1:660140/1‑90 (MQ=255)                                                           aaCCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACg                        <  2:660140/90‑1 (MQ=255)                                                           aaCCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACg                        <  1:192862/90‑1 (MQ=255)                                                             cccgaatcagcggaacgctttctgccgCGCCGTAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGc                      <  2:708117‑M1/63‑1 (MQ=255)                                                               cgaatcagcggaacgctttctgccgCGCCGCAGGCGTTCATGCGGTGATTGTGGCGAAAGATCGCTCCGCCCAGCTCAACGCTCCGGcgg                    <  1:362033‑M1/65‑2 (MQ=255)                                                                cGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGaa                   >  2:610826/1‑90 (MQ=255)                                                                     agcggaacgctttctgccgCGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAg              <  1:638788‑M1/71‑1 (MQ=255)                                                                      gTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGt             >  1:359378/1‑90 (MQ=255)                                                                      gTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGt             >  1:137939/1‑90 (MQ=255)                                                                                 gTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTgcggcg  >  1:78213/1‑90 (MQ=255)                                                                                 gTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTgcggcg  >  2:601859/1‑90 (MQ=255)                                                                                        |                                                                                   GATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCG  >  NZ_CP009273/4399355‑4399524

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