New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_002947	4237571 =	29 (0.640)	2 (0.050)	4/234	NT	6.3%	coding (723/1122 nt)	catB	muconate cycloisomerase 1
	?	NC_002947	4237600 =	32 (0.750)					coding (694/1122 nt)	catB	muconate cycloisomerase 1

CGGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTG‑G‑CCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGAC  >  NC_002947/4237454‑4237747                                                                                                                                                     |                                                                                                                                                    cgGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTG‑G‑CCAGCGCGGtt                                                                                                                                                  <  1:751975/150‑1 (MQ=255)        ccaTTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTG‑G‑CCAGCGCGGTTGATgcgc                                                                                                                                           <  2:519539/150‑1 (MQ=255)        ccaTTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTG‑G‑CCAGCGCGGTTGATgcgc                                                                                                                                           <  1:1010126/150‑1 (MQ=255)         caTTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTG‑G‑CCAGCGCGGTTGATGcgcg                                                                                                                                          >  2:36414/1‑150 (MQ=255)         caTTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTG‑G‑CCAGCGCGGTTGATGcgcg                                                                                                                                          >  1:690275/1‑150 (MQ=255)                 gCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTG‑G‑CCAGCGCGGTTGATGCGCGAAATCGgc                                                                                                                                  <  2:832270/150‑1 (MQ=255)                              cGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTG‑G‑CCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGt                                                                                                                     <  1:53776/150‑1 (MQ=255)                                tgggacgagtcccaggccatccgcgcctgccaggtattgggcgacaacggcatcgacctgatcgagcagccgatttcgcgcatcaaccgcgctggttcaggcgcacctggccagcgcGGTTg                                                                                                                                                 >  2:661837‑M2/118‑122 (MQ=255)                                tgggacgagtcccaggccatccgcgcctgccaggtattgggcgacaacggcatcgacctgatcgagcagccgatttcgcgcatcaaccgcgctggttcaggcgcacctggccagcgcGGTTg                                                                                                                                                 <  1:661837‑M2/5‑1 (MQ=255)                                  gATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTG‑G‑CCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGAt                                                                                                                 <  2:1054071/150‑1 (MQ=255)                                      cTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTG‑G‑CCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCg                                                                                                             <  1:890451/150‑1 (MQ=255)                                                cGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTG‑G‑CCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCa                                                                                                   >  1:45345/1‑150 (MQ=255)                                                  tCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTG‑G‑CCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAAt                                                                                                 <  1:895415/150‑1 (MQ=255)                                                    ggcggcCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTG‑G‑CCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATAc                                                                                               <  2:400844/150‑1 (MQ=255)                                                       ggcCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTG‑G‑CCAGCGCgg                                                                                                                                                    >  1:550449/1‑94 (MQ=255)                                                       ggcCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTG‑G‑CCAGCGCgg                                                                                                                                                    <  2:550449/94‑1 (MQ=255)                                                               tGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTG‑G‑CCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGgcgc                                                                                    >  1:562043/1‑150 (MQ=255)                                                                gAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTG‑G‑CCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGcgcg                                                                                   >  2:429068/1‑150 (MQ=255)                                                                                  ttCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTG‑G‑CCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTc                                                                 <  2:45345/150‑1 (MQ=255)                                                                                       tcgaCTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTG‑G‑CCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACGAGTCCCAgg                                                            <  2:690275/150‑2 (MQ=255)                                                                                            tCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTG‑G‑CCAGCGCGGTTGATGCGCGAAATCGgctgct                                                                                                                              <  2:1185081/79‑1 (MQ=255)                                                                                            tCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTG‑G‑CCAGCGCGGTTGATGCGCGAAATCGgctgct                                                                                                                              >  1:1185081/1‑79 (MQ=255)                                                                                              aTCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTG‑G‑CCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGt                                                                                                                     >  2:415507/1‑86 (MQ=255)                                                                                              aTCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTG‑G‑CCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGt                                                                                                                     <  1:415507/86‑1 (MQ=255)                                                                                              aTCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTG‑G‑CCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTgg                                                       >  1:1017605/1‑148 (MQ=255)                                                                                              aTCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTG‑G‑CCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTgg                                                       <  2:1017605/148‑1 (MQ=255)                                                                                                    cATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTG‑G‑CCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTc                                               >  1:1087464/1‑150 (MQ=255)                                                                                                     aTGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTG‑G‑CCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCg                                              >  2:154151/1‑150 (MQ=255)                                                                                                           ggAGCCGGACTGCGCTGGTTCAGGCGCACCTG‑G‑CCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGc                                        <  1:36414/150‑1 (MQ=255)                                                                                                                 caaccGCGCTGGC‑CAGGTGCGCCTGAA‑CCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCa                                                              >  2:661837‑M1/6‑122 (MQ=255)                                                                                                                 caaccGCGCTGGC‑CAGGTGCGCCTG‑AACCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCa                                                              <  1:661837‑M1/117‑1 (MQ=255)                                                                                                                  gACTGCGCTGGTTCAGGCGCACCTG‑G‑CCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTgg                                 >  1:533747/1‑150 (MQ=255)                                                                                                                   aCTGCGCTGGTTCAGGCGCACCTG‑G‑CCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGgc                                <  2:562043/150‑1 (MQ=255)                                                                                                                     tGCGCTGGTTCAGGCGCACCTG‑G‑CCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGgcgc                              >  1:35858/1‑150 (MQ=255)                                                                                                                        gcgctggccaggtgcgcctgaaccagcgcGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTgg                                                                         <  2:1191495‑M2/77‑1 (MQ=255)                                                                                                                        gcgctggccaggtgcgcctgaaccagcgcGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTgg                                                                         >  1:1191495‑M2/30‑106 (MQ=255)                                                                                                                                    gcACCTG‑G‑CCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGt                                                <  2:492305/117‑1 (MQ=255)                                                                                                                                    gcACCTG‑G‑CCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGt                                                >  1:492305/1‑117 (MQ=255)                                                                                                                                    gcACCTG‑G‑CCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCa                         <  2:224454/140‑1 (MQ=255)                                                                                                                                    gcACCTG‑G‑CCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCa                         >  1:224454/1‑140 (MQ=255)                                                                                                                                      aCCTG‑G‑CCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCgcgc             <  1:429068/150‑1 (MQ=255)                                                                                                                                       ccTG‑G‑CCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCt            >  1:227313/1‑150 (MQ=255)                                                                                                                                                agcgcGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAg                                                                                        >  2:58466‑M2/6‑67 (MQ=255)                                                                                                                                                agcgcGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAg                                                                                        <  1:58466‑M2/62‑1 (MQ=255)                                                                                                                                                agcgcGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATcgc     <  2:533747‑M2/145‑1 (MQ=255)                                                                                                                                                   gcGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTg                                  <  2:161196‑M2/116‑1 (MQ=255)                                                                                                                                                   gcGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTg                                  >  1:161196‑M2/3‑118 (MQ=255)                                                                                                                                                   gcGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGAc  <  1:154151‑M2/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CGGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTG‑G‑CCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGCTGTCACCCAGCTCGCGCTTGATCGCGAC  >  NC_002947/4237454‑4237747

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