New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4631855181 (1.150)8 (0.050) 5/166 NT 4.4% coding (1123/1938 nt) slt lytic murein transglycosylase, soluble
?NC_000913 = 4631875 180 (1.220)coding (1143/1938 nt) slt lytic murein transglycosylase, soluble
Rejected: Frequency below cutoff threshold.

GGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGATTTATTGCTGGAACGCGGACGTGAAGCTGAAGCAAAAGAGATTTTGCATCAACTCATGCAAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4631752‑4631855
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgcaacCATCGGGTAGAAACCACGCTGTTGCATGAGTTGATGCAAAATCTCTTTTGCTTCAGCTTCACGTCCGCGTTCCAGCAATAAATCCGCCTGCCAGTAACG  <  NC_000913/4631875‑4631777
                                                                                                                                                                                                           
GGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGATTTATTGCTGGAACGCGGACGTGAAGCTGAAGCAAAAGAGATTTTGCATCAACTCATG                                                                                                         >  1:3850670/1‑100
GGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGATTTATTGCTGGAACGCGGACGTGAAGCTGAAGAAAAAGAGATTTTGCATCAACTCATG                                                                                                         <  1:1387850/100‑1
GGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGAGATATTGCTGGAACGCGGACGTGAAGCTGAAGCAAAAGAGATCTTGCATCAACTCATG                                                                                                         >  1:3581946/1‑100
GGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGAAGTATTGCTGGAACGCGGACGTGAAGCTGAAGCAAAAGAGATCTTGCATTAACTCATG                                                                                                         >  1:1652455/1‑100
 GAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGATTTATTGCTGGAACGCGGACGTGAAGCTGAAGCAAAAGAGATTTTGCATCAACTCATGC                                                                                                        >  1:1352729/1‑100
 GAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGATTTATTGCTGGAACGCGGACGTGAAGCTGAAGCAAAAGAGATTTTGCATCAACTCATGC                                                                                                        <  1:1134739/100‑1
 GAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGATTTATTGCTGGAACGCGGACGTGAAGCTGAAGCAAAAGAGATTTTGCATCAACTCATGC                                                                                                        >  1:931777/1‑100
 GAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGATTTATTGCTGGAACGCGGACGTGAAGCTGAAGCAAAAGAGATTTTGCATCAACTCATGC                                                                                                        <  1:1849686/100‑1
  AAGCGAAAGAGAAAGATGAATGGCGTTACTGGTAGGTGGATTTATTGCTGGAACGCGGACGTGAAGCTGAAGCAAAATATTTTTTGCACCAACTCATGCA                                                                                                       <  1:1827912/100‑1
  AAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGATTTATTGCTGGAACGCGGACGTGAAGCTGAAGCAAAAGAGATTTTGCATCAACTCATGCA                                                                                                       <  2:941826/100‑1
  AAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGATTTATTGCTGGAACGCGGACGTGAAGCTGAAGCAAAAGAGATTTTGCATCAACTCATGCA                                                                                                       <  2:1451158/100‑1
  AAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGATTTATTGCTGGAACGCGGACGTGAAGCTGAAGCAAAAGAGATTTTGCATCAACTCATGCA                                                                                                       <  1:2856464/100‑1
  AAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGATTTATTGCTGGAACGCGGACGTGAAGCTGAAGCAAAAGAGATTTTGCATCAACTCATGCA                                                                                                       >  1:1426518/1‑100
  AAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGATTTATTGCTGGAACGCGGACGTGAAGCTGAAGCAAAAGAGATTTTGCATCAACTCATGCA                                                                                                       <  1:2277436/100‑1
  AAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGATTTATTGCTGGAACGCGGACGTGAAGCTGAAGCAAAAGAGATTTTGCATCAACTCATGCA                                                                                                       <  1:1829090/100‑1
    GCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGATTTATTGCTGGAACGCGGACGTGAAGCTGAAGCAAAAGAGATTTTGCATCAACTCATGCAAC                                                                                                     <  2:1249638/100‑1
    GCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGATTTATTGCTGGAACGCGAAAGTGAAGCTGAAGCAAAAGAGAGTTTGCATCAACTCAAGCAA                                                                                                      >  1:691922/1‑99
     CGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGATTTATTGCTGGAACGCGGACGTGAAGCTGAAGCAAAAGAGATTTTGCATCAACTCATGCAACC                                                                                                    >  2:2518175/1‑100
                                GGCAGGCGGATTTATTGCTGGAACGCGGACGTGAAGCTGAAGCAAAAGAGATTTTGCATCAACTCATGCAACCATCGGGTAGAAACCACGCTGTTGCATG                                                                         <  1:2937553/100‑1
                                GGCAGGCGGATTTATTGCTGGAACGCGGACGTGAAGCTGAAGCAAAAGAGATTTTGCATCAACTCATGCAACCATCGGGTAGAAACCACGCTGTTGCATG                                                                         <  1:1773834/100‑1
                                  CAGGCGGATTTATTGCTGGAACGCGGACGTGAAGCTGAAGCAAAAGAGATTTTGCATCAACTCATGCAACCATCGGGTAGAAACCACGCTGTTGCATGAG                                                                       <  2:2863268/100‑1
                                  CAGGCGGATTTATTGCTGGAACGCGGACGTGAAGCTGAAGCAAAAGAGATTTTGCATCAACTCATGCAACCATCGGGTAGAAACCACGCTGTTGCATGAG                                                                       <  1:2518163/100‑1
                                   acGCGGATTTATTGCTGGAACGCGGACGTGAAGCTTAAGCAAAACAGATTTTGCATCAACTCATGCAACCATCGGGTAGAAACCACGCTGTTGCATGAGT                                                                      <  2:3816480/98‑1
                                   AGGAGGATTTATTGCTGGAACGCGGACGTGAAGCTGAAGCAAAAGAGATTTTGCATCAACTCATGCAACCATCGGGTAGAAACCACGCTGTTGCATGAG                                                                       <  1:1652386/99‑1
                                    GGCGGATTGATTGCTGGAACGCGGACGTGAAGCTGAAGCAAAAGAGATTTTGCATCAACTCATGCAACCAGCGGGTAGAAACCACGCTGTTGCATGAGTT                                                                     <  1:856038/100‑1
                                      CGGATTTATTGCTGGAACGCGGACGTGAAGCTGAAGCAAAAGAGATTTTGCATCAACTCATGCAACCATCGGGTAGAAACCACGCTGTTGCATGAGTTGA                                                                   <  1:2780816/100‑1
                                                                                                   GCAACCATCGGGTAGAAACCACGCTGTTGCATGAGTTGATGCAAAATCTCTTTTGCTTCAGCTTCACGTCCGCGTTCCAGCAATAAATCCGCCTGCCAGT      >  2:1199087/1‑100
                                                                                                    CAACCATCGGGTAGAAACCACGCTGTTGCATGAGTTGATGCAAAATCTCTTTTGCTTCAGCTTCACGTCCGCGTTCCAGCAATAAATCCGCCTGCCAGTA     >  2:1065178/1‑100
                                                                                                    CAACCATCGGGTAGAAACCACGCTGTTGCATGAGTTGATGCAAAATCTCTTTTGCTTCAGCTTCACGTCCGCGTTCCAGCAATAAATCCGCCTGCCAGTA     >  1:1434115/1‑100
                                                                                                    CAACCATCGGGTAGAAACCACGCTGTTGCATGAGTTGATGCAAAATCTCTTTTGCTTCAGCTTCACGTCCGCGTTCCAGCAATAAATCCGCCTGCCAGT      >  1:3364782/1‑99
                                                                                                     AACCATCGGGTAGAAACCACGCTGTTGCATGAGTTGATGCAAAATCTCTTTTGCTTCAGCTTCACGTCCGCGTTCCAGCAATAAATCCGCCTGCCAGTAA    <  1:3654237/100‑1
                                                                                                       CCATCGGGTAGAAACCACGCTGTTGCATGAGTTGATGCAAAATCTCTTTTGCTTCAGCTTCACGTCCGCGTTCCAGCAATAAATCCGCCTGCCAGTAACG  >  2:2324960/1‑100
                                                                                                                                                                                                           
GGAAGCGAAAGAGAAAGATGAATGGCGTTACTGGCAGGCGGATTTATTGCTGGAACGCGGACGTGAAGCTGAAGCAAAAGAGATTTTGCATCAACTCATGCAAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4631752‑4631855
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tgcaacCATCGGGTAGAAACCACGCTGTTGCATGAGTTGATGCAAAATCTCTTTTGCTTCAGCTTCACGTCCGCGTTCCAGCAATAAATCCGCCTGCCAGTAACG  <  NC_000913/4631875‑4631777

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