New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =14 (0.420)4 (0.150) 3/128 2.7 25.7% coding (290/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 1209619 = 12 (0.460)pseudogene (37/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

GCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGA  >  NC_000913/1207714‑1207870
                                                                            |                                                                                
caataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGt                                                                              <  1:239872‑M1/5‑1 (MQ=255)
  ataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGttt                                                                            >  2:239917‑M1/75‑81 (MQ=255)
   ctgctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTG                                                                           <  1:239647‑M1/81‑1 (MQ=255)
   ctgctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTG                                                                           >  1:239680‑M1/1‑81 (MQ=255)
     agcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTggg                                                                         >  2:239894‑M1/72‑81 (MQ=255)
          gttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAg                                                                    <  1:239871‑M1/15‑1 (MQ=255)
           aaGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAG                                                                    <  1:239648‑M1/80‑1 (MQ=255)
           aaGCCGCAAGAATCCGTTTGGCGATATCAAATCAGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAGGG                                                                   >  2:239847‑M1/1‑81 (MQ=255)
             gCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGct                                                                 >  1:239852‑M1/1‑79 (MQ=255)
               cGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctct                                                               <  2:239649‑M1/81‑5 (MQ=255)
               cGCAAGAATCCGGTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCa                                                               <  1:239650/81‑1 (MQ=255)
                 cAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgc                                                             <  2:239651‑M1/81‑7 (MQ=255)
                 cAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgc                                                             >  2:1042591‑M1/1‑75 (MQ=255)
                    gAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACt                                                          >  2:239677/1‑81 (MQ=255)
                            aagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgcag                                                  >  2:239896‑M1/49‑81 (MQ=255)
                             tGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcgttacctgttg                                                  <  1:239652‑M1/80‑18 (MQ=255)
                             tGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgcag                                                  >  2:239672/1‑80 (MQ=255)
                               gCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcgttacctgttggtg                                               <  1:239653‑M1/81‑21 (MQ=255)
                               gCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcgttacctgttggtg                                               >  1:239868‑M1/1‑61 (MQ=255)
                               gCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgcagca                                                <  1:239654/80‑1 (MQ=255)
                                cGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAAc                                              >  1:239685/1‑81 (MQ=255)
                                aggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAAc                                              >  1:239891‑M1/45‑81 (MQ=255)
                                    atCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGc                                           >  2:239675/1‑80 (MQ=255)
                                     tCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCt                                          <  1:239655/80‑1 (MQ=255)
                                            ggATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcgttacctgttggtgtcccggttccgtg                                  >  2:239866‑M1/1‑48 (MQ=255)
                                            ggATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAgg                                  <  2:239656/81‑1 (MQ=255)
                                                 gCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGc                              <  2:239657/80‑1 (MQ=255)
                                                 gCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCa                             >  2:239681/1‑81 (MQ=255)
                                                   aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCAcc                           <  1:239659/81‑1 (MQ=255)
                                                   aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCAcc                           <  1:239658/81‑1 (MQ=255)
                                                    gttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGAcgcgc                            >  2:239881‑M1/25‑77 (MQ=255)
                                                       tgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGt                        <  1:239660/80‑1 (MQ=255)
                                                          aaaCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcgttacctgttggtgtcccggttccgtggccttcagccactc                    <  1:1042591‑M1/81‑48 (MQ=255)
                                                            aCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                  <  2:239661/81‑1 (MQ=255)
                                                                  ttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCtat             >  2:239889‑M1/11‑80 (MQ=255)
                                                                           cTTGGTTTGGGAGAAGGctctgcgttacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctg   <  1:239863‑M1/81‑65 (MQ=255)
                                                                            ttGGTTTGGGAGAAGGctctgcgttacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctga  <  1:239885‑M1/81‑66 (MQ=255)
                                                                            |                                                                                
GCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGA  >  NC_000913/1207714‑1207870

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