New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 520522146 (1.060)5 (0.120) 5/246 NT 10.8% intergenic (+64/+51) PP_4585/PP_4586 TrmH family RNA methyltransferase/hypothetical protein
?NC_002947 = 5205267 38 (0.910)intergenic (+110/+5) PP_4585/PP_4586 TrmH family RNA methyltransferase/hypothetical protein

CATGAACCTGGCGGCGACGGTGAATGTGGTGTTGTATGACCGCATGGCCAAAGGGCTGAATACCCGCTCAGGGCCCAAGTTCAAATAACAAGGCCGTGCGGCGAAGCTTGGGCTGTCAAAGGCCGGGACAGGCGCCAGGGATGCCTTGGGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/5205070‑5205221
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggccCGTATGAAAAGCCTTGCCTGACCTTCGTCTCAAGGCTTTTCATACAGGCCCAA  <  NC_002947/5205267‑5205215
                                                                                                                                                                                                             
CATGAACCTGGCGGCGACGGTGAATGTGGTGTTGTATGACCGCATGGCCAAAGGGCTGAATACCCGCTCAGGGCCCAAGTTCAAATAACAAGGCCGTGCGGCGAAGCTTGGGCTGTCAAAGGCCGGGACAGGCGCCAGGGATGCCTTGGG                                                         >  2:489955/1‑150
                                                  AAGGGCTGAATACCCGCTCAGGGCCCAAGTTCAAATAACAAGGCCGTGCGGCGAAGCTTGGGCTGTCAAAGGCCGGGACAGGCGCCAGGGATGCCTTGGGCCCGTATGAAAAGCCTTGCCTGACCTTCGTCTCAAGGCTTTTCATACAGG       <  1:489955/150‑1
                                                               CCGCTCAGGGACCAAGTTCAAATAACAAGGCCGTGCGGCGAAGCTTGGGCTGTCAAAGGCCGGGACAGGCGCCAGGGATGCCTTGGGCCTGTATGAAAAGCCTTGCCTGACCTTCGTCTCAAGGCTT                 >  1:389434/1‑127
                                                               CCGCTCAGGGACCAAGTTCAAATAACAAGGCCGTGCGGCGAAGCTTGGGCTGTCAAAGGCCGGGACAGGCGCCAGGGATGCCTTGGGCCTGTATGAAAAGCCTTGCCTGACCTTCGTCTCAAGGCTT                 <  2:389434/127‑1
                                                                                  AAATAACAAGGCCGTGCGGCGAAGCTTGGGCTGTCAAAGGCCGGGACAGGCGCCAGGGATGCCTTGGGCCCGTATGAAAAGCCTTGCCTGACCTTCGTCTCAAGGCTTTTCATACAGGCCCAA  >  1:434036/1‑123
                                                                                  AAATAACAAGGCCGTGCGGCGAAGCTTGGGCTGTCAAAGGCCGGGACAGGCGCCAGGGATGCCTTGGGCCCGTATGAAAAGCCTTGCCTGACCTTCGTCTCAAGGCTTTTCATACAGGCCCAA  <  2:434036/123‑1
                                                                                                                                                                                                             
CATGAACCTGGCGGCGACGGTGAATGTGGTGTTGTATGACCGCATGGCCAAAGGGCTGAATACCCGCTCAGGGCCCAAGTTCAAATAACAAGGCCGTGCGGCGAAGCTTGGGCTGTCAAAGGCCGGGACAGGCGCCAGGGATGCCTTGGGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/5205070‑5205221
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggccCGTATGAAAAGCCTTGCCTGACCTTCGTCTCAAGGCTTTTCATACAGGCCCAA  <  NC_002947/5205267‑5205215

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