New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =7 (0.160)33 (0.990) 23/96 0.1 74.8% coding (290/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 1209619 = 17 (0.510)pseudogene (37/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

GAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATT  >  NC_000913/1209556‑1209680
                                                               |                                                             
gAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACg                                                              <  1:3459737/65‑1 (MQ=255)
 gaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGt                                                             <  1:2838646‑M2/3‑1 (MQ=255)
 gaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGt                                                             >  1:1881500‑M2/63‑65 (MQ=255)
    accggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTa                                                          >  1:295610‑M2/60‑65 (MQ=255)
     ccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTAt                                                         >  1:2631505‑M2/59‑65 (MQ=255)
            gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                                  >  1:1046760/1‑65 (MQ=255)
            gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                                  >  1:2228342/1‑65 (MQ=255)
            gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                                  >  1:1566971/1‑65 (MQ=255)
              cctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaa                                                >  1:1920418‑M2/50‑65 (MQ=255)
                taatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaaga                                              <  1:1780948‑M2/18‑1 (MQ=255)
                taatctagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaaga                                              <  1:871527‑M2/18‑1 (MQ=255)
                  atgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGaaa                                            <  1:2939223‑M2/20‑1 (MQ=255)
                   tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                           >  1:1183816‑M2/45‑65 (MQ=255)
                    gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAatg                                          >  1:1602575‑M2/44‑64 (MQ=255)
                    gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                          >  1:2347006‑M2/44‑65 (MQ=255)
                    gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                          >  1:1228257‑M2/44‑65 (MQ=255)
                      agccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAc                                        >  1:582239‑M2/42‑65 (MQ=255)
                       gggACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACa                                       <  1:2530935/65‑1 (MQ=255)
                        ccgttgctgcattgagttttgctccttcttccaaaccaaCGTTTATGAAAATGAAGAAATAAcat                                      >  1:440312‑M2/40‑64 (MQ=255)
                        ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                      >  1:1843522‑M2/40‑65 (MQ=255)
                        ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                      >  1:1000518‑M2/40‑65 (MQ=255)
                         cgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAGCAAg                                     >  1:1713839‑M2/39‑65 (MQ=255)
                         cgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAg                                     >  1:1455441‑M2/39‑65 (MQ=255)
                            acCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCaa                                  >  1:769384/1‑65 (MQ=255)
                             gctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCaaa                                 <  1:1610014‑M2/31‑1 (MQ=255)
                             ccAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCaaa                                 <  1:3012429/65‑1 (MQ=255)
                               tgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATg                               >  1:384097‑M2/33‑65 (MQ=255)
                                gcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATgg                              <  1:190200‑M2/34‑1 (MQ=255)
                                 catttagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGc                             <  1:1779073‑M2/35‑1 (MQ=255)
                                 cAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGc                             <  1:1514973/65‑1 (MQ=255)
                                  attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                            >  1:448382‑M2/30‑65 (MQ=255)
                                    tgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                          >  1:533736‑M2/28‑65 (MQ=255)
                                     tAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                         <  1:2319318/65‑1 (MQ=255)
                                     gagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                         <  1:1620809‑M2/39‑1 (MQ=255)
                                       gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                       >  1:1437535‑M2/25‑65 (MQ=255)
                                       gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                       >  1:3178052‑M2/25‑65 (MQ=255)
                                       gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                       >  1:1807349‑M2/25‑65 (MQ=255)
                                       aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                       <  1:1490592/65‑1 (MQ=255)
                                       aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                       <  1:1079167/65‑1 (MQ=255)
                                         tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                     >  1:984216‑M2/23‑65 (MQ=255)
                                          cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                    >  1:949114/1‑65 (MQ=255)
                                          cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                    >  1:3364987/1‑65 (MQ=255)
                                          cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                    >  1:868677/1‑65 (MQ=255)
                                            gctccttctcccaagccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                  <  1:200690‑M2/46‑1 (MQ=255)
                                            gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                  <  1:1938398‑M2/46‑1 (MQ=255)
                                            gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                  >  1:2154290‑M2/20‑65 (MQ=255)
                                            gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                  <  1:1523227‑M2/46‑1 (MQ=255)
                                            gagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                  <  1:427901/65‑1 (MQ=255)
                                              gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                >  1:602001/1‑65 (MQ=255)
                                              gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                >  1:1673560/1‑65 (MQ=255)
                                               ccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttt               >  1:937626‑M2/17‑65 (MQ=255)
                                                         aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg     <  1:2349075‑M2/59‑1 (MQ=255)
                                                          accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa    <  1:1578320‑M2/60‑1 (MQ=255)
                                                          accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa    <  1:510620‑M2/60‑1 (MQ=255)
                                                          accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa    <  1:1466891‑M2/60‑1 (MQ=255)
                                                          accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa    <  1:785804‑M2/60‑1 (MQ=255)
                                                          accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa    <  1:863361‑M2/60‑1 (MQ=255)
                                                          accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa    <  1:2114827‑M2/60‑1 (MQ=255)
                                                          accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa    <  1:2684179‑M2/60‑1 (MQ=255)
                                                          accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa    <  1:1908569‑M2/60‑1 (MQ=255)
                                                           ccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAt   >  1:43207‑M2/5‑65 (MQ=255)
                                                            caaCGTTTATGAAAATGAAGAAATAACAAGCGAATGGCATCATTCCTGCTTTTACCAGGGGGAtt  >  1:1960132‑M2/4‑65 (MQ=255)
                                                               |                                                             
GAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATT  >  NC_000913/1209556‑1209680

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