New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 264548462 (0.680)3 (0.030) 3/252 NT 4.5% coding (1842/2313 nt) pbpC penicillin‑insensitive murein repair transglycosylase; inactive transpeptidase domain protein
?NC_000913 = 2645503 68 (0.770)coding (1823/2313 nt) pbpC penicillin‑insensitive murein repair transglycosylase; inactive transpeptidase domain protein
Rejected: Frequency below/above cutoff threshold.

AGCCAGGGTTCCAGCGGTAGCGGCCAGACATTTATCATTTCTTGTCGCGCTTGCGGGCAATCAGCGGCAACACGTTTGCCATTTTCATCCAGCCAGATGGGGAAACGAATGCCATTAATGCCTTCCTGCTCCGGCAGTAATAGAGTTGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/2645335‑2645484
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tggcTGCTGGACGGGAGTCAGCCGCCAACTCTATTACTGCCGGAGCAGGAAGGCATTAATGGCATTCGTTTCCCCATCTGGCTGGATGAAAATGGCAAACGTGTTGCCGCTGATTGCCCGCAAGCGCGACAAGAAATGATAAATGTCTGG  <  NC_000913/2645503‑2645358
                                                                                                                                                                                                                                                                                                         
AGCCAGGGTTCCAGCGGTAGCGGCCAGACAGTTATCATTTCTTGTCGCGCTTGCGGGCAATCAGCGGCAACACGTTTGCCATTTTCATCCAGCCAGATGGGGAAACGAATGGCATTAATTCCTTCCTGAGCCGGCAGTAATAGAGTTGG                                                                                                                                                      >  2:220041/1‑149
 GCCAGGGTTCCAGCGGTAGCGGCCAGACATTTATCATTTCTTGTCGCGCTTGCGGGCAATCAGCGGCAACACGTTTGCCATTTTCATCCAGCCAGATGGGGAAACGAATGCCATTAATGCCTTCCTGCTCCGGCAGTAATAGAGTTGGC                                                                                                                                                     >  7:52415/1‑149
    AGGGTTCCAGCGGTAGCGGCCAGACATTTATCATTTCTTGTCGCGCTTGCGGGCAATCAGCGGCAACACGTTTGCCATTTTCATCCAGCCAGATGGGGAAACGAATGCCATTAATGCCTTCCTGCTCCGGCAGTAATAGAGTTGGCTGC                                                                                                                                                  >  1:172180/1‑149
                                                                                                                CATTAATGCCTTCCTGCTCCGGCAGTAATAGAGTTGGCTGCTGGACGGGAGTCAGCCGCC                                                                                                                              <  1:172060/60‑1
                                                                                                                CATTAATGCCTTCCTGCTCCGGCAGTAATAGAGTTGGCTGCTGGACGGGAGTCAGCCGCC                                                                                                                              >  2:172060/1‑60
                                                                                                                          TTCCTGCTCCGGCAGTAATAGAGTTGGCGGCTGACTCCCGTCCAGCAGCCAACTCTATTACTGCCGGAGCAGGAAGGCATTAATGGCATTCGTTTCCCCATCTGGCTGGATGAAAATGGCAAACGTGTTGCCGCTGATTGCCCGCAAGC                           >  2:172180/1‑149
                                                                                                                                                    GCTGCTGGACGGGAGTCAGCCGCCAACTCTATTACTGCCGGAGCAGGAAGGCATTAATGGCATTCGTTTCCCCATCTGGCTGGATGAAAATGGCAAACGTGTTGCCGCTGATTGCCCGCAAGCGCGACAAGAAATGATAAATGTCTGG  <  4:352698/148‑1
                                                                                                                                                    GCTGCTGGACGGGAGTCAGCCGCCAACTCTATTACTGCCGGAGCAGGAAGGCATTAATGGCATTCGTTTCCCCATCTGGCTGGATGAAAATGGCAAACGTGTTGCCGCTGATTGCCCGCAAGCGCGACAAGAAATGATAAATGTCTGG  <  8:192854/148‑1
                                                                                                                                                                                                                                                                                                         
AGCCAGGGTTCCAGCGGTAGCGGCCAGACATTTATCATTTCTTGTCGCGCTTGCGGGCAATCAGCGGCAACACGTTTGCCATTTTCATCCAGCCAGATGGGGAAACGAATGCCATTAATGCCTTCCTGCTCCGGCAGTAATAGAGTTGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/2645335‑2645484
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tggcTGCTGGACGGGAGTCAGCCGCCAACTCTATTACTGCCGGAGCAGGAAGGCATTAATGGCATTCGTTTCCCCATCTGGCTGGATGAAAATGGCAAACGTGTTGCCGCTGATTGCCCGCAAGCGCGACAAGAAATGATAAATGTCTGG  <  NC_000913/2645503‑2645358

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.