New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =18 (0.330)32 (0.790) 18/96 0.3 61.2% coding (290/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 1209619 = 27 (0.670)pseudogene (37/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

TGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAA  >  NC_000913/1209559‑1209683
                                                            |                                                                
tGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGACTTCTCCCAAACCAACGttt                                                              >  1:2787429/1‑65 (MQ=255)
 accggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTa                                                             >  1:1417814‑M2/60‑65 (MQ=255)
      aaGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGaaa                                                        >  1:336133/1‑65 (MQ=255)
        ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAAt                                                      >  1:3718900/1‑65 (MQ=255)
         gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAACTg                                                     >  1:1212541/1‑65 (MQ=255)
         gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                                     >  1:3534202/1‑65 (MQ=255)
         gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                                     >  1:487369/1‑65 (MQ=255)
         gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                                     >  1:3507150/1‑65 (MQ=255)
         gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                                     >  1:316534/1‑65 (MQ=255)
         gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                                     >  1:323732/1‑65 (MQ=255)
               gAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGaaa                                               >  1:623747/1‑65 (MQ=255)
               atgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGaaa                                               <  1:4331100‑M2/20‑1 (MQ=255)
                tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGCAGAAAt                                              >  1:2009805‑M2/45‑65 (MQ=255)
                tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                              >  1:1384546‑M2/45‑65 (MQ=255)
                tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                              >  1:2186044‑M2/45‑65 (MQ=255)
                tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                              >  1:617915‑M2/45‑65 (MQ=255)
                 gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                             >  1:2872726‑M2/44‑65 (MQ=255)
                 gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                             >  1:2806303‑M2/44‑65 (MQ=255)
                 gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                             >  1:634625‑M2/44‑65 (MQ=255)
                    gccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACa                                          >  1:1761631‑M2/41‑65 (MQ=255)
                     ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                         >  1:3523587‑M2/40‑65 (MQ=255)
                      gACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAg                                        >  1:446843/1‑65 (MQ=255)
                       acacCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGc                                       >  1:3668653/1‑65 (MQ=255)
                        cacCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCa                                      <  1:1392881/65‑1 (MQ=255)
                         tgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCaa                                     <  1:1519940‑M2/30‑1 (MQ=255)
                          ccAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCaaa                                    <  1:2491624/65‑1 (MQ=255)
                             aCAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATgg                                 >  1:2151598/1‑65 (MQ=255)
                              cattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGc                                <  1:1431394‑M2/35‑1 (MQ=255)
                               attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                               >  1:4065248‑M2/30‑65 (MQ=255)
                               attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                               <  1:1829450‑M2/36‑1 (MQ=255)
                               attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                               >  1:2007951‑M2/30‑65 (MQ=255)
                               aGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                               >  1:2749580/1‑65 (MQ=255)
                                 gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                             >  1:2314801/1‑65 (MQ=255)
                                 gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                             >  1:1860399/1‑65 (MQ=255)
                                 gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                             >  1:2930274/1‑65 (MQ=255)
                                 gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                             >  1:2836903/1‑65 (MQ=255)
                                  gagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                            <  1:2558294‑M2/39‑1 (MQ=255)
                                   agttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatcat                           <  1:972103‑M2/40‑1 (MQ=255)
                                    gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:1004254‑M2/41‑1 (MQ=255)
                                    gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:711312‑M2/41‑1 (MQ=255)
                                    gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:2193238‑M2/41‑1 (MQ=255)
                                    gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:3446995‑M2/41‑1 (MQ=255)
                                    gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:3362368‑M2/41‑1 (MQ=255)
                                    gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          >  1:841420‑M2/25‑65 (MQ=255)
                                    gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:2677966‑M2/41‑1 (MQ=255)
                                    gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:2928142‑M2/41‑1 (MQ=255)
                                      tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                        >  1:1538336‑M2/23‑65 (MQ=255)
                                      gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                        >  1:852266/1‑65 (MQ=255)
                                      gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                        >  1:781514/1‑65 (MQ=255)
                                        agagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg                      >  1:217212/1‑65 (MQ=255)
                                         gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     <  1:2203390‑M2/46‑1 (MQ=255)
                                         gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     >  1:1628253‑M2/20‑65 (MQ=255)
                                         gagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     <  1:369957/65‑1 (MQ=255)
                                           tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                   >  1:965417‑M2/18‑65 (MQ=255)
                                           tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                   >  1:760148‑M2/18‑65 (MQ=255)
                                           tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                   >  1:2178107‑M2/18‑65 (MQ=255)
                                           gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                   <  1:2583804/65‑1 (MQ=255)
                                           gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                   >  1:1334259/1‑65 (MQ=255)
                                           gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGACATCATTCCTGCtt                   >  1:2923112/1‑65 (MQ=255)
                                            ccttctcccaaaccaaCGTTTATGAAAATGGAGAAATAACAAGCAAATGGCATCATTCCTGCttt                  >  1:2941539‑M2/17‑65 (MQ=255)
                                              ttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa                <  1:1897913‑M2/51‑1 (MQ=255)
                                                 tcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCa             <  1:1377921‑M2/54‑1 (MQ=255)
                                                  cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAg            >  1:71087‑M2/11‑65 (MQ=255)
                                                   ccaaaccaaCGTTTATGAAAATGAAGAAATAGCAAGCAAATGGCATCATTCCTGCTTTTACCAgg           <  1:4315222‑M2/56‑1 (MQ=255)
                                                   ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgg           <  1:3448552‑M2/56‑1 (MQ=255)
                                                      aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg        <  1:3165987‑M2/59‑1 (MQ=255)
                                                       accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:1102414‑M2/60‑1 (MQ=255)
                                                       accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:2581647‑M2/60‑1 (MQ=255)
                                                       accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:2537539‑M2/60‑1 (MQ=255)
                                                       accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:3543916‑M2/60‑1 (MQ=255)
                                                       accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:132362‑M2/60‑1 (MQ=255)
                                                       accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:4249822‑M2/60‑1 (MQ=255)
                                                       accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:4322499‑M2/60‑1 (MQ=255)
                                                        ccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAt      <  1:4265503‑M2/61‑1 (MQ=255)
                                                           aCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTa   <  1:449220‑M2/64‑1 (MQ=255)
                                                            cGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa  <  1:1785132/65‑1 (MQ=255)
                                                            cGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa  >  1:2804620/1‑65 (MQ=255)
                                                            |                                                                
TGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAA  >  NC_000913/1209559‑1209683

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