New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 2548123 =26 (1.180)3 (0.160) 3/114 1.8 12.4% coding (22/1305 nt) yfeW penicillin binding protein PBP4B; weak DD‑carboxypeptidase activity
?NC_000913 2548199 = 20 (1.050)coding (98/1305 nt) yfeW penicillin binding protein PBP4B; weak DD‑carboxypeptidase activity
Rejected: Frequency below/above cutoff threshold.

TTCACCACGCTCTTTGGCTGCCGTAACGTCAATCTGGACTGATATCATGAAACGGACAATGCTCTATCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAGCGCCGCAAAATACCCTGTTCTGACAGAAAGCTCGCC  >  NC_000913/2548056‑2548190
                                                                   |                                                                   
taaccgacatcaacttgctggctaatccagcgatccatctggttaagccgttcgactttaatccctgc                                                                     >  2:83657‑M1/68‑68 (MQ=255)
   ccgacatcaacttgctggctaatccagcgatccatctggttaagccgttcgacgttaaaccctgCttt                                                                  <  2:790893‑M1/4‑1 (MQ=255)
    ccacGCGCTTTGGCTGCCGTAACGTCAATCTGGACTGATATCATGAAACGGACAATGCTCTATCTTtc                                                                 <  2:71095/68‑1 (MQ=255)
      acGCTCTTTGGCTGCCGTAACGTCAATCTGGACTGATATCATGAAACGGACAATGCTCTATCTTtctc                                                               <  1:19576/68‑1 (MQ=255)
        gCTCTTTGGCTGCCGTAACGTCAATCTGGACTGATATCATGAAACGGACAATGCTCTATCTTTCTctg                                                             <  1:801089/68‑1 (MQ=255)
        atcaacttgctggctaatccagcgatccatctggttaagccgttcgacgttaaaccctgCTTTCTctg                                                             >  2:382449‑M1/60‑68 (MQ=255)
            acttgctggctaatccagcgatccatctggttaagccgttcgacgttaaaccctgc                                                                     <  1:542163‑M1/1‑1 (MQ=255)
            acttgctggctaatccagcgatccatctggttaagccgttcgacgttaaaccctgc                                                                     >  2:542163‑M1/56‑56 (MQ=255)
                gctggctaatccagcgatccatctggttaagccgttcgacgttaaaccctgCTTTCTCTGCTGGCTGt                                                     >  1:27594‑M1/52‑68 (MQ=255)
                gCTGCCGTAACGTCAATCTGGACTGATATCATGAAACGGACAATGCTCTATCTTTCTCTGCTGGCTGt                                                     >  1:148158/1‑68 (MQ=255)
                gCTGCCGTAACGTCAATCTGGACTGATATCATGAAACGGACAATGCTCTATCTTTCTCTGCTGGCTGt                                                     >  2:750649/1‑68 (MQ=255)
                    ccGTAACGTCAATCTGGACTGATATCATGAAACGGACAATGCTCTATCTTTCTCTGCTGGCTGTTAGt                                                 <  1:750649/68‑1 (MQ=255)
                        aaCGTCAATCTGGACTGATATCATGAAACGGACAATGCTCTATCTTTCTCTgctgg                                                         >  1:177318/1‑56 (MQ=255)
                        aaCGTCAATCTGGACTGATATCATGAAACGGACAATGCTCTATCTTTCTCTgctgg                                                         <  2:177318/56‑1 (MQ=255)
                                tCTGGACTGATATCATGAAACGGACAATGCTCTATCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAgc                                     >  2:402970/1‑68 (MQ=255)
                                tCTGGACTGATATCATGAAACGGACAATGCTCTATCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAgc                                     <  1:88673/68‑1 (MQ=255)
                                tCTGGACTGATATCATGAAACGGACAATGCTCTATCTTTCTCTGCTGGCTGTTAGTTGTAGCGGCAgc                                     <  1:154999/68‑1 (MQ=255)
                                 cTGGACTGATATCATGAAACGGACAATGCTCTATCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAgcg                                    >  2:526966/1‑68 (MQ=255)
                                     aCTGATATCATGAAACGGACAATGCTCTATCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAGcgccgc                                <  1:163324/68‑1 (MQ=255)
                                     aCTGATATCATGAAACGGACAATGCTCTATCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAGcgccgc                                >  2:74040/1‑68 (MQ=255)
                                     aCTGATATCATGAAACGGACAATGCTCTATCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAGcgccgc                                >  1:219104/1‑68 (MQ=255)
                                     aCTGATATCATGAAACGGACAATGCTCTATCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAGcgccgc                                <  1:526385/68‑1 (MQ=255)
                                      ctggttaagccgttcgacgttaaaccctgCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAGCGCCGCa                               <  2:167147‑M1/39‑1 (MQ=255)
                                       tggttaagccgttcgacgttaaaccctgCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAGCGCCGCaa                              >  1:167147‑M1/29‑68 (MQ=255)
                                       tGATATCATGAAACGGACAATGCTCTATCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAGCGCCGCaa                              >  2:761310/1‑68 (MQ=255)
                                       tGATATCATGAAACGGACAATGCTCTATCTTTCTCTGATGGCTGTTAGTTGTAGCGTCAGCGCCGCaa                              <  1:87001/68‑1 (MQ=255)
                                          tatCATGAAACGGACAATGCTCTATCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAGCGCCGCAAAAt                           >  2:511890/1‑68 (MQ=255)
                                          tatCATGAAACGGACAATGCTCTATCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAGCGCCGCAAAAt                           >  2:203857/1‑68 (MQ=255)
                                                  aaCGGACAATGCTCTATCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAGCGCCGCAAAATACCCTGtt                   >  2:303749/1‑68 (MQ=255)
                                                  aaCGGACAATGCTCTATCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAGCGCCGCAAAATACCCTGtt                   >  2:302744/1‑68 (MQ=255)
                                                    cGGACAATGCTCTATCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAGCGCCGCAAAATACCCTGTTCt                 <  1:402971/68‑1 (MQ=255)
                                                          aTGCTCTATCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAGCGCCGCAAAATACCCTGTTCTGACAGa           >  1:375203/1‑68 (MQ=255)
                                                             ctctATCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAGCGCCGCAAAAt                           <  2:554195/49‑1 (MQ=255)
                                                             ctctATCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAGCGCCGCAAAAt                           >  1:554195/1‑49 (MQ=255)
                                                                   CTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAGCGCCGCAAAATACCCTGTTCTGACAGAAAGCTCGcc  >  2:219079‑M1/1‑68 (MQ=255)
                                                                   CTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAGCGCCGCAAAATACCCTGTTCTGACAGAAAGCTCGcc  <  1:219080‑M1/68‑1 (MQ=255)
                                                                   |                                                                   
TTCACCACGCTCTTTGGCTGCCGTAACGTCAATCTGGACTGATATCATGAAACGGACAATGCTCTATCTTTCTCTGCTGGCTGTTAGTTGTAGCGTCAGCGCCGCAAAATACCCTGTTCTGACAGAAAGCTCGCC  >  NC_000913/2548056‑2548190

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