New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4524799	19 (0.680)	4 (0.160)	4/246	2.9	20.2%	coding (212/906 nt)	yjhH	putative lyase/synthase
	?	NC_000913	= 4524824	14 (0.550)					coding (187/906 nt)	yjhH	putative lyase/synthase

CTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAACAGCCCGTCGACCCCTTTATTAATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGA  >  NC_000913/4524677‑4524963                                                                                                                                                    |                                                                                                                                            cTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGATACAGCTTCTTCg                                                                                                                                            >  4:43147‑M2/1‑148 (MQ=255)                         cgacgaTACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTa                                                                                                                    <  5:115244/149‑1 (MQ=255)                                                               gcCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTgg                                                                                                                          <  8:7141/105‑1 (MQ=255)                                                               gcCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTgg                                                                                                                          >  7:7141/1‑105 (MQ=255)                                                                      tGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAACAGccc                                                                        <  8:53015/148‑1 (MQ=255)                                                                                   aGTTGAAGGGGATCCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGTTTACAGCTTCTTCGGCGAGTGCCATGCGCTGTGCTGTATTCAttt                                                                                                             <  8:7777/98‑1 (MQ=255)                                                                                   aGTGGAAGGGGAACCGACGCCAATCAAAACCGGCAATAGCCCGTCGAAAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTCTATTCAttt                                                                                                             >  7:7777/1‑98 (MQ=255)                                                                                         aGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAACAg                                                                           <  2:71494/126‑1 (MQ=255)                                                                                         aGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAACAg                                                                           >  1:71494/1‑126 (MQ=255)                                                                                                      ccAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAACAGCCCGTCGACCCCTTTATTAATCAGGAAGTCGGCAAc                                       >  1:72488/1‑149 (MQ=255)                                                                                                                   gCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAACAGCCCGTCGACCCCTTTATTAATCAGGAAGTCGGCAACTTCGCGCATTGcc                          <  2:72488/149‑1 (MQ=255)                                                                                                                     aCTCGCCCGTCGACAATGGTTACAGCTTCTTCggcgagtgccggtattgattggcgtcggttccccttccactgacgaagcggtcaaactggcgcagcatgcgcaagcctacggcgctgatggtatcgtcgccatcaacccctactact                        <  3:4052‑M2/149‑118 (MQ=255)                                                                                                                           ccGTCGACAATGGTTACAGCTTCTTCggcgagtgccggtattgattggcgtcggttccccttccactgacgaagcggtcaaactggcgcagcatgcgcaagcctacggcgctgatggtatcgtcgcca                                       <  7:46823‑M2/128‑103 (MQ=255)                                                                                                                           ccGTCGACAATGGTTACAGCTTCTTCggcgagtgccggtattgattggcgtcggttccccttccactgacgaagcggtcaaactggcgcagcatgcgcaagcctacggcgctgatggtatcgtcgcca                                       >  8:46823‑M2/1‑26 (MQ=255)                                                                                                                           ccGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAACAGCCCGTCGACCCCTTTATTAATCAGGaa                                                <  2:81781/119‑1 (MQ=255)                                                                                                                             gTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAACAGCCCGTCGACCCCTTTATTAATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCaa                <  1:68951/149‑1 (MQ=255)                                                                                                                                gACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAACAGCCCGTCGACCCCTTTATTAATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTa                    <  5:119036/142‑1 (MQ=255)                                                                                                                                gACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAACAGCCCGTCGACCCCTTTATTAATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTa                    >  6:119036/1‑142 (MQ=255)                                                                                                                                           cAGCTTCTTCggcgagtgccggtattgattggcgtcggttccccttccactgacgaagcggtcaaactggcgcagcatgcgcaagcctacggcgctgatggtatcgtcgccatcaacccctactactggaaagtcgcaccacgaaatct  >  4:4052‑M2/1‑10 (MQ=255)                                                                                                                                                    |                                                                                                                                            CTTTCCAGTAGTAGGGGTTGATGGCGACGATACCATCAGCGCCGTAGGCTTGCGCATGCTGCGCCAGTTTGACCGCTTCGTCAGTGGAAGGGGAACCGACGCCAATCAATACCGGCACTCGCCCGTCGACAATGGTTACAGCTTCTTCGGCGAGTGCCATGCGCTGGGCTGTATTCATTTGGCTAAATTCACCACCGGTACCCAGATAAAACAGCCCGTCGACCCCTTTATTAATCAGGAAGTCGGCAACTTCGCGCATTGCCTTTTTATCAAGGGTTCCGTCACGA  >  NC_000913/4524677‑4524963

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