New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =15 (0.380)25 (0.840) 17/96 0.2 55.4% coding (290/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 1209619 = 29 (0.970)pseudogene (37/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

AGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATA  >  NC_000913/1207726‑1207854
                                                                |                                                                
aGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTTGAAAACCT                                                                  >  1:3057666‑M1/1‑65 (MQ=255)
aGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCT                                                                  >  1:2662576‑M1/1‑65 (MQ=255)
 gCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTT                                                                 >  1:2740477‑M1/1‑65 (MQ=255)
 gCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTT                                                                 >  1:569329‑M1/1‑65 (MQ=255)
  ccGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTG                                                                <  1:1185816‑M1/65‑1 (MQ=255)
   cGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGG                                                               <  1:1214035‑M1/65‑1 (MQ=255)
    tccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGt                                                              <  1:624221‑M1/5‑1 (MQ=255)
    tccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGt                                                              <  1:3080814‑M1/5‑1 (MQ=255)
    tccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGt                                                              <  1:1322701‑M1/5‑1 (MQ=255)
    tccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGt                                                              <  1:2576652‑M1/5‑1 (MQ=255)
     ccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGtt                                                             <  1:380161‑M1/6‑1 (MQ=255)
     ccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGtt                                                             <  1:2014775‑M1/6‑1 (MQ=255)
     ccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGtt                                                             <  1:1808922‑M1/6‑1 (MQ=255)
     ccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGtt                                                             <  1:990062‑M1/6‑1 (MQ=255)
        gAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGG                                                          <  1:834694‑M1/65‑1 (MQ=255)
        cctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTgg                                                          >  1:111741‑M1/57‑65 (MQ=255)
         ctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTggg                                                         >  1:1996062‑M1/56‑65 (MQ=255)
         aaTCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGG                                                         >  1:3170011‑M1/1‑65 (MQ=255)
             taaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGaa                                                     <  1:424401‑M1/14‑1 (MQ=255)
             taaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGaa                                                     <  1:2109020‑M1/14‑1 (MQ=255)
                aagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGa                                                  >  1:36750‑M1/49‑65 (MQ=255)
                 agcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAg                                                 <  1:2481624‑M1/18‑1 (MQ=255)
                  gcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGaga                                                <  1:1701186‑M1/19‑2 (MQ=255)
                  gcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGc                                                <  1:2708353‑M1/19‑1 (MQ=255)
                  gcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGc                                                <  1:445434‑M1/19‑1 (MQ=255)
                  gcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGc                                                <  1:709043‑M1/19‑1 (MQ=255)
                  gcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGc                                                >  1:737452‑M1/47‑65 (MQ=255)
                    ggaTATTAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctg                                              <  1:3290920‑M1/64‑6 (MQ=255)
                    cGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctg                                              <  1:916528‑M1/65‑6 (MQ=255)
                     ggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCaaa                                             >  1:1028671‑M1/44‑65 (MQ=255)
                     gATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgc                                             <  1:2936105‑M1/65‑7 (MQ=255)
                     gATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgc                                             <  1:1047752‑M1/65‑7 (MQ=255)
                     gATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgc                                             <  1:1087228‑M1/65‑7 (MQ=255)
                       aatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAAc                                           <  1:206429‑M1/24‑1 (MQ=255)
                       aatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAAc                                           >  1:1908641‑M1/42‑65 (MQ=255)
                       aatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAAc                                           >  1:744316‑M1/42‑65 (MQ=255)
                          cAAATCGGATGGCACTGTGCAAACTGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattac                                        <  1:2812204‑M1/65‑12 (MQ=255)
                          cAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattac                                        <  1:891717‑M1/65‑12 (MQ=255)
                          cAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattac                                        <  1:789946‑M1/65‑12 (MQ=255)
                            tgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAAt                                      >  1:278468‑M1/37‑65 (MQ=255)
                              tCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgt                                    <  1:2479570‑M1/65‑16 (MQ=255)
                              tCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgc                                    >  1:2311895/1‑65 (MQ=255)
                                atttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgcag                                  >  1:342763‑M1/33‑65 (MQ=255)
                                 tttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgcagc                                 <  1:2053806‑M1/34‑1 (MQ=255)
                                     gCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtc                             <  1:2530081‑M1/65‑23 (MQ=255)
                                     gCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtc                             <  1:1437485‑M1/65‑23 (MQ=255)
                                     cttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACg                             >  1:1192392‑M1/28‑65 (MQ=255)
                                      ttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                            >  1:1365462‑M1/27‑65 (MQ=255)
                                      ttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                            >  1:236504‑M1/27‑65 (MQ=255)
                                      ttgttatttcttcattctcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                            >  1:3171871‑M1/27‑65 (MQ=255)
                                      cACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                            <  1:2677424/65‑1 (MQ=255)
                                      cACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                            <  1:353260/65‑1 (MQ=255)
                                       tgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGc                           <  1:805639‑M1/40‑1 (MQ=255)
                                       aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctttgcattaccggttggggtccc                           >  1:3206595‑M1/1‑41 (MQ=255)
                                       aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgccttacctgttggtgtccc                           >  1:3269393‑M1/1‑41 (MQ=255)
                                       aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtccc                           >  1:1863204‑M1/1‑41 (MQ=255)
                                       aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtccc                           >  1:1308768‑M1/1‑41 (MQ=255)
                                       aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtccc                           >  1:3171676‑M1/1‑41 (MQ=255)
                                       aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcactacctgttggtgtccc                           >  1:2012939‑M1/1‑41 (MQ=255)
                                       aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGc                           >  1:1721801/1‑65 (MQ=255)
                                       aCTGTGCAAACGGCTCTCGAAAACCTTGGCTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGc                           >  1:91358/1‑65 (MQ=255)
                                          tatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                        >  1:1191621‑M1/23‑65 (MQ=255)
                                          tatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                        >  1:2367687‑M1/23‑65 (MQ=255)
                                          gtgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                        <  1:1055822/65‑1 (MQ=255)
                                           atttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACa                       >  1:1377319‑M1/22‑65 (MQ=255)
                                            gCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttc                      >  1:1266378‑M1/1‑36 (MQ=255)
                                            gCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAAAGCAGCAACGGCTACAt                      >  1:2947881/1‑65 (MQ=255)
                                             cAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttcc                     <  1:1907652‑M1/65‑31 (MQ=255)
                                               aaCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACAGTac                   >  1:2691372/1‑64 (MQ=255)
                                                 cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtgg                 >  1:3265621‑M1/1‑31 (MQ=255)
                                                 cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtgg                 >  1:1644113‑M1/1‑31 (MQ=255)
                                                 cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggtcccgtgg                 >  1:3186706‑M1/1‑31 (MQ=255)
                                                 cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGa                 >  1:2077251/1‑65 (MQ=255)
                                                 cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGa                 >  1:3216228/1‑65 (MQ=255)
                                                 cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGa                 >  1:1969599/1‑65 (MQ=255)
                                                 cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGa                 >  1:16375/1‑65 (MQ=255)
                                                  ggCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggc                <  1:2279338‑M1/65‑36 (MQ=255)
                                                  ggCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggc                <  1:805851‑M1/65‑36 (MQ=255)
                                                  ggCACTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggc                <  1:116076‑M1/65‑36 (MQ=255)
                                                     tttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCa             >  1:2097070‑M1/12‑65 (MQ=255)
                                                     tctcGAAAACCTTAGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCa             <  1:2906893/65‑1 (MQ=255)
                                                          taaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGt        >  1:3230902‑M1/7‑65 (MQ=255)
                                                          aaaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTATATTAGGACGCACCGGt        >  1:2002696‑M1/7‑65 (MQ=255)
                                                          aaaaCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGt        <  1:1639589/65‑1 (MQ=255)
                                                          aaaaCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGt        <  1:22484/65‑1 (MQ=255)
                                                             aCCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagcctgt     >  1:297922‑M1/1‑19 (MQ=255)
                                                                TTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa  >  1:2222959‑M1/1‑65 (MQ=255)
                                                                TTGGTGTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa  >  1:1563739‑M1/1‑65 (MQ=255)
                                                                |                                                                
AGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATA  >  NC_000913/1207726‑1207854

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