New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	4542682 =	22 (0.620)	4 (0.120)	4/246	2.9	11.2%	intergenic (+49/‑433)	fimE/fimA	regulator for fimA/type 1 fimbriae major subunit
	?	NC_000913	4542996 =	43 (1.310)					intergenic (+363/‑119)	fimE/fimA	regulator for fimA/type 1 fimbriae major subunit
Rejected: Frequency below/above cutoff threshold.

ACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATG‑‑‑GGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACC  >  NC_000913/4542550‑4542814                                                                                                                                        |                                                                                                                                     aCCGCCAGTAATGCTGCTCGTTTTGCCGGATTATG‑‑‑GGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT                                                                                                                                      <  2:364576‑M1/133‑1 (MQ=255) aCCGCCAGTAATGCTGCTCGTTTTGCCGGATTATG‑‑‑GGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT                                                                                                                                      <  1:471125‑M1/133‑1 (MQ=255)   cGCCAGTAATGCTGCTCGTTTTGCCGGATTATG‑‑‑GGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTAATGATAATAAAGTTAGAAAAACAAATCAATAGCAGACAATcg                                                                                                                                    >  2:281831‑M1/1‑131 (MQ=255)   cGCCAGTAATGCTGCTCGTTTTGCCGGATTATG‑‑‑GGAAAGAAATAATCGCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTG                                                                                                                                    >  2:399670‑M1/1‑133 (MQ=255)     ccAGTAATGCTGCTCGTTTTGCCGGATTATG‑‑‑GGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGG                                                                                                                                  >  1:585952‑M1/1‑133 (MQ=255)     ccAGTAATGCTGCTCGTTTTGCCGGATTATG‑‑‑GGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT                                                                                                                                      >  2:389669‑M1/1‑129 (MQ=255)     ccAGTAATGCTGCTCGTTTTGCCGGATTATG‑‑‑GGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAT                                                                                                                                      <  1:389669‑M1/129‑1 (MQ=255)       taattttcatgctgctttcctttcaaaaaactatttctaaatcgacatgggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTgg                                                                                                                                   <  1:420456‑M1/4‑1 (MQ=255)        gTAATGCTGCGCGTTCTGCCGGATCATG‑‑‑GGAACGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCC                                                                                                                               <  2:247605‑M1/133‑1 (MQ=255)        gTAACGCTGCTCGTTTTGCCGGATTATG‑‑‑GGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAATAAAACAAGTAATTAGAAGACATTTGGGGCC                                                                                                                               >  1:164552‑M1/1‑133 (MQ=255)           aTGCTGCTCGTTTTGCCGGATTATG‑‑‑GGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCaaa                                                                                                                            <  2:585952/133‑1 (MQ=255)              ctgctCGTTTTGCCGGATTATGGCTGGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATAGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAtt                                                                                                                            <  2:350549‑M1/133‑3 (MQ=255)                gctCGTTTTGCCGGACTAAG‑‑‑GGGAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTc                                                                                                                       <  1:281831/133‑1 (MQ=255)                 ctCGTTTTGCCGGATTATG‑‑‑GGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTc                                                                                                                       <  1:556225/132‑1 (MQ=255)                 ctCGTTTTGCCGGATTATG‑‑‑GGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTc                                                                                                                       >  2:556225/1‑132 (MQ=255)                   cGTTTTGCTGGATTACG‑‑‑GCAAACAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactca                                                                                                                    <  1:286469‑M1/133‑11 (MQ=255)                      tttGCCGGATTATG‑‑‑GGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATc                                                                                                                 <  2:23205/133‑1 (MQ=255)                       ttcctttcaaaaaactatttctaaatcgacatgggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCata                                                                                                                   >  2:456987‑M1/114‑133 (MQ=255)                            ggATTATG‑‑‑GGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCataaat                                                                                                           <  2:118772/133‑1 (MQ=255)                                  tG‑‑‑GGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTa                                                                                                     >  2:521016/1‑133 (MQ=255)                                                  gacatgggcagtcgttctgtacactttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTtctc                                                                                        >  1:508101‑M1/87‑133 (MQ=255)                                                  atCTCATAAATGAAATATTAAAAACAGACGAGTTTTGATTTCACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTtctc                                                                                        <  2:164552/133‑1 (MQ=255)                                                       aTAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCa                                                                                   >  1:603089/1‑133 (MQ=255)                                                       aTAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCa                                                                                   >  1:5093/1‑133 (MQ=255)                                                                         ctttgttttgtaaacgagtttgactgccaacactgcacagttttcccccaaaagatgaaacatTTGGGGCCAAACTGTCTATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACg                                                                 <  2:591977‑M1/70‑1 (MQ=255)                                                                           agaagagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAAATGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGAc                                                               >  1:441344/1‑133 (MQ=255)                                                                             aagagGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattct                                                             >  2:16918‑M1/1‑68 (MQ=255)                                                                                agGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACaaaaa                                                          <  2:603089/133‑1 (MQ=255)                                                                                     tGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATc                                                     >  1:394047/1‑133 (MQ=255)                                                                                      gATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCt                                                    <  2:5093/133‑1 (MQ=255)                                                                                       aTTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatct                                                   <  1:16918‑M1/133‑76 (MQ=255)                                                                                        tttAACTTACTGATATTAAAGTTAAACAGACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTaa                                                  <  2:278189/133‑1 (MQ=255)                                                                                            aCTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGTCTCTAACttt                                              >  2:304899/1‑131 (MQ=255)                                                                                             cTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGtt                                             >  2:536996/1‑133 (MQ=255)                                                                                                aTTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatttata                                          <  1:597675‑M1/133‑85 (MQ=255)                                                                                                 ttGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGat                                         >  1:127898/1‑133 (MQ=255)                                                                                                              tAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTAttt                            <  1:521016/133‑1 (MQ=255)                                                                                                                   ttcccccaaaagatgaaacatTTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTAtt                       <  2:531074‑M1/112‑1 (MQ=255)                                                                                                                              taCAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTGTATTGTAaattaatt            >  2:247576/1‑133 (MQ=255)                                                                                                                                   gACAATTGGGGCCAACCTGTTCATATCAGAAAAAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTccca       <  2:253470/133‑1 (MQ=255)                                                                                                                                   gACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTcaca       >  1:62277/1‑133 (MQ=255)                                                                                                                                        ttGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatttatagatgtttatattgcatgaggtggtttttggagag        >  1:603417‑M1/1‑9 (MQ=255)                                                                                                                                        ttGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatttatagatgtttatattgcatgaggtggtttttggagag        <  2:603417‑M1/127‑119 (MQ=255)                                                                                                                                        TTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCAcc  <  1:419489‑M1/133‑1 (MQ=255)                                                                                                                                        TTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTAGTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATCATTTCTATTGTAAACTAATTTGACATCAcc  <  1:253470‑M1/133‑1 (MQ=255)                                                                                                                                        TTGGGGCCAAACTCTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCAcc  <  2:382185‑M1/133‑1 (MQ=255)                                                                                                                                        |                                                                                                                                     ACCGCCAGTAATGCTGCTCGTTTTGCCGGATTATG‑‑‑GGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACC  >  NC_000913/4542550‑4542814

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