New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1979270	NA (NA)	18 (0.500)	17/486	1.0	100%	noncoding (1/768 nt)	IS1	repeat region
	?	NC_000913	= 4541023	0 (0.000)					coding (67/603 nt)	fimB	Type 1 fimbriae regulatory protein FimB

GTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACGTTAGTTATTAAGCAAACGTTCGCTTTTATAAACATAATCAGGATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTA  >  NC_000913/4540786‑4541245                                                                                                                                                                                                                                              |                                                                                                                                                                                                                               gTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACGTTAGTTATTAAGCAAACGTTCGCTTTTATAAACATAATCAGGATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAggtgatgctgcca                                                                                                                                                                                                                   >  2:15738‑M2/1‑238 (MQ=255)             aaCAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACGTTAGTTATTAAGCAAACGTTCGCTTTTATAAACATAATCAGGATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAggtgatgctgccaacttactgattt                                                                                                                                                                                                       <  2:19045‑M2/251‑26 (MQ=255)                                                                            aaaCGTTCGCTTTTATAAACATAATCAGGATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgt                                                                                                                                        <  1:252690‑M2/251‑89 (MQ=255)                                                                                                   aTCAGGATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcg                                                                                                                 <  1:232974‑M2/251‑112 (MQ=255)                                                                                                             aaaTGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgtgcagctactgacggggtggtgcgtaacggcaaa                                                                                                       >  2:38742‑M2/1‑130 (MQ=255)                                                                                                                            gCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacat                                                                                        <  1:204470‑M2/251‑137 (MQ=255)                                                                                                                                     gCACAGCTAGTGCGCGTCTGTAATTATAAGGGCAAAACGATGAAGAATAAGGCTGATAACAAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctat                                                                               <  2:30470‑M2/251‑146 (MQ=255)                                                                                                                                                       tGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatca                                                                                      >  1:317894‑M2/1‑88 (MQ=255)                                                                                                                                                       tGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatca                                                                                      <  2:317894‑M2/226‑139 (MQ=255)                                                                                                                                                        gTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgt                                                             >  2:183993‑M2/1‑87 (MQ=255)                                                                                                                                                         tAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCTCTCCTTAggtgatgctgccaacttcctgatttagtgtatgatggtgtttttgaggtgctccggtggcttctgtttctatcagctggccctcctgttcagctactgactgggtggtacttaacggcaaaagcgccgccggacatcagcgctatctctgctctcactgccgta                                                            >  2:61221‑M2/1‑86 (MQ=255)                                                                                                                                                          aaTTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaa                                                           <  1:338193‑M2/250‑166 (MQ=255)                                                                                                                                                                   ggAAAAACGATGAAGAATAAGGCTGATAACAAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggca                                                  <  2:301815‑M2/250‑175 (MQ=255)                                                                                                                                                                                agaaTAAGGCTGATAACAAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcactt                                    <  2:294544‑M2/251‑189 (MQ=255)                                                                                                                                                                                     aaGGCTGATAACAAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcacttacac                                >  2:38298‑M2/1‑58 (MQ=255)                                                                                                                                                                                                               ccTGACCCATAGTGAAATCGAATCACTCCTTAggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcacttacaccgcttctcaacccggtacgcaccaga      <  2:7070‑M2/250‑219 (MQ=255)                                                                                                                                                                                                                  gACCCATAGTGAAATCGAATCACTCCTTAggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcacttacaccgcttctcaacccggtacgcaccagaaaat  >  1:57211‑M2/1‑29 (MQ=255)                                                                                                                                                                                                                  gACCCATAGTGAAATCGAATCACTCCTTAggtgatgctgccaacctactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcacttacaccgcttctcaacccggtacgcaccagaaaat  >  1:148252‑M2/1‑29 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                               GTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACGTTAGTTATTAAGCAAACGTTCGCTTTTATAAACATAATCAGGATAAAAATGTTGGATTATTGCTAACCCAGCACAGCTAGTGCGCGTCTGTAATTATAAGGGAAAAACGATGAAGAATAAGGCTGATAACAAAAAAAGGAACTTCCTGACCCATAGTGAAATCGAATCACTCCTTAAAGCAGCAAATACCGGGCCTCATGCAGCACGTAATTATTGTCTGACTTTGCTTTGTTTTATTCATGGTTTCCGGGCGAGTGAAATTTGTCGATTGAGGATTTCGGATATTGATCTTAAGGCAAAGTGTATATATATCCATCGATTAAAAAAAGGCTTTTCAACAACGCACCCGCTATTGAATAAAGAAGTTCAGGCTTTAAAAAACTGGTTGAGTATCCGTA  >  NC_000913/4540786‑4541245

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