New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4542690	22 (0.620)	5 (0.150)	5/246	2.6	13.6%	intergenic (+57/‑425)	fimE/fimA	regulator for fimA/type 1 fimbriae major subunit
	?	NC_000913	= 4542986	43 (1.310)					intergenic (+353/‑129)	fimE/fimA	regulator for fimA/type 1 fimbriae major subunit
Rejected: Frequency below/above cutoff threshold.

ATGCTGCTCGTTTTGCCGGATTATG‑‑‑GGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTAT  >  NC_000913/4542560‑4542822                                                                                                                                      |                                                                                                                                     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                <  2:603417‑M1/127‑119 (MQ=255)                                                                                                                              ttGGGGCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatttatagatgtttatattgcatgaggtggtttttggagag                >  1:603417‑M1/1‑9 (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)                                                                                                                               TGGGGCCAAACTGTCTATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCt         >  1:591977‑M1/1‑133 (MQ=255)                                                                                                                               TGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCt         >  1:531074‑M1/1‑133 (MQ=255)                                                                                                                                 GGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGACATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTcc       <  1:307468‑M1/133‑1 (MQ=255)                                                                                                                                   gCCAttttgactcatagaggaaagcatcgcggacaaactttttcagtttatttgttggcttaatattctattgttatctttatttatagatgtttatattgcatgaggtggtttttggagagaagaatgagga     >  2:141730‑M1/1‑4 (MQ=255)                                                                                                                                   GCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGc     <  1:349561‑M1/133‑1 (MQ=255)                                                                                                                                    CCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCt    <  1:504615‑M1/133‑1 (MQ=255)                                                                                                                                      AAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTAGTGTAAATAAATGTCGCATCACCTCCGCtat  <  1:169036‑M1/133‑1 (MQ=255)                                                                                                                                      |                                                                                                                                     ATGCTGCTCGTTTTGCCGGATTATG‑‑‑GGAAAGAAATAATCTCATAAACGAAAAATTAAAAAGAGAAGAGGTTTGATTTAACTTATTGATAATAAAGTTAAAAAAACAAATAAATACAAGACAATTGGGGCCAAACTGTCCATATCATAAATAAGTTACGTATTTTTTCTCAAGCATAAAAATATTAAAAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTAT  >  NC_000913/4542560‑4542822

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