New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 45520173 (0.190)4 (0.330) 4/190 NT 60.5% intergenic (+67/+114) pydA/rutR NADP‑dependent dihydropyrimidine dehydrogenase subunit PreA/transcriptional regulator, TetR family
?NC_002947 = 4552068 3 (0.250)intergenic (+118/+63) pydA/rutR NADP‑dependent dihydropyrimidine dehydrogenase subunit PreA/transcriptional regulator, TetR family

GCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCTGGTCACCGTCTGCACCGCATTCTCGAAG  >  NC_002947/4551923‑4552194
                                                                                                                                                 |                                                                                                                              
gcgTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcgg                                                                                                                            >  1:205615‑M2/1‑146 (MQ=255)
              gTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggcccgatcgcgacac                                                                                                              >  2:331939‑M2/1‑132 (MQ=255)
                gAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggcccgatcgcgacacaa                                                                                                            >  2:40920‑M2/1‑130 (MQ=255)
                                                                                tCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggcccgatcgcgacacaaggccgctcccacaggggccgcgtcagctgttcaggcggtgacccgcttcaactggcctcacggt                                            <  1:8789‑M2/150‑85 (MQ=255)
                                                                                   cgATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggcccgatcgcgacacaaggccgctcccacaggggccgcgtcagctgttcaggcggtgacccgcttcaactggcctcacggtagg                                         >  2:8789‑M2/1‑63 (MQ=255)
                                                                                        gggCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCgcggcccgatcgcgacacaaggccgctcccacaggggccgcgtcagctgttcaggcggtgacccgcttcaactggcctcacggtaggggtta                                    <  1:218164‑M2/150‑93 (MQ=255)
                                                                                         ggCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCa                                                                                      >  1:298486/1‑99 (MQ=255)
                                                                                         ggCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCa                                                                                      <  2:298486/99‑1 (MQ=255)
                                                                                                       cccAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACaaaaa                                                                                           <  1:292630/80‑1 (MQ=255)
                                                                                                       cccAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACaaaaa                                                                                           >  2:292630/1‑80 (MQ=255)
                                                                                                            cagTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCTGGTCACCGTCAGCa                >  2:88119/1‑150 (MQ=255)
                                                                                                                          aTACAGACATCGAGGGGCCGCTTCgcggcccgatcgcgacacaaggccgctcccacaggggccgcgtcagctgttcaggcggtgacccgcttcaactggcctcacggtaggggttacgcggatcctgcgtccagttcaggtacggcttgc  >  2:218164‑M2/1‑24 (MQ=255)
                                                                                                                                                 |                                                                                                                              
GCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCCCCCGCAAGATCACGCTGGTCACCGTCTGCACCGCATTCTCGAAG  >  NC_002947/4551923‑4552194

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