New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 3342514 =38 (1.110)3 (0.090) 3/234 2.4 8.0% coding (655/729 nt) mtgA monofunctional biosynthetic peptidoglycan transglycosylase
?NZ_CP009273 3342530 = 33 (1.020)coding (639/729 nt) mtgA monofunctional biosynthetic peptidoglycan transglycosylase
Rejected: Frequency below/above cutoff threshold.

AACCCGCGAGCAAACTTACCCGGTCGGAAGCTGCATTACTGGCAGCTGTATTACCTAATCCACTTCGTTTCAAAGTCTCCTCGCCATCGGGCTACGTGCGTAGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/3342618‑3342514
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgtagccCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCAGCTTCGACGCCAAACACACCGTC  >  NZ_CP009273/3342530‑3342664
                                                                                                                                                                                                                                                
AACCCGCGAGCAAACTTACCCGGTCGGAAGCTGCATTACTGGCAGCTGTATTACCTAATCCACTTCGTTTCAAAGTCTCCTCGCCATCGGGCTACGTGCGTAGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTA                                                                                                       <  1:156010/139‑1
                       TCGGAAGCTGCATTACTGGCAGCTGTATTACCTAATCCACTTCGTTTCAAAGTCTCCTCGCCATCGGGCTACGTGCGTAGCCCGATGGCGAGGAGACTTTGAAACGA                                                                                                                <  1:263095/107‑1
                       TCGGAAGCTGCATTACTGGCAGCTGTATTACCTAATCCACTTCGTTTCAAAGTCTCCTCGCCATCGGGCTACGTGCGTAGCCCGATGGCGAGGAGACTTTGAAACGA                                                                                                                >  2:263095/1‑107
                                                        AATCCACTTCGTTTCAAAGTCTCCTCGCCATCGGGCTACGTGCG                                                                                                                                              <  1:330753/44‑1
                                                        AATCCACTTCGTTTCAAAGTCTCCTCGCCATCGGGCTACGTGCG                                                                                                                                              >  2:330753/1‑44
                                                                                                     AGCCCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCAGCTTCGACGCCAAACACACCGTC  >  2:276885/1‑139
                                                                                                                                                                                                                                                
AACCCGCGAGCAAACTTACCCGGTCGGAAGCTGCATTACTGGCAGCTGTATTACCTAATCCACTTCGTTTCAAAGTCTCCTCGCCATCGGGCTACGTGCGTAGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/3342618‑3342514
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgtagccCGATGGCGAGGAGACTTTGAAACGAAGTGGATTAGGTAATACAGCTGCCAGTAATGCAGCTTCCGACCGGGTAAGTTTGCTCGCGGGTTTGTGGAAATAACGTTGTGCCGCAGCTTCGACGCCAAACACACCGTC  >  NZ_CP009273/3342530‑3342664

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