New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4074181	62 (0.960)	3 (0.050)	3/252	NT	4.9%	coding (443/897 nt)	yihV	6‑deoxy‑6‑sulphofructose kinase
	?	NC_000913	= 4074200	56 (0.890)					coding (462/897 nt)	yihV	6‑deoxy‑6‑sulphofructose kinase
Rejected: Frequency below/above cutoff threshold.

TGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4074031‑4074181 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcacACCCGCCTGACGGGCCAGGGTGAAGGCTTTTTTAGCGCCGTCGTGCCAGCGTACATCTGCCAGCACAACATCCCACTGAGAGAAATCAATTTCCTCCAACCACTCTGCGTCAGGCAGCAGATCCGGGCTGGGGTAGTTAATGATTATC  <  NC_000913/4074200‑4074053                                                                                                                                                                                                                                                                                                             TGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTC                                                                                                                                                        <  3:177040/149‑1                          AATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTT‑‑CACCCGCCTGACGGGCCAGGGTGAAG                                                                                                                             <  4:375794/149‑1                                                                                            TCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACACCCGCCTGACGGGCCAGGGT                                                                                                                                 >  1:345030/1‑81                                                                                            TCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACACCCGCCTGACGGGCCAGGGT                                                                                                                                 <  2:345030/81‑1                                                                                                                                                     CACACCCGCCTGACGGGCCAGGGTGAAGGCTTTTTTAGCGCCGTCGTGCCAGCGTACATCTGCCAGCACAACATCCCACTGAGAGAAATCAATTTCCTCCAACCACTCTGCGTCAGGCAGCAGATCCGGGCTGGGGTAGTTAATGATTA    >  2:261561/1‑149                                                                                                                                                     CACACCCGCCTGACGGGCCAGGGTGAAGGCTTTTTTAGCGCCGTCGTGCCAGCGTACATCTGCCAGCACAACATCCCACTGAGAGAAATCAATTTCCTCCAACCACTCTGCGTCAGGCAGCAGATCCGGGCTGGGG                 >  1:439503/1‑136                                                                                                                                                     CACACCCGCCTGACGGGCCAGGGTGAAGGCTTTTTTAGCGCCGTCGTGCCAGCGTACATCTGCCAGCACAACATCCCACTGAGAGAAATCAATTTCCTCCAACCACTCTGCGTCAGGCAGCAGATCCGGGCTGGGG                 <  2:439503/136‑1                                                                                                                                                     CACACCCGCCTGACGGGCCAGGGTGAAGGCTTTTTTAGCGCCGTCGTGCCAGCGTACATCTGCCAGCACAACATCCCACTGAGAGAAAACAATTTCCTCCAACCACTCTGCGTCAGGCAGCAGATCCGGGCTGGGGTAGTTAATGATTA    <  1:19724/149‑1                                                                                                                                                      ACACCCGCCTGACGGGCCAGGGTGAAGGCTTTTTTAGCGCCGTCGTGCCAGCGTACATCTGCCAGCACAACATCCCACTGAGAGAAATCAATTTCCTCCAACCACTCTGCGTCAGGCAGCAGATCCGGGCTGGGGTAGTTAATGATTAT   >  3:394257/1‑149                                                                                                                                                       CACCCGCCTGACGGGCCAGGGTGAAGGCCTTTTTAGCGCCGTCGTGCCAGCGTACATCTGCCAGCACAACATCCCACTGAGAGAAAGCAATTTCCTCCAACCACTCTGCGTCAGGCAGCAGATCCGGGCTGGGGTAGTTAATGATTATC  >  4:359565/1‑149                                                                                                                                                                                                                                                                                                             TGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4074031‑4074181 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcacACCCGCCTGACGGGCCAGGGTGAAGGCTTTTTTAGCGCCGTCGTGCCAGCGTACATCTGCCAGCACAACATCCCACTGAGAGAAATCAATTTCCTCCAACCACTCTGCGTCAGGCAGCAGATCCGGGCTGGGGTAGTTAATGATTATC  <  NC_000913/4074200‑4074053

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