New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =6 (0.140)40 (1.280) 21/96 0.1 82.9% coding (290/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 1209619 = 12 (0.380)pseudogene (37/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

GAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAA  >  NC_000913/1209556‑1209683
                                                               |                                                                
gAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGGAATGCCGAGCCTTCTCCCAAACCAACg                                                                 >  1:2009316/1‑65 (MQ=255)
    accggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTa                                                             >  1:1111604‑M2/60‑65 (MQ=255)
    accggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTa                                                             >  1:1549530‑M2/60‑65 (MQ=255)
         tgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaa                                                        >  1:919607‑M2/55‑65 (MQ=255)
            gtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATg                                                     >  1:2512115‑M2/52‑65 (MQ=255)
            gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                                     >  1:209132/1‑65 (MQ=255)
            gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                                     >  1:172297/1‑65 (MQ=255)
            gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCCTCTCCCAAACAAACGTTTATGAAAATg                                                     >  1:2427509/1‑65 (MQ=255)
               ctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaag                                                  >  1:862889‑M2/49‑65 (MQ=255)
               ctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaag                                                  <  1:704716‑M2/17‑1 (MQ=255)
                taatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaaga                                                 <  1:1109424‑M2/18‑1 (MQ=255)
                   tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                              >  1:1119380‑M2/45‑65 (MQ=255)
                   tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                              >  1:53844‑M2/45‑65 (MQ=255)
                   tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                              >  1:1158679‑M2/45‑65 (MQ=255)
                   tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                              >  1:1199010‑M2/45‑65 (MQ=255)
                   tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTGATGAAAATGAAGAAAt                                              >  1:2957562‑M2/45‑65 (MQ=255)
                    gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                             >  1:42486‑M2/44‑65 (MQ=255)
                    gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                             >  1:161161‑M2/44‑65 (MQ=255)
                        ccgttgctgcattgagttttgctccttctcccaaaccgaCGTTTATGAAAATGAAGAAATAACaa                                         >  1:3135869‑M2/40‑65 (MQ=255)
                        ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                         >  1:3320272‑M2/40‑65 (MQ=255)
                        ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                         >  1:1715431‑M2/40‑65 (MQ=255)
                        ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                         >  1:2665451‑M2/40‑65 (MQ=255)
                         gACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAg                                        >  1:566019/1‑65 (MQ=255)
                         cgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAg                                        >  1:1126302‑M2/39‑65 (MQ=255)
                             gctgcattgagttttgctccttctccccaaccaaCGTTTATGAAAATGAAGAAATAACAAGCaaa                                    >  1:2651243‑M2/35‑65 (MQ=255)
                             gctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCaaa                                    <  1:631286‑M2/31‑1 (MQ=255)
                             gctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCaaa                                    >  1:3318047‑M2/35‑65 (MQ=255)
                             ccAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCaaa                                    >  1:1245601/1‑65 (MQ=255)
                              ctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAAt                                   >  1:1238194‑M2/34‑65 (MQ=255)
                               tgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAGGAAATAACAAGCAAATg                                  >  1:1148200‑M2/33‑65 (MQ=255)
                               aaCAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATg                                  <  1:3124948/65‑1 (MQ=255)
                                 cattgagttttgctccttctcccaaaccagCGTTTATGAAAATGAAGAAATAACAAGCAAATGGc                                <  1:1519820‑M2/35‑1 (MQ=255)
                                  attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                               <  1:603685‑M2/36‑1 (MQ=255)
                                    gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                             >  1:227719/1‑65 (MQ=255)
                                     gagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                            <  1:2752647‑M2/39‑1 (MQ=255)
                                       gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:3065946‑M2/41‑1 (MQ=255)
                                       gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:109007‑M2/41‑1 (MQ=255)
                                       gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:51312‑M2/41‑1 (MQ=255)
                                       gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          >  1:35357‑M2/25‑65 (MQ=255)
                                       gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          >  1:931829‑M2/25‑65 (MQ=255)
                                       gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:701533‑M2/41‑1 (MQ=255)
                                       gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:2737547‑M2/41‑1 (MQ=255)
                                       gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          >  1:2245877‑M2/25‑65 (MQ=255)
                                       cttttgctccttctcccaaaccaaCGTTTATGAAAAGGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:17864‑M2/41‑1 (MQ=255)
                                       aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          <  1:2543502/65‑1 (MQ=255)
                                         tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                        >  1:1705866‑M2/23‑65 (MQ=255)
                                          cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                       >  1:2811523/1‑65 (MQ=255)
                                            gctccttctcccagaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     <  1:104959‑M2/46‑1 (MQ=255)
                                            gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     <  1:2481611‑M2/46‑1 (MQ=255)
                                            gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     <  1:2486191‑M2/46‑1 (MQ=255)
                                            gccccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     <  1:1503155‑M2/46‑1 (MQ=255)
                                              gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCCAATGGCATCATTCCTGCtt                   >  1:1689408/1‑65 (MQ=255)
                                              gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                   >  1:1835006/1‑65 (MQ=255)
                                                      ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgg           >  1:257210‑M2/10‑65 (MQ=255)
                                                      ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgg           <  1:788251‑M2/56‑1 (MQ=255)
                                                         aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg        <  1:2768388‑M2/59‑1 (MQ=255)
                                                          accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:2929412‑M2/60‑1 (MQ=255)
                                                          accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:3072010‑M2/60‑1 (MQ=255)
                                                          accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:1428576‑M2/60‑1 (MQ=255)
                                                          accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:1227884‑M2/60‑1 (MQ=255)
                                                            caaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAtt     <  1:2514773‑M2/62‑1 (MQ=255)
                                                             aaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt    <  1:3051557‑M2/63‑1 (MQ=255)
                                                               cGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa  <  1:396564/65‑1 (MQ=255)
                                                               |                                                                
GAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAA  >  NC_000913/1209556‑1209683

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