New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =25 (0.780)12 (0.500) 9/100 0.8 42.2% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 1209619 = 14 (0.580)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

CGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTA  >  NC_000913/1209554‑1209682
                                                                 |                                                               
tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGt                                                               >  1:981976‑M2/66‑68 (MQ=255)
           tgcgtcctaatgtagccgatgctgcattgagttttgctccttctcacaaaccaaCGTTTATGAAAATg                                                    >  2:136758‑M2/55‑68 (MQ=255)
              gcaaCGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaag                                                 >  1:102046/4‑68 (MQ=255)
                  cGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAAt                                             <  2:391623/68‑1 (MQ=255)
                  cGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTGTATGAAAATGAAGAAAt                                             <  2:625757/68‑1 (MQ=255)
                     tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAc                                          >  2:1049365‑M2/45‑68 (MQ=255)
                     tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAc                                          >  2:1032169‑M2/45‑68 (MQ=255)
                      gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACa                                         >  2:944131‑M2/44‑68 (MQ=255)
                      gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACa                                         >  2:1181025‑M2/44‑68 (MQ=255)
                      gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACa                                         >  1:526511‑M2/44‑68 (MQ=255)
                             cacCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAAt                                  <  2:259973/68‑1 (MQ=255)
                              acCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATg                                 >  1:937412/1‑68 (MQ=255)
                               gctgcattgagttttgctccttctcccaaaccaaCg                                                                >  1:626024‑M2/35‑36 (MQ=255)
                               gctgcattgagttttgctccttctcccaaaccaaCg                                                                <  2:626024‑M2/2‑1 (MQ=255)
                                ctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGc                               <  1:1049366‑M2/35‑1 (MQ=255)
                                 tgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                              <  1:1032170‑M2/36‑1 (MQ=255)
                                    aGGTAATGCAGAGCATTCTACCAAAACAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                           >  1:1180065/1‑68 (MQ=255)
                                      gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                         >  1:782531/1‑68 (MQ=255)
                                        aaTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                       >  1:382567/1‑68 (MQ=255)
                                          tGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg                     >  1:529851/1‑68 (MQ=255)
                                           gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                    >  1:1173259/1‑68 (MQ=255)
                                            cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCt                   >  2:6988/1‑68 (MQ=255)
                                             agagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                  >  1:636551/1‑68 (MQ=255)
                                             agagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                  >  1:785548/1‑68 (MQ=255)
                                                tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa               <  1:1181027‑M2/51‑1 (MQ=255)
                                                gacTTATCCCAGACATACGTATATGAAAATGAAGAAATAACAAGCTAATGGCCTCATTCCTGCtattc               >  1:483910/3‑67 (MQ=255)
                                                gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa               >  2:1144778/1‑68 (MQ=255)
                                                cccttctcccaaaccaaCGTTTATGAAAATGATGAAATAACAAGCAAATTGCATCATTCCTGCTTTTa               <  2:54956‑M2/51‑1 (MQ=255)
                                                       cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAACAgggg        >  2:265139‑M2/11‑68 (MQ=255)
                                                        ccaaaccaaCGTATATTAAAATGAAGAAATAACAAGCAAATTGCATCATTCCTGCTTTTACCaggggt       <  1:670513‑M2/59‑2 (MQ=255)
                                                         caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa      <  1:224087‑M2/60‑1 (MQ=255)
                                                         caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa      <  2:730488‑M2/60‑1 (MQ=255)
                                                         caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa      <  1:693399‑M2/60‑1 (MQ=255)
                                                            accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt   >  2:1107382‑M2/6‑68 (MQ=255)
                                                             ccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTa  <  2:529851‑M2/64‑1 (MQ=255)
                                                                 |                                                               
CGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTA  >  NC_000913/1209554‑1209682

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.