New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 4322733 =10 (0.600)4 (0.250) 4/270 1.5 29.1% intergenic (‑3/‑415) PP_3791/PP_5629 site‑specific recombinase/protease/amidase
?NC_002947 4322760 = 10 (0.640)intergenic (‑30/‑388) PP_3791/PP_5629 site‑specific recombinase/protease/amidase

CTCCGCGATTTGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGCCCGACCAATTTTCCTTTGTTCAAAGGCTGATGGCTGCTGGTCGTAATGGCGTTCATGGCTTGTCCTCCATGTGAGGTAGGACAAGGATGATCAACCCAATCCAGAGGCGCTAACCGACCCATTGTTGTCCTTCGTGGAGAACAGTTATCGGCCAAACGCGGTCACTTACTGAAGACAAAGTCTGAGGAGCTCACTAAATGCCCCCGGCGGTCTGGAGTCAATTGCGT  >  NC_002947/4322611‑4322901
                                                                                                                                                     |                                                                                                                                             
ctccGCGATTTGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGCCCGACCAATTTTCCTTTGTTCAAAGGCTGATGGCTGCTGGTCGTAATGGCGTTCATGGCTTGTCCTCCATGTGAGGTAGGACAAgg                                                                                                                                               <  2:187560/150‑1 (MQ=255)
 tccGCGATTTGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGCCCGACCAATTTTCCTTTGTTCAAAGGCTGATGGCTGCTGGTCGTAATGGCGTTCATGGCTTGTCCTCCATGTGAGGTAGGACAAGga                                                                                                                                              >  1:263265/1‑150 (MQ=255)
                    cccGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGCCCGACCAATTTTCCTTTGTTCAAAGGCTGATGGCTGCTGGTCGTAATGGCGTTCATGGCTTGTCCTCCATGTGAGGTAGGACAAGGATGATCAACCCAATCCagag                                                                                                                           >  2:122674/1‑150 (MQ=255)
                                   tatCTTTCAGTCGGAGCGGAGCTTTCTGCCCGACCAATTTTCCTTTGTTCAAAGGCTGATGGCTGCTGGTCGTAATGGCGTTCATGGCTTGTCCTCCATGTGAGGTAGGACAAGGATGATCAACCCAATCCAGAGGCGCTAACCGACCCa                                                                                                            <  1:233807/150‑1 (MQ=255)
                                       tttCAGTCGGAGCGGAGCTTTCTGCCCGACCAATTTTCCTTTGTTCAAAGGCTGATGGCTGCTGGTCGTAATGGCGTTCATGGCTTGTCCTCCATGTGAGGTAGGACAAGGATGATCAACCCAATCCAGAGGCGCTAACCGACCCAttgt                                                                                                        >  2:314619/1‑150 (MQ=255)
                                           tgaccgcgtttggccgataactgttctccacgaaggacaacaatgggtcggttagcgcctctggattgggttgatcatcCTTGTCCTCCATGTGAGGTAGGACAAGgatg                                                                                                                                            >  1:269055‑M1/80‑110 (MQ=255)
                                           tgaccgcgtttggccgataactgttctccacgaaggacaacaatgggtcggttagcgcctctggattgggttgatcatcCTTGTCCTCCATGTGAGGTAGGACAAGgatg                                                                                                                                            <  2:269055‑M1/31‑1 (MQ=255)
                                                                   aCCAATTTTCCTTTGTTCAAAGGCTGATGGCTGCTGGTCGTAATGGCGTTCATGGCTTGTCCTCCATGTGAGGTAGGACAAGGATGATCAACCCAATCCAGAGGCGCTAACCGACCCATTGTTGTCCTTCGTGGAGAACAGTTATCGGcc                                                                            >  1:241039/1‑150 (MQ=255)
                                                                               ttGTTCAAAGGCTGATGGCTGCTGGTCGTAATGGCGTTCATGGCTTGTCCTCCATGTGAGGTAGGACAAGGATGATCAACCCAATCCAGAGGCGCTAACCGACCCATTGTTGTCCTTCGTGGAGAACAGTTATCGGCCAAACGCGGTCAc                                                                <  2:228999/150‑1 (MQ=255)
                                                                                   tCAAAGGCTGATGGCTGCTGGTCGTAATGGCGTTCATGGCTTGTCCTCCATGTGAGGTAGGACAAGGATGATCAACCCAATCCAGAGGCGCTAACCGACCCATTGTTGTCCTTCGTGGAGAACAGTTATCGGCCAAGCGCGGTCACTTAc                                                            >  1:193019/1‑150 (MQ=255)
                                                                                                                      catccttgtcctacctcacatggaggacaagGATGATCAACCCAATCCAGAGGCGCTAACCGACCCATTGTTGTCCTTCGTGGAGAACAGTTATCGGCCAAACGCGGTCa                                                                 >  1:269055‑M2/32‑110 (MQ=255)
                                                                                                                      catccttgtcctacctcacatggaggacaagGATGATCAACCCAATCCAGAGGCGCTAACCGACCCATTGTTGTCCTTCGTGGAGAACAGTTATCGGCCAAACGCGGTCa                                                                 <  2:269055‑M2/79‑1 (MQ=255)
                                                                                                                           ttgtcctacctcacatggaggacaagGATGATCAACCCAATCCAGAGGCGCTAACCGACCCATTGTTGTCCTTCGTGGAGAACAGTTATCGGCCAAACGCGGTCACTTACTGAAGACAAAGTCTGAGGAGCTCACTAAATGCCCCcggcg                    >  2:310815‑M2/27‑150 (MQ=255)
                                                                                                                                CTCCATGTGAGGTAGGACAAGGATGATCAACCCAATCCAGAGGCGCTAACCGACCCATTGTTGTCCTTCGTGGAGAACAGTTATCGGCCAAACGCGGTCACTTACTGAAGACaaa                                                  <  1:69709‑M1/115‑1 (MQ=255)
                                                                                                                                CTCCATGTGAGGTAGGACAAGGATGATCAACCCAATCCAGAGGCGCTAACCGACCCATTGTTGTCCTTCGTGGAGAACAGTTATCGGCCAAACGCGGTCACTTACTGAAGACaaa                                                  >  2:69709‑M1/1‑115 (MQ=255)
                                                                                                                                CTCCATGTGAGGTAGGACAAGGATGATCAACCCAATCCAGAGGCGCTAACCGACCCATTGTTGTCCTTCGTGGAGAACAGTTATCGGCCAAACGCGGTCACTTACTGAAGACAAAGTCTGAGGAGCTCACTAAATGCCCCCGGCGGTCTg               <  1:78611‑M1/150‑1 (MQ=255)
                                                                                                                                CTCCATGTGAGGTAGGACAAGGATGATCAACCCAATCCAGAGGCGCTAACCGACCCATTGTTGTCCTTCGTGGAGAACAGTTATCGGCCAAACGCGGTCACTTACTGAAGACAAAGTCTGAGGAGCTCACTAAATGCCCCCGGCGGTCTg               >  1:314619‑M1/1‑150 (MQ=255)
                                                                                                                                CTCCATGTGAGGTAGGACAAGGATGATCAACCCAATCCAGAGGCGCTAACCGACCCATTGTTGTCCTTCGTGGAGAACAGTTATCGGCCAAACGCGGTCACTTACTGAAGACAAAGTCTGAGGAGCTCACTAAATGCCCCCGGCGGTCTg               >  2:263265‑M1/1‑150 (MQ=255)
                                                                                                                                   cctcacatggaggacaagGATGATCAACCCAATCCAGAGGCGCTAACCGACCCATTGTTGTCCTTCGTGGAGAACAGTTATCGGCCAAACGCGGTCACTTACTGAAGACAAAGTCTGAGGAGCTCACTAAATGCCCCCGGCGGTCTGGAg            <  1:310815‑M2/132‑1 (MQ=255)
                                                                                                                                   cATGTGAGGTAGGACAAGGATGATCAACCCAATCCAGAGGCGCTAACCGACCCATTGTTGTCCTTCGTGGAGAACAGTTATCGGCCAAGCGCGGTCACTTACTGAAGACAAAGTCTGAGGAGCTCACTAAATGCCCCCGGCGGTCTGGAg            <  2:193019/150‑1 (MQ=255)
                                                                                                                                     tgtgAGGTAGGACAAGGATGATCAACCCAATCCAGAGGCGCTAACCGACCCATTGTTGTCCTTCGTGGAGAACAGTTATCGGCCAAACGCGGTCACTTACTGAAGACAAAGTCTGAGGAGCTCACTAAATGCCCCCGGCGGTCTGGAGTc          >  2:193820/1‑150 (MQ=255)
                                                                                                                                             aggacaagGATGATCAACCCAATCCAGAGGCGCTAACCGACCCATTGTTGTCCTTCGTGGAGAACAGTTATCGGCCAAACGCGGTCACTTACTGAAGACAAAGTCTGAGGAGCTCACTAAATGCCCCCGGCGGTCTGGAGTCAATTGCgt  <  2:359776‑M2/142‑1 (MQ=255)
                                                                                                                                                    gGATGATCAACCCAATCCAGAGGCGCTAACCGACCCATTGTTGTCCTTCGTGGAGAACAGTTATCGGCCAAACGCGGTCACTTACTGAAGACAAAGTCTGAGGAGCTCACTAAATGCCCCcggc                     >  2:10987‑M2/2‑124 (MQ=255)
                                                                                                                                                    gGATGATCAACCCAATCCAGAGGCGCTAACCGACCCATTGTTGTCCTTCGTGGAGAACAGTTATCGGCCAAACGCGGTCACTTACTGAAGACAAAGTCTGAGGAGCTCACTAAATGCCCCcggc                     <  1:10987‑M2/123‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
CTCCGCGATTTGGAGCCTTACCCGAATGGCCCAGATATCTTTCAGTCGGAGCGGAGCTTTCTGCCCGACCAATTTTCCTTTGTTCAAAGGCTGATGGCTGCTGGTCGTAATGGCGTTCATGGCTTGTCCTCCATGTGAGGTAGGACAAGGATGATCAACCCAATCCAGAGGCGCTAACCGACCCATTGTTGTCCTTCGTGGAGAACAGTTATCGGCCAAACGCGGTCACTTACTGAAGACAAAGTCTGAGGAGCTCACTAAATGCCCCCGGCGGTCTGGAGTCAATTGCGT  >  NC_002947/4322611‑4322901

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.