New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 358813815 (0.520)6 (0.220) 6/254 2.7 22.5% intergenic (+16/‑59) catA‑II/benE‑II catechol 1,2‑dioxygenase/benzoate transport protein
?NC_002947 = 3588208 27 (0.970)coding (12/1272 nt) benE‑II benzoate transport protein

GCGGTGCCGATGAGCAGCGCAGCGGGCGGCCGCGAGCTTTGCAGGAGGCCTGATCCACCGAGTGGGCTGTTTCGCGGGCATGCCCGCTCCCACAGGTGAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGCTCTGTGTCGGCCATCGTCGCCGGCAGCATTGCCACGACGATCTCCTACGCCGGCCCGTTGGTGA  >  NC_002947/3588070‑3588356
                                                                                                                                          |                                                                                                                                                    
gCGGTGCCGATGAGCAGCGCAGCGGGCGGCCGCGAGCTTTGCAGGAGGCCTGATCCACCGAGTGGGCTGTTTCGCGGGCATGCCCGCTCCCACAGGTGAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACAcc                                                                                                                                           <  2:200995/150‑1 (MQ=255)
                     gcgggcggCCGCGAGCTTTGCAGGAGGCCTGATCCACCGAGTGGGCTGTTTCGCGGGCATGCCCGCTCCCACAGGTGAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACaaa                                                                                                                      <  2:39382/150‑1 (MQ=255)
                                    cTTTGCAGGAGGCCTGATCCACCGAGTGGGCTGTTTCGCGGGCATGCCCGCTCCCACAGGTGAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCCGCAAACCCACATACGAGAGCCCGCCACAAAGGCCGGCCGCTGCCt                                                                                                       >  1:523813/1‑150 (MQ=255)
                                                ccTGATCCACCGAGTGGGCTGTTTCGCGGGCATGCCCGCTCCCACAGGTGAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGcagccaactcggtggatcaggcctcctgcaaagctcgcggccgcccgctgcgctgctca                                                                                           <  2:249694‑M2/150‑60 (MQ=255)
                                                        aCCGAGTGGGCTGTTTCGCGGGCATGCCCGCTCCCACAGGTGAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGaaaa                                                                                   >  1:440642/1‑150 (MQ=255)
                                                                 gCTGTTTCGCGGGCATGCCCGCTCCCACAGGTGAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAggcc                                                                                                                  <  2:484920/110‑1 (MQ=255)
                                                                 gCTGTTTCGCGGGCATGCCCGCTCCCACAGGTGAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAggcc                                                                                                                  >  1:484920/1‑110 (MQ=255)
                                                                     tttCGCGGGCATGCCCGCTCCCACAGGTGAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGa                                                                      <  2:210202/150‑1 (MQ=255)
                                                                     tttCGCGGGCATGCCCGCTCCCACAGGTGAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGa                                                                      >  2:297399/1‑150 (MQ=255)
                                                                          cgGGCATGCCCGCTCCCACAGGTGAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGct                                                                 <  2:525687/150‑1 (MQ=255)
                                                                             gCATGCCCGCTCCCACAGGTGAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGCTCtg                                                              >  2:184922/1‑150 (MQ=255)
                                                                               aTGCCCGCTCCCACAGGTGAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGcagcccactcggtggatcaggcctcctgcaaagctcgcggccgcccgctgcgctgctcatcggcaccgccctgggcctgctgcaactgga                                                            <  1:399066‑M2/150‑91 (MQ=255)
                                                                                tGCCCGCTCCCACAGGTGAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGACGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTCCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGCTCTgtgt                                                           >  2:223190/1‑150 (MQ=255)
                                                                                   ccGCTCCCACAGGTGAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGcagcccactcggtggatcaggcctcctgcaaagctcgcggccgcccgctgcgctgctcatc                                                                                         >  1:379108‑M2/1‑56 (MQ=255)
                                                                                   ccGCTCCCACAGGTGAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGcagcccactcggtggatcaggcctcctgcaaagctcgcggccgcccgctgcgctgctcatc                                                                                         <  2:379108‑M2/117‑62 (MQ=255)
                                                                                   ccGCTCCCACAGGTGAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGCTCTGTGTCgg                                                        >  1:579037/1‑150 (MQ=255)
                                                                                      cTCCCACAGGTGAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGCTCTGTGTCGGCCa                                                     >  1:531422/1‑150 (MQ=255)
                                                                                       tCCCACAGGTGAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGcagcccactcggtggatcaggcctcctgcaaagctcgcggccgcccgctgcgctgctcatcggcaccgccctgggcctgctgcaactggaagtcgaac                                                    >  2:399066‑M2/1‑52 (MQ=255)
                                                                                             aGGTGAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCCTGGCAACGCTGATCAAGGACTGCTCTGTGTCGGCCATCGTCGc                                              >  1:371871/1‑150 (MQ=255)
                                                                                              ggTGAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGACACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGCTCTGTGTCGGCCATCGTCGcc                                             >  2:512490/1‑150 (MQ=255)
                                                                                                    ggcTGTGCAGTCGTTGGTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGCTCTGTGTCGGCCATCGTCGCCCgcagc                                       >  2:122221/3‑150 (MQ=255)
                                                                                                        gtgCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGCTCtgt                                                             <  2:693585/124‑1 (MQ=255)
                                                                                                        gtgCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGCTCtgt                                                             >  1:693585/1‑124 (MQ=255)
                                                                                                           cAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATc                                                                           >  2:337891/1‑107 (MQ=255)
                                                                                                           cAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATc                                                                           <  1:337891/107‑1 (MQ=255)
                                                                                                               cGTTGTTGGAGCGGGCATGCCCGCGAAGcagccaactcggtggatcaggcctcctgcaaagctcgcggccgcccgctgcgctgctcatcggcaccgccctgggcctgctgcaactggaagtcgaacttcagctcggcatgccgaccatcc                            >  1:249694‑M2/1‑28 (MQ=255)
                                                                                                                      ttcgCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGCTCTGTGTCGGCCATCGTCGCCGGCAGCATTGCCACGACGATCTcc                     >  1:131957/4‑150 (MQ=255)
                                                                                                                         gCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGCTCTGTGTCGGCCATCGTCGCCGGCAGCATTGCCACGACGATCTCCTAc                  >  1:646151/1‑150 (MQ=255)
                                                                                                                             gaaTGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGCTCTGTGTCGGCCATCGTCGCCGGCAGCATTGCCACGACGATCTCCTACGCCg              >  1:11629/3‑150 (MQ=255)
                                                                                                                                    cgcgAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGCTCTGTGTCGGCCATCGTCGCCGGCAGCATTGCCACGACGATCTCCTACGCCGGCCCGtt       >  1:239427/1‑150 (MQ=255)
                                                                                                                                       gAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGCTCTGTGTCGGCCATCGTCGCCGGCAGCATTGCCACGACGATCTCCTACGCCGGCCCGTTGGt    <  2:440642/150‑1 (MQ=255)
                                                                                                                                         agcagcCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGCTCTGTGTCGGCCATCGTCGCCGGCAGCATTGCCACGACGATCTCCTACGCCGGCCCGTTGGTGa  <  1:184922/150‑1 (MQ=255)
                                                                                                                                         agcagcCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGCTCTGTGTCGGCCATCGTCGCCGGCAGCATTGCCACGACGATCTCCTACGCCGGCCCGTTGGTGa  >  2:631332/1‑150 (MQ=255)
                                                                                                                                          |                                                                                                                                                    
GCGGTGCCGATGAGCAGCGCAGCGGGCGGCCGCGAGCTTTGCAGGAGGCCTGATCCACCGAGTGGGCTGTTTCGCGGGCATGCCCGCTCCCACAGGTGAAGCCTGTGCAGTCGTTGTTGGAGCGGGCATGCCCGCGAAGCAGCCAACACCCAAACAAGAGCCCGCCACAAAGGCCGGCCGCTGCCTTGCGAGAAGCTTCATGAAAACGCTGATCAAGGACTGCTCTGTGTCGGCCATCGTCGCCGGCAGCATTGCCACGACGATCTCCTACGCCGGCCCGTTGGTGA  >  NC_002947/3588070‑3588356

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