New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =2 (0.080)17 (0.870) 11/100 0.4 90.6% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 1209619 = 2 (0.100)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

GGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTT  >  NC_000913/1209552‑1209681
                                                                   |                                                              
gctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaac                                                                <  1:79323‑M2/1‑1 (MQ=255)
gctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaac                                                                >  1:89645‑M2/68‑68 (MQ=255)
gctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaac                                                                <  1:166608‑M2/1‑1 (MQ=255)
  tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGt                                                              >  1:716228‑M2/66‑68 (MQ=255)
       aaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATg                                                         >  1:522512‑M2/61‑68 (MQ=255)
       aaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATg                                                         >  1:559532‑M2/61‑68 (MQ=255)
        accggtgcgtcctagtgtagccgttgctgcattgagttttgctccttctcccaaaccagCGTATATGa                                                        >  2:271449‑M2/60‑68 (MQ=255)
        accggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGa                                                        >  2:959556‑M2/60‑68 (MQ=255)
        accggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGa                                                        >  2:382300‑M2/60‑68 (MQ=255)
                       tgtagccgttgctgcattgagtgttgctccttctcccaagccaaCGTTTATGAAAATGAAGAAATAAc                                         >  2:363405‑M2/45‑68 (MQ=255)
                        gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACa                                        >  2:888533‑M2/44‑68 (MQ=255)
                        gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACa                                        >  2:364059‑M2/44‑68 (MQ=255)
                             cgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAagcta                                   >  1:760305‑M2/39‑66 (MQ=255)
                                  ctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGc                              <  2:774351‑M2/35‑1 (MQ=255)
                                   aaCAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                             <  2:278217/68‑1 (MQ=255)
                                           gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                     >  2:779485‑M2/25‑68 (MQ=255)
                                             tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                   >  1:791000‑M2/23‑68 (MQ=255)
                                             tttgctccttctcccaaaccaaCGTATATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                   <  1:413091‑M2/46‑1 (MQ=255)
                                                  gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa              >  2:763506/1‑68 (MQ=255)
                                                           caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa     <  2:444747‑M2/60‑1 (MQ=255)
                                                           caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa     <  2:522512‑M2/60‑1 (MQ=255)
                                                              accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt  <  1:888536‑M2/63‑1 (MQ=255)
                                                              accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt  >  2:857200‑M2/6‑68 (MQ=255)
                                                                   |                                                              
GGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTT  >  NC_000913/1209552‑1209681

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