New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 120780525 (0.780)8 (0.330) 8/100 1.0 32.8% coding (305/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 = 1209602 14 (0.580)pseudogene (18/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

ATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTA  >  NC_000913/1209535‑1209669
                                                                   |                                                                   
aTTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACCCCAACAGGTAATGCagag                                                                     >  1:259973/1‑68 (MQ=255)
aTTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCagag                                                                     >  2:1117227/1‑68 (MQ=255)
 tttCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                    >  1:1117229‑M2/1‑67 (MQ=255)
 tttCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                    >  1:929837‑M2/1‑67 (MQ=255)
     agccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttc                                                                <  1:113564‑M2/68‑6 (MQ=255)
               ttGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaacca                                                      <  1:112880‑M2/68‑16 (MQ=255)
                 ggCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaag                                                    <  1:734800‑M2/68‑18 (MQ=255)
                               aGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                      >  2:187729‑M2/1‑37 (MQ=255)
                                 gcaaCGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaag                                    >  1:102046/4‑68 (MQ=255)
                                    aCGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgc                                 <  1:804416‑M2/68‑37 (MQ=255)
                                     cGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAAt                                <  2:391623/68‑1 (MQ=255)
                                     cGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTGTATGAAAATGAAGAAAt                                <  2:625757/68‑1 (MQ=255)
                                         aCCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                            >  1:444574‑M2/1‑27 (MQ=255)
                                                cacCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAAt                     <  2:259973/68‑1 (MQ=255)
                                                 acCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATg                    >  1:937412/1‑68 (MQ=255)
                                                    aaCAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatt                 <  2:79671‑M2/68‑53 (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)
                                                              gCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgcc       <  2:252371‑M2/68‑63 (MQ=255)
                                                              gCAGAGccttctcccaaaccaaggctttcgagagccgtttgcacagtgccatccgatttgatatcgcc       >  2:130682‑M2/1‑6 (MQ=255)
                                                              gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc       >  1:1173259/1‑68 (MQ=255)
                                                               cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCt      >  2:6988/1‑68 (MQ=255)
                                                                agagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt     >  1:785548/1‑68 (MQ=255)
                                                                agagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt     >  1:636551/1‑68 (MQ=255)
                                                                   gccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacg  <  1:226293‑M2/68‑68 (MQ=255)
                                                                   gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa  >  2:1144778/1‑68 (MQ=255)
                                                                   |                                                                   
ATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTA  >  NC_000913/1209535‑1209669

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.