New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3557740	19 (0.580)	5 (0.160)	5/230	3.3	22.3%	coding (340/1227 nt)	rtcB	RNA‑splicing ligase
	?	NC_000913	= 3557761	17 (0.550)					coding (319/1227 nt)	rtcB	RNA‑splicing ligase
Rejected: Coverage evenness skew score above cutoff.

ATAGTTATTGGTATTCAGGAAACGGGGATATTTTTGCGTTAACCACTGATAACCGGCTTCAAGCTCAGCCCATTTAGCATCGACGTTAACAGGTGGATTTTCCCAGGCACCTTTATCACGTTTACAACGGCCAGTGGTACGCCCGTGCGGCACGGCCGTTTCAATCGCCTGACGCAGCTCTGCCAGGTTTTCAGGCAGGTCTTCCGCCGTTAACGCGGTACGCAGCGCGTTCATTCCACAGCCAATATCCACGC  >  NC_000913/3557609‑3557862                                                                                                                                    |                                                                                                                           ataGTTATTGGTATTCAGGAAACGGGGATATTTTTGCGTTAACCACTGATAACCGGCTTCAAGCTCAGCCCATTTAGCATCGACGTTAACAGGTGGATTTTCCCAGGCACCTTTATCACGTTTACAACGGCCAGTGGTACGCCCGTGc                                                                                                            <  8:113294/148‑1 (MQ=255)            tATTCAGGAAACGGGGATATTTTTGCGTTAACCACTGATAACCGGCTTCAAGCTCAGCCCATTTAGCATCGACGTTAACAGGTGGATTTTCCCAGGCACCTTTATCACGTTTACAACGGCCAGTGGTACGCCCGTGCGGCACGGCCGtt                                                                                                >  2:106184/1‑149 (MQ=255)                          ggATATTTTTGCGTTAACCACTGATAACCGGCTTCAAGCTCAGCCCATTTAGCATCGACGTTAACAGGTGGATTTTCCCAGGCACCTTTATCACGTTTACAACGGCCAGTGGTACGCCCGTGc                                                                                                            >  8:46012/1‑123 (MQ=255)                          ggATATTTTTGCGTTAACCACTGATAACCGGCTTCAAGCTCAGCCCATTTAGCATCGACGTTAACAGGTGGATTTTCCCAGGCACCTTTATCACGTTTACAACGGCCAGTGGTACGCCCGTGc                                                                                                            <  7:46012/123‑1 (MQ=255)                                     cGTTAACCACTGATAACCGGCTTCAAGCTCAGCCCATTTAGCATCGACGTTAACAGGTGGATTTTCCCAGGCACCTTTATCACGTTTACAACGGCCAGTGGTACGCCCGTGCGGCACGGCCGTTTCAATCGCCTGACGCAGCTCTGCCa                                                                       <  1:106184/149‑1 (MQ=255)                                                         cTTCAAGCTCAGCCCATTTAGCATCGACGTTAACAGGTGGATTTTCCCAGGCACCTTTATCACGTTTACAACGGCCAGTGGTACGCCCGTGCGGCACGGCCGTTTCAATCGCCTGACGCAGCTCTGCCAGGTTTTCAGGCAGGTCTTcc                                                   <  4:73752/149‑1 (MQ=255)                                                         cTTCAAGCTCAGCCCATTTAGCATCGACGTTAACAGGTGGATTTTCCCAGGCACCTTTATCACGTTTACAACGGCCAGTGGTACGCCCGTGCGGCACGGCCGTTTCAATCGCCTGACGCAGCTCTGCCAGGTTTTCAGGCAGGTCTTcc                                                   >  3:73752/1‑149 (MQ=255)                                                             aaGCTCAGCCCATTTAGCATCGACGTTAACAGGTGGATTTTCCCAGGCACCTTTATCACGTTTACAACGGCCAGTGGTACGCCCGTGCGGCACGGCCGTTTCAATCGCCTGACGCAGCTCTGCCAGGTTTTCAGGCAGGTCTTccgccg                                               >  1:59974/1‑149 (MQ=255)                                                                                               ggATTTTCCCAGGCACCTTTATCACGTTTACAACGGCCAGTGGTACGCCCGTGCGGCACGGCCGTTTCAATCGCCTGACGCAGCTCTGCCAGGTTTTcaggcagg                                                         >  2:31184/1‑105 (MQ=255)                                                                                               ggATTTTCCCAGGCACCTTTATCACGTTTACAACGGCCAGTGGTACGCCCGTGCGGCACGGCCGTTTCAATCGCCTGACGCAGCTCTGCCAGGTTTTcaggcagg                                                         <  1:31184/105‑1 (MQ=255)                                                                                                        cAGGCACCTTTATCACGTTTACAACGGCCAGTGGTACGCCCGTGCGGCACGGCCGTTTCAATCGCCTGACGCAGCTCTGCCAGGTTTTCAGGCAGGTCTTCCGCCGTTAACGCGGTACGCa                                <  5:115045/121‑1 (MQ=255)                                                                                                        cAGGCACCTTTATCACGTTTACAACGGCCAGTGGTACGCCCGTGCGGCACGGCCGTTTCAATCGCCTGACGCAGCTCTGCCAGGTTTTCAGGCAGGTCTTCCGCCGTTAACGCGGTACGCa                                >  6:115045/1‑121 (MQ=255)                                                                                                          ggCACCTTTATCACGTTTACAACGGCCAGTGGTACGCCCGTGCGGCACGGCCGTTTCAATCGCCTGACGCAGCTCTGCCAGGTTTTCAGGCAGGTCTTCCGCCGTTAACGCGGTACGCAGCGCGTTAATTCCACAGCCAATATCCACGc  <  2:59974/149‑1 (MQ=255)                                                                                                                    tCACGTTTACAACGGCCAGTGGTACGCCCGTGCGGCACGGCCGTTTCAATCGCCTGACGCAGCTCTGCCAGGTTTTcaggcag                                                          >  3:59933/1‑83 (MQ=255)                                                                                                                    tCACGTTTACAACGGCCAGTGGTACGCCCGTGCGGCACGGCCGTTTCAATCGCCTGACGCAGCTCTGCCAGGTTTTcaggcag                                                          <  4:59933/83‑1 (MQ=255)                                                                                                                             ccccGGCCAGTGGTACGCCCGTGCGGCACGGCCGTTTCAATCGCCTGACGCAGCTCTGCCAGGTTTTCAGGCAGGGCTTCCGCCGTTAACGCGGTACGCAGCGCGtt                         <  6:9702/104‑1 (MQ=255)                                                                                                                             cAACGGCCAGTGGTACGCCCGTGCGGCACGGCCGTTTCAATCGCCTGACGCAGCTCTGCCAGGTTTTCAGGCAGGTCTTCCGCCGTTAACGCGGTACGCAGCGCGtt                         >  5:7406/1‑107 (MQ=255)                                                                                                                             cAACGGCCAGTGGTACGCCCGTGCGGCACGGCCGTTTCAATCGCCTGACGCAGCTCTGCCAGGTTTTCAGGCAGGTCTTCCGCCGTTAACGCGGTACGCAGCGCGtt                         <  6:7406/107‑1 (MQ=255)                                                                                                                             cAACGGCCAGTGGTACGCCCGTGCGGCACGGCCGTGTCAATCGCCTGACGCAGCTCTGCCAGGTTTTCAGGCAGGTCTTCCGCCGTTAACTCGGTACGCAGCGCGtt                         >  5:9702/1‑107 (MQ=255)                                                                                                                                    |                                                                                                                           ATAGTTATTGGTATTCAGGAAACGGGGATATTTTTGCGTTAACCACTGATAACCGGCTTCAAGCTCAGCCCATTTAGCATCGACGTTAACAGGTGGATTTTCCCAGGCACCTTTATCACGTTTACAACGGCCAGTGGTACGCCCGTGCGGCACGGCCGTTTCAATCGCCTGACGCAGCTCTGCCAGGTTTTCAGGCAGGTCTTCCGCCGTTAACGCGGTACGCAGCGCGTTCATTCCACAGCCAATATCCACGC  >  NC_000913/3557609‑3557862

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