New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 12078053 (0.100)13 (0.540) 13/128 0.7 80.3% coding (305/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 = 1209602 4 (0.170)pseudogene (54/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

GCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACC  >  NC_000913/1209522‑1209671
                                                                                |                                                                     
gCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCagag                                                                       <  2:221260/81‑1 (MQ=255)
 cGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGTTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                      >  2:221263‑M2/1‑80 (MQ=255)
 cGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                      <  1:221264‑M2/81‑2 (MQ=255)
 cGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                      >  2:221262‑M2/1‑80 (MQ=255)
 cGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                      >  2:221261‑M2/1‑80 (MQ=255)
    cACCATCGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGcctt                                                                   <  2:221265‑M2/81‑5 (MQ=255)
      ccATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttct                                                                 >  2:221266‑M2/1‑75 (MQ=255)
        aTTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcc                                                               >  1:221267‑M2/1‑73 (MQ=255)
         ttGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctccc                                                              >  2:221268‑M2/1‑72 (MQ=255)
          tGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctccc                                                              <  2:221269‑M2/80‑10 (MQ=255)
            cATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaa                                                           >  1:221270‑M2/1‑69 (MQ=255)
               ttCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaacc                                                         <  2:221271‑M2/80‑15 (MQ=255)
                   gccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAAAGCAGAGccttctcccaaaccaaggt                                                    >  2:221272‑M2/1‑62 (MQ=255)
                    ccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggtt                                                   <  1:221273‑M2/81‑21 (MQ=255)
                            ttGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgag                                             >  1:221274‑M2/1‑53 (MQ=255)
                             tGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagc                                          <  1:221075‑M2/81‑30 (MQ=255)
                               gCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagcc                                         >  1:221275‑M2/1‑50 (MQ=255)
                                       gCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATaa                                >  2:221290/1‑81 (MQ=255)
                                          gAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAAc                               >  2:221294/1‑79 (MQ=255)
                                                 aCGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccg                      >  2:221277‑M2/1‑32 (MQ=255)
                                                 aCGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATg                      <  2:221276/81‑1 (MQ=255)
                                                  cGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATgg                     <  1:221278/81‑1 (MQ=255)
                                                   ggAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgat                    <  2:221074‑M2/81‑52 (MQ=255)
                                                    gAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatt                   <  1:221073‑M2/81‑53 (MQ=255)
                                                      aCCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgattt                  <  2:221072‑M2/80‑54 (MQ=255)
                                                            acacCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcg           <  1:221071‑M2/81‑61 (MQ=255)
                                                             cacCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgc          >  1:221279‑M2/1‑20 (MQ=255)
                                                                cAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaa       >  1:221280‑M2/1‑17 (MQ=255)
                                                                     ggTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccgtccgatttgatatcgccaaacgga  >  2:221281‑M2/1‑12 (MQ=255)
                                                                                |                                                                     
GCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACC  >  NC_000913/1209522‑1209671

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