New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	4237457 =	57 (1.050)	3 (0.060)	3/250	NT	6.5%	coding (793/1122 nt)	catB	muconate cycloisomerase 1
	?	NC_002947	4237489 =	32 (0.620)					coding (761/1122 nt)	catB	muconate cycloisomerase 1

CATGGTGCCGCCGTACAAGGCGATGCCAGCGGCCTCGGCGATCTGTGCAGTGCGCAGAACCGCGCGCGGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTC  >  NC_002947/4237344‑4237637                                                                                                                                                  |                                                                                                                                                     gacaacggcatcgaccttatcgagcagccgatttcgcgcatcaaccgcgctggccaggtgcgcctgaaccagctcagtccggctccgatcatggccgatgagtcgatcgaaagCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCttt                                                                                                                                                   >  8:6547‑M1/114‑149 (MQ=255) cATGGTGCCGCCGTACAAGGCGATGCCAGCGGCCTCGGCGATCTGTGCAGTGCGCAGAACCGCGCGCGGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCttt                                                                                                                                                   <  7:302721/149‑1 (MQ=255) cATGGTGCCGCCGTACAAGGCGATGCCAGCGGCCTCGGCGATCTGTGCAGTGCGCAGAACCGCGCGCGGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCttt                                                                                                                                                   <  5:58956/149‑1 (MQ=255) cATGGTGCCGCCGTACAAGGCGATGCCAGCGGCCTCGGCGATCTGTGCAGTGCGCAGAACCGCGCGCGGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCttt                                                                                                                                                   >  2:311380/1‑149 (MQ=255) cATGGTGCCGCCGTACAAGGCGATGCCAGCGGCCTCGGCGATCTGTGCAGTGCGCAGAACCGCGCGCGGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCtt                                                                                                                                                    >  2:307365/1‑148 (MQ=255)    ggTGCCGCCGTACAAGGCGATGCCAGCGGCCTCGGCGATCTGTGCAGTGCGCAGAACCGCGCGCGGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTtcga                                                                                                                                                >  7:32180/1‑149 (MQ=255)     gTGCCGCCGTACAAGGCGATGCCAGCGGCCTCGGCGATCTGTGCAGTGCGCAGAACCGCGCGCGGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTtcgat                                                                                                                                               <  6:181065/149‑1 (MQ=255)                 aaGGCGATGCCAGCGGCCTCGGCGATCTGTGCAGTGCGCAGAACCGCGCGCGGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGc                                                                                                                                   >  5:328564/1‑149 (MQ=255)                  aGGCGATGCCAGCGGCCTCGGCGATCTGTGCAGTGCGCAGAACCGCGCGCGGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGc                                                                                                                                   >  3:166713/1‑148 (MQ=255)                   ggCGATGCCAGCGGCCTCGGCGATCTGTGCAGTGCGCAGAACCGCGCGCGGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCATCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCa                                                                                                                                 >  8:191768/1‑149 (MQ=255)                   atcgagcagccgatttcgcgcatcaaccgcgctggccaggtgcgcctgaaccagcgcagtccggctccgatcatggccgatgagtcgatcgaaagCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTtcgatc                                                                                                                                              >  4:90271‑M1/96‑136 (MQ=255)                   atcgagcagccgatttcgcgcatcaaccgcgctggccaggtgcgcctgaaccagcgcagtccggctccgatcatggccgatgagtcgatcgaaagCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTtcgatc                                                                                                                                              <  3:90271‑M1/41‑1 (MQ=255)                       aTGCCAGCGGCCTCGGCGATCTGTGCAGTGCGCAGAACCGCGCGCGGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGAt                                                                                                                             >  6:72933/1‑149 (MQ=255)                       aTGCCAGCGGCCTCGGCGATCTGTGCAGTGCGCAGAACCGCGCGCGGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGAt                                                                                                                             >  7:25399/1‑149 (MQ=255)                        tGCCAGCGGCCTCGGCGATCTGTGCAGTGCGCAGAACCGCGCGCGGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATc                                                                                                                            <  7:202509/149‑1 (MQ=255)                                       cGATCTGTGCAGTGCGCAGAACCGCGCGCGGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCAt                                                                                                                                <  4:5984/130‑1 (MQ=255)                                       cGATCTGTGCAGTGCGCAGAACCGCGCGCGGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCAt                                                                                                                                >  3:5984/1‑130 (MQ=255)                                              tgCAGTGCGCAGAACCGCGCGCGGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCt                                                                                                            <  8:32180/143‑1 (MQ=255)                                              tgCAGTGCGCAGAACCGCGCGCGGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTg                                                                                                           <  5:72933/144‑1 (MQ=255)                                              tgCAGTGCGCAGAACCGCGCGCGGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCa                                                                                                      <  3:300908/149‑1 (MQ=255)                                                cAGTGCGCAGAACCGCGCGCGGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAgg                                                                                                    >  2:94926/1‑149 (MQ=255)                                                 aGTGCGCAGAACCGCGCGCGGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGgc                                                                                                   <  6:328564/149‑1 (MQ=255)                                                  gTGCGCAGAACCGCGCGCGGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGgcg                                                                                                  >  6:178229/1‑149 (MQ=255)                                                   tGCGCAGAACCGCGCGCGGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAAGCGTCGAGGACGCCTtca                                                                                                                                                 >  6:240876/1‑100 (MQ=255)                                                        aGAACCGCGCGCGGGCCACCATTCTTGGCGATTTTCAGGGCGAAGAGGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTg                                                                                            <  5:268075/149‑1 (MQ=255)                                                                   cgGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGt                                                                                 <  5:178229/149‑1 (MQ=255)                                                                      gCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCAGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGa                                                                              >  3:232682/1‑149 (MQ=255)                                                                       ccaccaTTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGTCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCTGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGAt                                                                             <  2:232652/149‑1 (MQ=255)                                                                       ccaccaTTCTTGGCGATTTTCAGGGCGAAGATGCGGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGAt                                                                             <  4:232682/149‑1 (MQ=255)                                                                        caccaTTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATg                                                                            >  2:62346/1‑149 (MQ=255)                                                                         accaTTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATgc                                                                           <  3:226028/149‑1 (MQ=255)                                                                                ttGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAAt                                                                    >  7:324439/1‑149 (MQ=255)                                                                                  ggCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATc                                                                   <  2:316046/148‑1 (MQ=255)                                                                                     gATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATgc                                                                           >  2:108585/1‑137 (MQ=255)                                                                                     gATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATgc                                                                           <  1:108585/137‑1 (MQ=255)                                                                                     gATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGgct                                                               >  1:137130/1‑149 (MQ=255)                                                                                          tCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATAGGCTGCTc                                                           >  4:189588/1‑148 (MQ=255)                                                                                                         gatcgaaagcgtcgaggacgccttcagcctggccgccgacgCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGAt                                                        >  4:90271‑M2/42‑136 (MQ=255)                                                                                                         gatcgaaagcgtcgaggacgccttcagcctggccgccgacgCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGAt                                                        <  3:90271‑M2/95‑1 (MQ=255)                                                                                                              aaagcgtcgaggacgccttcagcctggccgccgacgCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGAGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATAAGGTCGATGCCGTTGTc                                      >  8:6547‑M2/37‑149 (MQ=255)                                                                                                                            cAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTgg                         >  4:157120/1‑148 (MQ=255)                                                                                                                                    aGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGcgcg                 >  4:218158/1‑148 (MQ=255)                                                                                                                                     ggCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGa               >  3:331521/1‑149 (MQ=255)                                                                                                                                      gccaagacgacgCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGa               <  3:180020‑M2/136‑1 (MQ=255)                                                                                                                                       cggccTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATg             <  7:259632/146‑1 (MQ=255)                                                                                                                                           cTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCt         >  4:177154/1‑149 (MQ=255)                                                                                                                                            tCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGcc                                            >  6:123626/1‑113 (MQ=255)                                                                                                                                            tCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGcc                                            <  5:123626/113‑1 (MQ=255)                                                                                                                                              gacgCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCTGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTgg       <  1:62346‑M2/144‑1 (MQ=255)                                                                                                                                                  cTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTc  <  7:191768/149‑1 (MQ=255)                                                                                                                                                  cTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTc  <  3:157120/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     CATGGTGCCGCCGTACAAGGCGATGCCAGCGGCCTCGGCGATCTGTGCAGTGCGCAGAACCGCGCGCGGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTC  >  NC_002947/4237344‑4237637

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