New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =5 (0.090)11 (0.230) 11/230 NT 72.4% coding (290/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 1209619 = 4 (0.080)pseudogene (37/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

GAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTG  >  NC_000913/1207654‑1207930
                                                                                                                                        |                                                                                                                                            
gtggctgccggtgcgaGCCTGGAGAAGGGGCGGCACGACGAGTTACTGACCGGACTTCGCGTGCTGCTGTTTAGCTGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTT                                                                                                                                        <  1:515835‑M1/136‑1 (MQ=255)
   tcagcgcattacggtgtagatcggtgctctggtcatttgttgatacgcgtacagagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTggg                                                                                                                                     >  1:231303‑M1/134‑143 (MQ=255)
          attacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGa                                                                                                                              <  2:599175‑M1/17‑1 (MQ=255)
             ccAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctct                                                                                                                           >  1:985570‑M1/1‑139 (MQ=255)
              cAGCCTGGAGAAGGGGCGGCACGACCAGGTCCTTACCGCACTTCGCGCGCTGCTGTTAAGCCACAAGACTCCGCTTGGCGATCGCAAATCGGATGGCACGGCGCAAACGGCTCTCGAACACCTCGGGCTGGCACAAGGAGCac                                                                                                                          >  2:270959/1‑142 (MQ=14)
               aGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATTGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgc                                                                                                                         <  1:548414‑M1/143‑7 (MQ=255)
                     gatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCaa                                                                                                                   >  1:183192‑M1/116‑143 (MQ=255)
                        aaGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgt                                                                                                                <  1:779483‑M1/143‑16 (MQ=255)
                                  aCGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccg                                                                                                      <  1:734270‑M1/143‑26 (MQ=255)
                                          ttacttacCGCACTTCGCGCGCTGCCGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtg                                                                                              <  2:785423‑M1/143‑34 (MQ=255)
                                                 ccGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTCGGTTTGGGCGAAGGctctgcattacctggtggtgtcccggttccgtggccctca                                                                                       >  2:972716‑M1/1‑103 (MQ=255)
                                                  cGCACTTCGCGCGCTGCTGTTAAGCCGCCAGAATCCGATGAGCGATAGCAAATCGGATGGCACTGAGCAAACGGCTCTCGAAAACAATGGCGTGGGAGAAAGAGCAAAACTCAATGAGGCAAAGGAGACAATAGGGAGCACCg                                                                                      >  2:826714/1‑143 (MQ=255)
                                                   gcaggtcgcgcCCGGCCGTTCGGCCGCAAGAATCCGTTTGGCGATAGCAAATCGGATGGCACTGTGCAAACGGCTCTCGGAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagc                                                                                     <  2:815810‑M1/138‑43 (MQ=255)
                                                                   catgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCg                                                                     <  1:955000‑M1/74‑1 (MQ=255)
                                                                      tgagcCGCAAGAATCCGTTTGGCGATATCAGATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctgg                                                                  <  2:387817‑M1/141‑62 (MQ=255)
                                                                       aaGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGCTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGCTTa                                                                 >  2:184646/1‑143 (MQ=255)
                                                                            tccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAGAGTAGGACGCTCCGGTTTCATAGCTATACCGGTtatggtt                                                            <  2:506266‑M1/83‑4 (MQ=255)
                                                                            gCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggcgataccgcgcctggtccttgagcccctccgccaacaggctggctgaaa                                                            >  1:236131‑M1/1‑76 (MQ=255)
                                                                               ccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGCGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggt                                                         >  2:569819‑M1/58‑143 (MQ=255)
                                                                                    gtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTAt                                                    <  1:450382‑M1/91‑1 (MQ=255)
                                                                                        aagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAg                                                >  1:599175‑M1/49‑143 (MQ=255)
                                                                                                             cttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTgg                           <  1:495825‑M1/116‑1 (MQ=255)
                                                                                                              ttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTggg                          <  1:104668‑M1/117‑1 (MQ=255)
                                                                                                                  tatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTg                      <  1:619856‑M1/121‑1 (MQ=255)
                                                                                                                    tttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCCATCAGTAATCATTCAGTGGGGGTGGa                    <  2:414509‑M1/123‑1 (MQ=255)
                                                                                                                       aaCGGCTCTCGAGAACCTTGGTTTGGGGGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGGACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATg                 <  2:1039450/143‑1 (MQ=255)
                                                                                                                        ttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGGGAATGc                <  2:854920‑M1/127‑1 (MQ=255)
                                                                                                                             tctcGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCGTTCAGTGGGGGTGGAATGCCGCaa           >  1:171392/1‑143 (MQ=255)
                                                                                                                                      ccTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccga  <  1:326024‑M1/143‑126 (MQ=255)
                                                                                                                                        |                                                                                                                                            
GAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTG  >  NC_000913/1207654‑1207930

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