New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 788831	0 (0.000)	23 (0.640)	19/228	0.7	100%	coding (7/1041 nt) coding (1149/1149 nt)	galM galK	aldose 1‑epimerase, type‑1 mutarotase, galactose mutarotase galactokinase
	?	NC_000913	789981 =	0 (0.000)					intergenic (‑2/+2)	galK/galT	galactokinase/galactose‑1‑phosphate uridylyltransferase

TCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGCCTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGTTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTA  >  NC_000913/789870‑790104                                                                                                                |                                                                                                                            ggaaagtaaagtcgcaccccagtccatcagcgtgactaccatccctgcgttgttacgcaaagttaacagtcggtacggctgaccatcgggtgccagtgcgggagtttcgttTCTTACACTCCGGATTCGCGAAAATGGa                                                                                                  <  1:455046‑M2/28‑1 (MQ=255)  gaaagtaaagtcgcaccccagtccatcagcgtgactaccatccctgcgttgttacgcaaagttaacagtcggtacggctgaccatcgggtgccagtgcgggagtttcgttTCTTACACTCCGGATTCGCGAAAATGGat                                                                                                 <  1:311435‑M2/29‑1 (MQ=255)  gaaagtaaagtcgcaccccagtccatcagcgtgactaccatccctgcgttgttacgcaaagttaacagtcggtacggctgaccatcgggtgccagtgcgggagtttcgttTCTTACACTCCGGATTCGCGAAAATGGat                                                                                                 >  2:311435‑M2/111‑139 (MQ=255)                            cagcgtgactaccatccctgcgttgttacgcaaagttaacagtcggtacggctgaccatcgggtgccagtgcgggagtttcgttTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTctg                                                                       >  1:177521‑M2/85‑139 (MQ=255)                            cagcgtgactaccatccctgcgttgttacgcaaagttaacagtcggtacggctgaccatcgggtgccagtgcgggagtttcgttTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTctg                                                                       >  2:645566‑M2/85‑139 (MQ=255)                                 tgactaccatccctgcgttgttacgcaaagttaacagtcggtacggctgaccatcgggtgccagtgcgggagtttcgttTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGcc                                                                  >  2:684622‑M2/80‑139 (MQ=255)                                            cctgcgttgttacgcaaagttaacagtcggtacggctgaccatcgggtgccagtgcgggagtttcgttTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGt                                                       <  2:197294‑M2/71‑1 (MQ=255)                                                 gttgttacgcaaagttaacagtcggtacggctgaccatcgggtgccagtgcgggagtttcgttTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGt                                                  <  2:497943‑M2/76‑1 (MQ=255)                                                      taactttgcgtaacagtcggtacggctgaccatcgggtgccagtgcgggagtttcgttTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGc                                             >  1:9583‑M2/59‑139 (MQ=255)                                                       acgcaaagttaacagtcggtacggctgaccatcgggtgccagtgcgggagtttcgttTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCt                                            <  1:652054‑M2/82‑1 (MQ=255)                                                             agttaacagtcggtacggctgaccatcgggtgccagtgcgggagtttcgttTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGtct                                      <  1:218037‑M2/88‑1 (MQ=255)                                                             agttaacagtcggtacgactgaccatcgggtgccagtgcgggagtttcgttTCTTACACTCCGGATTCGCGACAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGtct                                      >  2:440536‑M2/52‑139 (MQ=255)                                                               ttaacagtcggtacggctgaccatcgggtgccagtgcgggagtttcgttTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGt                                                       >  1:413313‑M2/50‑120 (MQ=255)                                                               ttaacagtcggtacggctgaccatcgggtgccagtgcgggagtttcgttTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGt                                                       <  2:413313‑M2/71‑1 (MQ=255)                                                                    agtcggtacggctgaccatcgggtgccagtgcgggagtttcgttTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGTCTCTGCCAg                               <  2:702536‑M2/95‑1 (MQ=255)                                                                         gtacggctgaccatcgggtgccagtgcgggagtttcgttTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGTCTCTGCCAGCAttt                          >  1:21230‑M2/40‑139 (MQ=255)                                                                           acggctgaccatcgggtgccagtgcgggagtttcgttTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGTCTCTGCCAg                               >  2:26814‑M2/38‑132 (MQ=255)                                                                           acggctgaccatcgggtgccagtgcgggagtttcgttTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGTCTCTGCCAg                               <  1:26814‑M2/95‑1 (MQ=255)                                                                                  accatcgggtgccagtgcgggagtttcgttTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGTCTCTGCCAGCATTTCATaaccaa                 >  2:76287‑M2/31‑139 (MQ=255)                                                                                   ccatcgggtgccagtgcgggagtttcgttTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGTCTCTGCCAGCATTTCATaaccaac                >  2:572858‑M2/30‑139 (MQ=255)                                                                                                 tgcgggagtttcgttTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTtacgta  <  1:684622‑M2/124‑1 (MQ=255)                                                                                                   cgggagtttcgttTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGTCTCTGc                                  >  2:35555‑M2/14‑105 (MQ=255)                                                                                                   cgggagtttcgttTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGTCTCTGc                                  <  1:35555‑M2/92‑1 (MQ=255)                                                                                                                |                                                                                                                            TCGGTGTGTTCACCAATCAAATTCACGCGGCCAGGCGCCTGAATGGTGTGAGTGGCAGGGTAGCCAAATGCGTTGGCAAACAGAGATTGTGTTTTTTCTTTCAGACTCATTTCTTACACTCCGGATTCGCGAAAATGGATATCGCTGACTGCGCGCAAACGCTCTGCTGCCTGTTCTGCGGTCAGGTCTCGCTGGGTCTCTGCCAGCATTTCATAACCAACCATAAATTTACGTA  >  NC_000913/789870‑790104

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