New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	4238005 =	47 (1.110)	4 (0.100)	4/244	NT	8.3%	coding (245/1122 nt)	catB	muconate cycloisomerase 1
	?	NC_002947	4238048 =	44 (1.100)					coding (202/1122 nt)	catB	muconate cycloisomerase 1

CAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGATGGTGGTGGCTTCACCGATGCCTTCCACGCCATCGCTGCAGCGCAGTCGCAATACCACCAGG  >  NC_002947/4237866‑4238151                                                                                                                                            |                                                                                                                                                    cagcagTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATc                                                                                                                                             >  2:12/1‑148 (MQ=255)           cTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGc                                                                                                                                  >  5:258023/1‑149 (MQ=255)              aCCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAg                                                                                                                               >  7:49991/1‑149 (MQ=255)              aCCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAg                                                                                                                               <  5:253325/149‑1 (MQ=255)                cGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGt                                                                                                                             <  5:129754/149‑1 (MQ=255)                  gCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑c                                                                                                                          <  2:138010/149‑1 (MQ=255)                    aGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑cgc                                                                                                                        <  8:112419/149‑1 (MQ=255)                     ggCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑cgcg                                                                                                                       <  7:96912/149‑1 (MQ=255)                      gCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGt                                                                                                                      <  6:52493/149‑1 (MQ=255)                          aGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGAt                                                                                                                  >  5:218730/1‑149 (MQ=255)                           ggCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATg                                                                                                                 >  2:86618/1‑149 (MQ=255)                                   ccctggtgtgattgcgactgcgctgcagcgatggcgtggacggcatcggtgaagccaacaccatctgtggcctggcgtatggctacgaaagcccctaagggatcaAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCtt                                                                                                         >  6:62480‑M1/106‑149 (MQ=255)                                    ccTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCtt                                                                                                         >  7:245080/1‑148 (MQ=255)                                        gtggtattgcgactgcgctgcatcgatggcgtggaaggcatcggtgaagccaccaccatcggtggcctggcgtatggctacgaaagccccgaagggatcaAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTGAt                                                                                                      >  3:72427‑M1/101‑147 (MQ=255)                                              gAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTGATCCCTTCgg                                                                                              <  7:249467/149‑1 (MQ=255)                                               aGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTGATCCCTTCggg                                                                                             <  8:38955/149‑1 (MQ=255)                                                     gcgcTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTc                                                                                       >  7:200294/1‑149 (MQ=255)                                                      cgcTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCg                                                                                      >  5:41045/1‑149 (MQ=255)                                                          tgcagcgatggcgtggaaggcatcggtgaagccaccaccatcggtggcctggcgtatggctacgaaagccccgaagggatcaAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTg                                                                                                        >  6:87777‑M1/83‑127 (MQ=255)                                                          tgcagcgatggcgtggaaggcatcggtgaagccaccaccatcggtggcctggcgtatggctacgaaagccccgaagggatcaAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTg                                                                                                        <  5:87777‑M1/45‑1 (MQ=255)                                                           tCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAg                                                                                   <  4:201015/147‑1 (MQ=255)                                                           tCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGc                                                                                  <  8:200294/148‑1 (MQ=255)                                                             gATGCCGGACTTGGCGAAGGTGTTTCCCTTGGGCAGCTTGTCCAGGTTGAGCATGGGGTCATTGAGGTTGTCTGCCGGGAGGCCAATCAAAGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCGTTGATCCCTTCGGGGCTTTCGTAGCCa                                                                                <  6:213104/148‑1 (MQ=255)                                                             gATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTGATCCCTTc                                                                                                <  4:24382/132‑1 (MQ=255)                                                             gATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTGATCCCTTc                                                                                                >  3:24382/1‑132 (MQ=255)                                                             gATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCa                                                                                <  5:47352/148‑1 (MQ=255)                                                                  cGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCAt                                                                               <  6:31040/144‑1 (MQ=255)                                                                  cGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCAt                                                                               >  5:31040/1‑144 (MQ=255)                                                                    gACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAg                                                                        <  8:49991/149‑1 (MQ=255)                                                                         ggCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCAc                                                                   >  4:142694/1‑149 (MQ=255)                                                                         ggCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCAc                                                                   >  6:264231/1‑149 (MQ=255)                                                                          gCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCAcc                                                                  >  7:31962/1‑149 (MQ=255)                                                                          gCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCAc                                                                   <  6:41045/148‑1 (MQ=255)                                                                              aGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGAtg                                                              <  8:245080/149‑1 (MQ=255)                                                                              aGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGAtg                                                              <  5:276619/149‑1 (MQ=255)                                                                              aGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGAtg                                                              <  6:218730/149‑1 (MQ=255)                                                                               ggTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGAtgg                                                             <  1:12/149‑1 (MQ=255)                                                                                 tgtTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGAtggtg                                                           <  1:142897/149‑1 (MQ=255)                                                                                 tgtTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGAtgg                                                             <  3:57152/147‑1 (MQ=255)                                                                                 tgtTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGAtgg                                                             >  4:57152/1‑147 (MQ=255)                                                                                       ccTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGATGGTGGTGGCt                                                     >  3:41055/1‑149 (MQ=255)                                                                                              cAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGAtggtggtgg                                                       >  2:33648/1‑140 (MQ=255)                                                                                              cAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGAtggtggtgg                                                       <  1:33648/140‑1 (MQ=255)                                                                                                gCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTGATCCCTTCg                                                                                               >  1:116729/1‑98 (MQ=255)                                                                                                gCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTGATCCCTTCg                                                                                               <  2:116729/98‑1 (MQ=255)                                                                                                               tcATGGCGGCATTGATGGTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGTCTTATCGCGTCGAAGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGATGGTGGTGGCTTCACCGATGCCTTCCACGCCATc                              <  1:86618/148‑1 (MQ=255)                                                                                                                   ggcggcATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGATGGTGGTGGCTTCACCGATGCCtt                                        <  5:136596/134‑1 (MQ=255)                                                                                                                   ggcggcATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGATGGTGGTGGCTTCACCGATGCCtt                                        >  6:136596/1‑134 (MQ=255)                                                                                                                    gcggcATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGATGGTGGTGGCTTCACCGATGCCTTCCACGCCATCGCTGc                         <  3:142694/148‑1 (MQ=255)                                                                                                                         aTTGATGTTGGCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGATGGTGGTGGCTTCACCGATGCCTTCCACGCCATCGCTGCAGCGCa                   <  4:202642/149‑1 (MQ=255)                                                                                                                             aTGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCTGGCCACCGATGGTGGTGGCTTCACCGATGCCTTCCACGCCATCGCTGCAGCGCAGTCg               >  2:38477/1‑149 (MQ=255)                                                                                                                                       ccGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGATGGTGGTGGCTTCACCGATGCCTTCCACGCCATCGCTGCAGCGCAGTCGCAATaccac      >  6:78958/1‑148 (MQ=255)                                                                                                                                          gCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGATGGTGGTGGCTTCACCGATGCCTTCCACGCCATCGCTGCAGCGCAGTCGCAATACCACCAgg  >  6:205641/1‑149 (MQ=255)                                                                                                                                            |                                                                                                                                                    CAGCAGTTCGCTGACCGGCAGGCCCAGGCGTTTGCCCTGGGCGTCGAGCAAGGCGCTTTCGATGCCGGACTTGGCGAAGGTGTTGCCCTTGGCCAGCTTGTCCAGCTTGAGCATGGCGGCATTGATGTTGTCTGCCGGCAGGCCAATCAACGCTGGGGCGAGGTA‑CGCGTCGATGTTGGCCTTGATCCCTTCGGGGCTTTCGTAGCCATACGCCAGGCCACCGATGGTGGTGGCTTCACCGATGCCTTCCACGCCATCGCTGCAGCGCAGTCGCAATACCACCAGG  >  NC_002947/4237866‑4238151

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.