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.

TATTCAGGAAACGGGGATATTTTTGCGTTAACCACTGATAACCGGCTTCAAGCTCAGCCCATTTAGCATCGACGTTAACAGGTGGATTTTCCCAGGCACCTTTATCACGTTTACAACGGCCAGTGGTACGCCCGTGCGGCACGGCCGTTTCAATCGCCTGACGCAGCTCTGCCAGGTTTTCAGGCAGGTCTTCCGCCGTTAACGCGGTACGCAGCGCGTTCATTCCACAGCCAATATCCACGCCCACCGCCGCCGGAATAATCGCCCCTTTGGTCGGGA  >  NC_000913/3557620‑3557898                                                                                                                                              |                                                                                                                                          tATTCAGGAAACGGGGATATTTTTGCGTTAACCACTGATAACCGGCTTCAAGCTCAGCCCATTTAGCATCGACGTTAACAGGTGGATTTTCCCAGGCACCTTTATCACGTTTACAACGGCCAGTGGTACGCCCGTGCGGCACGGCCGtt                                                                                                                                    >  2:106184/1‑149 (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                                                                                             <  1:31184/105‑1 (MQ=255)                                                                                    ggATTTTCCCAGGCACCTTTATCACGTTTACAACGGCCAGTGGTACGCCCGTGCGGCACGGCCGTTTCAATCGCCTGACGCAGCTCTGCCAGGTTTTcaggcagg                                                                                             >  2:31184/1‑105 (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)                                                                                                                                             aCGGCCGTTTCAATCGCCTGACGCAGCTCTGCCAGGTTTTCAGGCAGGTCTTCCGCCGTTAACGCGGTACGCAGCGCGTTCATTCCACAGCCAATATCCACGCCCACCGCCGCCGGAATAATCGCCCCTTTGGTCGGGa  <  5:63926/139‑1 (MQ=255)                                                                                                                                             aCGGCCGTTTCAATCGCCTGACGCAGCTCTGCCAGGTTTTCAGGCAGGTCTTCCGCCGTTAACGCGGTACGCAGCGCGTTCATTCCACAGCCAATATCCACGCCCACCGCCGCCGGAATAATCGCCCCTTTGGTCGGGa  >  6:63926/1‑139 (MQ=255)                                                                                                                                              |                                                                                                                                          TATTCAGGAAACGGGGATATTTTTGCGTTAACCACTGATAACCGGCTTCAAGCTCAGCCCATTTAGCATCGACGTTAACAGGTGGATTTTCCCAGGCACCTTTATCACGTTTACAACGGCCAGTGGTACGCCCGTGCGGCACGGCCGTTTCAATCGCCTGACGCAGCTCTGCCAGGTTTTCAGGCAGGTCTTCCGCCGTTAACGCGGTACGCAGCGCGTTCATTCCACAGCCAATATCCACGCCCACCGCCGCCGGAATAATCGCCCCTTTGGTCGGGA  >  NC_000913/3557620‑3557898

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