New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	CP009273	= 4532507	7 (0.390)	9 (0.540)	7/258	NT	53.7%	intergenic (+57/‑425)	fimE/fimA	tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
	?	CP009273	= 4532803	9 (0.540)					intergenic (+353/‑129)	fimE/fimA	tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)

TCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAA  >  CP009273/4532656‑4532943                                                                                                                                                    |                                                                                                                                             tCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAt                                                                                                                                             >  1:248241‑M2/1‑148 (MQ=255)                        gctcttcctctccaaaaacCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttat                                                                                                                      <  2:160084‑M2/142‑25 (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)                                                             cATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaataagttaaatcaa                                                                                 <  2:188781‑M2/149‑62 (MQ=255)                                                                     aataaaGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtgtttatttgtttttttaactttattatcaataagttaaatcaaacctcttc                                                                         <  2:103305‑M2/149‑70 (MQ=255)                                                                       taaaGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAACtggccccaattgtcttgtatttatttgtttttttaactttattatcaata                                                                                            >  1:284890‑M2/1‑78 (MQ=255)                                                                       taaaGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaata                                                                                            <  2:284890‑M2/128‑51 (MQ=255)                                                                                 aaTAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGaa                                                             <  1:148596/149‑1 (MQ=255)                                                                                          ttAAGCCAACAAATAAACTGAAAGAGTTTGTGCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaataagttaaatcaaacctcttctctttttaatttttcgtttat                                                    >  2:289582‑M2/1‑59 (MQ=255)                                                                                          ttAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTCCAGAACGACTGccc                                                    >  1:272464/1‑149 (MQ=255)                                                                                                          aCTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGACAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGACGATCTAGAAAt                                    >  1:244408/1‑149 (MQ=255)                                                                                                               aaaaGTTTGCCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGttt                               <  2:81902/149‑1 (MQ=255)                                                                                                                          cgcgATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttcttatcaataagttaaatcaaacctcttctctttttaatttttcgtttatgagattatttctttcccataatccggccaaac                    <  2:121281‑M2/149‑123 (MQ=255)                                                                                                                          cgcgATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaataagttaaatcaaacctcttctctttttaatttttcgtttatgagattatttctttcccataatccggcaaaac                    <  2:121282‑M2/149‑123 (MQ=255)                                                                                                                                 cTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaataagttaaatcaaacctcttctctttttaatttttcgtttatgagattatttctttcccataatccggcaaaacgagcagc             >  2:31833‑M2/1‑20 (MQ=255)                                                                                                                                            gAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTaaa  <  2:35182/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                             TCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAA  >  CP009273/4532656‑4532943

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