New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =3 (0.100)17 (0.710) 17/128 0.4 84.2% coding (290/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 1209619 = 4 (0.170)pseudogene (37/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

GCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTAT  >  NC_000913/1207714‑1207868
                                                                            |                                                                              
gcgcTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGT                                                                            >  1:221228‑M1/1‑81 (MQ=255)
  gctgctGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTT                                                                          >  2:221243‑M1/1‑81 (MQ=255)
     gctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGG                                                                        >  2:221246‑M1/1‑80 (MQ=255)
      ctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGG                                                                       <  2:221067‑M1/80‑1 (MQ=255)
      ctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGG                                                                        <  1:221068‑M1/79‑1 (MQ=255)
           ttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAg                                                                  >  1:221305‑M1/66‑80 (MQ=255)
           aaGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAG                                                                  <  2:221069‑M1/80‑1 (MQ=255)
             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)
                                                                            |                                                                              
GCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTAT  >  NC_000913/1207714‑1207868

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