New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =23 (0.550)11 (0.350) 9/96 0.8 27.4% coding (290/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 1209619 = 41 (1.310)pseudogene (37/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

GGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAA  >  NC_000913/1209555‑1209683
                                                                |                                                                
ggAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAAc                                                                  <  1:1473745/65‑1 (MQ=255)
        tGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGa                                                          >  1:2299092/1‑65 (MQ=255)
        tGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGa                                                          >  1:621090/1‑65 (MQ=255)
            ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAAt                                                      >  1:2863277/1‑65 (MQ=255)
             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                                     >  1:2037395/1‑65 (MQ=255)
             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                                     >  1:72706/1‑65 (MQ=255)
             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                                     >  1:1833586/1‑65 (MQ=255)
             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                                     >  1:1630555/1‑65 (MQ=255)
             gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                                     >  1:1671771/1‑65 (MQ=255)
                 taatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaaga                                                 <  1:118759‑M2/18‑1 (MQ=255)
                   atgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGaaa                                               >  1:3173815‑M2/46‑65 (MQ=255)
                    tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                              >  1:537673‑M2/45‑65 (MQ=255)
                     gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                             >  1:2107770‑M2/44‑65 (MQ=255)
                     gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                             >  1:2292389‑M2/44‑65 (MQ=255)
                     gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                             >  1:1879831‑M2/44‑65 (MQ=255)
                     aCCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATa                                             >  1:860162/1‑65 (MQ=255)
                          gACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAg                                        >  1:2656194/1‑65 (MQ=255)
                          gACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAg                                        >  1:2374380/1‑65 (MQ=255)
                          cgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAg                                        <  1:1824574‑M2/27‑1 (MQ=255)
                           acacCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGc                                       >  1:1495494/1‑65 (MQ=255)
                           acacCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGc                                       >  1:1083895/1‑65 (MQ=255)
                             acCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCaa                                     >  1:3021361/1‑65 (MQ=255)
                              gctgcattgagttttgctcattctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCaaa                                    >  1:2742533‑M2/35‑65 (MQ=255)
                                 aCAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATgg                                 <  1:1968614/65‑1 (MQ=255)
                                   attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                               >  1:987563‑M2/30‑65 (MQ=255)
                                   aGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                               >  1:868970/1‑65 (MQ=255)
                                     gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                             >  1:2878148/1‑65 (MQ=255)
                                     gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                             >  1:2761472/1‑65 (MQ=255)
                                     gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                             >  1:1073119/1‑65 (MQ=255)
                                     gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                             >  1:2134200/1‑65 (MQ=255)
                                      tAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                            <  1:2260433/65‑1 (MQ=255)
                                        gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:2015744‑M2/41‑1 (MQ=255)
                                        aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:1743336/65‑1 (MQ=255)
                                        aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:2600251/65‑1 (MQ=255)
                                         tGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc                         >  1:1285369/1‑65 (MQ=255)
                                         tGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc                         >  1:262828/1‑65 (MQ=255)
                                          gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                        >  1:2148940/1‑65 (MQ=255)
                                          gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                        >  1:2911199/1‑65 (MQ=255)
                                          gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                        >  1:473071/1‑65 (MQ=255)
                                          gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                        >  1:1789268/1‑65 (MQ=255)
                                          gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                        >  1:633184/1‑65 (MQ=255)
                                          gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                        >  1:13845/1‑65 (MQ=255)
                                           cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                       >  1:1792771/1‑65 (MQ=255)
                                           cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                       >  1:2878315/1‑65 (MQ=255)
                                           cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                       >  1:334186/1‑65 (MQ=255)
                                           cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                       >  1:2081904/1‑65 (MQ=255)
                                           cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                       >  1:1654311/1‑65 (MQ=255)
                                            agagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg                      >  1:2610355/1‑65 (MQ=255)
                                            agagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg                      >  1:2292862/1‑65 (MQ=255)
                                             gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     <  1:522455‑M2/46‑1 (MQ=255)
                                               gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                   >  1:3502/1‑65 (MQ=255)
                                               gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                   >  1:2958434/1‑65 (MQ=255)
                                                ccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttt                  >  1:3149680‑M2/17‑65 (MQ=255)
                                                 cttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtttt                 >  1:877740‑M2/16‑65 (MQ=255)
                                                  ttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa                <  1:118040‑M2/51‑1 (MQ=255)
                                                      cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAg            >  1:1581312‑M2/11‑65 (MQ=255)
                                                      cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAg            >  1:2893735‑M2/11‑65 (MQ=255)
                                                      cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAg            >  1:2257295‑M2/11‑65 (MQ=255)
                                                      cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAg            >  1:1395896‑M2/11‑65 (MQ=255)
                                                      cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAg            >  1:1923531‑M2/11‑65 (MQ=255)
                                                      cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGCCACCATTCCTGCTTTTACAAg            >  1:3140007‑M2/11‑65 (MQ=255)
                                                       ccaaagcaaCGTTTATGAAAATGAAGAAATAACGAGCAAATGGCATCATTCCTGCTTTTACCAgg           <  1:157798‑M2/56‑1 (MQ=255)
                                                         aaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgggg         <  1:1027654‑M2/58‑1 (MQ=255)
                                                          aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg        <  1:1935997‑M2/59‑1 (MQ=255)
                                                           accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:2777312‑M2/60‑1 (MQ=255)
                                                           accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:68454‑M2/60‑1 (MQ=255)
                                                           accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:2326378‑M2/60‑1 (MQ=255)
                                                           accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       >  1:768582‑M2/6‑65 (MQ=255)
                                                           accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:1778133‑M2/60‑1 (MQ=255)
                                                           accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:2966024‑M2/60‑1 (MQ=255)
                                                                cGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa  >  1:349333/1‑65 (MQ=255)
                                                                cGTTTACGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa  >  1:929889/1‑65 (MQ=255)
                                                                |                                                                
GGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAA  >  NC_000913/1209555‑1209683

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