New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1634457	NA (NA)	3 (0.050)	3/482	3.2	7.9%	coding (429/576 nt)	tfaQ	putative tail fiber assembly protein TfaQ
	?	NC_000913	4536248 =	35 (0.610)					intergenic (‑2195/‑349)	topAI/ythB	toxin of the TopAI‑YjhQ toxin‑antitoxin system, TopA inhibitor/protein YthB
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

TCCAGCAACGAGATTTCTTCCTCCGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATC  >  NC_000913/1634392‑1634705                                                                  |                                                                                                                                                                                                                                                         tCCAGCAACGAGATTTCTTCCTCCGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATGTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCg                                                                   <  2:599319/249‑1 (MQ=11)                                  accAGATCTGCAGCATCCTGAAGCGGCGCAATATGCGCACTGGCTACCTGCATCAGGTTGTGTTTTGTGTCTTCCGCCTCCCGGATCCGGAACAGTTTATCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCGTCCGGCGATAACCAGGTAACATTTTCCTGTAACGGACCGAGTTCAGAAATAAATAACGCATCGCCGGAATCCACTTCATAAACCGTTTTACACCgat                                 <  2:243487/249‑1 (MQ=11)                                    cAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCgatgat                              >  1:433300/1‑251 (MQ=31)                                       aTCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCtt                           >  1:108584/1‑251 (MQ=31)                                        tCTGCAACATACTGAAGCGGCGAACTATGCTCACTGGCTACCTGCATAAGGTTGTTTTTTGTTTCTTACGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAACTAAATACAGCGTCGCCGGAAGCCACGTCATACACCGGTTTACCCCGATGCTCTTc                          >  2:251732/1‑251 (MQ=11)                                            cagcaTCCTGAAGCGGCGCATTATGCTCCCTTCTTACCTGCGTCAGTTTGTTTTTTGTTGCTTCCGCCTCCCGGATCCGGAGCAGTTTTTCTGTTTCCGTATGCTGCACCCATGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCATGGAGCATTTTCCGGTAACGGGCCGAGTTGAGAAATAAGTAACGCTTCGCCGATAGGCACGGCAGGAACCGTTTGAcccccaggagctacaacg                      <  2:300101/251‑15 (MQ=2)                                                     gAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGcc               <  1:405333/249‑1 (MQ=31)                                                               aaTATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCATCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGACGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACtcat   >  1:39424/1‑251 (MQ=14)                                                                atatGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCAACGATGACtcatc  >  1:20971/1‑251 (MQ=14)                                                                  |                                                                                                                                                                                                                                                         TCCAGCAACGAGATTTCTTCCTCCGTTGCAATTTCCAGATCTGCAGCATCCTGAAGCGGCGCAATATGCTCACTGGCTACCTGCATCAGGTTGTTTTTTGTTTCTTCCGCCTCCCGGATCCGGAACAGTTTTTCTGCTTCCGTATCCTTCACCCAGGCTGTGCCGTTCCACTTCTGAAACTCCCCTTCCGGCGATAACCAGGTAACATTTTCCGGTAACGGACCGAGTTCAGAAATAAATAACGCGTCGCCGGAAGCCACGTCATAAACCGTTTTACCCCGATGATCTTCAACGAGATGCCACGATGACTCATC  >  NC_000913/1634392‑1634705

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