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

GAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTA  >  NC_002947/1526579‑1526863
                                                                                                                                             |                                                                                                                                               
gaacaggaacTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAgcctggtg                                                                                                                                         <  2:320700‑M2/150‑9 (MQ=255)
                aCCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTg                                                                                                                         >  2:43040/1‑150 (MQ=255)
                   tGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCgg                                                                                                                      >  2:187105/1‑150 (MQ=255)
                                cttcGATGCCAATGGCATCGCAGGCTGTGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTg                                                                                                              <  1:265057/142‑1 (MQ=255)
                                cttcGATGCCAATGGCATCGCAGGCTGTGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTg                                                                                                              >  2:265057/4‑145 (MQ=255)
                                          aTTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTg                                                                                               <  1:43040/150‑1 (MQ=255)
                                                aTCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCGGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTg                                                                                         <  1:277504/150‑1 (MQ=255)
                                                                                          ggCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGt                                                        <  1:307666/141‑1 (MQ=255)
                                                                                          ggCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGt                                                        >  2:307666/1‑141 (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)
                                                                                                                                cATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTAcc         <  1:48874/150‑1 (MQ=255)
                                                                                                                                    aCAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGcc     >  1:247942/1‑150 (MQ=255)
                                                                                                                                       gATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTa  <  1:317796/150‑1 (MQ=255)
                                                                                                                                       gATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTa  <  1:318344/150‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                               
GAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACTGGTGCTGGCAATCCTGGCCGAAAGCGTGGCGACGCGAGGTTAATAACCATGCAAGGCGCGATGATACGTCAGCAGCAACCCGTTACCGGCCGTA  >  NC_002947/1526579‑1526863

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.