New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =10 (0.210)16 (0.390) 15/232 NT 56.4% coding (290/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 1209619 = 16 (0.390)pseudogene (37/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

TGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGC  >  NC_000913/1207649‑1207931
                                                                                                                                             |                                                                                                                                             
tggtggAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTT                                                                                                                                              <  1:446756‑M1/143‑1 (MQ=255)
  gtggAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAACTCGGATGGCACTGTGCAAACGGCACTCGCACACCTTGG                                                                                                                                            >  2:75852‑M1/1‑143 (MQ=255)
       ttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgattgttatttcttcattttcatagacgTTGGTATgg                                                                                                                                       <  2:47569‑M1/9‑1 (MQ=255)
       ggCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGG                                                                                                                                       <  1:474925‑M1/143‑1 (MQ=255)
       ggCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGGTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGG                                                                                                                                       <  2:814591‑M1/143‑1 (MQ=255)
          agcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGgag                                                                                                                                    <  1:70118‑M1/12‑1 (MQ=255)
          agcgcagtacgttggagcgctgtgttctggtcatttgttgatacgcgtacatagccaattagcatgttaaatccccgtggtaaaagcaggaatgatgccatttgcgtgttatttcttcattttcataaacgTTGGTTTGGgag                                                                                                                                    >  2:357534‑M1/132‑143 (MQ=255)
           tCCGGTGCCAGCCTGGAGAAGGGGCGGATGGACCATGGACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCGAGAATCCGTCTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGCAAACCTAGGTTCGGGAGA                                                                                                                                   >  2:506085‑M1/1‑143 (MQ=255)
             cgcagccagCCGGGGGAGGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAG                                                                                                                                 <  2:679031‑M1/139‑1 (MQ=255)
                ttacgctgcagacctgtgttcttttcattcgttgatacgcgtacatagccaataagcatgttaaatcccccgggtaaaagcaggaatgatgccatgtgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAg                                                                                                                              >  1:59818‑M1/126‑143 (MQ=255)
                    gttgtagatctgtgttctggtaatttgttgatacgcgtacatagccaataagcatgtaaaatcatcctggtaaaacaaggaatgatgaaatttgcttgatatttcttcattttcataaacgAAGGTTTGGGAGAaggaggaat                                                                                                                          <  2:153834‑M1/22‑5 (MQ=255)
                     ttgtagttctgtgttctggtcttttgttgatacgagtacatagccaataagcatgttaaatcccccttgtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCaaaa                                                                                                                         >  2:284348‑M1/121‑143 (MQ=255)
                        tGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGTTGTTAAGCCGCAAGAAGCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcatta                                                                                                                      <  2:59160‑M1/143‑11 (MQ=255)
                        tGGAGAAGGGGCGGCACGACCAGTTACGTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTc                                                                                                                      <  1:287533/143‑1 (MQ=255)
                         agatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCa                                                                                                                     >  2:710596‑M1/117‑143 (MQ=255)
                                        cGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccgg                                                                                                      <  2:433617‑M1/143‑27 (MQ=255)
                                         catttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                                                                                     >  1:647493‑M1/101‑143 (MQ=255)
                                          aCCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggtt                                                                                                    <  1:837525‑M1/143‑29 (MQ=255)
                                           tttctggatacgcgtacacccccaacacgcatgttaactccccctgggaacagcaggaatgatgccatttgcttgttatttctgcattttcataaacgTTGGTTTGTGAGAAGGAGCAAACCTCAATGCAGCAACGGCTACAt                                                                                                   >  1:634486‑M1/99‑143 (MQ=255)
                                                tacttacCGCACTTCGCGAGCTGCTGTAAAGCCGCAAGAATCCGTTGGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacatgttggtgtcccggttccgtgg                                                                                              >  1:217182‑M1/1‑109 (MQ=255)
                                                 acttacCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggc                                                                                             >  1:13036‑M1/1‑108 (MQ=255)
                                                  atccgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACg                                                                                            >  1:541017‑M1/92‑143 (MQ=255)
                                                   ttacCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggcct                                                                                           >  1:492037‑M1/1‑106 (MQ=255)
                                                   ttacCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGAGGc                                                                                           >  2:734551/1‑143 (MQ=255)
                                                   tacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGCTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACCTTaggagga                                                                                           <  2:182042‑M1/53‑4 (MQ=255)
                                                       cggacatagccaatacgcatgttaaatgcccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTGGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCg                                                                                       >  2:58163‑M1/87‑143 (MQ=255)
                                                       cGCACTTCGCGCGCTGCTGTTAAGCCGCAAGACTCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAGAACTCAATGCAGCAAGGGCTACATTAGGACGCACCg                                                                                       >  1:103361/1‑143 (MQ=255)
                                                        gCATTTCGCGCGCTGCTGTTATGCCGCAAGAATCCGGTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagc                                                                                      <  2:354194‑M1/143‑43 (MQ=255)
                                                               gcgcgcTGCTGTTAAGCCGCAAGAATCCGTGTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                                                                               <  1:315368/143‑1 (MQ=255)
                                                                cctataaccatgttaaatccccctgttaaaagcaggaatattgcctttttcttgttatttcttcattttaataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAg                                                                              >  2:44918‑M1/78‑143 (MQ=255)
                                                                       ctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCACTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacagg                                                                       <  2:451067‑M1/143‑58 (MQ=255)
                                                                            ttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTa                                                                  <  1:351515‑M1/78‑1 (MQ=255)
                                                                                     cctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGCAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggtg                                                         <  1:29579‑M1/87‑1 (MQ=255)
                                                                                      ctggtaaaagcagtaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggtgg                                                        >  2:290595‑M1/56‑143 (MQ=255)
                                                                                      aaTCCGTTTGGCGATATGAACGCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaa                                                        <  2:629508‑M1/143‑73 (MQ=255)
                                                                                      aaTCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtgcccgtcagccactccgccaacaggctggctgaaatgcaa                                                        <  1:629381‑M1/143‑73 (MQ=255)
                                                                                         ccGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTAt                                                     <  1:794078/143‑1 (MQ=255)
                                                                                          cGTTTGGCGAGTGAAAATCGGATGGCACTGTGCAAGCGGCTCTCGAAAGCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTAtt                                                    <  1:887428/143‑1 (MQ=255)
                                                                                            tttGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgc                                                  <  1:918346‑M1/143‑79 (MQ=255)
                                                                                            aaagcaggaatgatgctatttgcttctcatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAACGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGa                                                  >  2:358984‑M1/50‑143 (MQ=255)
                                                                                            aaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGa                                                  >  1:611541‑M1/50‑143 (MQ=255)
                                                                                                    tatCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCACTCAg                                          >  2:824604/1‑143 (MQ=255)
                                                                                                     atgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGAGTGGTGGTATTGAGCACTCAGt                                         <  2:4254‑M1/103‑1 (MQ=255)
                                                                                                      tCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTa                                        >  2:691672/1‑143 (MQ=255)
                                                                                                         actcgGACGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccgggtccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgct                                     <  1:511464‑M1/141‑92 (MQ=255)
                                                                                                                      tgtgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggcgccgctttttctgctgaagcagacccgg                        >  2:34205‑M1/1‑39 (MQ=255)
                                                                                                                             ttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGc                 <  2:623141‑M1/127‑1 (MQ=255)
                                                                                                                                     cGAAAACCGTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggctt         >  2:701014‑M1/1‑24 (MQ=255)
                                                                                                                                      gAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAATAGc        >  1:538210/1‑143 (MQ=255)
                                                                                                                                      gAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGc        >  1:729448/1‑143 (MQ=255)
                                                                                                                                      ataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAACAGc        <  2:801311‑M1/136‑1 (MQ=255)
                                                                                                                                        ccaccTTGGTTCGGGAGAAGGctctgcattccctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatc      <  1:219284‑M1/141‑123 (MQ=255)
                                                                                                                                         aaCGTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGGAACAGCTACATTAGGACGCACAGGTTTCATAACTATACCGGTTATGAGAGGTGGTATTGAGCAATAAGTAATGAGTGCGGTGCGGTGGAATGTGGAAACAGCATc     >  2:571890/1‑143 (MQ=255)
                                                                                                                                           ccTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccga   >  1:522288‑M1/1‑18 (MQ=255)
                                                                                                                                            cTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGc  >  2:873930/1‑143 (MQ=255)
                                                                                                                                             |                                                                                                                                             
TGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGC  >  NC_000913/1207649‑1207931

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