New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	= 4532507	3 (0.140)	8 (0.420)	6/148	1.0	57.8%	intergenic (+57/‑425)	fimE/fimA	type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA
	?	NZ_CP009273	= 4532803	9 (0.480)					intergenic (+353/‑129)	fimE/fimA	type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA

AATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAA  >  NZ_CP009273/4532724‑4532873                                                                                |                                                                       aataaaGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaat                                                              >  1:131807‑M2/1‑80 (MQ=255) aataaaGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaat                                                              >  1:259704‑M2/1‑80 (MQ=255) aataaaGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCaaa                                                              >  1:483368/1‑90 (MQ=255) aataaaGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCaaa                                                              <  1:541386/90‑1 (MQ=255)              aTAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCtttt                                                 >  2:191304/1‑90 (MQ=255)                 gAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggg                                              >  1:282606/1‑90 (MQ=255)                      ttAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggcc                                                                   >  2:363566‑M2/1‑59 (MQ=255)                      ttAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggcc                                                                   <  1:363566‑M2/64‑6 (MQ=255)                      ttAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaa                                         >  1:367095/1‑90 (MQ=255)                        aaGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttt                                       <  2:494755‑M2/90‑34 (MQ=255)                        aaGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttt                                       <  1:81535‑M2/90‑34 (MQ=255)                           ccAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaa                                    <  2:352377‑M2/90‑37 (MQ=255)                             aacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGc                                  <  2:367095/90‑1 (MQ=255)                             aacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGc                                  >  1:147100/1‑90 (MQ=255)                              acaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCa                                 <  1:191304/90‑1 (MQ=255)                                   tAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttatta                            >  1:352377‑M2/1‑46 (MQ=255)                                                 ttGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaataagtt                  >  1:392887‑M2/1‑32 (MQ=255)                                                 ttGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaataagtt                  <  2:392887‑M2/86‑55 (MQ=255)                                                             cTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaacaa  <  2:147100/90‑1 (MQ=255)                                                                                |                                                                       AATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAA  >  NZ_CP009273/4532724‑4532873

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