New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	4532499 =	3 (0.330)	7 (0.850)	6/154	0.3	64.6%	intergenic (+49/‑433)	fimE/fimA	type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA
	?	NZ_CP009273	4532813 =	5 (0.610)					intergenic (+363/‑119)	fimE/fimA	type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA

AATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTG  >  NZ_CP009273/4532724‑4532890                                                                                          |                                                                              acagttagatgctttttgtcgttttttaatatttttatgcttgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaa                                                                               <  2:111485‑M2/1‑1 (MQ=255) acagttagatgctttttgtagttttttaatatttttatgcttgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaa                                                                               <  1:183892‑M2/1‑1 (MQ=255)     aaGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAATTGtc                                                                           >  2:254223/1‑89 (MQ=255)        aTAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAATTGtcttg                                                                        <  1:64068/90‑5 (MQ=255)                     atTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaa                                                           >  1:109446/1‑90 (MQ=255)                             aacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGc                                                   >  1:224806/1‑90 (MQ=255)                                  tttatgcttgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAgtgt                                              >  2:132257‑M2/57‑90 (MQ=255)                                                 ttGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCg                               >  1:212242/1‑90 (MQ=255)                                                     acgtaacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTa                           >  1:237619‑M2/38‑90 (MQ=255)                                                           cttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaac                     <  2:218519‑M2/59‑1 (MQ=255)                                                                           acagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGa     <  2:237619‑M2/75‑1 (MQ=255)                                                                              gtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTg  >  2:37734‑M2/13‑90 (MQ=255)                                                                                          |                                                                              AATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTG  >  NZ_CP009273/4532724‑4532890

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