New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 120780523 (0.550)9 (0.290) 8/96 0.9 23.6% coding (305/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 = 1209602 41 (1.310)pseudogene (54/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

TTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATG  >  NC_000913/1207741‑1207869
                                                                |                                                                
aaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAgg                                                                  >  1:3149680‑M1/50‑65 (MQ=255)
   gcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGc                                                               <  1:522455‑M1/19‑1 (MQ=255)
    gCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctct                                                              <  1:1568749‑M1/65‑5 (MQ=255)
    gCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCa                                                              >  1:2105564/1‑65 (MQ=255)
        gatCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAAc                                                          <  1:1048655/64‑1 (MQ=255)
        aatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAAc                                                          <  1:2015744‑M1/24‑1 (MQ=255)
           cAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattac                                                       <  1:106126‑M1/65‑12 (MQ=255)
             tgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAAt                                                     >  1:987563‑M1/37‑65 (MQ=255)
              aTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGGAAAACTAAATg                                                    >  1:2151042/1‑65 (MQ=255)
                  tttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAATGAGCAAAACTCAATgcagc                                                >  1:2742533‑M1/32‑65 (MQ=255)
                      gCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACg                                            <  1:2662405/65‑1 (MQ=255)
                      gCACTATGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACg                                            <  1:111590/65‑1 (MQ=255)
                      cttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACg                                            <  1:1824574‑M1/38‑1 (MQ=255)
                       cACTTTGCAAACGGCGCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                                           <  1:1412334/65‑1 (MQ=255)
                       cACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                                           <  1:1457113/65‑1 (MQ=255)
                        aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtccc                                          >  1:3032498‑M1/1‑41 (MQ=255)
                        aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGc                                          >  1:913172/1‑65 (MQ=255)
                        aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGc                                          >  1:1564354/1‑65 (MQ=255)
                        aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGc                                          >  1:125487/1‑65 (MQ=255)
                          tgtgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccgg                                        >  1:1203739‑M1/1‑39 (MQ=255)
                           tatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                       >  1:2107770‑M1/23‑65 (MQ=255)
                           tatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                       >  1:2292389‑M1/23‑65 (MQ=255)
                           tatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                       >  1:1879831‑M1/23‑65 (MQ=255)
                            tgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACa                                      <  1:1021216/65‑1 (MQ=255)
                            atttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACa                                      >  1:537673‑M1/22‑65 (MQ=255)
                             tttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACAt                                     >  1:3173815‑M1/21‑65 (MQ=255)
                             gCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttc                                     >  1:179313‑M1/1‑36 (MQ=255)
                             gCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTaaaa                                     >  1:421045/1‑62 (MQ=255)
                             gCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACAt                                     >  1:2409700/1‑65 (MQ=255)
                               tcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTa                                   <  1:118759‑M1/47‑1 (MQ=255)
                                  cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggtttcgtgg                                >  1:509289‑M1/1‑31 (MQ=255)
                                  cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgagg                                >  1:497123‑M1/1‑31 (MQ=255)
                                  cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACCTTaacg                                >  1:2636474/1‑62 (MQ=255)
                                  cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGa                                >  1:460540/1‑65 (MQ=255)
                                  cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGa                                >  1:529248/1‑65 (MQ=255)
                                  cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGa                                >  1:837945/1‑65 (MQ=255)
                                  cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAAGa                                >  1:3240771/1‑65 (MQ=255)
                                  cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACAATAGGa                                >  1:2296044/1‑65 (MQ=255)
                                  cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACAATAAGa                                >  1:2793105/1‑65 (MQ=255)
                                  cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAAATTAGGa                                >  1:391912/1‑65 (MQ=255)
                                   ggCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggc                               <  1:453020‑M1/65‑36 (MQ=255)
                                         cGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcag                         >  1:854620‑M1/1‑24 (MQ=255)
                                              aCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGCCGCACCGGTTTc                    >  1:3139671/1‑65 (MQ=255)
                                              aCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCCGTTTc                    >  1:3067359/1‑65 (MQ=255)
                                                ggtggTTTGGGAGAAGGctctgcattacctgttggtgtcccggtgccgtggccttcagccactcc                  <  1:1783374‑M1/63‑49 (MQ=255)
                                                 TTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                 >  1:57480‑M1/1‑65 (MQ=255)
                                                  tGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgc                <  1:152167‑M1/65‑51 (MQ=255)
                                                   GGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGc               >  1:2408273‑M1/1‑65 (MQ=255)
                                                   GGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGc               >  1:323407‑M1/1‑65 (MQ=255)
                                                    GTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCt              >  1:2626588‑M1/1‑65 (MQ=255)
                                                    GTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCt              >  1:1914168‑M1/1‑65 (MQ=255)
                                                     tttGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaa             >  1:2757441‑M1/1‑12 (MQ=255)
                                                     tttGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaa             >  1:1638337‑M1/1‑12 (MQ=255)
                                                     TTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCta             <  1:6312‑M1/65‑1 (MQ=255)
                                                     TTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGGTTCATAGCta             <  1:1557410‑M1/65‑1 (MQ=255)
                                                      TTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCtat            >  1:1392502‑M1/1‑65 (MQ=255)
                                                         ggAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacagg         >  1:217338‑M1/1‑8 (MQ=255)
                                                          GAGAAGGAGCAAAACTCAATGCGGCAACGGCTACATTAGGACGCCCCGGTTTCCTAGCTATACCg        <  1:1326646‑M1/65‑1 (MQ=255)
                                                          GAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCg        <  1:649152‑M1/65‑1 (MQ=255)
                                                          GAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGAAGCACCGGGTTCATAGCTATACCg        <  1:1029218‑M1/65‑1 (MQ=255)
                                                           agaAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggct       >  1:1129020‑M1/1‑6 (MQ=255)
                                                           AGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCgg       >  1:2192888‑M1/1‑65 (MQ=255)
                                                           AGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCgg       >  1:1493188‑M1/1‑65 (MQ=255)
                                                           AGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCgg       >  1:2435810‑M1/1‑65 (MQ=255)
                                                           AGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCgg       >  1:835857‑M1/1‑65 (MQ=255)
                                                            GAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGt      >  1:906622‑M1/1‑65 (MQ=255)
                                                                gctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctg  >  1:2134750‑M1/1‑1 (MQ=255)
                                                                gctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggacg  >  1:909323‑M1/1‑1 (MQ=255)
                                                                GAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCAACGGTTTCATAGCTATACCGGTTATg  <  1:1836222‑M1/65‑1 (MQ=255)
                                                                |                                                                
TTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATG  >  NC_000913/1207741‑1207869

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