New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 152670711 (0.660)4 (0.260) 3/268 1.8 24.8% coding (852/957 nt) ddlB D‑alanine‑‑D‑alanine ligase B
?NC_002947 = 1526720 14 (0.900)coding (865/957 nt) ddlB D‑alanine‑‑D‑alanine ligase B

GAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACC  >  NC_002947/1526579‑1526856
                                                                                                                                |                                                                                                                                                     
gaacaggaacTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAgcctggtg                                                                                                                                  <  2:320700‑M2/150‑9 (MQ=255)
       aacTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCtatgatctgtcatgcctggtgcggtgtt                                                                                                                           <  1:310298‑M1/150‑29 (MQ=255)
       aacTGATCGACCTGACCGCGCGCGCCTGCGATGCCAGTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCtatgatctgtcatgcctggtgcggtgtt                                                                                                                           <  1:310304‑M1/150‑29 (MQ=255)
                aCCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTg                                                                                                                  >  2:43040/1‑150 (MQ=255)
                   tGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCgg                                                                                                               >  2:187105/1‑150 (MQ=255)
                                cttcGATGCCAATGGCATCGCAGGCTGTGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTg                                                                                                       >  2:265057/4‑145 (MQ=255)
                                cttcGATGCCAATGGCATCGCAGGCTGTGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTg                                                                                                       <  1:265057/142‑1 (MQ=255)
                                          aTTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTg                                                                                        <  1:43040/150‑1 (MQ=255)
                                                aTCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCGGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTg                                                                                  <  1:277504/150‑1 (MQ=255)
                                                                             gCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCtatgatctgtcatgcctggtgcgg                                                                                                                               <  1:349279‑M1/76‑25 (MQ=255)
                                                                             gCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCtatgatctgtcatgcctggtgcgg                                                                                                                               >  2:349279‑M1/1‑52 (MQ=255)
                                                                                          ggCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGt                                                 >  2:307666/1‑141 (MQ=255)
                                                                                          ggCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGt                                                 <  1:307666/141‑1 (MQ=255)
                                                                                                       tgctCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCAAGGTTAATAACCATGCAAGGCGCgat                           <  1:323439/150‑1 (MQ=255)
                                                                                                              aGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGt                    <  1:73919/150‑1 (MQ=255)
                                                                                                                  aaCACCGCACCAGGCATGACAGATCATAgcctggtgcggtgttgacttcgagcagccagaaccggccctgctcgtcctgcatcacgtccaaccgcccccagccttcgatgccactggcatcgcaggcgcgcgcggtcaggtcgatcagtt                <  1:310304‑M2/150‑123 (MQ=255)
                                                                                                                  aaCACCGCACCAGGCATGACAGATCATAgcctggtgcggtgttgacttcgagcagccagaaccggccctgctcgtcctgcatcacgtccaaccgcccccagccttcgatgccaatggcatcgcaggcgcgcgcggtcaggtcgatcagtt                <  1:310298‑M2/150‑123 (MQ=255)
                                                                                                                      ccGCACCAGGCATGACAGATCATAgcctggtgcggtgttgacttcgagcagccagaaccggccctgctcgtcctgc                                                                                      <  1:349279‑M2/76‑53 (MQ=255)
                                                                                                                      ccGCACCAGGCATGACAGATCATAgcctggtgcggtgttgacttcgagcagccagaaccggccctgctcgtcctgc                                                                                      >  2:349279‑M2/1‑24 (MQ=255)
                                                                                                                         cACCAGGCtatgatctgtcatgcctggtgcggtgttgacttcgagcagccagaaccggccctgctcgtcctgcatcacgtccaaccgcccccagccttcgatgccaatggcatcgcaggcgcgcgcggtcaggtcgatcagttcctgttc         <  2:320700‑M1/150‑143 (MQ=255)
                                                                                                                                cATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTAcc  <  1:48874/150‑1 (MQ=255)
                                                                                                                                |                                                                                                                                                     
GAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACC  >  NC_002947/1526579‑1526856

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