New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =34 (0.820)13 (0.410) 9/102 0.9 26.8% coding (290/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 1209619 = 45 (1.420)pseudogene (37/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

TAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATA  >  NC_000913/1207724‑1207854
                                                                  |                                                                
tAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTT                                                                 >  1:1325298‑M1/1‑68 (MQ=255)
tAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTT                                                                 >  1:2764843‑M1/1‑68 (MQ=255)
 ttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTg                                                                >  1:1585870‑M1/66‑68 (MQ=255)
  aGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGCAAACCTTGG                                                               >  1:1548711‑M1/1‑68 (MQ=255)
  aGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGG                                                               >  1:2054097‑M1/1‑68 (MQ=255)
  aGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGG                                                               >  1:1756851‑M1/1‑68 (MQ=255)
  aGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGG                                                               >  1:2623497‑M1/1‑68 (MQ=255)
   gCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGT                                                              >  1:2368912‑M1/1‑68 (MQ=255)
   aaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGt                                                              >  1:329907‑M1/64‑68 (MQ=255)
     cGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTT                                                            >  1:2634673‑M1/1‑68 (MQ=255)
      tccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTg                                                           <  1:1275736‑M1/8‑1 (MQ=255)
      tccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTg                                                           <  1:3359645‑M1/8‑1 (MQ=255)
      tccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTg                                                           <  1:738824‑M1/8‑1 (MQ=255)
      tccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTg                                                           <  1:776469‑M1/8‑1 (MQ=255)
       ccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTggtttag                                                          <  1:374550‑M1/9‑3 (MQ=255)
        cccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTggg                                                         <  1:328033‑M1/10‑1 (MQ=255)
         ccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGga                                                        >  1:3218361‑M1/58‑68 (MQ=255)
          cctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaccgTTGGTTTGGgag                                                       <  1:2582785‑M1/12‑1 (MQ=255)
          cctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaaggTTGGTTTGGgag                                                       >  1:427893‑M1/57‑68 (MQ=255)
            tggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGaa                                                     <  1:1176465‑M1/14‑1 (MQ=255)
             tCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAG                                                    >  1:329248‑M1/1‑68 (MQ=255)
               cGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGa                                                  <  1:965731/68‑1 (MQ=255)
                 tttGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGc                                                >  1:197191/1‑68 (MQ=255)
                 tttGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGc                                                <  1:2436038/68‑1 (MQ=255)
                 aaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGc                                                >  1:1515277‑M1/50‑68 (MQ=255)
                    gcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCaaa                                             >  1:2796000‑M1/47‑68 (MQ=255)
                    gcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCaaa                                             <  1:2331553‑M1/22‑1 (MQ=255)
                    gcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCaaa                                             >  1:3269219‑M1/47‑68 (MQ=255)
                    gcacgaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCaaa                                             >  1:1773364‑M1/47‑68 (MQ=255)
                     gCGATATCAAATCGGATTGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgca                                            <  1:320992‑M1/68‑8 (MQ=255)
                      cGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcat                                           >  1:141324‑M1/1‑60 (MQ=255)
                      cGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAAc                                           >  1:723638/1‑68 (MQ=255)
                      cGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAAc                                           <  1:2355983/68‑1 (MQ=255)
                      aggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAAc                                           >  1:2013320‑M1/45‑68 (MQ=255)
                         tatCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCa                                        <  1:2714843/68‑1 (MQ=255)
                         aatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCa                                        >  1:2663567‑M1/42‑68 (MQ=255)
                           tgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAAt                                      >  1:1442997‑M1/40‑68 (MQ=255)
                           tCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacct                                      <  1:1099998‑M1/68‑14 (MQ=255)
                           tCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAAt                                      <  1:799145/68‑1 (MQ=255)
                                   gATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAAc                              >  1:2420110/1‑68 (MQ=255)
                                    ttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACg                             >  1:2795696‑M1/31‑68 (MQ=255)
                                    aTGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtc                             <  1:177341‑M1/68‑23 (MQ=255)
                                     tgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                            >  1:1656825‑M1/30‑68 (MQ=255)
                                     tgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                            >  1:929459‑M1/30‑68 (MQ=255)
                                     tGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                            <  1:544711/68‑1 (MQ=255)
                                     tGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                            <  1:77518/68‑1 (MQ=255)
                                     aggCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                            <  1:2717798/67‑1 (MQ=255)
                                         caggtgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                        <  1:1807489/65‑1 (MQ=255)
                                         aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggt                        >  1:2016828‑M1/1‑41 (MQ=255)
                                         aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAGCTCAATGCAGCAACGGCTAc                        >  1:3037374/1‑68 (MQ=255)
                                         aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTaa                        >  1:1055559/1‑67 (MQ=255)
                                         aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                        <  1:1411835/68‑1 (MQ=255)
                                         aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                        >  1:1895977/1‑68 (MQ=255)
                                         aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                        >  1:798747/1‑68 (MQ=255)
                                         aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                        >  1:2408110/1‑68 (MQ=255)
                                         aCTGTGCAAACGGCTCTCGAAAACCTTGGTGTGGGAGAAGGAGCAAAACTCAATGCAGCAACGTCTaa                        >  1:2165995/1‑67 (MQ=37)
                                         aCTGTGCAAACGCCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTaa                        >  1:1374723/1‑67 (MQ=255)
                                          gttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACa                       >  1:679757‑M1/25‑68 (MQ=255)
                                          gttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACa                       >  1:1935830‑M1/25‑68 (MQ=255)
                                              gCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgt                   >  1:2431696‑M1/1‑36 (MQ=255)
                                              gCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAg                   >  1:1094793/1‑68 (MQ=255)
                                                 aaCGGCTCTCGAAAGCCTTGGTTTGGGAGAAGGctctgcattccctgttggtgtcccggttccgtggc                <  1:1147287‑M1/68‑36 (MQ=255)
                                                 aaCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggc                <  1:2291533‑M1/68‑36 (MQ=255)
                                                   cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGAccc              >  1:1747419/1‑66 (MQ=255)
                                                   cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGc              <  1:363230/68‑1 (MQ=255)
                                                   cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGc              >  1:3285671/1‑68 (MQ=255)
                                                   cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGc              <  1:557109/68‑1 (MQ=255)
                                                   cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGc              >  1:772274/1‑68 (MQ=255)
                                                   cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGc              >  1:1941579/1‑68 (MQ=255)
                                                   cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAAATTAGGAAGc              >  1:2103108/1‑68 (MQ=255)
                                                   cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAAATTAGGAAGc              >  1:2075100/1‑68 (MQ=255)
                                                   cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAAAATAGGaacc              >  1:2669173/1‑65 (MQ=37)
                                                          cGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGtt       >  1:3290258/1‑68 (MQ=255)
                                                          cGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGtt       <  1:1150882/68‑1 (MQ=255)
                                                               aCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa  >  1:2997453/1‑68 (MQ=255)
                                                               aCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa  <  1:1941146/68‑1 (MQ=255)
                                                                  |                                                                
TAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATA  >  NC_000913/1207724‑1207854

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