New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =15 (0.770)6 (0.410) 5/100 0.8 36.0% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 1209619 = 10 (0.680)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

GGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACAT  >  NC_000913/1209552‑1209686
                                                                   |                                                                   
ggCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAAc                                                                     >  2:29871/1‑68 (MQ=255)
 gCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACg                                                                    >  2:458158/1‑68 (MQ=255)
  tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGt                                                                   >  1:167732‑M2/66‑68 (MQ=255)
     aGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTa                                                                >  1:104669/1‑68 (MQ=255)
       tGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGGGCCTTCTCCCAAGCCAACGTTTATg                                                              >  2:612575/1‑68 (MQ=255)
             aaGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                                        >  2:43144/1‑68 (MQ=255)
             aaGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                                        >  2:362556/1‑68 (MQ=255)
                        gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACa                                             >  1:133329‑M2/44‑68 (MQ=255)
                            ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCa                                         >  2:221943‑M2/40‑68 (MQ=255)
                            acgttgctgaattgagttttgctccttctcacaaaccaaCGTTTATGAAAATGAAGAAATAACAAgct                                         >  2:4955‑M2/40‑67 (MQ=255)
                                tgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATg                                     <  1:638021‑M2/33‑1 (MQ=255)
                                    aCAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcat                                 >  1:66639/1‑68 (MQ=255)
                                      attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                               >  1:300473‑M2/30‑68 (MQ=255)
                                        gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                             <  2:662910/68‑1 (MQ=255)
                                           aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                          >  2:282401/1‑68 (MQ=255)
                                             gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                        <  1:29871/68‑1 (MQ=255)
                                                   ccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAc                  >  1:204948‑M2/17‑68 (MQ=255)
                                                             aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                        >  2:103288‑M2/7‑52 (MQ=255)
                                                             aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                        <  1:103288‑M2/46‑1 (MQ=255)
                                                              accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt       >  1:409540‑M2/6‑68 (MQ=255)
                                                              accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt       <  2:311387‑M2/63‑1 (MQ=255)
                                                               ccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTa      <  1:221943‑M2/64‑1 (MQ=255)
                                                                caaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa     <  2:663976‑M2/65‑1 (MQ=255)
                                                                   cGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACAt  <  2:409540/68‑1 (MQ=255)
                                                                   |                                                                   
GGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACAT  >  NC_000913/1209552‑1209686

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