New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4409595	35 (0.860)	6 (0.150)	6/284	2.6	14.6%	coding (321/732 nt)	rlmB	23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent
	?	NC_000913	= 4409617	37 (0.950)					coding (343/732 nt)	rlmB	23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent
Rejected: Frequency below cutoff threshold.

GAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTT  >  NC_000913/4409454‑4409745                                                                                                                                              |                                                                                                                                                       gaAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCgg                                                                                                                                                >  1:80059/1‑150 (MQ=255)  aaaaGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACtacgcaggca                                                                                                                                              <  1:97162‑M1/151‑11 (MQ=255)    aaGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGcgcg                                                                                                                                            <  1:65063/151‑1 (MQ=255)       cGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACtacgcaggcacgcgc                                                                                                                                         <  1:439202‑M1/151‑16 (MQ=255)         aCGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGcc                                                                                                                                        <  2:62464/150‑1 (MQ=255)         aCGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGcc                                                                                                                                        <  2:491868/150‑1 (MQ=255)            gTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCg                                                                                                                                    >  2:478055/1‑151 (MQ=255)             tGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGt                                                                                                                                   <  1:607212/151‑1 (MQ=255)               ccGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAg                                                                                                                                 >  2:553474‑M2/1‑150 (MQ=255)               ccGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAg                                                                                                                                 >  2:604493‑M2/1‑150 (MQ=255)                         gggCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAgtgcggatcgg                                                                                                                       <  1:438920‑M2/151‑12 (MQ=255)                         gggCATTAGCCCCCGCGTGACGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGAcgc                                                                                                                        >  2:19475/1‑150 (MQ=255)                            cATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAgtgcggatcggtta                                                                                                                    >  2:620496‑M2/1‑137 (MQ=255)                                     ccGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCAt                                                                                                           <  2:203015/151‑1 (MQ=255)                                          tGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGt                                                                                                      >  1:173900/1‑151 (MQ=255)                                            aaGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTg                                                                                                     <  2:524669/150‑1 (MQ=255)                                               ccAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGAt                                                                                                   <  1:421793/149‑1 (MQ=255)                                                   gACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATtgtg                                                                                               <  2:363713/149‑1 (MQ=255)                                                     cGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGa                                                                                           >  2:78493/1‑151 (MQ=255)                                                        cAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTA‑CCGATCCGCACTACGCAGGCACGCGCCGAGGTTgtgcggatcggttacgccgtcgaggatcagcaggaacggttga                                                                                       <  2:438911‑M2/151‑44 (MQ=255)                                                                ggAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTc                                                                                 >  2:605558/1‑150 (MQ=255)                                                                 gAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTc                                                                                 <  1:460872/149‑1 (MQ=255)                                                                    aaCGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGcac                                                                            <  1:477492/151‑1 (MQ=255)                                                                    aaCGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGcac                                                                            <  2:259530/151‑1 (MQ=255)                                                                    aaCGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGcac                                                                            <  2:345837/151‑1 (MQ=255)                                                                      cGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAg                                                                          >  2:366335/1‑151 (MQ=255)                                                                         tCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCt                                                                        >  1:19410/1‑150 (MQ=255)                                                                         tCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTc                                                                       <  1:54134/151‑1 (MQ=255)                                                                            gCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCaa                                                                     >  1:368/1‑150 (MQ=255)                                                                             ccGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACTACGCAGGCACGCGCCGAGGTTgtgcggatcggttacgccgtcgaggatcagcaggaacggttgatcgagcgaagcgatcagatc                                                                   >  1:333227‑M2/1‑88 (MQ=255)                                                                             ccGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAAc                                                                    <  2:529662/150‑1 (MQ=255)                                                                             ccGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAAc                                                                    <  2:532141/150‑1 (MQ=255)                                                                              cGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGc                                                                  <  1:466053/151‑1 (MQ=255)                                                                              cGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGc                                                                  <  2:166753/151‑1 (MQ=255)                                                                                gATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACaacctcggcgcgtgcctgcgtagtgcggatcggttacgccgtcgaggatcagcaggaacggttgatcgagcgaagcgatcagatccgg                                                                >  1:333227‑M1/1‑63 (MQ=255)                                                                                  tCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACg                                                              >  1:506893/1‑151 (MQ=255)                                                                                             cgctCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGt                                                   <  2:536703/151‑1 (MQ=255)                                                                                                    tCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACaacctcggcgcgtgcctgcgtagtgcggatcggtacgccgtcgaggatcagcaggaacggttgatcgagcgaagcgatcagatccggcagatcgttttcctgatactg                                            <  2:438911‑M1/151‑109 (MQ=255)                                                                                                                 ctgATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTgcggcgtccgcact                                <  1:454131/150‑9 (MQ=255)                                                                                                                    aTCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCg                            >  2:283446/1‑151 (MQ=255)                                                                                                                    aTCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCg                            >  1:153392/1‑151 (MQ=255)                                                                                                                            cGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCg                       >  1:108864/1‑148 (MQ=255)                                                                                                                                 tAACCGATCCGCACtacgcaggcacgcgccgaggttgtgcggatcggttacgccgtcgaggatcagcaggaacggttgatcgagcgaagcgatcagatccggcagatcgttttcctgatactgacgtcctggcttcacgcgggcgataatg               >  2:620496‑M1/1‑14 (MQ=255)                                                                                                                                 tAACCGATCCGCACAACCTCGGCGCGTGCCTGCGCAGTGCGGACCCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGAGAGAAAAATAGCCTGCGGCGCGGCAGAAAGCGATCCGCTGATTCg               >  2:533692/1‑151 (MQ=255)                                                                                                                                   aCCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGAGGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCggg             <  1:177253/151‑1 (MQ=255)                                                                                                                                    ccGATCCGCACtacgcaggcacgcgccgaggttgtgcggatcggttacgccgtcgaggatcagcaggaacggttgatcgagcgaagcgatcagatccggcagatcgttttcctgatactgacgtcctggcttcacgcgggcgataatgccc            <  1:438920‑M1/151‑141 (MQ=255)                                                                                                                                            cacaACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAAc     <  1:90354/150‑1 (MQ=255)                                                                                                                                              ctacgcaggcacgcgccgaggttgtgcggatcggttacgccgtcgaggatcagcaggaacggttgatcgagcgaagcgatcagatccggcagatcgttttcctgatactgacgtcctggcttcacgcgggcgataatgccctgatgcacgg  >  2:553474‑M1/1‑1 (MQ=255)                                                                                                                                              ctacgcaggcacgcgccgaggttgtgcggatcggttacgccgtcgaggatcagcaggaacggttgatcgagcgaagcgatcagatccggcagatcgttttcctgatactgacgtcctggcttcacgcgggcgataatgccctgatgcacgg  >  2:604493‑M1/1‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                       GAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTCGATCAACCGTTCCTGCTGATCCTCGACGGCGTAACCGATCCGCACAACCTCGGCGCGTGCCTGCGTAGTGCGGACGCCGCAGGCGTTCATGCGGTGATTGTGCCGAAAGATCGCTCCGCACAGCTCAACGCTACGGCGAAAAAAGTAGCCTGCGGCGCGGCAGAAAGCGTTCCGCTGATTCGGGTGACTAACCTT  >  NC_000913/4409454‑4409745

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