New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 11971360 (0.000)18 (0.340)
+CAGATAATT
9/50 1.8 49.5% coding (397/666 nt) ymfD e14 prophage; putative SAM‑dependent methyltransferase
?NC_000913 3979137 = 50 (0.690)coding (537/1353 nt) wzyE enterobacterial common antigen polymerase

CGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCG  >  NC_000913/3979103‑3979171
                                  |                                  
cGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGc                                    <  1:7754144/35‑1 (MQ=255)
cGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGc                                    >  1:3577778/1‑35 (MQ=255)
 gATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCa                                   >  1:3680707/1‑35 (MQ=255)
 gATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCa                                   <  1:5199869/35‑1 (MQ=255)
 gATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCa                                   <  1:5460062/35‑1 (MQ=255)
 gATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCa                                   <  1:2049566/35‑1 (MQ=255)
 gATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCa                                   <  1:7314119/35‑1 (MQ=255)
  aTGCTGGTGGTCTACTTTCTGCGCCAGGACAGCaa                                  >  1:3838801/1‑35 (MQ=255)
  aTGCTGGTGGTCTACTTTCTGCGCCAGGACAGCaa                                  >  1:121923/1‑35 (MQ=255)
  aTGCTGGTGGTCTACTTTCTGCGCCAGGACAGCaa                                  <  1:4999769/35‑1 (MQ=255)
  aTGCTGGTGGTCTACTTTCTGCGCCAGGACAGCaa                                  <  1:5084836/35‑1 (MQ=255)
  aTGCTGGTGGTCTACTTTCTGCGCCAGGACAGCaa                                  <  1:2482979/35‑1 (MQ=255)
   tGCTGGTGGTCTACTTTCTGCGCCAGGACAGCaaa                                 <  1:854723/35‑1 (MQ=255)
   tGCTGGTGGTCTACTTTCTGCGCCAGGACAGCaaa                                 >  1:1570101/1‑35 (MQ=255)
    gCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAg                                >  1:466424/1‑35 (MQ=255)
     cTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGc                               >  1:4591897/1‑35 (MQ=255)
      tggtggTCTACTTTCTGCGCCAGGACAGCAAAGCg                              <  1:7461749/35‑1 (MQ=255)
      tggtggTCTACTTTCTGCGCCAGGACAGCAAAGCg                              <  1:478278/35‑1 (MQ=255)
      tcttgatgctaagaactatcagataattCAAAGCg                              <  1:153616‑M2/7‑1 (MQ=255)
       ggtggtCTACTTTCTGCGCCAGGACAGCAACGCGt                             >  1:7675047/1‑35 (MQ=255)
       ggtggtCTACTTTCTGCGCCAGGACAGCAAAGCGt                             >  1:7507164/1‑35 (MQ=255)
       ggtggtCTACTTTCTGCGCCAGGACAGCAAAGCGt                             >  1:1067782/1‑35 (MQ=255)
       ggtggtCTACTTTCTGCGCCAGGACAGCAAAGCGt                             >  1:4172134/1‑35 (MQ=255)
       cttgatgctaagaactatcagataattCAAAGCGt                             >  1:2125030‑M2/28‑35 (MQ=255)
        ttgatgctaagaactatcagataattCAAAGCGTg                            >  1:5188728‑M2/27‑35 (MQ=255)
        ttgatgctaagaactatcagataattCAAAGCGTg                            >  1:4844387‑M2/27‑35 (MQ=255)
        ttgatgctaagaactatcagataattCAAAGCGTg                            >  1:3298595‑M2/27‑35 (MQ=255)
        gtggtCTACTTTCTGCGCCAGGACAGCAAAGCGTg                            >  1:633424/1‑35 (MQ=255)
        gtggtCTACTTTCTGCGCCAGGACAGCAAAGCGTg                            >  1:2834595/1‑35 (MQ=255)
          ggtCTACTTTCTGCGCCAGGACAGCAAAGCGTGGc                          <  1:4876928/35‑1 (MQ=255)
          gatgctaagaactatcagataattCAAAGCGTGGc                          >  1:6303147‑M2/25‑35 (MQ=255)
          gatgctaagaactatcagataattCAAAGCGTGGc                          <  1:1060763‑M2/11‑1 (MQ=255)
          gatgctaagaactatcagataattCAAAGCGTGGc                          >  1:6225558‑M2/25‑35 (MQ=255)
          gatgctaagaactatcagataattCAAAGCGTGGc                          >  1:3214146‑M2/25‑35 (MQ=255)
          gatgctaagaactatcagataattCAAAGCGTGGc                          >  1:4749639‑M2/25‑35 (MQ=255)
          gatgctaagaactatcagataattCAAAGCGTGGc                          <  1:407697‑M2/11‑1 (MQ=255)
           gtCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCt                         >  1:682743/1‑35 (MQ=255)
           gtCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCt                         >  1:6616659/1‑35 (MQ=255)
           gtCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCt                         >  1:5774302/1‑35 (MQ=255)
           atgctaagaactatcagataattCAAAGCGTGGCt                         >  1:1701366‑M2/24‑35 (MQ=255)
           atgctaagaactatcagataattCAAAGCGTGGCt                         >  1:2648770‑M2/24‑35 (MQ=255)
           atgctaagaactatcagataattCAAAGCGTGGCt                         >  1:2505700‑M2/24‑35 (MQ=255)
            tgctaagaactatcagataattCAAAGCGTGGCTg                        <  1:4142612‑M2/13‑1 (MQ=255)
            tgctaagaactatcagataattCAAAGCGTGGCTg                        >  1:598199‑M2/23‑35 (MQ=255)
            tCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTg                        >  1:1463993/1‑35 (MQ=255)
            tCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTg                        >  1:1740127/1‑35 (MQ=255)
             gctaagaactatcagataattCAAAGCGTGGCTGt                       <  1:3781179‑M2/14‑1 (MQ=255)
              tACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGtt                      <  1:1021110/35‑1 (MQ=255)
              tACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGtt                      <  1:383263/35‑1 (MQ=255)
              tACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGtt                      >  1:2126446/1‑35 (MQ=255)
              tACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGtt                      <  1:7015172/35‑1 (MQ=255)
              ctaagaactatcagataattCAAAGCGTGGCTgtg                      >  1:7039226‑M2/21‑34 (MQ=255)
               aCTTTCTGCGCCAGGACAGCAAAGCGTGGCTGttt                     <  1:7337200/35‑1 (MQ=255)
                cTTTCTGCGCCAGGACAGCAAAGCGTGGCTGtttt                    <  1:708014/35‑1 (MQ=255)
                  ttCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTc                  <  1:6560648/35‑1 (MQ=255)
                   tCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTcc                 <  1:5471606/35‑1 (MQ=255)
                    cTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCt                <  1:6833944/35‑1 (MQ=255)
                    cTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCt                <  1:7572019/35‑1 (MQ=255)
                     tGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTc               >  1:4093730/1‑35 (MQ=255)
                     tGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTc               >  1:5879527/1‑35 (MQ=255)
                       cgcCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGt             >  1:6687647/1‑35 (MQ=255)
                       cgcCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGt             <  1:3322892/35‑1 (MQ=255)
                          cAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAg          <  1:4631356/35‑1 (MQ=255)
                            ggACAGCAAAGCGTGGCTGTTTTTCCTTGTCAGCa        >  1:6243997/1‑35 (MQ=255)
                            ggACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCa        >  1:6442728/1‑35 (MQ=255)
                             gACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCAc       >  1:7228911/1‑35 (MQ=255)
                             gACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCAc       <  1:458125/35‑1 (MQ=255)
                             gACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCAc       >  1:1086845/1‑35 (MQ=255)
                                  cAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCg  >  1:670390/1‑35 (MQ=255)
                                  cAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCg  >  1:7296842/1‑35 (MQ=255)
                                  |                                  
CGATGCTGGTGGTCTACTTTCTGCGCCAGGACAGCAAAGCGTGGCTGTTTTTCCTCGTCAGCACGGTCG  >  NC_000913/3979103‑3979171

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