New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 315290	12 (0.740)	4 (0.330) +TTAGCCGCGTGTCCACTATTGCTGGGTAAGATCA	3/196	NT	26.5%	noncoding (62/1255 nt)	IS3	repeat region
	?	NC_000913	= 572471	18 (1.110)					coding (6/456 nt)	ybcN	DLP12 prophage; DNA base‑flipping protein

GTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCA  >  NC_000913/315152‑315422                                                                                                                                           |                                                                                                                                     gTCGGGCTGTTTGGTAAAGATAAGCGCCAGGGAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTCACCTCTGACTGAGAGTTTACTCACttag                                                                                                                                  <  1:285777‑M1/143‑5 (MQ=255)          tttGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAg                                                                                                                         >  2:92058/1‑143 (MQ=255)             ggTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACACAACCAGTAt                                                                                                                      >  1:22324/1‑143 (MQ=255)                aaaGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGTCAAAAACAGTATCaa                                                                                                                   >  2:141392/1‑143 (MQ=255)                              gAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAAccc                                                                                                     <  1:96633/143‑1 (MQ=255)                                          tGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATAAACCAGTAAAAAACCCCGTAAACAGCAt                                                                                         >  2:175086/1‑143 (MQ=255)                                              aTTACTGCCGAGGAGAGCTATACGCCAGTGCCACTTACCCAGCAATAGTGGACACGCGGCTAAGTGAATAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGc                                                                                     <  2:56561/143‑1 (MQ=14)                                                               cTATACGCCAGTGCCTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGaa                                                                    >  1:252683/1‑143 (MQ=255)                                                                   aCGCCAGTGCCTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGccc                                                                <  2:39559/143‑1 (MQ=38)                                                                       cAGTGCCTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGTCAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGaa                                                            <  1:141392/143‑1 (MQ=25)                                                                       cAGTGCCTCTTACCCAGCAATACTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGaa                                                            <  1:245469/143‑1 (MQ=25)                                                                         gTGCCTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGc                                                          >  1:160268/1‑143 (MQ=33)                                                                             ctctTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCGACCAGTAAAAAACCTCGTAAACCGCCTCACCCTGAATTTCGCAGTGAAGACCTGAAGCTTGc                                                      >  1:161855/1‑143 (MQ=11)                                                                                                       ggCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACttagccgcgtgtccactattgctgggtaagatcagttcacgcttcccctccgcagaggtcaaacgctagatgcaacgaattgcaggtgcatttctgcatctgtgaa                            >  1:170123‑M1/1‑37 (MQ=255)                                                                                                          tAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACttagccgcgtgtccactattgctgggtaagatcagttcacgcttcccctccgcagaggtcaaacgctagatgcaaagaattgcaggtgcatttctgcatctgtgaaggg                         >  1:73322‑M1/1‑34 (MQ=255)                                                                                                                      cgctcagTCAGAGGTGACTCACttagccgcgtgtccactattgctgggtaagatcagttcacgcttcccctccgcagaggtcaaacgctagatgcaaagaattgcaggtgcatttctgcatctgtgaagggagaagagagttt             >  1:305931‑M1/3‑22 (MQ=255)                                                                                                                                 agttgacGCACtctgccgcctgcccgctattgctgggtaagatcagttcacgcttcccctccccagaggtcaaacgctagatgcaaagaattgcaggtgcatttctgcatctgtgaagggagaagagagtttggattgtatgt  <  1:22122‑M1/140‑133 (MQ=255)                                                                                                                                           |                                                                                                                                     GTCGGGCTGTTTGGTAAAGATAAGCGCCAGAAAGACCCGCATGCTATTTCTGCCGAGGTGACCTATACGCCAGTGCCTCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCA  >  NC_000913/315152‑315422

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