New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_012967	717135 =	27 (0.090)	18 (0.070)	12/124	14.4	32.3%	coding (62/312 nt)	ECB_RS23290	DUF6258 family protein
	?	NC_012967	717390 =	52 (0.210)					intergenic (+5/‑328)	ECB_RS23290/ECB_RS03505	DUF6258 family protein/ISAs1 family transposase
Rejected: Coverage evenness skew score above cutoff.

ATTTAATCTCAAGGATAGCCTATGGTACAACAGAATGGATTTTTTATGATAATGAAGACTTAGAAGATGGTTATTTGGTATTTTCAGATGTTAGTTTCTTCAATATAACTCCTAATGGCTCTATTCCAGATGATTATA  >  NC_012967/717069‑717206                                                                   |                                                                        atttaatCTCAAGGATAGCCTATGGTACAACAGAATGGATTTTTTATGATAATGAAGACTTAGAAGATGGtt                                                                    >  1:2862639/1‑72 (MQ=255)       tctcAAGGATAGCCTATGGTACAACAGAATGGATTTTTTATGATAATGAAGAGTTAAAAGATGGTTATTTgg                                                              <  2:6042871/72‑1 (MQ=255)        ctcAAGGATAGCCTATGGTACAACAGAATGGATTTTTTATGATAATGAAGACTTAGAAGATGGTTATTTGGt                                                             >  2:6494843/1‑72 (MQ=255)          cAAGGATAGCCTATGGTACAACAGAATGGATTTTTTATGATAATGAAGACTTAGAAGATGGTTATTTGGTAt                                                           >  1:637207/1‑72 (MQ=255)           aaGGATAGCCTATGGTACAACAGAATGGATTTTTTATGATAATGAAGACTTAGAAGATGGTTATTTGGTAtt                                                          >  2:8453551/1‑72 (MQ=255)            aGGATAGCCTATGGTACAACAGAATGGATTTTTTATGATAATGAAGACTTAGAAGATGGTTATTTGGTAttt                                                         >  2:5675162/1‑72 (MQ=255)             ggATATCCTATGGCACAACAGAATGGATTTTTTATGATAATGAAGACTTAGAAGATGGTTATTTGGTAtttt                                                        >  1:3131329/1‑72 (MQ=255)             ggATAGCCTATGGTACAACAGAATGGATTTTTTATGATAATGAAGACTTAGAAGATGGTTATTTGGTAtttt                                                        <  2:5157757/72‑1 (MQ=255)                 aGCCTATGGTACAACAGAATGGATTTTTTATGATAATGAAGACTTAGAAGATGGTTATTTGGTATTTTccga                                                    <  2:2378826/72‑4 (MQ=255)                 aGCCTATGGTACAACAGAATGGATTTTTTATGATAATGAAGACTTAGAAGATGGTTATTTGGTATTTTCAGa                                                    >  1:5815478/1‑72 (MQ=255)                    cTATGGTACAACAGAATGGATTTTTTATGATAATGAAGACTTAGAAGATGGTTATTTGGTATTTTCCGAtac                                                 <  2:4754649/72‑3 (MQ=255)                    cTATGGTACAACAGAATGGATTTTTTATGATAATGAAGACTTAGAAGATGGTTATTTGGTATTTTCAGATGt                                                 >  2:4880728/1‑72 (MQ=255)                       tGGTATAACAGAATGGATTTTTTATGATAATGAAGACTTAGAAGATGGTTATTTGGTATTTTCAGATGTTAg                                              <  2:1905683/72‑1 (MQ=255)                               caGAATGGGTTTTTTATGATAATGAAGACTTAGAAGATGGTTATTTGGTATTTTCAGATGTTAGTTTCTTCa                                      >  1:5951314/1‑72 (MQ=255)                                aGAATGGATTTTTTATGATAATGAAGACTTAGAAGATGGTTATTTGGTATTTTCAGATGTTAGTTTCTTCaa                                     >  1:7643376/1‑72 (MQ=255)                                     ggATTTTTTATGATAATGAAGACTTAGAAGATGGTTATTTGGTATTTTCAGATGTTAGTTTCTTCAATATaa                                >  1:5419838/1‑72 (MQ=255)                                     ggATTTTTTATGATAATGAAGACTTAGAAGATGGTTATTTGGTATTTTCAGATGTTAGTTTCTTCAATATaa                                >  2:8498833/1‑72 (MQ=255)                                              aTGATAATGAAGACTTAGAAGATGGTTATTTGGTATTTTCAGATGTTAGTTTCTTCAATATAACTCCTAATg                       <  1:339020/72‑1 (MQ=255)                                                gATAATGAAGACTTAGAAGATGGTTATTTGGTATTTTCAGATGTTAGTTTCTTCAATATAACTCCTAATGGc                     <  1:2049202/72‑1 (MQ=255)                                                      gAAGACTTAGAAGATGGTTATTTGGTATTTTCAGATGTTAGTTTCTTCAATATAACTCCTAATGGCTCTAtt               <  2:5469298/72‑1 (MQ=255)                                                       aaGACTTAGAAGATGGTTATTTGGTATTTTCAGATGTTAGTTTCTTCAATATAACTCCTAATGGCTCTATTc              >  1:4999109/1‑72 (MQ=255)                                                        aGACTTAGAAGATGGTTATTTGGTATTTTCAGATGTTAGTTTCTTCAATATAACTCCTAATGGCTCTATTcc             >  1:4584880/1‑72 (MQ=255)                                                          aCTTAGAAGATGGTTATTTGGTATTTTCAGATGTTAGTTTCTTCAATATAACTCCTAATGGCTCTATTCCAg           >  1:3265323/1‑72 (MQ=255)                                                              agaagaTGGTTATTTGGTATTTTCAGATGTTAGTTTCTTCAATATAACTCCTAATGGCTCTATTCCAgatga       >  1:8847194/1‑72 (MQ=255)                                                                aagaTGGTTATTTGGTATTTTCAGATGTTAGTTTCTTCAATATAACTCCTAATGGCTCTATTCCAGATGAtt     >  1:6852429/1‑72 (MQ=255)                                                                aagaTGGTTATTTGGTATTTTCAGATGTTAGTTTCTTCAATATAACTCCTAATGGCTCTATTCCAGATGAtt     <  1:3996266/72‑1 (MQ=255)                                                                aacaTGGTTATTTGGTATTTTCAGATGTTAGTTTCTTCAATATAACTCCTAATGGCTCTATTCCAGATGAtt     >  2:9454321/4‑72 (MQ=255)                                                                aacaTGGTTATTTGGTATTTTCAGATGTTAGTTTCTTCAATATAACTCCTAATGGCTCTATTCCAGATGAtt     >  2:9934696/4‑72 (MQ=255)                                                                aacaTGGTTATTTGGTATTTTCAGATGTTAGTTTCTTCAATATAACTCCTAATGGCTCTATTCCAGATGAtt     >  2:109014/4‑72 (MQ=255)                                                                 agaTGGTTATTTGGTATTTTCAGATGTTAGTTTCTTCAATATAACTCCTAATGGCTCTATTCCAGATGATta    >  1:6083571/1‑72 (MQ=255)                                                                 acaTGGTTATTTGGTATTTTCAGATGTTAGTTTCTTCAATATAACTCCTAATGGCTCTATTCCAGATGATta    >  2:2986728/3‑72 (MQ=255)                                                                 acaTGGTTATTTGGTATTTTCAGATGTTAGTTTCTTCAATATAACTCCTAATGGCTCTATTCCAGATGATta    >  2:9480674/3‑72 (MQ=255)                                                                   aTGGTTATTTGGTATTTTCAGATGTTAGTTTCTTCAATATAACTCCTAATGGCTCTATTCCAGATGATtata  >  1:9635407/1‑72 (MQ=255)                                                                   |                                                                        ATTTAATCTCAAGGATAGCCTATGGTACAACAGAATGGATTTTTTATGATAATGAAGACTTAGAAGATGGTTATTTGGTATTTTCAGATGTTAGTTTCTTCAATATAACTCCTAATGGCTCTATTCCAGATGATTATA  >  NC_012967/717069‑717206

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