New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =2 (0.040)59 (1.550) 25/96 0.1 90.4% coding (290/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 1209619 = 11 (0.290)pseudogene (37/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

AGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATT  >  NC_000913/1209557‑1209680
                                                              |                                                             
gaaaccggagcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaacgaaCGt                                                             >  1:1831366‑M2/63‑65 (MQ=255)
  aaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGttt                                                           >  1:3345990‑M2/61‑65 (MQ=255)
   accggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTa                                                          >  1:1530846‑M2/60‑65 (MQ=255)
   accggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTa                                                          >  1:882479‑M2/60‑65 (MQ=255)
    ccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTAt                                                         >  1:3864717‑M2/59‑65 (MQ=255)
        tgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTAtgaca                                                     >  1:2125170‑M2/55‑63 (MQ=255)
        aaGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGaaa                                                     >  1:3589668/1‑65 (MQ=255)
           gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                                  >  1:3889196/1‑65 (MQ=255)
               taatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaaga                                              <  1:1794793‑M2/18‑1 (MQ=255)
                  tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAACATGAAGAAAt                                           >  1:1238045‑M2/45‑65 (MQ=255)
                  tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                           >  1:432026‑M2/45‑65 (MQ=255)
                  tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                           <  1:162497‑M2/21‑1 (MQ=255)
                  tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                           >  1:976198‑M2/45‑65 (MQ=255)
                  tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                           >  1:3615421‑M2/45‑65 (MQ=255)
                  tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                           >  1:2865657‑M2/45‑65 (MQ=255)
                  tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                           >  1:3698471‑M2/45‑65 (MQ=255)
                   gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                          >  1:3697511‑M2/44‑65 (MQ=255)
                   gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                          <  1:3485049‑M2/22‑1 (MQ=255)
                     agccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAc                                        >  1:2282644‑M2/42‑65 (MQ=255)
                       ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATaactg                                      >  1:1701449‑M2/40‑63 (MQ=255)
                       ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                      >  1:881672‑M2/40‑65 (MQ=255)
                       ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                      >  1:3912238‑M2/40‑65 (MQ=255)
                       ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                      >  1:2820831‑M2/40‑65 (MQ=255)
                       ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                      >  1:2888690‑M2/40‑65 (MQ=255)
                       ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                      >  1:1812277‑M2/40‑65 (MQ=255)
                       ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                      >  1:1594793‑M2/40‑65 (MQ=255)
                       ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                      >  1:3559537‑M2/40‑65 (MQ=255)
                       ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                      >  1:2461393‑M2/40‑65 (MQ=255)
                       ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                      >  1:1394187‑M2/40‑65 (MQ=255)
                        gACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAg                                     >  1:364228/1‑65 (MQ=255)
                        cgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAg                                     <  1:1835039‑M2/27‑1 (MQ=255)
                          ttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCa                                   >  1:3333031‑M2/37‑65 (MQ=255)
                          ttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCa                                   >  1:2529445‑M2/37‑65 (MQ=255)
                            gctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCaaa                                 <  1:3980310‑M2/31‑1 (MQ=255)
                            gctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCaaa                                 >  1:256351‑M2/35‑65 (MQ=255)
                              aaCAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATg                               >  1:3350054/1‑65 (MQ=255)
                               gcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATgg                              <  1:177289‑M2/34‑1 (MQ=255)
                                cattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACACGCAAATGGc                             >  1:752209‑M2/31‑65 (MQ=255)
                                   gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                          >  1:1885884/1‑65 (MQ=255)
                                    gagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                         <  1:1156673‑M2/39‑1 (MQ=255)
                                     agttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatcat                        <  1:891146‑M2/40‑1 (MQ=255)
                                      gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                       <  1:1930626‑M2/41‑1 (MQ=255)
                                      gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                       <  1:3304987‑M2/41‑1 (MQ=255)
                                      gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                       <  1:3303006‑M2/41‑1 (MQ=255)
                                      gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                       >  1:1774397‑M2/25‑65 (MQ=255)
                                      gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                       <  1:764140‑M2/41‑1 (MQ=255)
                                      gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                       >  1:791063‑M2/25‑65 (MQ=255)
                                      gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                       >  1:131762‑M2/25‑65 (MQ=255)
                                      gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCACCAtt                       >  1:1220914‑M2/25‑65 (MQ=255)
                                      gttttgctccttctcccaaaccaaCGTTTATGAAAATGAACAAATAACAAGCAAATGGCATCAtt                       >  1:724492‑M2/25‑65 (MQ=255)
                                      aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                       <  1:980932/65‑1 (MQ=255)
                                      aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                       <  1:3252761/65‑1 (MQ=255)
                                       tGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc                      >  1:1278946/1‑65 (MQ=255)
                                       tGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc                      >  1:2362191/1‑65 (MQ=255)
                                        tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                     >  1:1368338‑M2/23‑65 (MQ=255)
                                        tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                     >  1:2219573‑M2/23‑65 (MQ=255)
                                        tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                     >  1:2274572‑M2/23‑65 (MQ=255)
                                           gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                  <  1:2025474‑M2/46‑1 (MQ=255)
                                           gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                  <  1:943071‑M2/46‑1 (MQ=255)
                                           gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                  <  1:1716769‑M2/46‑1 (MQ=255)
                                           gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                  <  1:3986739‑M2/46‑1 (MQ=255)
                                           gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                  <  1:2353577‑M2/46‑1 (MQ=255)
                                           gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                  >  1:2516478‑M2/20‑65 (MQ=255)
                                           gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                  <  1:2924492‑M2/46‑1 (MQ=255)
                                           gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                  <  1:2322912‑M2/46‑1 (MQ=255)
                                           gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAACCAAATGGCATCATTCCTGc                  <  1:19163‑M2/46‑1 (MQ=255)
                                           gagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                  <  1:2972821/65‑1 (MQ=255)
                                             tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                >  1:3043319‑M2/18‑65 (MQ=255)
                                             tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                >  1:3963361‑M2/18‑65 (MQ=255)
                                             gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                >  1:479985/1‑65 (MQ=255)
                                                ttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa             <  1:2287755‑M2/51‑1 (MQ=255)
                                                   tcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCa          <  1:1522460‑M2/54‑1 (MQ=255)
                                                    cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAg         >  1:482763‑M2/11‑65 (MQ=255)
                                                    cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAg         >  1:596733‑M2/11‑65 (MQ=255)
                                                      caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg       <  1:306569‑M2/57‑1 (MQ=255)
                                                        aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTGCCAggggg     <  1:1412119‑M2/59‑1 (MQ=255)
                                                        aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg     <  1:2589746‑M2/59‑1 (MQ=255)
                                                        aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg     <  1:1282040‑M2/59‑1 (MQ=255)
                                                         accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa    <  1:4018122‑M2/60‑1 (MQ=255)
                                                         accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa    <  1:174323‑M2/60‑1 (MQ=255)
                                                         accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa    <  1:2602146‑M2/60‑1 (MQ=255)
                                                         accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa    <  1:1393175‑M2/60‑1 (MQ=255)
                                                         accaaCGTTTACGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa    <  1:1969385‑M2/60‑1 (MQ=255)
                                                          ccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAt   >  1:3708416‑M2/5‑65 (MQ=255)
                                                           caaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAtt  >  1:2171581‑M2/4‑65 (MQ=255)
                                                              |                                                             
AGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATT  >  NC_000913/1209557‑1209680

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