New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1432951	NA (NA)	7 (0.170)	7/276	1.8	NA	coding (541/576 nt)	tfaR	Rac prophage; putative tail fiber assembly protein TfaR
	?	NC_000913	= 1433168	NA (NA)					coding (507/591 nt)	pinR	Rac prophage; putative site‑specific recombinase

ATTTTTTATATGTCAGGTGTCTGAAGTTTTGCTTTGGCTCTTAAAATGGTTTGCCGCGAGGTTTTGAATTCCCGGGCAATGGCACTTATACTTACACCTGACTTAATTCGTTCGAATACCACCTGTTTCTGTTCTTCATTTAACACAGGAGGTCGACCAAAACGTTTCCCTGCGCCGCGGGCTCTTACTATCCCGGAAT  >  NC_000913/1433079‑1433277                                                                                          |                                                                                                              attTTTTATATGTCAGGTGTCTGAAGTTTTGCTTTGGCTCTTAAAATGGTTTGCCGCGAGGTTTTGAATTCCCGGGCAATGGCACTTATACTTACACCTGACTTAATTCGTTCGAATACCACCTGTTTCTGTTCTTCAtt                                                             <  1:596862/140‑1 (MQ=32) attTTTTATATGTCAGGTGTCTGAAGTTTTGCTTTGGCTCTTAAAATGGTTTGCCGCGAGGTTTTGAATTCCCGGGCAATGGCACTTATACTTACACCTGACTTAATTCGTTCGAATACCACCTGTTTCTGTTCTTCAtt                                                             <  2:620744/140‑1 (MQ=32)                           ttttGCTTTGGCTCTTAAAATGGTTTGCCGCGAGGTTTTGAATTCCCGGGCAATGGCACTTATACTTACACCTGACTTAATTCGTTCGAATACCACCTGTTTCTGTTCTTCATTTAACACAGGAGGTCGACCAAAACGtt                                   >  1:663250/1‑140 (MQ=34)                                     gCTCTTAAAATGGTTTGCCGCGAGGTTTTGAATTCCCGGGCAATGGCACTTATACTTACACCTGACTTAATTCGTTCGAATACCACCTGTTTCTGTTCTTCATTTAACACAGGAGGTCGACCAAAACGTTTCCCTGcgcc                         >  2:200585/1‑140 (MQ=34)                                     gCTCTTAAAATGGTTTGCCGCGAGGTTTTGAATTCCCGGGCAATGGCACTTATACTTACACCTGACTTAATTCGTTCGAATACCACCTGTTTCTGTTCTTCATTTAACACAGGAGGTCGACCAAAACGTTTCCCTGcgcc                         >  2:397327/1‑140 (MQ=34)                                                ggTTTGCCGCGAGGTTTTGAATTCCCGGGCAATGGCACTTATACTTACACCTGACTTAATTCGTTCGAATACCACCTGTTTCTGTTCTTCATTTAACACAGGAGGTCGACCAAAACGTTTCCCTGCGCCGCGGGCTCTTa              <  2:663250/140‑1 (MQ=34)                                                         cgAGGTTTTGAATTCCCGGGCAATGGCACTTATACTTACACCTGACTTAATTCGTTCGAATACCACCTGTTTCTGTTCTTCATTTAACACAGGAGGTCGACCAAAACGTTTCCCTGCGCCGCGGGCTCTTACTATCCCgg     >  2:253433/1‑140 (MQ=34)                                                            ggTTTTGAATTCCCGGGCAATGGCACTTATACTTACACCTGACTTAATTCGTTCGAATACCACCTGTTTCTGTTCTTCATTTAACACAGGAGGTCGACCAAAACGTTTCCCTGCGCCGCGGGCTCTTACTATCCCGGAAt  >  1:558759/1‑140 (MQ=34)                                                            ggTTTTGAATTCCCGGGCAATGGCACTTATACTTACACCTGACTTAATTCGTTCGAATACCACCTGTTTCTGTTCTTCATTTAACACAGGAGGTCGACCAAAACGTTTCCCTGCGCCGCGGGCTCTTACTATCCCGGAAt  <  2:223625/140‑1 (MQ=34)                                                                                          |                                                                                                              ATTTTTTATATGTCAGGTGTCTGAAGTTTTGCTTTGGCTCTTAAAATGGTTTGCCGCGAGGTTTTGAATTCCCGGGCAATGGCACTTATACTTACACCTGACTTAATTCGTTCGAATACCACCTGTTTCTGTTCTTCATTTAACACAGGAGGTCGACCAAAACGTTTCCCTGCGCCGCGGGCTCTTACTATCCCGGAAT  >  NC_000913/1433079‑1433277

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