New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 395867 =25 (1.130)4 (0.200) 4/118 1.6 14.5% coding (421/1158 nt) ampH D‑alanyl‑D‑alanine‑ carboxypeptidase/endopeptidase; penicillin‑binding protein; weak beta‑lactamase
?NC_000913 395897 = 25 (1.270)coding (391/1158 nt) ampH D‑alanyl‑D‑alanine‑ carboxypeptidase/endopeptidase; penicillin‑binding protein; weak beta‑lactamase
Rejected: Frequency below/above cutoff threshold.

TCCGCCAGGCGCACGCGTGCCAACCTACAATGGGACGCCCATTACGCTGGTCAATCTGGCAACCCATACCAGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/395940‑395867
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑accagcgTAATGGGCGTCCCATTGTAGGTTGGCACGCGTGCGCCTGGCGGAGCATATTTACTTAGCGGATCGTT  >  NC_000913/395897‑395963
                                                                                                                                             
TCCGCCAGGCGCACGCGTGCCAACCTACAATGGGACGCCCATTACGCTGGTCAATCTGGCAACCCATA                                                                           >  1:542569/1‑68
TCCGCCAGGCGCACGCGTGCCAACCTACAATGGGACGCCCATTACGCTGGTCAATCTGGCAACCCATA                                                                           <  1:86940/68‑1
        GCGCACGCGTGCCAACCTACAATGGGACGCCCATTACGCTGGTCAATCTGGCAACCCATACCAGCGTA                                                                   <  1:411063/68‑1
                              TGGGACGCCCATTACGCTGGTAAATCTGGAAACCCATACAAGCGTAATGGGCGTCCCAGTGTAGGTTG                                             >  2:581339/1‑68
                                                  TCAATCTGGCAACCCATACCAGCGTAATGGGCGTCCCATTGTAGGTTGGCACGCGTGCGCCTGGC                            <  1:738066/65‑1
                                                  TCAATCTGGCAACCCATACCAGCGTAATGGGCGTCCCATTGTAGGTTGGCACGCGTGCGCCTGGC                            >  2:738066/1‑65
                                                                    CCAGCGTAATGGGCGTCCCATTGTAGGTTGGCACGCGTGCGCCTGGCGGAGCATATTTACTTAGCGGA       >  2:141473/1‑68
                                                                      AGCGTAATTGGCGTCCCATTGTAGGTTTGCACGCGTGCGCCTGGCGGAGCATATTTACTTAGCGGATC     <  2:662846/68‑1
                                                                      AGCGTAATGGGCGTCCCATTGTAGGTTGGCACGCGTGCGCCTGGCGGAGCATATTTACTTAGCGGATC     <  2:51267/68‑1
                                                                      AGCGTAATGGGCGTCCCATTGTAGGTTGGCACGCGTGCGCCTGGCGGAGCATATTTACTTAGCGGATC     <  2:86940/68‑1
                                                                        CGTAATTGGCGTCCCATTGTAGGTTGGCACGCGTGCGCCTGGCGGAGCATATTTACTTAGCGGATCGT   <  2:533494/68‑1
                                                                         GTAATGGGCGTCCCATTGTAGGTTGGCACGCGTGCGCCTGGCGGAGCATATTTACTTAGCGGATCGTT  >  2:91525/1‑68
                                                                                                                                             
TCCGCCAGGCGCACGCGTGCCAACCTACAATGGGACGCCCATTACGCTGGTCAATCTGGCAACCCATACCAGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/395940‑395867
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑accagcgTAATGGGCGTCCCATTGTAGGTTGGCACGCGTGCGCCTGGCGGAGCATATTTACTTAGCGGATCGTT  >  NC_000913/395897‑395963

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