New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 4647267 =18 (1.080)4 (0.240) 4/280 1.6 19.7% coding (420/483 nt) PP_5658 hypothetical protein
?NC_002947 4647285 = 15 (0.920)coding (438/483 nt) PP_5658 hypothetical protein
Rejected: Frequency below/above cutoff threshold.

CCCAGGTGTTGAAAAAAGCCTGGCTGTCGTACACCAAGGCTGATCTCAACCTGGTCGGTGCGCTGGACCTGACCATCCTCAAGACCTTCACCGAGGCCGCCCTGCACACGCGCCTGCTGGCCTTGCTGCCGGCCAGCGGCCTGGCCCGTGAAGTGGATGGCCAGTGGATGCTGCAAGCCCAGACCGACTGACGGCCCGGTGGCAGGGCAGGATTCCACCGGTCTGCTAACTTGGTGGAATCCGTCCACCCACAGGTAGTCCGCCATGGCCCGCACAACGCCCATCGAGCTGTACC  >  NC_002947/4647140‑4647434
                                                                                                                                                 |                                                                                                                                                     
cccAGGTGTTGAAAAAAGCCTGGCTGTCGTACACCAAGGCTGATCTCAACCTGGTCGGTGCGCTGGACCTGACCATCCTCAAGACCTTCACCGAGGCCGCCCTGCACACGCGCCTGCTGGCCTTGCTGCCGGCCAGCGGCCTGGCCCGTg                                                                                                                                                   >  2:216666/1‑150 (MQ=255)
                aGCCTGGCTGTCGTACACCAAGGCTGATCTCAACCTGGTCGGTGCGCTGGACCTGACCATCCTCAAGACCTTCACCGAGGCCGCCCTGTACACGCGCCTGCTGGCCTTGCTGCCGGCCAGCGGCCTGGCCCGTGAAGTGGATGGCCAGTg                                                                                                                                   <  1:286701/150‑1 (MQ=255)
                    gtgggtggacggattccaccaagttagcagaccggtggaatcctgccctgccaccgggccgtcagtcggtctgggcttgcagcatccactggccatccacttcacgggccaggccgctggccggcCCGTGa                                                                                                                                                  >  2:257440‑M2/126‑131 (MQ=255)
                    gtgggtggacggattccaccaagttagcagaccggtggaatcctgccctgccaccgggccgtcagtcggtctgggcttgcagcatccactggccatccacttcacgggccaggccgctggccggcCCGTGa                                                                                                                                                  <  1:257440‑M2/6‑1 (MQ=255)
                                 ccAAGGCTGATCTCAACCTGGTCGGTGCGCTGGACCTGACCATCCTCAAGACCTTCACCGAGGCCGCCCTGCACACGCGCCTGCTGGCCTTGCTGCCGGCCAGCGGCCTGGCCCGTGAAGTGGATGGCCAGTGGATGCTGCAAGCCCAGa                                                                                                                  <  2:282276/150‑1 (MQ=255)
                                  cAAGGCTGATCTCAACCTGGTCGGTGCGCTGGACCTGACCATCCTCAAGACCTTCACCGAGGCCGCCCTGCACACGCGCCTGCTGGCCTTGCTGCCGGCCAGCGGCCTGGCCCgca                                                                                                                                                   >  2:250168/1‑114 (MQ=255)
                                  cAAGGCTGATCTCAACCTGGTCGGTGCGCTGGACCTGACCATCCTCAAGACCTTCACCGAGGCCGCCCTGCACACGCGCCTGCTGGCCTTGCTGCCGGCCAGCGGCCTGGCCCgca                                                                                                                                                   <  1:250168/116‑3 (MQ=255)
                                                           gcgcTGGACCTGACCATCCTCCAGACCTTCACCGAGGCCGCCCTGCACACGCGCCTGCTGGCCTTGCTGGCGGCCAGCGGCCTGGCCCGTGAAGTGGATGGCCAGTGGCTGCTGCAAGCCCAGACCGACTGACGGCCCGGTggcagggca                                                                                        <  1:314182/150‑1 (MQ=255)
                                                            cgcTGGACCTGACCATCCTCAAGACCTTCACCGAGGCCGCCCTGCACACGCGCCTGCTGGCCTTGCTGCCGGCCAGCGGCCTGGCCCGTGAAGTGGATGGCCAGTGGATGCTGCAAGCCCAGACCGACTGACGGCCCGGTggcagggcag                                                                                       <  2:243312/150‑1 (MQ=255)
                                                                 ccctgccaccgggccgtcagtcggtctgggcttgcagcatccactggccatccacttcacgggccggccagcggcctggcCCGTGAAGt                                                                                                                                               >  1:250555‑M2/81‑89 (MQ=255)
                                                                 ccctgccaccgggccgtcagtcggtctgggcttgcagcatccactggccatccacttcacgggccggccagcggcctggcCCGTGAAGt                                                                                                                                               <  2:250555‑M2/9‑1 (MQ=255)
                                                                          aTCCTCAAGACCTTCACCGAGGCCGCCCTGCACACGCGCCTGCTGGCCTTGCTGCCGGCCAGCGGCCTGGCCCGTGAAGTGGATGGCCAGTGGATGCTGCAAGCCCAGACCGACTGACGGCCCGGTGGCAGGGCAGGATTCCACCGGTCt                                                                         <  1:202973/150‑1 (MQ=255)
                                                                                      ttCACCGAGGCCGCCCTGCACACGCGCCTGCTGGCCTTGCTGCCGGCCAGCGGCCTGGCCCGTGAAGTGGATGGCCAGTGGATGCTGCAAGCCCAGACCGACTGACGGCCCGGTGGCAGGGCAGGATTCCACCGGTCTGCTAACTtggtg                                                             <  1:340259/150‑1 (MQ=255)
                                                                                      ttCACCGAGGCCGCCCTGCACACGCGCCTGCTGGCCTTGCTGCCGGCCAGCGGCCTGGCCCGTGAAGTGGATGGCCAGTGGATGCTGCAAGCCCAGACCGACTGACGGCCCGGTGGCAGGGCAGGATTCCACCGGTCTGCTAACTtggtg                                                             >  1:352697/1‑150 (MQ=255)
                                                                                      ttCACCGAGGCCGCCCTGCACACGCGCCTGCTGGCCTTGCTGCCGGCCAGCGGCCTGGCCCGTGAAGTGGATGGCCAGTGGATGCTGCAAGCCCAGACCGACTGACGGCCCGGTGGCAGGGCAGGATTCCACCGGTCTGCTAACTtggtg                                                             <  1:384595/150‑1 (MQ=255)
                                                                                              ggCCGCCCTGCACACGCGCCTGCTGGCCTTGCTGCCGGCCAGCGGCCTGGCCCGTGAAGTGGATGGCCAGTGGATGCTGCAAGCCCAGACCGACTGACGGCCCGGTGGCAGGGCAGGATTCCACCGGTCTGCTAACTTGGTGGAATCCGt                                                     <  1:146896/150‑1 (MQ=255)
                                                                                                 cgccCTGCACACGCGCCTGCTGGCCTTGCTGCCGGCCAGCGGCCTGGCCCGTGAAGTGGATGGCCAGTGGATGCTGCAAGCCCAGACCGACTGACGGCCCGGTGGCAGGGCAGGATTCCACCGGTCTGCTAACTTGGTGGAATCCGTcca                                                  <  2:89483/150‑1 (MQ=255)
                                                                                                     cTGCACACGCGCCTGCTGGCCTTGCTGCCGGCCAGCGGCCTGGCCCGTGAAGTGGATGGCCAGTGGATGCTGCAAGCCCAGACCGACTGACGGCCCGGTGGCAGGGCAGGATTCCACCGGTCTGCTAACTTGGTGGAATCCGTccaccca                                              >  1:99566/1‑150 (MQ=255)
                                                                                                           acGCGCCTGCTGGCCTTGCTGCCGGCCAGCGGCCTGGCCCGTGAAGTGGATGGCCAGTGGATGCTGCAAGCCCAGACCGACTGACGGCCCGGTGGCAGGGCAGGATTCCACCGGTCTGCTAACTTGGTGGAATCCGTCCACCCACAGGTa                                        <  2:352697/150‑1 (MQ=255)
                                                                                                                      acttcacggGCCAGGCCGCTGGCCGGCCCGTGAAGTGGATGGCCAGTGGATGCTGCAAGCCCAGACCGACTGACGGCCCGGTggcaggg                                                                                          <  2:250555‑M1/80‑1 (MQ=255)
                                                                                                                      acttcacggGCCAGGCCGCTGGCCGGCCCGTGAAGTGGATGGCCAGTGGATGCTGCAAGCCCAGACCGACTGACGGCCCGGTggcaggg                                                                                          >  1:250555‑M1/10‑89 (MQ=255)
                                                                                                                         tcacggGCCGGCCAGCGGCCTGGCCCGTGAAGTGGATGGCCAGTGGATGCTGCAAGCCCAGACCGACTGACGGCCCGGTGGCAGGGCAGGATTCCACCGGTCTGCTAACTTGGTGGAATCCGTccacccac                                             >  2:257440‑M1/7‑131 (MQ=255)
                                                                                                                         tcacggGCCGGCCAGCGGCCTGGCCCGTGAAGTGGATGGCCAGTGGATGCTGCAAGCCCAGACCGACTGACGGCCCGGTGGCAGGGCAGGATTCCACCGGTCTGCTAACTTGGTGGAATCCGTccacccac                                             <  1:257440‑M1/125‑1 (MQ=255)
                                                                                                                                                 ccGTGAAGTGGATGGCCAGTGGATGCTGCAAGCCCAGACCGACTGACGGCCCGGTGGCAGGGCAGGATTCCACCGGTCTGCTAACTTGGTGGAATCCGTCCACCCACAGGTa                                        <  1:51734/112‑1 (MQ=255)
                                                                                                                                                 ccGTGAAGTGGATGGCCAGTGGATGCTGCAAGCCCAGACCGACTGACGGCCCGGTGGCAGGGCAGGATTCCACCGGTCTGCTAACTTGGTGGAATCCGTCCACCCACAGGTa                                        >  2:51734/1‑112 (MQ=255)
                                                                                                                                                 ccGTGAAGTGGATGGCCAGTGGATGCTGCAAGCCCAGACCGACTGACGGCCCGGTGGCAGGGCAGGATTCCACCGGTCTGCTAACTTGGTGGAATCCGTCCACCCACAGGTAGTCCGCCATGGCCCGCACAACGCCCATCGAGCTGTAcc  <  1:216666/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CCCAGGTGTTGAAAAAAGCCTGGCTGTCGTACACCAAGGCTGATCTCAACCTGGTCGGTGCGCTGGACCTGACCATCCTCAAGACCTTCACCGAGGCCGCCCTGCACACGCGCCTGCTGGCCTTGCTGCCGGCCAGCGGCCTGGCCCGTGAAGTGGATGGCCAGTGGATGCTGCAAGCCCAGACCGACTGACGGCCCGGTGGCAGGGCAGGATTCCACCGGTCTGCTAACTTGGTGGAATCCGTCCACCCACAGGTAGTCCGCCATGGCCCGCACAACGCCCATCGAGCTGTACC  >  NC_002947/4647140‑4647434

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.