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

CGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTAT  >  NC_000913/1207711‑1207868
                                                                               |                                                                              
cgcgcgCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTT                                                                               >  1:219019‑M1/1‑81 (MQ=255)
cgccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgtt                                                                               <  2:219066‑M1/2‑1 (MQ=255)
 gcgcgcTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTG                                                                              >  2:219030‑M1/1‑81 (MQ=255)
    aataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGtt                                                                           >  1:219087‑M1/76‑81 (MQ=255)
        agcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTgg                                                                        <  1:219065‑M1/9‑1 (MQ=255)
          tgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGG                                                                        >  1:219020‑M1/1‑78 (MQ=255)
            ttAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGA                                                                    <  2:218814‑M1/80‑1 (MQ=255)
            ttAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAG                                                                     <  2:218815‑M1/79‑1 (MQ=255)
                aaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGa                                                                >  2:219090‑M1/64‑80 (MQ=255)
                   gCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctg                                                            >  1:219064‑M1/1‑76 (MQ=255)
                   gCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctg                                                            >  1:219063‑M1/1‑76 (MQ=255)
                        ctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTc                                                       >  2:219080‑M1/56‑81 (MQ=255)
                          ggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCaa                                                     <  2:218830‑M1/28‑1 (MQ=255)
                           gtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCaa                                                     >  2:219083‑M1/53‑80 (MQ=255)
                           gtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCa                                                      >  1:219089‑M1/53‑79 (MQ=255)
                           gtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAAt                                                    >  2:219082‑M1/53‑81 (MQ=255)
                            cGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctg                                                   <  1:218817‑M1/81‑15 (MQ=255)
                            cGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctg                                                   <  1:218816‑M1/81‑15 (MQ=255)
                             aaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATg                                                   <  2:218832‑M1/30‑1 (MQ=255)
                              aaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgc                                                  >  1:219088‑M1/50‑80 (MQ=255)
                                tGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttgg                                               >  1:219062‑M1/1‑63 (MQ=255)
                                 ggCGATATCAAATCGGATGGCACTGTGCAAGCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttgg                                               <  2:218818‑M1/80‑19 (MQ=255)
                                  gCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtg                                             >  1:219061‑M1/1‑61 (MQ=255)
                                    gATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtc                                           <  2:218819‑M1/81‑23 (MQ=255)
                                     atatCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcc                                          <  1:218820‑M1/81‑24 (MQ=255)
                                      tatCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtccc                                         <  1:218821‑M1/81‑25 (MQ=255)
                                      tatCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtccc                                         <  1:218822‑M1/81‑25 (MQ=255)
                                          atgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACa                                     <  2:218831‑M1/44‑1 (MQ=255)
                                                tttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGa                               >  1:219085‑M1/32‑81 (MQ=255)
                                                   ggCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggcct                            >  1:219060‑M1/1‑44 (MQ=255)
                                                     cACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttc                          <  1:218823‑M1/81‑40 (MQ=255)
                                                       gttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCg                        <  2:218833‑M1/57‑1 (MQ=255)
                                                        tgtgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagc                       >  2:219059‑M1/1‑39 (MQ=255)
                                                         tatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGt                      <  1:218835‑M1/59‑1 (MQ=255)
                                                            cAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccact                   >  1:219058‑M1/1‑35 (MQ=255)
                                                            cAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccac                    <  1:218825‑M1/80‑46 (MQ=255)
                                                            cAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccac                    <  1:218824‑M1/80‑46 (MQ=255)
                                                             aaaCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCa                  <  2:218826/81‑1 (MQ=255)
                                                               aCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccg                >  1:219057‑M1/1‑32 (MQ=255)
                                                                   ctctCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagcactccgccaac            <  2:218827‑M1/81‑54 (MQ=255)
                                                                     ctcGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaa            >  1:219056‑M1/1‑26 (MQ=255)
                                                                        ataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATAcc        >  1:219081‑M1/8‑80 (MQ=255)
                                                                          aaaCCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctg     <  1:219055‑M1/81‑61 (MQ=255)
                                                                              gTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTAt  <  1:218834‑M1/79‑1 (MQ=255)
                                                                               |                                                                              
CGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTAT  >  NC_000913/1207711‑1207868

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