New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 2300506 =14 (0.610)4 (0.180) 4/236 NT 24.8% coding (894/1239 nt) sbcD exonuclease subunit SbcD
?NC_002947 2300548 = 11 (0.500)coding (852/1239 nt) sbcD exonuclease subunit SbcD

CTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGTTACCGTCCAACTCCACTTCCAGGACCTGGTGCGGGTAGTTGATTTCGGCGAACGACAGCGGGATTGGCGAGCCGCTGTAGCGGATGCGCTCTTCGCGGTTG  >  NC_002947/2300443‑2300695
                                                                                                         |                                                                                                                                                   
cTCCAGCCAGGGCAGGCGGTGGGGATCTTCAAGCATGTCGACCACGGGCGGGGCGGGCAGTTGCTGCAACAGCACGCCCAGCGGTTCCCGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACCGGGCGTGGCTCGACGCCGACCAGGcc                                                                                                         >  2:11218/1‑150 (MQ=255)
    agccagGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGTTa                                                                                                     <  2:82491/150‑1 (MQ=255)
                       ggTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGTTACCGTCCAACTCCACTTCCa                                                                                  >  1:198929/1‑150 (MQ=255)
                         tCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGTTACCGTCCAACTCCACTTCCAgg                                                                                <  1:158133/150‑1 (MQ=255)
                                          cACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGTTACCGTCCAACTCCACTTCCAGGACCTGGTGCGGGTAGtt                                                               <  1:151250/150‑1 (MQ=255)
                                             gtaccacgcgccgtggccCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTg                                                                                                                                        >  1:598592‑M1/19‑74 (MQ=255)
                                             gtaccacgcgccgtggccCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTg                                                                                                                                        <  2:598592‑M1/56‑1 (MQ=255)
                                                 cacgcgccgtggccctgcaacgggtgggcccggcaccgctgggcgagctgttgcagGGCCACGGCGCGTGGTAc                                                                                                                                    >  1:598592‑M2/57‑74 (MQ=255)
                                                 cacgcgccgtggccctgcaacgggtgggcccggcaccgctgggcgagctgttgcagGGCCACGGCGCGTGGTAc                                                                                                                                    <  2:598592‑M2/18‑1 (MQ=255)
                                                      ggagaGTTGCTGCAACAGCACGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGTTACCGTCCAACTCCACTTCCAGGACCTGGTGCGGGGAGTTGATTTCGGCGaa                                                   >  2:127547/5‑150 (MQ=255)
                                                        ccAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGTTACCGTCCAACTCCACTTCCAGGACCTGGTGCGGGTAGTTGATTTCGGCGAACg                                                 >  2:158225/1‑150 (MQ=255)
                                                               ctgcaacgggtgggcccggcaccgctgggcgagctgttgcagGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGTTACCGTCCAACTCCACTTCCAGGACCTGGTGCggg                                                                    >  1:565794‑M2/43‑124 (MQ=255)
                                                               ctgcaacgggtgggcccggcaccgctgggcgagctgttgcagGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGTTACCGTCCAACTCCACTTCCAGGACCTGGTGCggg                                                                    <  2:565794‑M2/82‑1 (MQ=255)
                                                                                 cGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACACTAACCAGGCCGTTACCGTCCAACTCCACTTCCAGGACCTGGTGCGGGTAGTTGATTTCGGCGAACGACAGCGGGATTGGCGAGCCGCTGTa                        <  2:391646/150‑1 (MQ=255)
                                                                                          gCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGTTACCGTCCAACTCCACTTCCAGGACCTGGTGCGGGTAGTTGATTTCGGCGAACGACAGCGGGATTGGCGAGCCGCTGTAGCGGATgcg               >  1:139784/1‑150 (MQ=255)
                                                                                           cccACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGTTACCGTCCAACTCCACTTCCAGGAc                                                                              >  1:240326/1‑86 (MQ=255)
                                                                                           cccACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGTTACCGTCCAACTCCACTTCCAGGAc                                                                              <  2:240326/86‑1 (MQ=255)
                                                                                              aCCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGTTACCGTCCAACTCCACTTCCAGGACCTGGTGCGGGTAGTTGATTTCGGCGAACGACAGCGGGATTGGCGAGCCGCTGTAGCGGATGCGctct           >  2:294002/1‑150 (MQ=255)
                                                                                                       agGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGTTACCGTCCAACTCCACTTCCAGGACCTGGTGCGGGTAGTTGATTTCGGCGAACGACAGCGGGATTGGCGAGCCGCTGTAGCGGATGCGCTCTTCGCGGTTg  >  2:319665‑M2/3‑150 (MQ=255)
                                                                                                         ggCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGTTACCGTCCAACTCCa                                                                                        <  1:266961/62‑1 (MQ=255)
                                                                                                         ggCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGTTACCGTCCAACTCCa                                                                                        >  2:266961/1‑62 (MQ=255)
                                                                                                         |                                                                                                                                                   
CTCCAGCCAGGGCTGGCGGTTGGGGTCTTCAAGCAGGTCGACCACGGGCAGGTCGGCCAGTTGCTGCAACAGCTCGCCCAGCGGTGCCGGGCCCACCCGTTGCAGGGCCACGGCGCGTGGTACGGGGCGTGGCTCGACGCTGACCAGGCCGTTACCGTCCAACTCCACTTCCAGGACCTGGTGCGGGTAGTTGATTTCGGCGAACGACAGCGGGATTGGCGAGCCGCTGTAGCGGATGCGCTCTTCGCGGTTG  >  NC_002947/2300443‑2300695

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