New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	1203246 =	15 (0.480)	19 (0.760)	15/136	0.4	71.5%	coding (290/630 nt)	BW25113_RS05990	hypothetical protein
	?	NZ_CP009273	1205075 =	3 (0.120)					coding (40/540 nt)	BW25113_RS06005	phage tail protein

GCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAC  >  NZ_CP009273/1204989‑1205159                                                                                       |                                                                                     accaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGtt                                                                                   >  2:397475‑M2/87‑90 (MQ=255) accaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGtt                                                                                   >  2:239889‑M2/87‑90 (MQ=255)                      tatgaaaccggtgcgtcataatgtagcggttgctgcattgagttttgctccttctcccaaaccgaCGTTTATGAAAATGAAGAAATAACa                                                              >  1:418594‑M2/66‑90 (MQ=255)                       atgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                                             >  1:198525‑M2/65‑90 (MQ=255)                         gaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGc                                                           <  1:880476‑M2/28‑1 (MQ=255)                         gaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGc                                                           >  1:739157‑M2/63‑90 (MQ=255)                             ccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAAt                                                       <  2:189035‑M2/32‑1 (MQ=255)                               ggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATgg                                                     >  1:415408‑M2/57‑90 (MQ=255)                                       aCGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc                                             >  2:604946/1‑90 (MQ=255)                                        taatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                                            <  1:239889‑M2/43‑1 (MQ=255)                                        taatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                                            <  2:415408‑M2/43‑1 (MQ=255)                                            gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaa                                                                            >  1:152216‑M2/44‑54 (MQ=255)                                            gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaa                                                                            <  2:152216‑M2/11‑1 (MQ=255)                                               gccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtttt                                     >  1:767639‑M2/41‑90 (MQ=255)                                               gccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtttt                                     >  1:105353‑M2/41‑90 (MQ=255)                                                  gttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAcc                                  <  2:767639‑M2/53‑1 (MQ=255)                                                     gctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgg                               <  2:739157‑M2/56‑1 (MQ=255)                                                      ctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                              <  1:753975‑M2/57‑1 (MQ=255)                                                      ctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                              <  1:65847‑M2/57‑1 (MQ=255)                                                                 tttgccccttctcccaaaccaaCGTTTATGTAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACAt                   <  1:455383‑M2/68‑1 (MQ=255)                                                                    gagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                                            >  1:551849/1‑62 (MQ=255)                                                                    gagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                                            <  2:551849/62‑1 (MQ=255)                                                                       ccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAt             <  1:518633‑M2/74‑1 (MQ=255)                                                                       ccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAt             <  1:264654‑M2/74‑1 (MQ=255)                                                                                  accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAc  <  1:289784‑M2/85‑1 (MQ=255)                                                                                       |                                                                                     GCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAC  >  NZ_CP009273/1204989‑1205159

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