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

TCGGCGATCTGTGCAGTGCGCAGAACCGCGCGCGGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGC  >  NC_002947/4237422‑4237719
                                                                                                                                                     |                                                                                                                                                      
tCGGCGATCTGTGCAGTGCGCAGAACCGCGCGCGGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTg                                                                                                                                                        <  1:1146454/150‑1 (MQ=255)
 gtcgcgatcaagcgcgagctgggtgacagcgccagcgtgcgggtcgacgtcaaccagtactgggacgagtcccaggccatccgcgcctgccaggtattgggcgacaacggcatcgacctgatcgagcagccgatttcgcgcatcaaccgc                                                                                                                                                       <  1:154151‑M1/2‑1 (MQ=255)
  ggCGATCTGTGCAGTGCGCAGAACCGCGCGCGGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTgcg                                                                                                                                                      >  1:832270/1‑150 (MQ=255)
   gCGATCTGTGCAGTGCGCAGAACCGCGCGCGGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTgcgc                                                                                                                                                     >  1:400844/1‑150 (MQ=255)
    gcgatcaagcgcgagctgggtgacagcgccagcgtgcgggtcgacgtcaaccagtactgggacgagtcccaggccatccgcgcctgccaggtattgggcgacaacggcatcgacctgatcgagcagccgatttcgcgcatcaaccGCGCt                                                                                                                                                    <  2:533747‑M1/5‑1 (MQ=255)
    cGATCTGTGCAGTGCGCAGAACCGCGCGCGGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCt                                                                                                                                                    <  1:739740/150‑1 (MQ=255)
             cAGTGCGCAGAACCGCGCGCGGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGgc                                                                                                                                           <  2:327889/150‑1 (MQ=255)
                        aCCGCGCGCGGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCa                                                                                                                              <  1:293764/150‑1 (MQ=255)
                                cgGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGtt                                                                                                                      <  1:751975/150‑1 (MQ=255)
                                 cagcgtgcgggtcgacgtcaaccagtactgggacgagtcccaggccatccgcgcctgccaggtattgggcgacaacggcatcgacctgatcgagcagccgatttcgcgcatcaaccgc                                                                                                                                                       >  1:161196‑M1/117‑118 (MQ=255)
                                 cagcgtgcgggtcgacgtcaaccagtactgggacgagtcccaggccatccgcgcctgccaggtattgggcgacaacggcatcgacctgatcgagcagccgatttcgcgcatcaaccgc                                                                                                                                                       <  2:161196‑M1/2‑1 (MQ=255)
                                       ccaTTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATgcgc                                                                                                               <  1:1010126/150‑1 (MQ=255)
                                       ccaTTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATgcgc                                                                                                               <  2:519539/150‑1 (MQ=255)
                                        caTTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGcgcg                                                                                                              >  2:36414/1‑150 (MQ=255)
                                        caTTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGcgcg                                                                                                              >  1:690275/1‑150 (MQ=255)
                                                gCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGgc                                                                                                      <  2:832270/150‑1 (MQ=255)
                                                             cGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGt                                                                                         <  1:53776/150‑1 (MQ=255)
                                                                 gATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGAt                                                                                     <  2:1054071/150‑1 (MQ=255)
                                                                     cTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCg                                                                                 <  1:890451/150‑1 (MQ=255)
                                                                        ccaggccatccgcgcctgccaggtattgggcgacaacggcatcgacctgatcgagcagccgatttcgcgcatcaaccGCGCTGGTTCAGGCGCACCTGGCCAgcgc                                                                                                                          >  1:1191495‑M1/78‑106 (MQ=255)
                                                                        ccaggccatccgcgcctgccaggtattgggcgacaacggcatcgacctgatcgagcagccgatttcgcgcatcaaccGCGCTGGTTCAGGCGCACCTGGCCAgcgc                                                                                                                          <  2:1191495‑M1/29‑1 (MQ=255)
                                                                               cGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCa                                                                       >  1:45345/1‑150 (MQ=255)
                                                                                 tCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAAt                                                                     <  1:895415/150‑1 (MQ=255)
                                                                                   ggcggcCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATAc                                                                   <  2:400844/150‑1 (MQ=255)
                                                                                      ggcCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCgg                                                                                                                        >  1:550449/1‑94 (MQ=255)
                                                                                      ggcCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCgg                                                                                                                        <  2:550449/94‑1 (MQ=255)
                                                                                       ctgccaggtattgggcgacaacggcatcgacctgatcgagcagccgatttcgcgcatcaaccGCGCt                                                                                                                                                    <  1:58466‑M1/5‑1 (MQ=255)
                                                                                       ctgccaggtattgggcgacaacggcatcgacctgatcgagcagccgatttcgcgcatcaaccGCGCt                                                                                                                                                    >  2:58466‑M1/63‑67 (MQ=255)
                                                                                              tGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGgcgc                                                        >  1:562043/1‑150 (MQ=255)
                                                                                               gAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGcgcg                                                       >  2:429068/1‑150 (MQ=255)
                                                                                                                 ttCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTc                                     <  2:45345/150‑1 (MQ=255)
                                                                                                                      tcgaCTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACGAGTCCCAgg                                <  2:690275/150‑2 (MQ=255)
                                                                                                                           tCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGgctgct                                                                                                  <  2:1185081/79‑1 (MQ=255)
                                                                                                                           tCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGgctgct                                                                                                  >  1:1185081/1‑79 (MQ=255)
                                                                                                                             aTCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGt                                                                                         <  1:415507/86‑1 (MQ=255)
                                                                                                                             aTCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGt                                                                                         >  2:415507/1‑86 (MQ=255)
                                                                                                                             aTCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTgg                           >  1:1017605/1‑148 (MQ=255)
                                                                                                                             aTCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTgg                           <  2:1017605/148‑1 (MQ=255)
                                                                                                                                   cATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTc                   >  1:1087464/1‑150 (MQ=255)
                                                                                                                                    aTGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCg                  >  2:154151/1‑150 (MQ=255)
                                                                                                                                          ggAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGc            <  1:36414/150‑1 (MQ=255)
                                                                                                                                                caaccGCGCTGGCCAGGTGCGCCTGAACCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCa                                  >  2:661837‑M1/6‑122 (MQ=255)
                                                                                                                                                caaccGCGCTGGCCAGGTGCGCCTGAACCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCa                                  <  1:661837‑M1/117‑1 (MQ=255)
                                                                                                                                                 gACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTgg     >  1:533747/1‑150 (MQ=255)
                                                                                                                                                  aCTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGgc    <  2:562043/150‑1 (MQ=255)
                                                                                                                                                    tGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGgcgc  >  1:35858/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TCGGCGATCTGTGCAGTGCGCAGAACCGCGCGCGGGCCACCATTCTTGGCGATTTTCAGGGCGAAGATGCTGGCGGCGCCGTCGGCGGCCAGGCTGAAGGCGTCCTCGACGCTTTCGATCGACTCATCGGCCATGATCGGAGCCGGACTGCGCTGGTTCAGGCGCACCTGGCCAGCGCGGTTGATGCGCGAAATCGGCTGCTCGATCAGGTCGATGCCGTTGTCGCCCAATACCTGGCAGGCGCGGATGGCCTGGGACTCGTCCCAGTACTGGTTGACGTCGACCCGCACGCTGGCGC  >  NC_002947/4237422‑4237719

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.