New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =23 (0.550)11 (0.350) 9/96 0.8 27.4% coding (290/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 1209619 = 41 (1.310)pseudogene (37/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

AGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATA  >  NC_000913/1207726‑1207854
                                                                |                                                                
aGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCT                                                                  >  1:2851237‑M1/1‑65 (MQ=255)
aGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCT                                                                  >  1:2912637‑M1/1‑65 (MQ=255)
aGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCT                                                                  >  1:761420‑M1/1‑65 (MQ=255)
aGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCT                                                                  >  1:3153075‑M1/1‑65 (MQ=255)
 gCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTT                                                                 >  1:139090‑M1/1‑65 (MQ=255)
    tccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGt                                                              <  1:2326378‑M1/5‑1 (MQ=255)
    tccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGt                                                              <  1:2777312‑M1/5‑1 (MQ=255)
    tccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGt                                                              <  1:2966024‑M1/5‑1 (MQ=255)
    tccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGt                                                              >  1:768582‑M1/61‑65 (MQ=255)
    tccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGt                                                              <  1:1778133‑M1/5‑1 (MQ=255)
    tccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGt                                                              <  1:68454‑M1/5‑1 (MQ=255)
     ccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGtt                                                             <  1:1935997‑M1/6‑1 (MQ=255)
      cccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGttt                                                            <  1:1027654‑M1/7‑1 (MQ=255)
        cctggtaaaagcaggaatgatgccatttgctcgttatttcttcattttcataaacgTTGCTTTgg                                                          <  1:157798‑M1/9‑1 (MQ=255)
         cttgtaaaagcaggaatggtggcatttgcttgttatttcttcattttcataaacgTTGGTTTggg                                                         >  1:3140007‑M1/56‑65 (MQ=255)
         ctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTggg                                                         >  1:1581312‑M1/56‑65 (MQ=255)
         ctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTggg                                                         >  1:1395896‑M1/56‑65 (MQ=255)
         ctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTggg                                                         >  1:1923531‑M1/56‑65 (MQ=255)
         ctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTggg                                                         >  1:2257295‑M1/56‑65 (MQ=255)
         ctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTggg                                                         >  1:2893735‑M1/56‑65 (MQ=255)
           tCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAG                                                       >  1:84215‑M1/1‑65 (MQ=255)
             taaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGaa                                                     <  1:118040‑M1/14‑1 (MQ=255)
              aaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAg                                                    >  1:877740‑M1/51‑65 (MQ=255)
               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:2292389‑M1/23‑65 (MQ=255)
                                          tatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                        >  1:2107770‑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:529248/1‑65 (MQ=255)
                                                 cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGa                 >  1:460540/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)
                                                                TTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa  >  1:57480‑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 < 15 ≤ ATCG/ATCG < 23 ≤ 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.