New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =30 (1.060)16 (0.750) 13/100 0.3 48.8% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 1209619 = 11 (0.510)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

GCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATT  >  NC_000913/1209553‑1209680
                                                                  |                                                             
gCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACg                                                              <  1:664880/68‑1 (MQ=255)
ctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCg                                                              <  2:566166‑M2/2‑1 (MQ=255)
 tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaac                                                               >  1:312586‑M2/66‑66 (MQ=255)
 tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaac                                                               <  2:312586‑M2/1‑1 (MQ=255)
 tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGt                                                             >  1:364878‑M2/66‑68 (MQ=255)
 tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGt                                                             >  2:391783‑M2/66‑68 (MQ=255)
 tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGt                                                             >  1:665763‑M2/66‑68 (MQ=255)
   gAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGttt                                                           <  2:188799/68‑1 (MQ=255)
          ggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaaa                                                    <  1:242841‑M2/12‑1 (MQ=255)
             agtccACGGAACCGTGACACCATCAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATga                                                 <  2:176879/65‑1 (MQ=255)
              ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaa                                                >  1:432005/1‑68 (MQ=255)
               gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaag                                               >  2:208614/1‑68 (MQ=255)
                tcctaatgtagccgtttctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaaga                                              <  2:665763‑M2/18‑1 (MQ=255)
                      tgtagccgttgctgcattgagttttgctccttctcccacaccaaCGTTTATGAAAATGAAGAAATAAc                                        >  1:62502‑M2/45‑68 (MQ=255)
                         agccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAg                                     >  2:11637‑M2/42‑68 (MQ=255)
                           ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCa                                   >  2:313637‑M2/40‑68 (MQ=255)
                            gACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCaa                                  >  1:488196/1‑68 (MQ=255)
                            cgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCaa                                  >  2:899982‑M2/39‑68 (MQ=255)
                                 ctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGc                             <  2:290212‑M2/35‑1 (MQ=255)
                                  tgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACa                                       <  2:473440‑M2/25‑1 (MQ=255)
                                  tgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACa                                       >  1:473440‑M2/33‑57 (MQ=255)
                                     aGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                         >  1:997163/1‑68 (MQ=255)
                                       gTAATGCAGAGCCTTCTCCCAATCCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                       <  1:837368/68‑1 (MQ=255)
                                           tGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg                   >  1:734483/1‑68 (MQ=255)
                                            tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                  <  1:271007‑M2/46‑1 (MQ=255)
                                            gCAGAGCCGTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                  <  2:591702/68‑1 (MQ=255)
                                             cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCt                 >  2:929712/1‑68 (MQ=255)
                                                         ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg     <  2:337518‑M2/59‑1 (MQ=255)
                                                         ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg     <  2:200765‑M2/59‑1 (MQ=255)
                                                          caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa    <  1:836278‑M2/60‑1 (MQ=255)
                                                          caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa    <  1:480169‑M2/60‑1 (MQ=255)
                                                            aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                  <  1:94891‑M2/46‑1 (MQ=255)
                                                            aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                  >  2:94891‑M2/7‑52 (MQ=255)
                                                            aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAtt  <  2:62502‑M2/62‑1 (MQ=255)
                                                            aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAtt  >  1:437318‑M2/7‑68 (MQ=255)
                                                                  |                                                             
GCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATT  >  NC_000913/1209553‑1209680

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.