New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	CP009273	1413126 =	10 (0.580)	6 (0.360)	3/292	NT	40.5%	coding (199/489 nt)	sieB	Rac prophage; phage superinfection exclusion protein
	?	CP009273	1413206 =	8 (0.480)					coding (279/489 nt)	sieB	Rac prophage; phage superinfection exclusion protein

AACTTTTTGCTAACGAGCCACTCGAAAGACTTATGTATACGATTATCATTTTTGGTCTCACTCTCTGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATACTGAAATACCTTGGCTCTTTCAGATTATCGTTTTTGCCTTTTCTTTCGTGGTCGCCATTTCCTTCTCAAGATTGCGAGCACATATTCAAAAGCATTATTCATTACTACCAGAGCAACGAGTATTGCTTCGTTTATCTGAGAAAGAAATCGCTGTATTTAAAGATTTCCTTAAAACAGGAAATCTTATTATCACTTCTCCTTGCCGTAACCCGGTTATGAAAAAATTAGAACGGAAGGGCATCATTCAACATCAGAGTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAAT  >  CP009273/1412956‑1413410                                                                                                                                                                                                                                                           |                                                                                                                                                                                                             aaCTTTTTGCTAACGAGCCACTCGAAAGACTTATGTATACGATTATCATTTTTGGTCTCACTCTCTGGCTGATACCGAAAGAGTTTACTGTNGCATTCAATGCTTATACTGTAATACCTTGGCTCTTTCAGATTATCGTTTTTGCCTTTTCTTTCGTGGTCGCCATTTCCTTCTCAAGATTGCGAGCACATATTCAAAAGCATTATTCATTACTACCAGAGCAACGAGTATTGCTTCGTTTATCTGAGAAAGAAATCGCTGTATTTAAAGATTTc                                                                                                                                                                                      <  1:255143/275‑1 (MQ=255)                                           ttATCATTTTTGGTCTCACTCTCTGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATACTGAAATACCTTGGCTCTTTCAGATTATCGTTTTTGCCTTTTCTTTCGTGGTCGCCATTTCCTTCTCAAGATTGCGAGCACATATTCAACAGCATTATTCATTACTACCAGAGCCACGAGTATTGCTTCGTTTATCTGAGAAAGAAATCGCTGTATTTAAATCTTTCCTTAAAACAGGAAATCTTATTATCACTTCTCCTTTCCGtca                                                                                                                                             >  1:82378/1‑272 (MQ=255)                                                                          aCCGAAAGAGTTTACTGTCGCATTCAATGCTTATACTGAAATACCTTGGCTCTTTCAGATTATCGTTTTTGCCTTTTCTTTCGTGGTCGCCATTTCCTTCTCAAGATTGCGAGCACATATTCAAAAGCATTATTCATTACTACCAGAGCAACGAGTATTGCTTCGTTTATCTGAGAAAGAAATCGCTGTATTTAAAGATTTCCTTAAAACAGGAAATCTTATTATCACTTCTACTTGCCGTAACCCGGTTATGAAAAAATTAGAACGAAAGGGca                                                                                                             >  1:196771/1‑275 (MQ=255)                                                                                                                                                                                      gagcgCACATATTAAAAAGCAGTATTAATTACTACCAGAGCAACGAGTATTGCGTCGTTTATCTGAGAAAGAAATCGCTGTATTTAAAGATTTCCTTAAAAAAGGAAATATTATTATCACTTCTCCTTGCCGTAACCCGGTTATGAAAAAATTAGAACGGAAGGGCATCATTCAACATCAGAGTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAAt  <  2:82376/270‑1 (MQ=255)                                                                                                                                                                                            caatactcgttgctctggtagtaatgaataatgcttttgaatatgtgctcgcaatcttgagaaAGAAATCGCTGTATTTAAAGa                                                                                                                                                                                          >  2:38653‑M2/64‑84 (MQ=255)                                                                                                                                                                                            caatactcgttgctctggtagtaatgaataatgcttttgaatatgtgctcgcaatcttgagaaAGAAATCGCTGTATTTAAAGATTTCCTTAAAACAGGAAATCTTATTATCACTTCTCCTTGCCGTAACCCGGTTATGAAAAAAt                                                                                                                            >  2:36697‑M2/64‑146 (MQ=255)                                                                                                                                                                                            caatactcgttgctctggtagtaatgaataatgcttttgaatatgtgctcgcaatcttgagaaAGAAATCGCTGTATTTAAAGATTTCCTTAAAACAGGAAATCTTATTATCACTTCTCCTTGCCGTAACCCGGTTATGAAAAAAt                                                                                                                            >  2:49730‑M2/64‑146 (MQ=255)                                                                                                                                                                                                actcgttgctctggtagtaatgaataatgcttttgaatatgtgctcgcaatcttgagaaAGAAATCGCTGTATTTAAAGATTTc                                                                                                                                                                                      <  1:38653‑M2/25‑1 (MQ=255)                                                                                                                                                                                                actcgttgctctggtagtaatgaataatgcttttgaatatgtgctcgcaatcttgagaaAGAAATCGCTGTATTTAAAGATTTCCTTAAAACAGGAAATCTTATTATCACTTCTCCTTGCCGTAACCCGGTTATGAAAAAATTAGa                                                                                                                        <  1:36697‑M2/87‑1 (MQ=255)                                                                                                                                                                                                actcgttgctctggtagtaatgaataatgcttttgaatatgtgctcgcaatcttgagaaAGAAATCGCTGTATTTAAAGATTTCCTTAAAACAGGAAATCTTATTATCACTTCTCCTTGCCGTAACCCGGTTATGAAAAAATTAGa                                                                                                                        <  1:49731‑M2/87‑1 (MQ=255)                                                                                                                                                                                                            tattCATTACTACCAGAGCAACGAGTATTGCTTCGTTTATCTGAGAAAGAAATCGCTGTATTTAAAGATTTCCTTAAAACAgg                                                                                                                                                                           >  1:202965/1‑83 (MQ=255)                                                                                                                                                                                                                cATTACTACCAGAGCAACGAGTATTGCTTCGTTTATCTGAGAAAGAAATCGCTGTATTTAAAGATTTCCTTAAAACAGGAAAt                                                                                                                                                                       <  2:202959/83‑1 (MQ=255)                                                                                                                                                                                                                                             cGTTTATCTGAGAAAGAAATCGCTGTATTTAAAGATTTCCTTAAAAc                                                                                                                                                                              >  1:74802/1‑47 (MQ=255)                                                                                                                                                                                                                                             cGTTTATCTGAGAAAGAAATCGCTGTATTTAAAGATTTCCTTAAAACAGGAAATc                                                                                                                                                                      >  1:179764/1‑55 (MQ=255)                                                                                                                                                                                                                                                 tATCTGAGAAAGAAATCGCTGTATTTAAAGATTTCCTTAAAACAGGa                                                                                                                                                                          <  2:74800/47‑1 (MQ=255)                                                                                                                                                                                                                                                 tATCTGAGAAAGAAATCGCTGTATTTAAAGATTTCCTTAAAACAGGAAATCttat                                                                                                                                                                  <  2:179760/55‑1 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                             AACTTTTTGCTAACGAGCCACTCGAAAGACTTATGTATACGATTATCATTTTTGGTCTCACTCTCTGGCTGATACCGAAAGAGTTTACTGTCGCATTCAATGCTTATACTGAAATACCTTGGCTCTTTCAGATTATCGTTTTTGCCTTTTCTTTCGTGGTCGCCATTTCCTTCTCAAGATTGCGAGCACATATTCAAAAGCATTATTCATTACTACCAGAGCAACGAGTATTGCTTCGTTTATCTGAGAAAGAAATCGCTGTATTTAAAGATTTCCTTAAAACAGGAAATCTTATTATCACTTCTCCTTGCCGTAACCCGGTTATGAAAAAATTAGAACGGAAGGGCATCATTCAACATCAGAGTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAAT  >  CP009273/1412956‑1413410

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