New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	5731921 =	98 (1.110)	9 (0.110)	7/234	NT	8.9%	coding (691/1206 nt)	hutI	imidazolonepropionase
	?	NC_002947	5731965 =	91 (1.090)					coding (647/1206 nt)	hutI	imidazolonepropionase

GTCGGGCCGACGACTACATCGCGCACATCTGTGACGAGATGTTGCCGGCCTTGGCCGGCGAAGGCCTGGTGGATGCGGTGGATGCCTTCTGCGAACACCTCGCGTTCTCGCCGGCCCAGGTCGAGCGGCTGTTCATCAAGGCGCGCGAACTGGGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/5732076‑5731921 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ctgggccGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGC  >  NC_002947/5731965‑5732112                                                                                                                                                                                                                                                                                                                  GTCGGGCCGACGACTACATCGCGCACATCTGTGACGAGATGTTGCCGGCCTTGGCCGGCGAAGGCCTGGTGGATGCGGTGGATGCCTTCTGCGAACACCTCGCGTTCTCGCCGGCCCAGGTCGAGCGGCTGTTCATCAAGGCGCGCGAAC                                                                                                                                                            <  1:468084/150‑1     GGCCGACGACTACATCGCGCACATCTGTGACGAGATGTTGCCGGCCTTGGCCGGCGAAGGCCTGGTGGATGCGGTGGATGCCTTCTGCGAACACCTCGCGTTCTCGCCGGCCCAGGTCGAGCGGCTGTTCATCAAGTCGCGCGAACTGGG                                                                                                                                                        >  1:2295155/1‑150     GGCCGACGACTACATCGCGCACATCTGTGACGAGATGTTGCCGGCCTTGGCCGGCGAAGGCCTGGTGGATGCGGTGGATGCCTTCTGCGAACACCTCGCGTTCTCGCCGGCCCAGGTCGAGCGGCTGTTCATCAAGGCGCGCGAACTGGG                                                                                                                                                        >  2:1325506/1‑150      GCCGACGACTACATCGCGCACATCTGTGACGAGATGTTGCCGGCCTTGGCCGGCGAAGGCCTGGTGGATGCGGTGGATGCCTTCTGCGAACGCCTCGCGTTCTCGCCGGCCCAGGTCGAGCGGCTGTTCATCAAGGCGCGCGAACTGGGC                                                                                                                                                       <  2:1486040/150‑1            GACTACATCGCGCACATCTGTGACGAGATGTTGCCGGCCTTGGCCGGCGAAGGCCTGGTGGATGCGGTGGATGCCTTCTGCGAACACCTCGCGTTCTCGCCGGCCCAGGTCGAGCGGCTGTTCATCAAGGCGCGCGAACTGG                                                                                                                                                         <  1:833481/142‑1            GACTACATCGCGCACATCTGTGACGAGATGTTGCCGGCCTTGGCCGGCGAAGGCCTGGTGGATGCGGTGGATGCCTTCTGCGAACACCTCGCGTTCTCGCCGGCCCAGGTCGAGCGGCTGTTCATCAAGGCGCGCGAACTGG                                                                                                                                                         >  2:833481/1‑142                                                                                                                    CCAGGTCGAGCGGCTGTTCATCAAGGCGCGCGAACTGGGCCTGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCC                                         >  1:1468075/1‑150                                                                                                                    CCAGGTCGAGCGGCTGTTCATCAAGGCGCGCGAACTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCC                                         >  2:1283716/1‑150                                                                                                                    CCAGGTCGAGCGGCTGTTCATCAAGGCGCGCGAACTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCG                                                                                   <  1:668464/108‑1                                                                                                                    CCAGGTCGAGCGGCTGTTCATCAAGGCGCGCGAACTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCG                                                                                   >  2:668464/1‑108                                                                                                                       GGTCGAGCGGCTGTTCATCAAGGCGCGCGAACTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACA                                                                   <  1:1679763/121‑1                                                                                                                       GGTCGAGCGGCTGTTCATCAAGGCGCGCGAACTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACA                                                                   >  2:1679763/1‑121                                                                                                                             GCGGCTGTTCATCAAGGCGCGCGAACTGGGCCTGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGTCAACATCTCGTCACAGATGTGCGCGATGTA                                                    >  2:262851/1‑130                                                                                                                             GCGGCTGTTCATCAAGGCGCGCGAACTGGGCCTGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTA                                                    <  1:262851/130‑1                                                                                                                                    TTCATCAAGGCGCGCGAACTGGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGG                         <  1:1283716/150‑1                                                                                                                                                        GGGCCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCA     >  1:533976/1‑150                                                                                                                                                           CCGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGC  <  2:351975/150‑1                                                                                                                                                            CGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGG                      <  1:1958409/129‑1                                                                                                                                                            CGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGG                      >  2:1958409/1‑129                                                                                                                                                                                                                                                                                                                  GTCGGGCCGACGACTACATCGCGCACATCTGTGACGAGATGTTGCCGGCCTTGGCCGGCGAAGGCCTGGTGGATGCGGTGGATGCCTTCTGCGAACACCTCGCGTTCTCGCCGGCCCAGGTCGAGCGGCTGTTCATCAAGGCGCGCGAACTGGGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/5732076‑5731921 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ctgggccGGCGAGAACGCGAGGTGTTCGCAGAAGGCATCCACCGCATCCACCAGGCCTTCGCCGGCCAAGGCCGGCAACATCTCGTCACAGATGTGCGCGATGTAGTCGTCGGCCCGACCCGCGTACTCCGGGGGCAAAGCGTGCGCCGCCAGGC  >  NC_002947/5731965‑5732112

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