New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	CP009273	= 4532507	4 (0.150)	4 (0.170)	4/250	NT	26.0%	intergenic (+57/‑425)	fimE/fimA	tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)
	?	CP009273	= 4532803	19 (0.790)					intergenic (+353/‑129)	fimE/fimA	tyrosine recombinase/inversion of on/off regulator of fimA/major type 1 subunit fimbrin (pilin)

CTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGCTCTGTCCCTCAGTTCTACAGCGGCTCTGGCCGCTGCCACGACGGTTAATGGTGGGACCGTTCACTTTAAAGGGGAAGTTGTTAACGCCGCTTGCGCAGTTGATGCAG  >  CP009273/4532677‑4533076                                                                                                                               |                                                                                                                                                                                                                                                                                  cTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTCGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaataagttaaatcaaacctcttctctttttaatttttcgtttatgagattatttctttcccataatccggcaaaacgagcagcattactggcggtataacg                                                                                                                                >  1:277998‑M2/1‑127 (MQ=255)                                             ataaataaAGATAACAATAGAATAGTAAGACAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTNGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGCTCTGTCCCTCAGTTCTACAGCGGctc                                                                                   >  1:367051/1‑275 (MQ=255)                                                     aGATAACAATAGAATATTAAGTCGCCAAAAAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCAGCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACACGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGNTTTTTGAAAGGCAAGCAGCATGAAAATTAAAACTCTGGCAATCGATGTTCTGTCGGCTCTGGTCAACAGTTCTACAGAGGCACTggacaag                                                                           >  1:365018/1‑270 (MQ=255)                                                     aGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGCTCTGTCCCTCAGTTCTACAGCGGCTCTGGCCGCt                                                                           >  1:124307/1‑275 (MQ=255)                                                           aaaTAGCATATTATGCCAACAAAAAACCTGAAAAAGTTTGTCCGCGATTCTTTCCTCTATGAGTCAAAATGGACCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGCTCTGTCCCTCAGTTCTACAGCGGCTCTGGCCGCTGCCacg                                                                     <  2:367045/274‑1 (MQ=255)                                                                   tatagaGCCAAAAAATAAACTGAAAAAGGTTGTCCGCGATGCTTTCCTATATGAGTCAAAATGGCCCCAAATTTTTCATCTTGTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAAGTAAAACTCTGGCAATCGTTGTTCTGTCGGCTCTGTCCCTCAGTTCTACAGCGGCTCTGGCCGCTGCCACGACGGTTaa                                                             <  2:124306/270‑1 (MQ=255)                                                                      tAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCgttgtt                                                                                                                     >  2:504136/1‑216 (MQ=255)                                                                        aGCCAACCAATACACTGTAAAAGTTGGTCTGCGTTACTAAATTTTACGAGACACCATGGGCCCACCTGTTTCATCTTCACGGCGAACACTGTGCAGTGTTGGCCGTCAACCTCGTTTACAAATCAACGTGTACAGAACGACTACCCATGTCGATTTAGACATAGTTTTTTGACACGAAAGCAGCATGCAAATTCAAACTATCGCAATCGTTGTTCTGTCGGCTCTGTCCCTCAGTTATACACCGGCTCTGGCCGCTGCCACGCNCGTTaacgttc                                                        >  2:49114/1‑270 (MQ=255)                                                                        aGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGaaaa                                                                                                                                            >  1:87738/1‑191 (MQ=255)                                                                          ccAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAATTGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGt                                                                                                                 <  1:504145/216‑1 (MQ=255)                                                                            aacaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTaa                                                                                                                                        <  2:87738/191‑1 (MQ=255)                                                                             acaaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaataagttaaatcaaacctcttctctttttaatttttcgtttatgagattatttctttcccataatccggcaaaacgagcagcattactggcggtataacgcacagtatggcgaatatttcgatgcccgagataatcctgaat                                                                                      <  2:4048‑M2/240‑190 (MQ=255)                                                                              caaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaa                                                                                                                                                                                                                >  1:87306/1‑117 (MQ=255)                                                                                 aTAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaataagttaaatcaaacctcttctctttttaatttttcgtttatgagattatttctttcccataatccggcaaaacgagcagcattactggcggtataacgcacagtatggcgaatatttcgatgcccgagataatcctgaattaaa                                                                                  >  1:4048‑M2/1‑47 (MQ=255)                                                                                  tAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaacaaa                                                                                                                                                                                                            <  2:87306/117‑1 (MQ=255)                                                                                         aaaaaGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtg                                                                                                                                                                                                                                                                                 <  1:152465‑M2/41‑3 (MQ=255)                                                                                         aaaaaGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAActc                                                                                                                                   >  2:36792/1‑183 (MQ=255)                                                                                             aGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccc                                                                                                                                                                                                                                                                             >  2:152464‑M2/1‑35 (MQ=255)                                                                                             aGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGc                                                                                                                               <  1:36792/183‑1 (MQ=255)                                                                                                 tGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaa                                                                                                                                                                                                                >  1:335499/1‑98 (MQ=255)                                                                                                     cgcgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaacaaa                                                                                                                                                                                                            <  2:335494/98‑1 (MQ=255)                                                                                                                    tATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGTGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTaca                                                                                                                                                                                                                   >  2:259123/1‑76 (MQ=255)                                                                                                                    tATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGTCTCTGTCCCTCAGTTCTACAGCGGCTCTGGCCGCTGCCACGACGGTTAATGGTGGGACAGTTCACTTTAAAGGGGAAGTTGTTAACGCGGCTTgcg              >  1:231152/1‑273 (MQ=255)                                                                                                                        aGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTacaaaac                                                                                                                                                                                                               <  1:259127/76‑1 (MQ=255)                                                                                                                         gTCAAAAtggccccaattgtcttgtatttatttgtttttttaactt                                                                                                                                                                                                                                            <  1:100779‑M2/46‑40 (MQ=255)                                                                                                                             aaatggccccaattgtcttgtatttatttgtttttttaactttatt                                                                                                                                                                                                                                        >  2:100779‑M2/1‑3 (MQ=255)                                                                                                                              aaTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGCTCTGTCCCTCAGTTCTCCAGCGGCTCTGGCCGCTGCCACGACGGTTAATTGTGGGACCGGTCACTTTAAAGGGGAAGTAGTTAACGCCGACTGAGCAGTTTATGcga  >  1:262752/1‑273 (MQ=255)                                                                                                                               |                                                                                                                                                                                                                                                                                  CTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGCTCTGTCCCTCAGTTCTACAGCGGCTCTGGCCGCTGCCACGACGGTTAATGGTGGGACCGTTCACTTTAAAGGGGAAGTTGTTAACGCCGCTTGCGCAGTTGATGCAG  >  CP009273/4532677‑4533076

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