New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1530096 =	36 (0.470)	3 (0.040)	3/248	NT	8.2%	coding (175/483 nt)	ydcD	putative immunity protein for RhsE
	?	NC_000913	1530130 =	33 (0.450)					coding (209/483 nt)	ydcD	putative immunity protein for RhsE

CATATTCTCCAGTTTACTAATGTTTTTTAATCGCACGCTATAAAGAATATCATTGAAAAAATTAAGATTTATATCAACCACTCCACAGTGCGTCCTGTAATTATTAACCATTATACTTATAAAATGACCATTCCCTCTACTAATTTTTTGGTCATT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1530251‑1530096 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tggtcattTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTTATAGCGTGCGATTAAAAAACATTAGTAAACTGGAGAATATGGAGTTCTGCGCTACAAAACAACGG  >  NC_000913/1530130‑1530275                                                                                                                                                                                                                                                                                                                CATATTCTCCAGTTTACTAATGTTTTTTAATCGCACGCTATAAAGAATATCATTGAAAAAATTAAGATTTATATCAACCACTCCACAGTGCGTCCTGTAATTATTAACCATTATACTTATAAAATGACCATTCCCTCTACTAATTTTTT                                                                                                                                                           >  3:568167/1‑149    ATTCTCCAGTTTACTAATGTTTTTTAATCGCACGCTATAAAGAATATCATTGAAAAAATTAAGATTTATATCAACCACTCCACAGTGCGTCCTGTAATTATTAACCATTATACTTATAAAATGACCATTCCCTCTACTAATTTTTTGGT                                                                                                                                                        <  4:280277/149‑1       CTCCAGTTTACTAATGTTTTTTAATCGCACGCTATAAAGAATATCATTGAAAAAATTAAGATTTATATCAACCACTCCACAGTGCGTCCTGTAATTATTAACCATTATACTTATAAAATGACCATTCCCTCTACTAATTTTTTGGTCAT                                                                                                                                                     <  2:492094/149‑1        TCCAGTTTACTAATGTTTTTTAATCGCACGCTATAAAGAATATCATTGAAAAAATTAAGATTTATATGAACCACTCCACAGTGCGTCCTGTAATTATTAACCATTATACTTATAAAATGACCATTCCCTCTACTAATTTTTTGGTCATT                                                                                                                                                    >  4:48401/1‑149        TCCAGTTTACTAATGTTTTTTAATCGCACGCTATAAAGAATATCATTGAAAAAATTAAGATTTATATCAACCACTCCACAGTGCGTCCTGTAATTATTAACCATTATACTTATAAAATGACCATTCCCTCTACTAATTTTTTGGTCATT                                                                                                                                                    >  3:500453/1‑149                                                                                                            CCATTATACTTATAAAATGACCATTCCCTCTACTAATTTTTTGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAA                                                                              >  3:329136/1‑119                                                                                                            CCATTATACTTATAAAATGACCATTCCCTCTACTAATTTTTTGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAA                                                                              <  4:329136/119‑1                                                                                                                           AATGACCATTCCCTCTACTAATTTTTTGGTCATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTTATAGCGTGCGATTAAAAAACATTAGTAAACTGGA                                 >  1:338892/1‑149                                                                                                                                                          ATTTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTTATAGCGTGCGATTAAAAAACATTAGTAAACTGGAGAATATGGAGTTCTGCGCTACAAAACAACGG  <  2:523715/149‑1                                                                                                                                                                                                                                                                                                                CATATTCTCCAGTTTACTAATGTTTTTTAATCGCACGCTATAAAGAATATCATTGAAAAAATTAAGATTTATATCAACCACTCCACAGTGCGTCCTGTAATTATTAACCATTATACTTATAAAATGACCATTCCCTCTACTAATTTTTTGGTCATT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1530251‑1530096 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tggtcattTTATAAGTATAATGGTTAATAATTACAGGACGCACTGTGGAGTGGTTGATATAAATCTTAATTTTTTCAATGATATTCTTTATAGCGTGCGATTAAAAAACATTAGTAAACTGGAGAATATGGAGTTCTGCGCTACAAAACAACGG  >  NC_000913/1530130‑1530275

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