New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 2640634 =31 (0.580)3 (0.060) 3/150 4.3 8.7% coding (51/2313 nt) pbpC peptidoglycan glycosyltransferase PbpC
?NZ_CP009273 2640664 = 35 (0.730)coding (21/2313 nt) pbpC peptidoglycan glycosyltransferase PbpC
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

GGGCGACCGGCGCGGCTGAAGATCTGCTGATTGTCAGACCGTAAATGCCTCGCTTGTTAACCAAACGCGGCTGCTGGATAACGTTGGCAGCCGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/2640728‑2640634
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcagccgcgTTTGGTTAACAAGCGAGGCATTTACGGTCTGACAATCAGCAGATCTTCAGCCGCGCCGGTCGCCCGCCATTGGGGAACATACAT  >  NZ_CP009273/2640664‑2640747
                                                                                                                                                                                   
gcGCGACCGGCGCGGCTGAAGATCTGCTGATTGTCAGACCGTAAATGCCTCGCTTGTTAACCAAACGCGGCTGCTGGATAACGTTGGCAG                                                                                           >  2:139046/3‑90
GGGCGACCGGCGCGGCTGAAGATCTGCTGATTGTCAGACCGTAAATGCCTCGCTTGTTAACCAAACGCGGCTGCTGGATAACGTTGGCAG                                                                                           >  2:1456379/1‑90
  GCGACCGGCGCGGCTGAAGATCTGCTGATTGTCAGACCGTAAATGCCTCGCTTGTTAACCAAACGCGGCTGCTGGATAACGTTGGCAGCC                                                                                         <  1:196284/90‑1
  GCGACCGGCGCGGCTGAAGATCTGCTGATTGTCAGACCGTAAATGCCTCGCTTGTTAACCAAACGCGGCTGCTGGATAACGTTGGCAGCC                                                                                         <  2:191124/90‑1
           GCGGCTGAAGATCTGCTGATTGTCAGACCGTAAATGCCTCGCTTGTTACCCAAGCGCGGCTGCTGGATAACGTGGGAAGCCGCGTTTGGT                                                                                <  1:1577509/90‑1
                        TGCTGATTGTCAGACCGTAAATGCCTCGCTTGTTAACCAAACGCGGCTGCTGGATAACGTTGGCAGCCGCGTTTGGTTAACAAGCGAGGC                                                                   <  1:1020988/90‑1
                                              GCCTCGCTTGTTAACCAAACGCGGCTGCTGGATAACGTTGGCAGCCGCGTTTGGTTAACAAGCGAGGCATTTACGGTCTGACAATCAGCA                                             <  1:603812/90‑1
                                                                                         GCCGCGTTTGGTTAACAAGCGCGGCATTTACGGTCTGAAAAACAGCAGATCTTCACCCGCCCCGGTCGCCCGCCATTGGGGAACATACAT  <  1:796536/90‑1
                                                                                                                                                                                   
GGGCGACCGGCGCGGCTGAAGATCTGCTGATTGTCAGACCGTAAATGCCTCGCTTGTTAACCAAACGCGGCTGCTGGATAACGTTGGCAGCCGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NZ_CP009273/2640728‑2640634
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcagccgcgTTTGGTTAACAAGCGAGGCATTTACGGTCTGACAATCAGCAGATCTTCAGCCGCGCCGGTCGCCCGCCATTGGGGAACATACAT  >  NZ_CP009273/2640664‑2640747

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