New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3583428 =	NA (NA)	20 (0.620)	16/486	0.7	100%	noncoding (1/768 nt)	IS1	repeat region
	?	NC_000913	= 4542049	0 (0.000)					coding (13/597 nt)	fimE	regulator for fimA

CCACTGCAATCTATTGTTATATTGAATCAAATCAATGAAAATAGATGTTGTCACATCAGTGATATTTTATTTTTGTATGATATTTAATGTAATTGACTGATAGCCACATCACTCCGTGTGTGGTTATCTTTTTATCTATTGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGA  >  NC_000913/4541824‑4542287                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                               ccACTGCAATCTATTGTTATATTGAATCAAATCAATGAAAATAGATGTTGTCACATCAGTGATATTTTATTTTTGTATGATATTTAATGTAATTGACTGATAGCCACATCACTCCGTGTGTGGTTATCTTTTTATCTATTGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGGAGTACTACTATTGTGAGTAAACGTCggtgatgctgccaacttactgattt                                                                                                                                                                                                                       <  1:150023‑M2/251‑26 (MQ=255)            tATTGTTATATTGAATCAAATCAATGAAAATAGATGTTGTCACATCAGTGATATTTTATTTTTGTATGATATTTAATGTAATTGACTGATAGCCACATCACTCCGTGTGTGGTTATCTTTTTATCTATTGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGGAGTACTACTATTGTGAGTAAACGTCggtgatgctgccaacttactgatttagtgtatgatg                                                                                                                                                                                                            <  1:317550‑M2/251‑37 (MQ=255)                                      aaaaTAGATGTTGTCACATCAGTGATATTTTATTTTTGTATGATATTTAATGTAATTGACTGATAGCCACATCACTCCGTGTGTGGTTATCTTTTTATCTATTGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGGAGTACTACTAGTGTGAGTAAACGTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggctt                                                                                                                                                                                  >  1:55997‑M2/1‑189 (MQ=255)                                              tgttgtCACATCAGTGATATTTTATTTTTGTATGATATTTAATGTAATTGACTGATAGCCACATCACTCCGTGTGTGGTTATCTTTTTATCTATTGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGGAGTACTACTATTGTGAGTAAACGTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttct                                                                                                                                                                          <  1:20966‑M2/251‑71 (MQ=255)                                                                   ttattttTGTATGATATTTAAGGTAATTGACTGATAGCCACATCACTCCGTGTGTGGTTATCTTTTTATCTATTGGGCTAAGTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGGAGTACTACTATTGTGAGTAAACGTCggtgatgctgccaacttactgatttagtgtatgatggtgttgttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttc                                                                                                                                                      >  1:54132‑M2/1‑160 (MQ=255)                                                                                    ttAATGTAATTGACTGATAGCCACATCACTCCGTGTGTGGTTATCTTTTTATCTATTGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGGAGTACTACTAGTGTGAGTAAACGTCggtgacgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctggttctatcagctgtcccgcctgttcagctactgacgaggtgg                                                                                                                                     >  2:317550‑M2/1‑143 (MQ=255)                                                                                           aaTTGACTGATAGCCACATCACTCCGTGTGTGGTTATCTTTTTATCTATTGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGGAGTACTACTATTGTGAGTAAACGTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaac                                                                                                                             >  1:160497‑M2/1‑136 (MQ=255)                                                                                                         cacatatctccgagtGTGGTTATCTTTTTCTCTCTTGGGCTAATTTTGCAAGCTTGTGGTTTCCTATAGGTATTCATTCATCTATATCTCAGTTAGGAGTACTACTATTGTGAGTAAACGTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaccgg                                                                                                                           <  2:303249‑M2/232‑118 (MQ=255)                                                                                                         cacaTCACTCCGTGTGTGGTTATCTTTTTATCTATTGGGCTAATGTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCCCAGTGAGGAGTACTACTATTGTGAGTAAACGTCggtgatgctgccaacttactgatttagtgtatgagggtgtttttgaggtgctccagtggcttctgttgctatcagctggccctcctgttcagctactgacggggtggtgcgtaacgg                                                                                                                           >  1:303249‑M2/1‑122 (MQ=255)                                                                                                                                   tttATCTATTGAACTACTTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGGTGTACTACTATTGTGAGTAAACGTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtcactcctgttcagctactgacggggtggtgcgtaacggcaaacgcaccgccggacatcagcgctatctctgctctc                                                                                     <  2:278191‑M2/251‑156 (MQ=255)                                                                                                                                           ttGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGGAGTACTACTATTGTGAGTAAACGTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgt                                                                             <  2:166783‑M2/251‑164 (MQ=255)                                                                                                                                                            gATTGAGGTTTCTTATAGGTATTCATTCAAATATATCTCAGTTAGGAGTACTACTATTGTGAGTAAACGTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgc                                                             <  2:258938‑M2/250‑180 (MQ=255)                                                                                                                                                                                                ctcAGTTAGGAGTACTACTATTGTGAGTAAACGTCggtgatgctgccaacttactgctttagtgtatgatggtgtttttgaggtgatccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcacttacaccgctactcaccccggtacgcacc                         <  2:326962‑M2/250‑216 (MQ=255)                                                                                                                                                                                                ctcAGTTAGGAGTACTACTATTGTGAGTAAACGTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcacttacaccgcttctcaacccggtacgcacc                         <  1:3743‑M2/250‑216 (MQ=255)                                                                                                                                                                                                ctcAGTTAGGAGTACTACTATTGTGAGTAAACGTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaac                                                                        <  2:299882‑M2/203‑169 (MQ=255)                                                                                                                                                                                                ctcAGTTAGGAGTACTACTATTGTGAGTAAACGTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaac                                                                        >  1:299882‑M2/1‑35 (MQ=255)                                                                                                                                                                                                   aGTTAGGAGTCCTACTATTGTGAGTTAACGTCggtgatgctgccaacttactgatttagtttatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctactgttcagctactgacggggaggtgagtcacggcaaaagcaccgccggacatcagcgctctctcttctctccctgccgtaaaacctggcaacttcagttcacttaccccgcttctcaccccggtacgcaccaga                      <  2:248190‑M2/250‑219 (MQ=255)                                                                                                                                                                                                     ttAGGAGTACTACTATTGTGAGTAAACGTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcacttacaccgcttctcaacccggtacgcaccagaaa                    <  2:172520‑M2/250‑221 (MQ=255)                                                                                                                                                                                                                    ttGTGAGTAAACGTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctacttacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcacttacaccgcttctcaacccggtacgcaccagaaaatcattgatatggcc     <  2:141495‑M2/250‑236 (MQ=255)                                                                                                                                                                                                                      gagagTAAACGTCggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcacttacaccgcttctcaacccggtacgccccagaaaatcattgctatggccatg  <  2:207959‑M2/249‑239 (MQ=255)                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                               CCACTGCAATCTATTGTTATATTGAATCAAATCAATGAAAATAGATGTTGTCACATCAGTGATATTTTATTTTTGTATGATATTTAATGTAATTGACTGATAGCCACATCACTCCGTGTGTGGTTATCTTTTTATCTATTGGGCTAATTTTGACCGATTGAGGTTTCCTATAGGTATTCATTCAAATATATCTCAGTTAGGAGTACTACTATTGTGAGTAAACGTCGTTATCTTACCGGTAAAGAAGTTCAGGCCATGATGCAGGCGGTTTGTTACGGGGCAACGGGAGCCAGAGATTATTGTCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGA  >  NC_000913/4541824‑4542287

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