New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3950464 =71 (0.980)17 (0.240) 9/68 2.8 25.6% intergenic (+44/‑43) ilvL/ilvX ilvXGMEDA operon leader peptide/uncharacterized protein IlvX
?NC_000913 = 4631789 28 (0.390)coding (1057/1938 nt) slt soluble lytic murein transglycosylase

CTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGG  >  NC_000913/3950430‑3950498
                                  |                                  
cTTAACGAACTAAGACCCCCGCACCGAAAGGTCCg                                    <  1:2842519/35‑1 (MQ=255)
 ttAACGAACTAAGACCCCCGCACCGAAAGGTCCgg                                   <  1:3502644/35‑1 (MQ=255)
 cgaaagagaaagatgaatggcgttactggcagggg                                   <  1:7557459‑M1/2‑1 (MQ=255)
  tAACGAACTAAGACCCCCGCACCGAAAGGTCCggg                                  >  1:379175/1‑35 (MQ=255)
  gaaagagaaagatgaatggcgttactggcaggggg                                  <  1:449134‑M1/3‑1 (MQ=255)
   aaagagaaagatgaatggcgttactggcagggggg                                 <  1:6487079‑M1/4‑1 (MQ=255)
   aaagagaaagatgaatggcgttactggcagggggg                                 <  1:444287‑M1/4‑1 (MQ=255)
   aaagagaaagatgaatggcgttactggcagggggg                                 <  1:4191996‑M1/4‑1 (MQ=255)
   aaCGAACTAAGACCCCCGCACCGAAAGGTCCgggg                                 >  1:4806066/1‑35 (MQ=255)
    aagagaaagatgaatggcgttactggcaggggggg                                <  1:7116817‑M1/5‑1 (MQ=255)
    aagagaaagatgaatggcgttactggcaggggggg                                <  1:5644267‑M1/5‑1 (MQ=255)
    aCGAACTAAGACCCCCGCACCGAAAGGTCCggggg                                <  1:904486/35‑1 (MQ=255)
    aCGAACTAAGACCCCCGCACCGAAAGGTCCggggg                                >  1:7388702/1‑35 (MQ=255)
     agagaaagatggatggcgttactggcagggggggg                               <  1:5394456‑M1/6‑2 (MQ=255)
     agagaaagatgaatggcgttactggcagggggggg                               <  1:4381184‑M1/6‑2 (MQ=255)
      gAACTAAGACCCCCGCACCGAAAGGTCCGGGGGtt                              <  1:3948932/35‑1 (MQ=255)
      gAACTAAGACCCCCGCACCGAAAGGTCCGGGGGtt                              <  1:3484832/35‑1 (MQ=255)
       aaCTAAGACCCCCGCACCGAAAGGTCCGGGGGttt                             >  1:6695159/1‑35 (MQ=255)
       aaCTAAGACCCCCGCACCGAAAGGTCCGGGGGttt                             >  1:7343586/1‑35 (MQ=255)
        gaaagatgaatggcgttactggcaggGGGGGgttt                            <  1:2684545‑M1/9‑1 (MQ=255)
        aCTAAGCCCCCCGCGCCGAAAGGTCCGGGGGtttt                            <  1:4115610/35‑1 (MQ=25)
        aCTAAGACCCCCGCACCGAAAGGTCCGGGGGtttt                            >  1:6158227/1‑35 (MQ=255)
        aCTAAGACCCCCGCACCGAAAGGTCCGGGGGtttt                            >  1:2622819/1‑35 (MQ=255)
        aCTAAGACCCCCGCACCGAAAGGCCCGGGGGgttt                            >  1:3118793/1‑35 (MQ=38)
        aCTAAGACCCCCGCACCGAAAAGTCCGGGGGtttt                            >  1:3278286/1‑35 (MQ=255)
         cTAAGACCCCCGCACCGAAAGGTCCGGGGGttttt                           <  1:7620874/35‑1 (MQ=255)
         cTAAGACCCCCGCACCGAAAGGTCCGGGGGttttt                           <  1:4997765/35‑1 (MQ=255)
         cTAAGACCCCCGCACCGAAAGGTCCGGGGGttttt                           <  1:2099331/35‑1 (MQ=255)
         aaagatgaatggcgttactggcaggGGGGGttttt                           <  1:7226863‑M1/10‑1 (MQ=255)
          tAAGACCCCCGCACCGAAAGGTCCGGGGGtttttt                          <  1:2572232/35‑1 (MQ=255)
          aagatgaatggcgttactggcaggGGGGGtttttt                          <  1:6907900‑M1/11‑1 (MQ=255)
           agatgaatggcgttactggcaggGGGGtgtttttg                         <  1:7748630‑M1/12‑2 (MQ=255)
           aaGACCCCCGCACCGAAAGGTCCGGGGGttttttt                         >  1:2406233/1‑35 (MQ=255)
            gatgaatggcgttactggcaggGGGGTttattttt                        <  1:184535‑M1/13‑1 (MQ=255)
            gatgaatggcgttactggcaggGGGGGtttttttt                        <  1:5529506‑M1/13‑1 (MQ=255)
            aGACCCCCGCACCGAAAGGTCCGGGGGtttttttt                        >  1:3760253/1‑35 (MQ=255)
            aGACCCCCGCACCGAAAGGTCCGGGGGtttttttt                        <  1:508175/35‑1 (MQ=255)
            aGACCCCCGCACCGAAAGGTCCGGGGGgttttttt                        >  1:3106252/1‑35 (MQ=255)
             gACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTg                       >  1:6999807/1‑35 (MQ=255)
             atgaatggcgttactggcaggGGGGGTTTTTTTTg                       <  1:5557622‑M1/14‑1 (MQ=255)
             atgaatggcgttactggcaggCGGGGGTTTTTTTg                       <  1:5457279‑M1/14‑1 (MQ=255)
              aCCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGa                      <  1:7059295/35‑1 (MQ=255)
                 cccGCACCGAAAGGTCCGGGGGTTTTTTTTGACCt                   <  1:1909023/35‑1 (MQ=255)
                 cccGCACCGAAAGGTCCGGGGGTTTTTTTTGACCt                   <  1:7293928/35‑1 (MQ=255)
                  ccGCACCGAAAGGTCCGGGGGTTTTTTTTGACCtt                  <  1:6727945/35‑1 (MQ=255)
                   cGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTa                 <  1:4746799/35‑1 (MQ=255)
                   cGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTa                 <  1:1439067/35‑1 (MQ=255)
                   cGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTa                 <  1:2830750/35‑1 (MQ=255)
                   cGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTa                 <  1:6191657/35‑1 (MQ=255)
                   cGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTa                 <  1:3635332/35‑1 (MQ=255)
                    gCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTaa                <  1:2143369/35‑1 (MQ=255)
                    gCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTaa                <  1:2896721/35‑1 (MQ=255)
                    gCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTaa                <  1:3488482/35‑1 (MQ=255)
                     cACCGAAAGGTCCGGGGGTTTTTTTTGACCTTaaa               <  1:3718092/35‑1 (MQ=255)
                     cACCGAAAGGTCCGGGGGTTTTTTTTGACCTTaaa               <  1:2373823/35‑1 (MQ=255)
                     cACCGAAAGGTCCGGGGGTTTTTTTTGACCTTaaa               >  1:7029232/1‑35 (MQ=255)
                     cACCGAAAGGTCCGGGGGTTTTTTTTGACCTTaaa               <  1:3566064/35‑1 (MQ=255)
                      aCCGAAAGGTCCGGGGGTTTTTTTTGACCTTaaaa              <  1:1988842/35‑1 (MQ=255)
                       ccGAAAGGTCCGGGGGTTTTTTTTGACCTTaaaaa             >  1:7459805/1‑35 (MQ=255)
                       ccGAAAGGTCCGGGGGTTTTTTTTGACCTTaaaaa             >  1:6963197/1‑35 (MQ=255)
                       ccGAAAGGTCCGGGGGTTTTTTTTGACCTTaaaaa             >  1:5230374/1‑35 (MQ=255)
                       ccGAAAGGTCCGGGGGTTTTTTTTGACCTTaaaaa             <  1:5689443/35‑1 (MQ=255)
                       ccGAAAGGTCCGGGGGTTTTTTTTGACCTTaaaaa             >  1:5395876/1‑35 (MQ=255)
                       ccGAAAGGTCCGGGGGTTTTTTTTGACCTTaaaaa             >  1:3044194/1‑35 (MQ=255)
                       ccGAAAGGTCCGGGGGTTTTTTTTGACCTTaaaaa             >  1:4271727/1‑35 (MQ=255)
                          aaaGGTCCGGGGGTTTTTTTTGACCTTAAAAACAt          >  1:1039944/1‑35 (MQ=255)
                          aaaGGTCCGGGGGTTTTTTTTGACCTTAAAAACAt          <  1:872809/35‑1 (MQ=255)
                          aaaGGTCCGGGGGTTTTTTTTGACCTTAAAAACAt          <  1:7797785/35‑1 (MQ=255)
                          aaaGGTCCGGGGGTTTTTTTTGACCTTAAAAACAt          >  1:2432091/1‑35 (MQ=255)
                          aaaGGTCCGGGGGTTTTTTTTGACCTTAAAAACAt          <  1:1133241/35‑1 (MQ=255)
                          aaaGGTCCGGGGGTTTTTTTTGACCTTAAAAACAt          <  1:3092667/35‑1 (MQ=255)
                          aaaGGTCCGGGGGTTTTTTTTGACCTTAAAAACAt          <  1:1284489/35‑1 (MQ=255)
                          aaaGGTCCGGGGGTTTTTTTTGACCTTAAAAACAt          <  1:1906352/35‑1 (MQ=255)
                          aaaGGTCCGGGGGTTTTTTTTGACCTTAAAAACAt          <  1:7117213/35‑1 (MQ=255)
                           aaGGTCCGGGGGTTTTTTTTGACCTTAAAAACATa         >  1:448948/1‑35 (MQ=255)
                           aaGGTCCGGGGGTTTTTTTTGACCTTAAAAACATa         >  1:5099439/1‑35 (MQ=255)
                           aaGGTCCGGGGGTTTTTTTTGACCTTAAAAACATa         >  1:5845566/1‑35 (MQ=255)
                            aGGTCCGGGGGTTTTTTTTGACCTTAAAAACATaa        <  1:1347308/35‑1 (MQ=255)
                            aGGTCCGGGGGTTTTTTTTGACCTTAAAAACATaa        <  1:4636528/35‑1 (MQ=255)
                            aGGTCCGGGGGTTTTTTTTGACCTTAAAAACATaa        >  1:1478836/1‑35 (MQ=255)
                             ggTCCGGGGGTTTTTTTTGACCTTAAAAACATAAc       >  1:1878391/1‑35 (MQ=255)
                             ggTCCGGGGGTTTTTTTTGACCTTAAAAACATAAc       <  1:3803173/35‑1 (MQ=255)
                              gTCCGGGGGTTTTTTTTGACCTTAAAAACATAAcc      <  1:6944183/35‑1 (MQ=255)
                              gTCCGGGGGTTTTTTTTGACCTTAAAAACATAAcc      >  1:6820882/1‑35 (MQ=255)
                              gTCCCGGGGGTTTTTTTGACCTTAAAAACATAAcc      <  1:2183687/35‑1 (MQ=25)
                               tCCGGGGGTTTTTTTTGACCTTAAAAACATAACCg     >  1:6863040/1‑35 (MQ=255)
                                ccGGGGGTTTTTTTTGACCTTAAAAACATAACCga    >  1:5803575/1‑35 (MQ=255)
                                 cGGGGGTTTTTTTTGACCTTAAAAACATAACCgag   <  1:5447828/35‑1 (MQ=255)
                                  gggggTTTTTTTTGACCTTAAAAACATAACCgagg  >  1:5729455/1‑35 (MQ=255)
                                  |                                  
CTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGG  >  NC_000913/3950430‑3950498

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