New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4542690	4 (0.060)	19 (0.270)	19/474	1.6	86.7%	intergenic (+57/‑425)	fimE/fimA	regulator for fimA/type 1 fimbriae major subunit
	?	NC_000913	= 4542986	2 (0.030)					intergenic (+353/‑129)	fimE/fimA	regulator for fimA/type 1 fimbriae major subunit

AAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGCTCTGTCCCTCAGTTCTACAGCGGCTCTGGCCGCTGCCACGACGGTTAATGGTGGGACCGTTC  >  NC_000913/4542746‑4543214                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                     aaaaCGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaat                                                                                                                                                                                                                            >  1:587846‑M2/1‑241 (MQ=255)                 tCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTAAAAAAACAACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGATTGTCCGAGATGCTTTCCTCTATGAGTCAAAAtagacccaattgtcttgtatttattt                                                                                                                                                                                                            <  1:458775‑M2/251‑27 (MQ=255)                    aaCTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatt                                                                                                                                                                                                             >  1:4464‑M2/1‑222 (MQ=255)                                 gTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaacttta                                                                                                                                                                                             >  1:96081‑M2/1‑209 (MQ=255)                                   aaaTTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGGCAAAAtggccccaattgtcttgtatttatttgtttttttaactttatta                                                                                                                                                                                          >  1:110960‑M2/1‑207 (MQ=255)                                                                                         aCGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaataagttaaatcaaacctcttcgctttttaattttgcgtttatgagattat                                                                                                                                    >  2:616287‑M2/1‑153 (MQ=255)                                                                                                   ggCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCGCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaataagttaaatcaaacctcttctctttttaatttttcgtttatgagattatttctttccc                                                                                                                           >  2:221829‑M2/1‑143 (MQ=255)                                                                                                                     cATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaataagttaaatcaaacctcttctctttttaatttttcgtttatgagattatttctttcccataatccggcaaaacgagc                                                                                                        >  1:316173‑M2/1‑125 (MQ=255)                                                                                                                      aTTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaataagttaaatcaaacctcttctctttttaatttttcgtttatgagattatttctttcccataatccggcaaaacgagca                                                                                                       >  1:203871‑M2/1‑124 (MQ=255)                                                                                                                                ccAAAAACCACCTCATGCAATATACACATCTATACATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaataagttaaatcaaacctcttctctttttaatttttcgtttatgagattatttctttcccataatccggcaaaacgagcagcattactgg                                                                                             <  2:138430‑M2/251‑138 (MQ=255)                                                                                                                                             cATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaataagttaaatcaaacctcttctctttttaatttttcgtttatgagattatttctttcccataatccggcaaaacgagcagcattactggcggtataacgcac                                                                                <  1:27715‑M2/251‑151 (MQ=255)                                                                                                                                             cATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCa                                                                                 >  1:21968/1‑250 (MQ=255)                                                                                                                                                                               aTAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaataagttaaatcaaacctcttctctttttaatttttcgtttatgagattatttctttcccataatccggcaaaacgagcagcattactggcggtataacgcacagtatggcgaatatttcgatgcccgagataatcc                                              <  2:533216‑M2/251‑185 (MQ=255)                                                                                                                                                                                 aGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaataagttaaatcaaacctcttctctttttaatttttcgtttatgagattatttctttcccataatccggcaaaacgagcagcattactggcggtataaagcacagtattgcgaatatttcgatgacggagataatcctg                                            <  1:388969‑M2/251‑187 (MQ=255)                                                                                                                                                                                     tatTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGCTGCTTTCCTCTATGAGTCAAAAtggccacaattgtcttgtatttatttgtttttttaactttattatcaataagttaaatcaaaactcttctctttttaatttttcgtttatgagattatttctttcccataatccggcaaaacgagcagcattactggcggaattacgaccaatatgacgaatatttcggtgcccgagataatcctgaatt                                        <  2:194828‑M2/251‑191 (MQ=255)                                                                                                                                                                                        tAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaataagttaaatcaaacctcttctctttttaatttttcgtttatgagattatttctttcccataatccggcaaaacgagcagcattactggcggtataacgcacagtatggcgaatatttcgatgcccgagataatcctgaattaa                                      >  1:533216‑M2/1‑58 (MQ=255)                                                                                                                                                                                                caaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaataagttaaatcaaacctcttctctttttaatttttcgtttatgagattatttctttcccataatccggcaaaacgagcagcattactggcggtataacgcacagtatggcgaatatttcgatgcccgcgataataatgaattaaacgagattt                             <  2:142238‑M2/251‑202 (MQ=255)                                                                                                                                                                                                             aaaGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaataagttaaatcaaacctcttctctttttaatttttcgtttatgagattatttctttcccataatccggcaaaacgagcagcattactggcggtataacgcacagtatcgcgcatatttcgatgcacgggatcatcctgcattaaacgagtatctgcaccacgctcc                <  2:359202‑M2/251‑215 (MQ=255)                                                                                                                                                                                                                       cgcgATGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttacctttattatcaataagttaaatcaaacctcttctctttttaatttttcgtttatgagattatttctttcccataatccggcaaaacgagcagcattactggcggtataacgcacagtatggcgaatatttcgatgcccgagataatcctgaattaaacgagtatctgcaccacgctccgccaattcat      <  1:207829‑M2/251‑225 (MQ=255)                                                                                                                                                                                                                       cgcgATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGTGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTGTTTGAAAGGAAAGCAGCATGAAAATGAAAACTCTGGCAATCGTTGTTCTGTCGGCTCTGTCCCTCAGTTCTACAGCGGCTCTGGCCGCTGCCACGACGGTTAATGGTGGGAcc      <  1:417867/251‑1 (MQ=255)                                                                                                                                                                                                                           aTGCTTTCCTCTATGAGTCAAAAtggccccaattgtcttgtatttatttgtttttttaactttattatcaataagttaaatcaaacctcttctctttttaatttttcgtttatgagattatttcgttcccataatccggcaaaacgagcagcattactggcggtataacgcacagtatggcgaatatttcgatgcccgagataatcctgaattaaacgagtatctgcaccacgctccgccaattcataacc  >  1:466056‑M2/1‑23 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                     AAAACGACAAAAAGCATCTAACTGTTTGATATGTAAATTATTTCTATTGTAAATTAATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCAGTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGCTCTGTCCCTCAGTTCTACAGCGGCTCTGGCCGCTGCCACGACGGTTAATGGTGGGACCGTTC  >  NC_000913/4542746‑4543214

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