New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =14 (0.420)4 (0.150) 3/128 2.7 25.7% coding (290/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 1209619 = 12 (0.460)pseudogene (37/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

CAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTAT  >  NC_000913/1209539‑1209699
                                                                                |                                                                                
cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAACGCAGAGCCTTCTCCCAAACCAAc                                                                                  >  2:239892/1‑81 (MQ=255)
  gccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAACGCAGAGCCTTCTCCCAAACCAACGt                                                                                >  1:239882/1‑81 (MQ=255)
      gccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAACGCAGAGCCTTCTCCCAAACCAACGTTTAt                                                                            >  2:239880/1‑81 (MQ=255)
          atagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaa                                                                         >  2:239889‑M2/71‑80 (MQ=255)
           ttGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAACGCAGAGCCTTCTCCCAAACCAACGTTTATGaaaa                                                                       <  2:239864/81‑1 (MQ=255)
                    tGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAACGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAAt                                                              <  1:239865/81‑1 (MQ=255)
                         gcgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAc                                                           >  2:239881‑M2/56‑79 (MQ=255)
                              ccACGGAACCGGGACACCAACAGGTAACGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCaaa                                                    >  1:239890/1‑81 (MQ=255)
                               cACGGAACCGGGACACCAACAGGTAACGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAAt                                                   >  2:239902/1‑81 (MQ=255)
                                        gggACACCAACAGGTAACGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                                          >  2:239900/1‑81 (MQ=255)
                                          gACACCAACAGGTAACGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                                        <  2:239867/81‑1 (MQ=255)
                                           gttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                                       >  1:239891‑M2/38‑81 (MQ=255)
                                            cacCAACAGGTAACGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg                                      >  2:239893/1‑81 (MQ=255)
                                               ctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                                   >  2:239896‑M2/34‑81 (MQ=255)
                                                 aCAGGTAACGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtttt                                 <  1:239869/81‑1 (MQ=255)
                                                            agagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAt                      <  1:239870/81‑1 (MQ=255)
                                                                 cttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAAc                 <  1:239871‑M2/66‑1 (MQ=255)
                                                                      cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCt            >  2:239894‑M2/11‑81 (MQ=255)
                                                                         aaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAt         >  2:239917‑M2/8‑81 (MQ=255)
                                                                           accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTg       <  1:239872‑M2/76‑1 (MQ=255)
                                                                                cGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTAt  <  2:239873/81‑1 (MQ=255)
                                                                                cGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTAt  >  2:239899/1‑81 (MQ=255)
                                                                                |                                                                                
CAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTAT  >  NC_000913/1209539‑1209699

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