New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 5868966 =67 (1.290)3 (0.060) 3/258 NT 4.6% coding (397/783 nt) PP_5143 membrane protein
?NC_002947 5869013 = 59 (1.170)coding (350/783 nt) PP_5143 membrane protein
Rejected: Frequency below/above cutoff threshold.

GATGACCATCGGCATCCTCCTGGGGGCCGGTGTCGGTGCCAAGACCGCTTCGCTGATCCAGGGGCCACTGCTGCAAAAGATCATCGGTGTGTTCGCCCTGGTCATCGCCGCGCAGATGGCCCTTGACCTCAAACCCAAGGCCAGCCGTGGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/5869117‑5868966
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcgATGACCAGGGCGAACACACCGATGATCTTTTGCAGCAGTGGCCCCTGGATCAGCGAAGCGGTCTTGGCACCGACACCGGCCCCCAGGAGGATGCCGATGGT  >  NC_002947/5869013‑5869113
                                                                                                                                                                                                                                                             
GATGACCATCGGCATCCTCCTGGGGGCCGGTGTCGGTGCCAAGACCGCTTCGCTGATCCAGCGGCCACTGCTGCAAAAGATCATCGGTGTGTTCCCCCTGGTCATCGCCGCGCAGATCGCCCTTGACCTCAAACCCAAGGCCCCCCGTG                                                                                                          >  7:162982/1‑149
                                                                                                      CATCGCCGCGCAGATGGCCCTTGACCTCAAACCCAAGGCCAGCCGTGGCGATGACCAGGGCGAACACACCGATGATCTTTTGCAGCAGTGGCCCCTGGATCAGCGAAGCGGTCTTGGCACCGACACCGGCCCCCAGGAGGATGCCGAT     <  5:317183/148‑1
                                                                                                      CATCGCCGCGCAGATGGCCCTTGACCTCAAACCCAAGGCCAGCCGTGGCGATGACCAGGGCGAACACACCGATGATCTTTTGCAGCAGTGGCCCCTGGATCAGCGAAGCGGTCTTGGCACCGACACCGGCCCCCAGGAGGATGCCGAT     >  6:317183/1‑148
                                                                                                        TCGCCGCGCAGATGGCCCTTGACCTCAAACCCAAGGCCAGCCGTGGCGATGACCAGGGCGAACACACCGATGATCTTTTGCAGCAGTGGCCCCTGGATCAGCGAAGCGGTCTTGGCACCGACACCGGCCCCCAGGAGGATGCCGATGGT  >  6:89557/1‑149
                                                                                                                                                                                                                                                             
GATGACCATCGGCATCCTCCTGGGGGCCGGTGTCGGTGCCAAGACCGCTTCGCTGATCCAGGGGCCACTGCTGCAAAAGATCATCGGTGTGTTCGCCCTGGTCATCGCCGCGCAGATGGCCCTTGACCTCAAACCCAAGGCCAGCCGTGGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_002947/5869117‑5868966
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcgATGACCAGGGCGAACACACCGATGATCTTTTGCAGCAGTGGCCCCTGGATCAGCGAAGCGGTCTTGGCACCGACACCGGCCCCCAGGAGGATGCCGATGGT  >  NC_002947/5869013‑5869113

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