New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 12078058 (0.160)41 (1.080) 19/102 0.2 80.3% coding (305/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 = 1209602 14 (0.370)pseudogene (54/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

ATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTA  >  NC_000913/1209535‑1209669
                                                                   |                                                                   
aTTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCagag                                                                     >  1:1595684/1‑68 (MQ=255)
aTTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGGAATGCagag                                                                     >  1:2937332/1‑68 (MQ=255)
aTTTCAGCCAGCCCGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGCAATGCagag                                                                     >  1:3520377/1‑68 (MQ=255)
 tttCAGCCGGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                    >  1:1135205‑M2/1‑67 (MQ=255)
 tttCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGc                                                                    >  1:3837920‑M2/1‑67 (MQ=255)
 tttCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCAGGACACCAACAGGTAATGCAGAGc                                                                    >  1:1549514‑M2/1‑67 (MQ=255)
  ttCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGcc                                                                   >  1:3550159‑M2/1‑66 (MQ=255)
   tCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCCCGGAACCGGGACACCAACAGGGAATGCAGAGcct                                                                  <  1:1865840‑M2/68‑4 (MQ=255)
     agccagccTGTTGGCGGAGTGGCTGAAGGACCCGGAACCGGGACACAAACAGGAAATGCAGAGccttc                                                                <  1:3060159‑M2/68‑6 (MQ=255)
      gccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttct                                                               >  1:1039885‑M2/1‑62 (MQ=255)
       ccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctc                                                              <  1:4042475‑M2/68‑8 (MQ=255)
       ccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGGAATGCAGAGccttctc                                                              <  1:1085320‑M2/68‑8 (MQ=255)
           ccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaa                                                          <  1:2513290‑M2/68‑12 (MQ=255)
            cTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaa                                                         <  1:3118333‑M2/68‑13 (MQ=255)
                  gCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaagg                                                   <  1:3032003‑M2/68‑19 (MQ=255)
                   gggaGTGGCTGAAGGCCACGGATCCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggt                                                  >  1:3075103‑M2/2‑49 (MQ=255)
                   cGGAGTGGCTGGAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggt                                                  >  1:988917‑M2/1‑49 (MQ=255)
                   cGGAGTGGCTGAATGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggt                                                  >  1:3944828‑M2/1‑49 (MQ=255)
                    ggAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggtt                                                 <  1:463006‑M2/68‑21 (MQ=255)
                      aGTGGCTGAAGTCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggtttt                                               >  1:2893858‑M2/1‑46 (MQ=255)
                        tGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcg                                             >  1:1484473‑M2/1‑44 (MQ=255)
                            tGAAGGCCACGGAACCGGGACACCAAAAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGaaaa                                         >  1:4020387/1‑68 (MQ=255)
                             gAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagc                                        >  1:1874743‑M2/1‑39 (MQ=255)
                              aaGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagcc                                       <  1:3427944‑M2/68‑31 (MQ=255)
                              aaGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagcc                                       <  1:281461‑M2/68‑31 (MQ=255)
                               aGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                      >  1:372667‑M2/1‑37 (MQ=255)
                               aGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                      >  1:2995277‑M2/1‑37 (MQ=255)
                               aGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                      >  1:512372‑M2/1‑37 (MQ=255)
                               aGGCCACGGAACCGGGACACCAACAGGGAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                      >  1:3424776‑M2/1‑37 (MQ=255)
                                 gCCACGGAACCGGGACACCAACAGGTAATGCAGCGccttctcccaaaccaaggttttcgagagccgtt                                    <  1:1605418‑M2/68‑34 (MQ=255)
                                 gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtt                                    <  1:1298081‑M2/68‑34 (MQ=255)
                                 gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtt                                    <  1:2297035‑M2/68‑34 (MQ=255)
                                 gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtt                                    <  1:1829860‑M2/68‑34 (MQ=255)
                                 gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtt                                    <  1:1158579‑M2/68‑34 (MQ=255)
                                 gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaag                                    >  1:3878709/1‑68 (MQ=255)
                                 gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaag                                    >  1:853225/1‑68 (MQ=255)
                                       gAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcaca                              <  1:267576‑M2/68‑40 (MQ=255)
                                        aaCCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacag                             >  1:1182852‑M2/1‑28 (MQ=255)
                                         cccgGGGCACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                            >  1:1607599‑M2/2‑27 (MQ=255)
                                         cccgGGACACCAACAGGTAATGCAGAGcctgctcccaaaccaaggttttcgagagccgtttgcacagt                            >  1:2628921‑M2/2‑27 (MQ=255)
                                         aCCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                            >  1:633008‑M2/1‑27 (MQ=255)
                                         aCCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                            >  1:4077793‑M2/1‑27 (MQ=255)
                                         aCCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                            <  1:2588345‑M2/68‑42 (MQ=255)
                                         aCCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                            >  1:2580540‑M2/1‑27 (MQ=255)
                                         aCCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                            >  1:3419152‑M2/1‑27 (MQ=255)
                                         aCCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                            <  1:153680‑M2/68‑42 (MQ=255)
                                         aCCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACa                            >  1:59235/1‑68 (MQ=255)
                                              gACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCaa                       >  1:1105303/1‑68 (MQ=255)
                                               acacCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCaaa                      <  1:1120779/68‑1 (MQ=255)
                                                cacCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAAt                     >  1:2552484/1‑68 (MQ=255)
                                                 acCAACAGGTAATGCAGAGCCTTCTGCCCAACCAACGTTTATGAAAATGAAGAAATAACAAGCAACTg                    >  1:463261/1‑68 (MQ=255)
                                                    aaCAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                 >  1:2864027/1‑68 (MQ=255)
                                                       aGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttga              <  1:3821807‑M2/68‑56 (MQ=255)
                                                         gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgata            <  1:1308139‑M2/68‑58 (MQ=255)
                                                         gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgata            <  1:1169121‑M2/68‑58 (MQ=255)
                                                         gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt            >  1:1670929/1‑68 (MQ=255)
                                                            aTGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcg         >  1:1890144‑M2/1‑8 (MQ=255)
                                                             tGCAGAGCCTTCTCCCAAACCAACGTTTATGAAACTGAAGAAATAACAAGCAAATGGCATCATTCCTg        >  1:2707404/1‑68 (MQ=255)
                                                             tGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg        >  1:1952180/1‑68 (MQ=255)
                                                               cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgcca      <  1:1268817‑M2/68‑64 (MQ=255)
                                                               cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgcca      <  1:1427828‑M2/68‑64 (MQ=255)
                                                               cAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgcca      <  1:530161‑M2/68‑64 (MQ=255)
                                                                agagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatattgccaa     <  1:2887769‑M2/68‑65 (MQ=255)
                                                                agagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaa     <  1:2603127‑M2/68‑65 (MQ=255)
                                                                agagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt     >  1:7095/1‑68 (MQ=255)
                                                                   gccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacg  <  1:2410349‑M2/68‑68 (MQ=255)
                                                                   gccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacg  <  1:2149779‑M2/68‑68 (MQ=255)
                                                                   gccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacg  <  1:1702277‑M2/68‑68 (MQ=255)
                                                                   gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa  <  1:3715096/68‑1 (MQ=255)
                                                                   |                                                                   
ATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTA  >  NC_000913/1209535‑1209669

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