New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 1573491 =26 (1.120)4 (0.180) 4/242 NT 13.2% coding (207/792 nt) pcaD 3‑oxoadipate enol‑lactonase
?NC_002947 1573532 = 28 (1.260)coding (248/792 nt) pcaD 3‑oxoadipate enol‑lactonase

TGTCCAACTCGCTGGGTACCGACCTGGGCATGTGGGACACCCAGATTCCGCTCTGGAGTCAGCACTTCCGGGTGCTGCGCTATGACACCCGTGGTCACGGCGCATCGCTGGTCACTGAAGGCCCTTACAGCATCGAACAGCTGGGCCGCGACGTGCTGGCCCTGCTCGATGGCCTGGACATTCAAAAGGCTCACTTCGTCGGCCTGTCGATGGGCGGCCTGATCGGCCAGTGGCTGGGTATCCATGCAGGTGAGCGCCTGCACAGCCTGACCCTGTGC  >  NC_002947/1573361‑1573638
                                                                                                                                  |                                                                                                                                                   
tGTCCAACTCGCTGGGTACCGACCTGGGCATGTGGGACACCCAGATTCCGCTCTGGAGTCAGCACTTCCGGGTGCTGCGCTATGACACCCGTGGTCACGGCGCATCGCTGGTCACTGAAGGCCCTTACAGCATCGAACAGCTGGGCcgcg                                                                                                                                  >  2:352613/1‑150 (MQ=255)
             cggccgtgttgcacagggtcaggctgtgcaggcgctcacctgcatggatacccagccactggccgatcaggccgcccatcgacaggccgacgaattgagccttttgaatgtcccggcCATCGA‑‑‑GCTGTGCCGCGACGTGCTGGCCCTGCt                                                                                                                  >  1:321493‑M1/118‑150 (MQ=255)
                     ttgcacagggtcaggctgtgcaggcgctcacctgcatggatacccagccactggccgatcaggccgcccatcgacaggccgacgaagtgagccttttgaatgtccaggcCATCGAACAGCTGGGCCGCGACGTGCTGGCCCTGCTCGATg                                                                                                             >  1:382628‑M1/110‑150 (MQ=255)
                     aCCTGGGCATGTGGGACACCCAGATTCCGCTCTGGAGTCAGCACTTCCGGGTGCTGCGCTATGACACCCGTGGTCACGGCGCATCGCTGGTCACTGAAGGCCCTTACAGCATCGAACAGCTGGGCCGCGACGTGCTGGCCCTGCTCGATg                                                                                                             <  2:59781/150‑1 (MQ=255)
                       cTGGGCATGTGGGACACCCAGATTCCGCTCTGGAGTCAGCACTTCCGGGTGCTGCGCTATGACACCCGTGGTCACGGCGCATCGCTGGTCACTGAAGGCCCTTACAGCATCGAACAGCTGGGCCGCGACGTGCTGGCCCTGCTCGATGGc                                                                                                           >  2:166995/1‑150 (MQ=255)
                            cATGTGGGACACCCAGATTCCGCTCTGGAGTCAGCACTTCCGGGTGCTGCGCTATGACACCCGTGGTCACGGCGCATCGCTGGTCACTGAAGGCCCTTACAGCATCGAACAGCTGGGCCGCGACGTGCTGGCCCTGCTCGATGGCCTGGa                                                                                                      >  1:121898/1‑150 (MQ=255)
                               tcaggctgtgcaggcgctcacctgcatggatacccagccactggccgatcaggccgcccatcgacaggccgacgaagtgagccttttgaatgtccaggcCATCGAACAGCTGGGCCGCGACGTGCTGGCCCTGCTCGATgg                                                                                                            <  2:356461‑M1/42‑1 (MQ=255)
                               tcaggctgtgcaggcgctcacctgcatggatacccagccactggccgatcaggccgcccatcgacaggccgacgaagtgagccttttgaatgtccaggcCATCGAACAGCTGGGCCGCGACGTGCTGGCCCTGCTCGATgg                                                                                                            >  1:356461‑M1/100‑141 (MQ=255)
                                tgGGACACCCAGATTCCGCTCTGGAGTCAGCACTTCCGGGTGCTGCGCTATGACACCCGTGGTCACGGCGCATCGCTGGTCACTGAAGGCCCTTACAGCATCGAACAGCTGGGCCGCGACGTGCTGGCCCTGCTCGATGGCCTGGACAtt                                                                                                  >  2:5515/1‑150 (MQ=255)
                                  ggACACCCAGATTCCGCTCTGGAGTCAGCACTTCCGGGTGCTGCGCTATGACACCCGTGGTCACGGCGCATCGCTGGTCACTGAAGGCCCTTACAGCATCGAACAGCTGGGCCGCGACGTGCTGGCCCTGCTCGATGGCCTGGACATTCa                                                                                                <  1:5515/150‑1 (MQ=255)
                                        ccAGATTCCGCTCTGGAGTCAGCACTTCCGGGTGCTGCGCTATGACACCCGTGGTCACGGCGCATCGCTGGTCACTGAAGGCCCTTACAGCATCGAACAGCTGGGCCGCGACGTGCTGGCCCTGCTCGATGGCCTGGACATTCaaaa                                                                                             >  1:23124/1‑147 (MQ=255)
                                        ccAGATTCCGCTCTGGAGTCAGCACTTCCGGGTGCTGCGCTATGACACCCGTGGTCACGGCGCATCGCTGGTCACTGAAGGCCCTTACAGCATCGAACAGCTGGGCCGCGACGTGCTGGCCCTGCTCGATGGCCTGGACATTCaaaa                                                                                             <  2:23124/147‑1 (MQ=255)
                                                 gCTCTGGAGTCAGCACTTCCGGGTGCTGCGCTATGACACCCGTGGTCACGGCGCATCGCTGGTCACTGAAGGCCCTTACAGCATCGAACAGCTGGGCCGCGACGTGCTGGCCCTGCTCGATGGCCTGGACATTCAAAAGGCTCACTtcgc                                                                                 >  2:193081/1‑149 (MQ=255)
                                                    ctGGAGTCAGCACTTCCGGGTGCTGCGCTATGACACCCGTGGTCACGGCGCATCGCTGGTCACTGAAGGCCCTTACAGCATCGAACAGCTGGGCCGCGACGTGCTGGCCCTGCTCGATGGCCTGGACATTCAAAAGGCTCACTTCGTCgg                                                                              >  2:149456/1‑150 (MQ=255)
                                                              cACTTCCGGGTGCTGCGCTATGACACCCGTGGTCACGGCGCATCGCTGGTCACTGAAGGCCCTTACAGCATCGAACAGCTGGGCCGCGACGTGCTGGCCCTGCTCGATGGCCTGGACATTCAAAAGGCTCACTTCGTCGGCCTGTCGATg                                                                    <  1:149456/150‑1 (MQ=255)
                                                                      ggTGCTGCGCTATGACACCCGTGGTCACGGCGCATCGCTGGTCACTGAAGGCCCTTACAGCATCGAACAGCTGGGCCGCGACGTGCTGGCCCCGCTCGATGGCCTGGACATTCAAAAGGCTCACTTCGTCGGCCTGTCGATGGGCGGACt                                                            <  1:352613/150‑1 (MQ=255)
                                                                       gTGCTGCGCTATGACACCCGTGGTCACGGCGCATCGCTGGTCACTGAAGGCCCTTACAGCATCGAACAGCTGGGCCGCGACGTGCTGGCCCTGCTCGATGGCCTGGACATTCAAAAGGCTCACTTCGTCGGCCTGTCGATGGGCGGCCTg                                                           >  2:166152/1‑150 (MQ=255)
                                                                                       cccGTGGTCACGGCGCATCGCTGGTCACTGAAGGCCCTTACAGCATCGAACAGCTGGGCCGCGACGTGCTGGCCCTGCTCGATGGCCTGGACATTCAAAAGGCTCACTTCGTCGGCCTGTCGATGGGCGGCCTGATCGGCCAGTGGCTgg                                           <  1:166152/150‑1 (MQ=255)
                                                                                        ccGTGGTCACGGCGCATCGCTGGTCACTGAAGGCCCTTACAGCATCGAACAGCTGGGCCGCGACGTGCTGGCCCTGCTCGATGGCCTGGACATTCAAAAGGCTCAc                                                                                      <  2:83064/106‑1 (MQ=255)
                                                                                        ccGTGGTCACGGCGCATCGCTGGTCACTGAAGGCCCTTACAGCATCGAACAGCTGGGCCGCGACGTGCTGGCCCTGCTCGATGGCCTGGACATTCAAAAGGCTCAc                                                                                      >  1:83064/1‑106 (MQ=255)
                                                                                               cACGGCGCATCGCTGGTCACTGAAGGCCCTTACAGCATCGAACAGCTGGGCCGCGACGTGCTGGCCCTGCTCGATGGCCTGGACATTCAAAAGGCTCACTTCGTCGGCCTGTCGATGGGCGGCCTGATCGGCCAGTGGCTGGGTATCCAt                                   >  1:218572/1‑150 (MQ=255)
                                                                                                       aTCGCTGGTCACTGAAGGCCCTTACAGCATCGAACAGCTGGGCCGCGACGTGCTGGCCCTGCTCGATGGCCTGGACATTCAAAAGGCTCACTTCGTCGGCCTGTCGATGGGCGGCCTGATCGGCCAGTGGCTGGGTATCCATGCAGGTGa                           <  2:121898/150‑1 (MQ=255)
                                                                                                             ggTCACTGAAGGCCCTTACAGCATCGAACAGCTGGGCCGCGACGTGCTGGCCCTGCTCGATGGCCTGGACATTCAAAAGGCTCACTTCGTCGGCCTGTCGATGGGCGGCCTGATCGGCCAGTGGCTGGGTATCCATGCAGGTGAGCGCCt                     <  2:218572/150‑1 (MQ=255)
                                                                                                              gTCACTGAAGGCCCTTACAGCATCGAACAGCTGGGCCGAGACGTGCTGGATATGCTTGATGGCCTTGACATTCAAAAGGCTCACTTCGTCGGCCTGTCGATGGGCGGCCTGATCGGCCAGTGGCTGGGTATCCATGCAGGTGAGCGCCTg                    >  1:168162/1‑150 (MQ=255)
                                                                                                               tCACTGAAGGCCCTTACAGCATCGAACAGCTGGGCCGCGACGTGCTGGCCCTGCTCGATGGCCTGGACATTCAAAAGGCTCACTTCGTCGGCCTGTCGATGGGCGGCCTGATCGGCCAGTGGCTGGGTATCCATGCAGGTGAGCGCCTGc                   >  1:492124/1‑150 (MQ=255)
                                                                                                                     aaGGCCCTTACAGCATCGAACAGCTGGGCCGCGACGTGCTGGCCCTGCTCGATGGc                                                                                                           >  1:29271/1‑56 (MQ=255)
                                                                                                                     aaGGCCCTTACAGCATCGAACAGCTGGGCCGCGACGTGCTGGCCCTGCTCGATGGc                                                                                                           <  2:29271/56‑1 (MQ=255)
                                                                                                                        gCCCTTACAGCATCGAACAGCTGGGCCGCGACGTGCTGGCCCTGCTCGATGGCCTGGACATTCAAAAGGCTCACTTCGTCGGCCTGTCGATGGGCGGCCTGATCGGCCAGTGGCTGGGTATCCATGCAGGTGAGCGCCTGCACAGCCTGa          >  1:400249/1‑150 (MQ=255)
                                                                                                                          ccTTACAGCATCGAACAGCTGGGCCGCGACGTGCTGGCCCTGCTCGATGGCCTGGACATTCAAAAGGCTCACTTCGTCGGCCTGTCGATGGGCGGCCTGATCGGCCAGTGGCTGGGTATCCATGCAGGTGAGCGCCTGCACAGCCTGAcc        >  2:610139/1‑150 (MQ=255)
                                                                                                                           cTTACAGCATCGAACAGCTGGGCCGCGACGTGCTGGCCCTGCTCGATGGCCTGGACATTCAAAAGGCTCACTTCGTCGGCCTGTCGATGGGCGGCCTGATCGGCCAGTGGCTGGGTATCCATGCAGGTGAGCGCCTGCACAGCCTGAccc       >  1:366468/1‑150 (MQ=255)
                                                                                                                                aGCATCGAACAGCTGGGCCGCGACGTGCTGGCCCTGCTCGATGGCCTGGACATTCAAAAGGCTCACTTCGTCGGCCTGTCGATGGGCGGCCTGATCGGCCAGTGGCTgg                                           >  2:214671/1‑109 (MQ=255)
                                                                                                                                aGCATCGAACAGCTGGGCCGCGACGTGCTGGCCCTGCTCGATGGCCTGGACATTCAAAAGGCTCACTTCGTCGGCCTGTCGATGGGCGGCCTGATCGGCCAGTGGCTgg                                           <  1:214671/109‑1 (MQ=255)
                                                                                                                                aGCATCGAACAGCTGGGCCGCGACGTGCTGGCCCTGCTCGATGGCCTGGACATTCAAAAGGCTCACTTCGTCGGCCTGTCGATGGGCGGCCTGATCGGCCAGTGGCTGGGTATCCATGCAGGTGAGCGCCTGCACAGCCTGACCCTGTGc  >  1:383816/1‑150 (MQ=255)
                                                                                                                                  |                                                                                                                                                   
TGTCCAACTCGCTGGGTACCGACCTGGGCATGTGGGACACCCAGATTCCGCTCTGGAGTCAGCACTTCCGGGTGCTGCGCTATGACACCCGTGGTCACGGCGCATCGCTGGTCACTGAAGGCCCTTACAGCATCGAACAGCTGGGCCGCGACGTGCTGGCCCTGCTCGATGGCCTGGACATTCAAAAGGCTCACTTCGTCGGCCTGTCGATGGGCGGCCTGATCGGCCAGTGGCTGGGTATCCATGCAGGTGAGCGCCTGCACAGCCTGACCCTGTGC  >  NC_002947/1573361‑1573638

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