New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	CP009273	4532499 =	8 (0.450)	6 (0.360)	5/258	NT	42.1%	intergenic (+49/‑433)	fimE/fimA	tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
	?	CP009273	4532813 =	9 (0.540)					intergenic (+363/‑119)	fimE/fimA	tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)

TTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCG  >  CP009273/4532674‑4532956                                                                                                                                            |                                                                                                                                                agcggaggtgatgtgaaattaatttacaatagaaataatttacatatcaaacagttagatgctttttgtcgttttttaatatttttatgcttgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGTTTCATc                                                                                                                                        >  2:49902‑M2/140‑149 (MQ=255)                         caccTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAAAAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCATGCATGAGTAAAAATGGCACCAAATGTTTCATCTTTTGGGGAAAAACTGTGCAgtgt                                                                                                                >  1:35182/1‑149 (MQ=255)                             tCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCCGTGTTGGc                                                                                                            >  1:35676/1‑149 (MQ=255)                             aaaTGCAACATAACCATCTATAAATAAAGACAACAAGAGAACATTAAGCTAACAATTAAACTGCAAAAGTTTGTCCTCGATGCTTTACTCTATGAGTCACAATGGCCCCAAATGTGTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGc                                                                                                            <  1:83477/147‑1 (MQ=255)                                tGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGATTCATCTTTTGGGGGAAGACTGTGCAGTGTTGGCAGt                                                                                                         >  2:10187/1‑149 (MQ=255)                                                         agattctttttggcgttttttaaaatttttaagcttgaaaaaaaatacgtaacttatttatgatatggacagtttggccccaaATCTTTCATCTTTTGGCGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAgtgt                                                                                <  1:133673‑M2/66‑1 (MQ=255)                                                         agatgctttttgtcgttttttaatatttttatgcttgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAgtgt                                                                                <  1:253555‑M2/66‑1 (MQ=255)                                                               aaTAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGaa                                                                          <  1:148596/149‑1 (MQ=255)                                                                        ttAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTCCAGAACGACTGccc                                                                 >  1:272464/1‑149 (MQ=255)                                                                               atatttttatgcttgagaaaaaatacgtaacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAATTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGa                                                          <  1:49902‑M2/88‑1 (MQ=255)                                                                                        aCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGACAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGACGATCTAGAAAt                                                 >  1:244408/1‑149 (MQ=255)                                                                                             aaaaGTTTGCCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGttt                                            <  2:81902/149‑1 (MQ=255)                                                                                                          taacttatttatgatatggacagtttggccccaaATGTTTCATCTTTTTGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTgaaaggaaag                               <  2:140820‑M2/115‑1 (MQ=255)                                                                                                              ttatttatgatatggacagtttggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAagcagc                           >  2:96037‑M2/31‑149 (MQ=255)                                                                                                                          gAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTaaa               <  2:35182/149‑1 (MQ=255)                                                                                                                               aagtgggccccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGAGTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAActct          <  1:59617‑M2/136‑1 (MQ=255)                                                                                                                                       cccaaATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCg  >  2:113424‑M2/6‑149 (MQ=255)                                                                                                                                            |                                                                                                                                                TTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCG  >  CP009273/4532674‑4532956

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