New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =0 (0.000)15 (0.620) 14/128 0.6 100% 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

GCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTT  >  NC_000913/1207712‑1207866
                                                                              |                                                                            
gccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgtt                                                                             >  1:216902‑M1/79‑80 (MQ=255)
  caataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGt                                                                          <  2:216869‑M1/5‑1 (MQ=255)
  caataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGt                                                                          >  1:216896‑M1/77‑81 (MQ=255)
    gctgctGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTT                                                                        >  1:216648‑M1/1‑81 (MQ=255)
     ctgctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTG                                                                       <  1:216637‑M1/81‑1 (MQ=255)
        gcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGga                                                                    >  2:216895‑M1/71‑81 (MQ=255)
        ctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGG                                                                     <  1:216638‑M1/80‑1 (MQ=255)
             ttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAg                                                                >  2:216901‑M1/66‑80 (MQ=255)
                 cGCCAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctct                                                           <  2:448292‑M1/81‑5 (MQ=255)
                       ctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTc                                                     >  1:216893‑M1/56‑81 (MQ=255)
                        aTCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattac                                                    <  1:216639‑M1/81‑12 (MQ=255)
                          ccGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacct                                                  >  1:216868‑M1/1‑68 (MQ=255)
                             aaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgc                                                >  2:216887‑M1/50‑80 (MQ=255)
                                  aggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAAc                                          >  1:216881‑M1/45‑81 (MQ=255)
                                    gaatgatgccatttgcttgttatttcttcattttcataaacgTAGGTTTGGCATAAGGAGCTAAACTCAATGCAGCAACgg                                        <  1:902732‑M1/39‑1 (MQ=255)
                                    atatCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcc                                        <  2:216640‑M1/81‑24 (MQ=255)
                                      atgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCt                                      >  2:216906‑M1/41‑81 (MQ=255)
                                      atgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCt                                      >  2:216905‑M1/41‑81 (MQ=255)
                                         aaaTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggtt                                   <  1:216641‑M1/81‑29 (MQ=255)
                                            ccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACAtt                                 >  1:216890‑M1/35‑80 (MQ=255)
                                              ggATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgt                               <  2:216642‑M1/80‑33 (MQ=255)
                                                ttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGAc                            <  1:216649‑M1/51‑1 (MQ=255)
                                                   gCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggcct                          <  1:216643‑M1/80‑38 (MQ=255)
                                                      gttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCAcc                       >  2:216891‑M1/25‑80 (MQ=255)
                                                      cTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcag                      <  1:216644‑M1/81‑42 (MQ=255)
                                                       ttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCgg                     <  1:216653‑M1/58‑1 (MQ=255)
                                                         atttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGt                    >  2:216889‑M1/22‑80 (MQ=255)
                                                            aaaCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactc                >  2:216867‑M1/1‑34 (MQ=255)
                                                                ggCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgcc            <  1:216646‑M1/81‑52 (MQ=255)
                                                                ggCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgcc            <  2:216647‑M1/81‑52 (MQ=255)
                                                                ggCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgcc            <  1:216645‑M1/81‑52 (MQ=255)
                                                                cattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa              >  1:216907‑M1/15‑79 (MQ=255)
                                                                      cataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATAc       >  1:216888‑M1/9‑80 (MQ=255)
                                                                       gAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggc     >  1:216866‑M1/1‑23 (MQ=255)
                                                                       ataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCg     <  2:216657‑M1/74‑1 (MQ=255)
                                                                        aaaaCCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggct    >  1:216865‑M1/1‑22 (MQ=255)
                                                                         aaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGt   <  1:216655‑M1/76‑1 (MQ=255)
                                                                          aaCCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctgg  <  2:216864‑M1/81‑62 (MQ=255)
                                                                              |                                                                            
GCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTT  >  NC_000913/1207712‑1207866

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