New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 521524219 (0.700)4 (0.160) 4/274 2.8 17.1% intergenic (‑111/+67) rnd/PP_4592 ribonuclease D/metal dependent hydrolase
?NC_002947 = 5215279 21 (0.830)intergenic (‑148/+30) rnd/PP_4592 ribonuclease D/metal dependent hydrolase
Rejected: Frequency below/above cutoff threshold.

TCTGGTCGTCACGGATCCAGTGAATTTCGATGGCCACAAGGCTCTCCCACTAACATTGGCGCGCAGTATATACGGCGCGGGCACTGCGGGTGAAACTGTCCCGCTTTTCGGGGGCAGTCCACCCTGTGGGAGCGGGCGCGCCCGCGAATGAAGGCAACACGGTGTATGGCCCGGGTTTGCCCGGGTTCGCGGGCGTGCCCGCTCCCACAGGGCTTCAGCGCTTGGCCGCCAGCCAAGGCCGGCAACCGGCAAACATGTCCAGCGACTGCTGGTAGAGCTCGGTATGTACCTGC  >  NC_002947/5215095‑5215387
                                                                                                                                                   |                                                                                                                                                 
tctGGTCGTCACGGATCCAGTGAATTTCGATGGCCACAAGGCTCTCCCACTAACATTGGCGCGCAGTATATACGGCGCGGGCACTGCGGGTGAAACTGTCCCGCTTTTCGGGGGCAGTCCACCCTGTGGGAGCGGGCGCGCCCGCGAATg                                                                                                                                                 >  2:80923/1‑150 (MQ=255)
tctGGTCGTCACGGATCCAGTGAATTTCGATGGCCACAAGGCTCTCCCACTAACATTGGCGCGCAGTATATACGGCGCGGGCACTGCGGGTGAAACTGTCCCGCTTTTCGGGGGCAGTCCACCCTGTGGGAGCGGGCGCGCCCGCGAATg                                                                                                                                                 >  1:324146/1‑150 (MQ=255)
          aCGGATCCAGTGAATTTCGATGGCCACAAGGCTCTCCCACTAACATTGGCGCGCAGTATATACGGCGCGGGCACTGCGGGTGAAACTGTCCCGCTTTTCGGGGGCAGTCCACCCTGTGGGAGCGGGCGCGCCCGCGAATGAAGGCAacac                                                                                                                                       >  2:506314/1‑150 (MQ=255)
                       aTTTCGATGGCCACAAGGCTCTCCCACTAACATTGGCGCGCAGTATATACGGCGCGGGCACTGCGGGTGAAACTGTCCCGCTTTTCGGGGGCAGTCCACCCTGTGGGAGCGGGCGCGCCCGCGAAcccgggcaaacccgggccatacacc                                                                                                                          >  1:53327‑M1/1‑125 (MQ=255)
                         ttCGATGGCCACAAGGCTCTCCCACTAACATTGGCGCGCAGTATATACGGCGCGGGCACTGCGGGTGAAACTGTCCCGCTTTTCGGGGGCAGTCCACCCTGTGGGAGCGGGCGCGCCCGCGACTGAAGGCAACCCGGTGTATGGCCCggg                                                                                                                        >  2:214339/1‑150 (MQ=255)
                                  cacaAGACTCTCCCACTAACATTGGCGCGCAGTATATACGGCGCGGGCACTGCGGGTGAAACTGTCCCGCTTTTCGGGGGCAGTCCACCCTGTGGGAGCGGGCGCGCCCGCGAAcccgggcaaacccgggccatacaccgtgttgccttc                                                                                                               <  1:572410‑M1/150‑37 (MQ=255)
                                   acaAGGCTCTCCCACTAACATTGGCGCGCAGTATATACGGCGCGGGCACTGCGGGTGAAACTGTCCCGCTTTTCGGGGGCAGTCCACCCTGTGGGAGCGGGCGCGCCCGCGAATGAAGGCAACACGGTGTATGGCCCGGGTTTGCCCggg                                                                                                              <  2:91328/150‑1 (MQ=255)
                                   acaAGGCTCTCCCACTAACATTGGCGCGCAGTATATACGGCGCGGGCACTGCGGGTGAAACTGTCCCGCTTTTCGGGGGCAGTCCACCCTGTGGGAGCGGGCGCGCCCGCGAATGAAGGCAACACGGTGTATGGCCCGGGTTTGCCCggg                                                                                                              >  2:108323/1‑150 (MQ=255)
                                        gCTCTCCCACTAACATTGGCGCGCAGTATATACGGCGCGGGCACTGCGGGTGAAACTGTCCCGCTTTTCGGGGGCAGTCCACCCTGTGGGAGCGGGCGCGCCCGCGAAcccgggcaaacccgggccatacaccgggttgccttcattcgc                                                                                                         <  2:468994‑M1/150‑43 (MQ=255)
                                         ctctCCCACTAACATTGGCGCGCAGTATATACGGCGCGGGCACTGCGGGTGAAACTGTCCCGCTTTTCGGGGGCAGTCCACCCTGTGGGAGCGGGCGCGCCCGCGAAcccgggcaaacccgggccatacaccgtgttgccttcattcgcg                                                                                                        <  2:53327‑M1/150‑44 (MQ=255)
                                                  tAACATTGGCGCGCAGTATATACGGCGCGGGCACTGCGGGTGAAACTGTCCCGCTTTTCGGGGGCAGTCCACCCTGTGGGAGCGGGCGCGCCCGCGAATGAAGGCAACACGGTGTATGGCCCGGGTTTGCCCGGGTTCGCGGGCGTGccc                                                                                               >  2:434275/1‑150 (MQ=255)
                                                          gcgcgcAGTATATACGGCGCGGGCACTGCGGGTGAAACTGTCCCGCTTTTCGGGGGCAGTCCACCCTGTGGGAGCGGGCGCGCCCGCGAATGAAGGCAACACGGTGTATGGCCCGGGTTTGCCCGGGTTCGCGGGCGTGCCCGCTCCcac                                                                                       >  2:170707/1‑150 (MQ=255)
                                                                    tataCTGCGCGGGCACTGCGGGTGAAACTGTCCCGCTTTTCGGGGGCAGTCCACCCTGTGGGAGCGGGCGCGCCCGCGAATGAAGGCAACACGGTGTATGGCCCGGGGTTGCCCGGGTTCGCGGGCGTGCCCGCTCCCACAGGGCTTCAg                                                                             >  2:258138/1‑150 (MQ=255)
                                                                                gCACTGCGGGTGAAACTGTCCCGCTTTTCGGGGGCAGTCCACCCTGTGGGAGCGGGCGCGCCCGCGAATGAAGGCAACACGGTGTATTGCCCGGGTTTGCCCGGGTTCGCGGGCGTGCCCCCTCCCACACGGCTTCAGCGCGTGgccgcc                                                                 >  1:36541/1‑150 (MQ=255)
                                                                                                            cGGGGGCAGTCCACCCTGTGGGAGCGGGCGCGCCCGCGAATGAAGGCAACACGGGGTATGGCCCGGGTATGCCCGGGTTCGCGGGCGTGCCCGCTCCCACAGGGCTTCAGCGCTTGGCCGCCAGCCAAGGCCGGCAACCGGCCAACATGt                                     >  2:525145/1‑150 (MQ=255)
                                                                                                                          ccTGTGGGAGCGGGCGCGCCCGCGAATGAAGGCAACACGGTGTATGGCCCGGGTTTGCCCGGGTTCGCGGGCGTGCCCGCTCCCACAGGGCTTCAGCGCTTGGCCGCCAGCCAAGGc                                                        <  2:487842/117‑1 (MQ=255)
                                                                                                                          ccTGTGGGAGCGGGCGCGCCCGCGAATGAAGGCAACACGGTGTATGGCCCGGGTTTGCCCGGGTTCGCGGGCGTGCCCGCTCCCACAGGGCTTCAGCGCTTGGCCGCCAGCCAAGGc                                                        >  1:487842/1‑117 (MQ=255)
                                                                                                                            tgtgGGAGCGGGCGCGCCCGCGAATGAAGGCAACACGGTGTATGGCCCGGGTTTGCCCGGGTTCGCGGGCGTGCCCGCTCCCACAGGGCTTCAGCGCTTGGCCGCCAGCCAAGGCCGGCAACCGGCAAACATGTCCAGCGACTGCTGGTa                     <  2:235352/150‑1 (MQ=255)
                                                                                                                              tgGGAGCGGGCGCGCCCGCGAATGACGGCAACACGGTGTATGGCCCGGGTTTGCCCGGGTTCGCGGGCGTGCCCGCTCCCACAGGGCTTCAGCGCTTGGCCGCCAGCCAAGGCCGGCAACCGGCAAACATGTCCAGCGACTGCTGGTaga                   <  2:162601/150‑1 (MQ=255)
                                                                                                                              tgGGAGCGGGCGCGCCCGCGAATGAAGGCAACACGGTGTATGGCCCGGGTTTGCCCGGGTTCGCGGGCGTGCCCGCTCCCACAGGGCTTCAGCGCTTGGCCGCCAGCCAAGGCCGGCAACCGGCAAACATGTCCAGCGACTGCTGGTaga                   >  2:473176/1‑150 (MQ=255)
                                                                                                                                  tgggggCGCGCCCGCGAATGAAGGCCCCCCGGTGTTTGGCCGGGGTTGGCCGGGGGGCGGGGGCGGGCCCGCTCCCACAGGGCTTCAGCGCTTGGCCGCCAGCCAAGGCCGGCAACCGGCAAACATGTCCAGCGACTGCTGGTAGAGCTc               <  1:89565/147‑1 (MQ=255)
                                                                                                                                   gCGGGCGCGCCCGCGAATGCCGGCAACCCGGTGTATGGCCCGGGTTTGCCCGGGTACGCGGGCGTGCCCGCTCCCACAGGGCTTCAGCGCTTGGCCGCCAGCCAAGGCCGGCAACCGGCAAACATGTCCAGCGACTGCTGGTAGAGCTCg              <  2:566340/150‑1 (MQ=255)
                                                                                                                                          cgcCCGCTATTGACGGAACCAGGGTGTATGGCCGGTGTTGGCCGGGGTTGGGGGGCGTGCCCGCTCCCACAGGGCTTCAGCGCTTGGCCGCCAGCCAAGGCCGGCAACCGGCAAACATGTCCAGCGACTGCTGGTAGAGCTCGGTATGTa       <  1:60441/150‑1 (MQ=255)
                                                                                                                                              cgcgAATGAAGGCAACACGGTGTATGGCCCGGGTTTGCCCGGGttcgcgggcgcgcccgctcccacagggtggactgcccccgaaaagcgggacagtttcacccgcagtgcccgcgccgtatatactgcgcgccaatgttagtgggagag   <  2:53327‑M2/150‑108 (MQ=255)
                                                                                                                                               gcgAATGAAGGCAACCCGGTGTATGGCCCGGGTTTGCCCGGGttcgcgggcgcgcccgctcccacagggtggactgcccccgaaaagcgggacagtttcacccgcagtgcccgcgccgtatatactgcgcgccaatgttagtgggagagc  <  2:468994‑M2/150‑109 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TCTGGTCGTCACGGATCCAGTGAATTTCGATGGCCACAAGGCTCTCCCACTAACATTGGCGCGCAGTATATACGGCGCGGGCACTGCGGGTGAAACTGTCCCGCTTTTCGGGGGCAGTCCACCCTGTGGGAGCGGGCGCGCCCGCGAATGAAGGCAACACGGTGTATGGCCCGGGTTTGCCCGGGTTCGCGGGCGTGCCCGCTCCCACAGGGCTTCAGCGCTTGGCCGCCAGCCAAGGCCGGCAACCGGCAAACATGTCCAGCGACTGCTGGTAGAGCTCGGTATGTACCTGC  >  NC_002947/5215095‑5215387

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