New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3873741 =	111 (1.870)	3 (0.050)	3/262	NT	2.9%	coding (734/879 nt)	dgoK	2‑oxo‑3‑deoxygalactonate kinase
	?	NC_000913	3873793 =	96 (1.710)					coding (682/879 nt)	dgoK	2‑oxo‑3‑deoxygalactonate kinase
Rejected: Frequency below/above cutoff threshold.

GGTCTTAATACGCCCGCCATATTGCCGCAGCTTTTTGAAGTTCGCGCCTCGCATGTGCTGGGAACACTTCCCCGCGAACAGGTCAGCGAATTTCTCTCTGGTTTGTTGATTGGCGCAGAGGTCGCCAGTATGCGCGACTATGTGGCCCATCAACA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3873895‑3873741 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑atcaacaAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGTATTAAGACCACGCTCAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGG  >  NC_000913/3873793‑3873940                                                                                                                                                                                                                                                                                                                 tGTCTTAATACGCCCGCCATATTGCCGCAGCTCTATGAAGTGCGCGCCTCGCATGTGCTGGGAACACTTCCCCGCGAACAGGTCAGCGAATTTCTCTCTGGATTTTTGATTGGCGCAGAGGTCGCCAGTATGCGCGACTATGTGGCCCA                                                                                                                                                            <  1:702100/148‑1  GTCTAAATACGCCCGCCATATTGCCGCAGCTTTTTGAAGTTCGCGCCTCGCATGTGCTGGGAACACTTCACCGCGAACAGGTCAGCGAAATTCTCTCTGGTTTGTTGATTGGCGCAGAGGTCGCCCGTATGCGCGACTATGTGGCCCAT                                                                                                                                                           >  1:426732/1‑149    CTTAATACGCCCGCCATATTGCCGCAGCTTTTTGAAGTTCGCGCCTCGCATGTGCTGGGAACACTTCCCCGCGAACAGGTCAGCGAATTTCTCTCTGGTTTGTTGATTGGCGCAGAGGTCGCCAGTATGCGCGACTATGTGGCCCATCA                                                                                                                                                         >  1:873704/1‑149    CTTAATACGCCCGCCATATTGCCGCAGCTTTTTGAAGTTCGCGCCTCGCATGTGCTGGGAACACTTCCCCGCGAACAGGTCAGCGAATTTCTCTCAGGTTTGTTGATTGGCGCAGAGGTCGCCAGTATCCCCGACTATCTGGCCCATCA                                                                                                                                                         >  1:119822/1‑149     TTAATACGCCCGCCATATTGCCGCAGCTTTTTGAAGTTCGCGCCTCGCATGTGCTGGGAACACTTCCCCGCGACCAGGTCAGCGAATTTCTCTCTGGTTTGTTGATTGGCGCAGAGGTCGCCAGTATGCGCGACAATGTGGCCCATCAA                                                                                                                                                        >  2:294808/1‑149       AATACGCCCGCCATATTGCCGCAGCTTTTTGAAGGTCGCGACGCGCATGTGCTGGGAACACTTCCCCGCGAACAGGTCAGCGAATTTCTCTCTGGTTTGTTGATTGGCGCAGAGGTCGCCAGTATGCGCGACTATGTGGCCCATCAACA                                                                                                                                                      <  2:574388/149‑1                                                                                                  CTGGTTTGTTGATTGGCGCAGAGGTCGCCAGTATGCGCGACTATGTGGCCCATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGG                                                           >  2:226927/1‑149                                                                                                     GTTTGTTGATTGGCGCAGAGGTCGCCAGTATGCGCGACTATGTGGCCCATCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGT                                                        >  1:1025134/1‑149                                                                                                      TTTGTTCATTGGCGCAGAGGTCGCCAGTATGCGCGACTATGTGGCCCATCACCAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCAAGCACAAGGGAGGCGCGAACTGCAAAAAGCTGCGGCAATATGGCGGGCGTA                                                       >  2:748064/1‑149                                                                                                                                                      TCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGTATTAAGACCCCGCTCAAGGCCCGCTGTGAAGGCATCGGCAGAGTTTTCC       >  2:856521/1‑149                                                                                                                                                      TCAACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGTATTAAGACCACGCTCAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCC       <  1:3648/149‑1                                                                                                                                                      TCAACAAACCAGAGAGAAATTAGCTGACCTGTTCGCGTGGAAGTGTTACCAGCACATGCGAGGCGAGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGTATTAAGACCACGCTCAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCC       <  2:659537/149‑1                                                                                                                                                       CAACAAACCAGAGAGAAATTCGCTGACCTGTTCGAGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAACGCTGCGGCAATATGGCGGGCGTATTAAGACCACGCTCAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCT      <  1:741030/149‑1                                                                                                                                                        AACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGTATTAAGACCACGCTCAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTG     <  1:214338/149‑1                                                                                                                                                        AACAAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGGATTAAGACCACGCTCAAGGACAGATGTGACGGCATCGGCAGAGTTTgcggg     >  1:561321/1‑144                                                                                                                                                           AAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGTATTAAGACCACGCTCAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGG  <  2:722998/149‑1                                                                                                                                                                                                                                                                                                                 GGTCTTAATACGCCCGCCATATTGCCGCAGCTTTTTGAAGTTCGCGCCTCGCATGTGCTGGGAACACTTCCCCGCGAACAGGTCAGCGAATTTCTCTCTGGTTTGTTGATTGGCGCAGAGGTCGCCAGTATGCGCGACTATGTGGCCCATCAACA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3873895‑3873741 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑atcaacaAACCAGAGAGAAATTCGCTGACCTGTTCGCGGGGAAGTGTTCCCAGCACATGCGAGGCGCGAACTTCAAAAAGCTGCGGCAATATGGCGGGCGTATTAAGACCACGCTCAAGGCCAGCTGTGAAGGCATCGGCAGAGTTTTCCTGCGG  >  NC_000913/3873793‑3873940

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