New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =2 (0.060)15 (0.520) 14/128 1.0 92.0% coding (290/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 1209619 = 1 (0.030)pseudogene (37/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

GCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTAT  >  NC_000913/1207714‑1207868
                                                                            |                                                                              
gcgcTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGT                                                                            <  2:231251‑M1/81‑1 (MQ=255)
caataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGt                                                                            <  2:231526‑M1/5‑1 (MQ=255)
  ataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGttt                                                                          >  1:231546‑M1/75‑81 (MQ=255)
  ataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGttt                                                                          <  1:231525‑M1/7‑1 (MQ=255)
   ctgctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTG                                                                         <  2:231252‑M1/81‑1 (MQ=255)
      ctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGG                                                                        >  2:231482‑M1/1‑79 (MQ=255)
        gTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAG                                                                     >  1:231461‑M1/1‑80 (MQ=255)
               cGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctct                                                             <  2:231253‑M1/81‑5 (MQ=255)
                    gAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcatt                                                        <  2:231254‑M1/81‑10 (MQ=255)
                    gAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcatt                                                        <  2:231255‑M1/81‑10 (MQ=255)
                     ctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACt                                                        <  1:231272‑M1/25‑1 (MQ=255)
                         cGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAAt                                                    <  1:231256/80‑1 (MQ=255)
                          aaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGGGAAGGAGCAAAACTCAATg                                                   <  2:231271‑M1/30‑1 (MQ=255)
                          aaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATg                                                   <  2:231274‑M1/30‑1 (MQ=255)
                           tttGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgtt                                                 >  1:231524‑M1/1‑65 (MQ=255)
                           aaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATg                                                   <  2:231280‑M1/30‑1 (MQ=255)
                              ggCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggt                                              >  1:231523‑M1/1‑62 (MQ=255)
                                  gaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                                          <  1:231273‑M1/39‑1 (MQ=255)
                                  atatCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcc                                          >  2:231522‑M1/1‑58 (MQ=255)
                                  atatCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcc                                          <  1:231257‑M1/81‑24 (MQ=255)
                                     tgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTa                                       >  2:231549‑M1/40‑81 (MQ=255)
                                     tCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccg                                        <  2:231259‑M1/80‑26 (MQ=255)
                                     tCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccg                                        <  2:231258‑M1/80‑26 (MQ=255)
                                       atgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                      >  1:231547‑M1/38‑80 (MQ=255)
                                       aaaTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggtt                                     <  1:231261‑M1/81‑29 (MQ=255)
                                       aaaTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggtt                                     <  2:231260‑M1/81‑29 (MQ=255)
                                         aTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACAtt                                   <  2:231262/81‑1 (MQ=255)
                                           cGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgt                                 >  2:231521‑M1/1‑49 (MQ=255)
                                            ggATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgt                                 <  1:231263‑M1/80‑33 (MQ=255)
                                             tttgcttgttatttcttcattatcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGa                               <  2:231277‑M1/50‑1 (MQ=255)
                                                 cttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCa                           >  1:231544‑M1/28‑81 (MQ=255)
                                                  cACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttc                          <  1:231265‑M1/81‑40 (MQ=255)
                                                  cACTGTGCAAACGGCTCTCGAAAACCTTGGTATGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttc                          <  2:231264‑M1/81‑40 (MQ=255)
                                                     ttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCg                        >  1:231541‑M1/24‑80 (MQ=255)
                                                       tgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagcc                      <  1:231266‑M1/80‑44 (MQ=255)
                                                         cAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccac                    <  1:231267‑M1/80‑46 (MQ=255)
                                                          tcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTc                   >  1:231537‑M1/19‑80 (MQ=255)
                                                          aaaCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactc                  >  2:231519‑M1/1‑34 (MQ=255)
                                                            ttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCa                  <  1:231276‑M1/63‑1 (MQ=255)
                                                            aCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccg                >  1:231518‑M1/1‑32 (MQ=255)
                                                             cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgc               <  1:231268‑M1/81‑51 (MQ=255)
                                                             cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgc               >  1:231517‑M1/1‑31 (MQ=255)
                                                              cattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                >  1:231543‑M1/15‑79 (MQ=255)
                                                                ctctCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaa            <  1:231269‑M1/81‑54 (MQ=255)
                                                                   tcGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaac           <  2:231516‑M1/79‑55 (MQ=255)
                                                                     gAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggc       >  2:231515‑M1/1‑23 (MQ=255)
                                                                     ataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATAcc        >  2:231536‑M1/8‑80 (MQ=255)
                                                                      aaaaCCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggct      <  1:231514‑M1/81‑60 (MQ=255)
                                                                       aaaCCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctg     >  1:231513‑M1/1‑21 (MQ=255)
                                                                        aacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGtt    >  1:231539‑M1/5‑81 (MQ=255)
                                                                        aaCCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctgg    <  1:231512‑M1/81‑62 (MQ=255)
                                                                         aCCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggc   >  2:231511‑M1/1‑19 (MQ=255)
                                                                           cTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggct  >  1:231510‑M1/1‑17 (MQ=255)
                                                                            |                                                                              
GCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTAT  >  NC_000913/1207714‑1207868

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