New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NZ_CP009273	4532499 =	3 (0.110)	13 (0.530)	7/150	1.1	56.9%	intergenic (+49/‑433)	fimE/fimA	type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA
	?	NZ_CP009273	4532813 =	17 (0.700)					intergenic (+363/‑119)	fimE/fimA	type 1 fimbria switch DNA invertase FimE/type 1 fimbrial major subunit FimA

ATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTG  >  NZ_CP009273/4532731‑4532890                                                                                   |                                                                              aTAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCa                                                                        >  1:128874/1‑90 (MQ=255)       ttttgtcgttttttaatatttttatgcttgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGTTTCATCtttt                                                                  >  1:259238‑M2/77‑90 (MQ=255)          gAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTggg                                                               >  1:413534/1‑90 (MQ=255)                      aacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaaa                                                         >  1:350363/1‑84 (MQ=255)                      aacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGaaaa                                                         <  2:350363/84‑1 (MQ=255)                      aacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGc                                                   >  1:325790/1‑90 (MQ=255)                      aacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGc                                                   >  2:196167/1‑90 (MQ=255)                       tatttttatgcttgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCa                                                  <  2:760983‑M2/30‑1 (MQ=255)                       acaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGACCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCa                                                  >  1:286798/1‑90 (MQ=255)                             aaaCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTg                                            >  2:260965/1‑90 (MQ=255)                                  gAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGt                                       >  1:31509/1‑90 (MQ=255)                                  gAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGt                                       >  1:658082/1‑90 (MQ=255)                                  gAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGt                                       >  2:193424/1‑90 (MQ=255)                                   tgtgaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGTTTCATCtttt                                                                  <  1:764493‑M2/14‑1 (MQ=255)                                   tgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGTTTCATCtttt                                                                  >  2:764493‑M2/49‑62 (MQ=255)                                   tgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTc                                      >  2:150167‑M2/49‑90 (MQ=255)                                     aaaGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCaa                                    <  1:533546/90‑1 (MQ=255)                                                   aTGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaa                      <  2:325790/90‑1 (MQ=255)                                                      tatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGGAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaacaa                   <  2:529378‑M2/61‑1 (MQ=255)                                                        tttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAg                 <  2:259238‑M2/63‑1 (MQ=255)                                                           ctctATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAgtgt              <  2:658082/90‑1 (MQ=255)                                                               aTGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAg          <  2:128874/90‑1 (MQ=255)                                                                 tggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGaa        >  1:419307‑M2/19‑90 (MQ=255)                                                                 tggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGaa        >  2:65879‑M2/19‑90 (MQ=255)                                                                       gtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAg                 <  1:722803‑M2/63‑1 (MQ=255)                                                                       gtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAg                 >  2:722803‑M2/13‑75 (MQ=255)                                                                       gtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTg  >  2:737140‑M2/13‑90 (MQ=255)                                                                       gtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTg  >  1:340140‑M2/13‑90 (MQ=255)                                                                       aaaTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACtgtg                                                    <  2:391049/40‑1 (MQ=255)                                                                       aaaTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACtgtg                                                    >  1:391049/1‑40 (MQ=255)                                                                                   |                                                                              ATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTG  >  NZ_CP009273/4532731‑4532890

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