New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 12078053 (0.100)13 (0.540) 13/128 0.7 80.3% coding (305/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 = 1209602 4 (0.170)pseudogene (54/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

GCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGC  >  NC_000913/1207727‑1207885
                                                                              |                                                                                
aaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGa                                                                                 >  2:221303‑M1/64‑80 (MQ=255)
 ccGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGc                                                                               <  1:221070/81‑1 (MQ=255)
 aatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGc                                                                               >  1:221298‑M1/63‑81 (MQ=255)
        ctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACt                                                                         <  1:221085‑M1/25‑1 (MQ=255)
          tCCGTTTGGCGATATCAAATCGGACGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacc                                                                      >  2:221281‑M1/1‑69 (MQ=255)
           gtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAAt                                                                     >  2:221292‑M1/53‑81 (MQ=255)
               ttGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttg                                                                 >  1:221280‑M1/1‑64 (MQ=255)
               aagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgcag                                                                 <  2:221078‑M1/33‑1 (MQ=255)
                agcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgcag                                                                 <  1:221079‑M1/33‑1 (MQ=255)
                  gCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtg                                                              >  1:221279‑M1/1‑61 (MQ=255)
                   cGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgt                                                             <  1:221071‑M1/81‑22 (MQ=255)
                     gaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                                                           >  1:221295‑M1/43‑81 (MQ=255)
                       atgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCt                                                         >  1:221296‑M1/41‑81 (MQ=255)
                          aaaTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggt                                                       <  2:221072‑M1/80‑28 (MQ=255)
                           aaTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttc                                                     <  1:221073‑M1/81‑30 (MQ=255)
                            aTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttcc                                                    <  2:221074‑M1/81‑31 (MQ=255)
                             ccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACAtt                                                    >  1:221300‑M1/35‑80 (MQ=255)
                              cGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgt                                                  >  2:221277‑M1/1‑49 (MQ=255)
                                    cttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCa                                            <  2:221081‑M1/54‑1 (MQ=255)
                                        ttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCg                                         >  2:221289‑M1/24‑80 (MQ=255)
                                           tttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGttt                                     <  2:221080‑M1/61‑1 (MQ=255)
                                                 ggCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgc                                >  1:221275‑M1/1‑30 (MQ=255)
                                                  gCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgcca                              <  1:221075‑M1/81‑53 (MQ=255)
                                                  attttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGc                               >  2:221301‑M1/14‑80 (MQ=255)
                                                  attttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCt                              >  2:221307‑M1/14‑81 (MQ=255)
                                                     ttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCtata                           <  2:221084‑M1/71‑1 (MQ=255)
                                                     ctcGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaa                             >  1:221274‑M1/1‑26 (MQ=255)
                                                      tcGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCtat                            <  2:221076/79‑1 (MQ=255)
                                                          aaaCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGt                      >  2:221086/1‑81 (MQ=255)
                                                           aaCCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctgg                     <  1:221273‑M1/81‑62 (MQ=255)
                                                            acgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTa                    >  1:221293‑M1/4‑81 (MQ=255)
                                                            aCCTTGGTTTGGGAGAAGGctctgctttacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggc                    >  2:221272‑M1/1‑19 (MQ=255)
                                                             cgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTAt                   >  1:221291‑M1/3‑81 (MQ=255)
                                                                 ggTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaa                <  2:221271‑M1/80‑67 (MQ=255)
                                                                   tttGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatg             >  1:221270‑M1/1‑12 (MQ=255)
                                                                      gggAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaa          >  2:221268‑M1/1‑9 (MQ=255)
                                                                      gggAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgca           <  2:221269‑M1/80‑72 (MQ=255)
                                                                       ggAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaat         >  1:221267‑M1/1‑8 (MQ=255)
                                                                         agaAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatgg       >  2:221266‑M1/1‑6 (MQ=255)
                                                                         AGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTAt       >  2:221077‑M1/1‑81 (MQ=255)
                                                                           aaGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcgatggtg     <  2:221265‑M1/81‑78 (MQ=255)
                                                                           AAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTAtt      <  2:221087‑M1/80‑1 (MQ=255)
                                                                            AGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGa    >  2:221288‑M1/1‑81 (MQ=255)
                                                                              gctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccg  >  2:221262‑M1/1‑1 (MQ=255)
                                                                              gctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccg  >  2:221261‑M1/1‑1 (MQ=255)
                                                                              gctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccg  <  1:221264‑M1/81‑81 (MQ=255)
                                                                              gctctgcattacctgttggtgtcccggttccgtggccttcaaccactccgccaacaggctggctgaaatgcaatggtgccg  >  2:221263‑M1/1‑1 (MQ=255)
                                                                              |                                                                                
GCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGC  >  NC_000913/1207727‑1207885

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.