New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =34 (0.820)13 (0.410) 9/102 0.9 26.8% coding (290/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 1209619 = 45 (1.420)pseudogene (37/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

GAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAA  >  NC_000913/1209556‑1209683
                                                               |                                                                
gAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGttt                                                              >  1:2845056/1‑68 (MQ=255)
gAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGGAATGCAGAGCCTTCTCCCAAACCAACGttt                                                              >  1:383592/1‑68 (MQ=255)
         aaGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                                     >  1:3161701/1‑68 (MQ=255)
         aaGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGACTTCTCCCAAACCAACGTTTATGAAAATg                                                     >  1:299983/1‑68 (MQ=255)
            gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaag                                                  >  1:752690/1‑68 (MQ=255)
            gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaag                                                  >  1:750868/1‑68 (MQ=255)
            gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaag                                                  >  1:2771249/1‑68 (MQ=255)
            gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaag                                                  >  1:2633246/1‑68 (MQ=255)
            gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaag                                                  >  1:1890832/1‑68 (MQ=255)
             ccACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaaga                                                 >  1:782220/1‑68 (MQ=255)
               aCGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGaaa                                               <  1:2291131/68‑1 (MQ=255)
                   tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAc                                           >  1:679757‑M2/45‑68 (MQ=255)
                   tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAc                                           >  1:1935830‑M2/45‑68 (MQ=255)
                        ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCa                                      >  1:929459‑M2/40‑68 (MQ=255)
                        ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCa                                      >  1:1656825‑M2/40‑68 (MQ=255)
                         gACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCaa                                     >  1:2770601/1‑68 (MQ=255)
                         cgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCaa                                     >  1:2795696‑M2/39‑68 (MQ=255)
                            acCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATg                                  >  1:2737733/1‑68 (MQ=255)
                             ccAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATgg                                 <  1:91937/68‑1 (MQ=255)
                              cAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGc                                <  1:717430/68‑1 (MQ=255)
                               aaCAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                               <  1:3109504/68‑1 (MQ=255)
                               aaCAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                               >  1:1575714/1‑68 (MQ=255)
                                  attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                            >  1:1442997‑M2/30‑68 (MQ=255)
                                  aGGTTATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                            <  1:2613670/68‑1 (MQ=255)
                                  aGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatna                            >  1:2354247/1‑68 (MQ=255)
                                  aGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                            >  1:2022981/1‑68 (MQ=255)
                                    tgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          >  1:2663567‑M2/28‑68 (MQ=255)
                                    gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          >  1:1486490/1‑68 (MQ=255)
                                    gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          >  1:1939127/1‑68 (MQ=255)
                                    gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          >  1:2392796/1‑68 (MQ=255)
                                    gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          >  1:1236635/1‑68 (MQ=255)
                                    gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          >  1:668449/1‑68 (MQ=255)
                                    gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                          >  1:665211/1‑68 (MQ=255)
                                       gttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                       >  1:2013320‑M2/25‑68 (MQ=255)
                                        tGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg                      >  1:1149920/1‑68 (MQ=255)
                                        tGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg                      >  1:1061699/1‑68 (MQ=255)
                                         tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCGTGc                     >  1:1773364‑M2/23‑68 (MQ=255)
                                         tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     >  1:2796000‑M2/23‑68 (MQ=255)
                                         tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     <  1:2331553‑M2/46‑1 (MQ=255)
                                         tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     >  1:3269219‑M2/23‑68 (MQ=255)
                                         gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     <  1:461558/68‑1 (MQ=255)
                                         gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     >  1:333336/1‑68 (MQ=255)
                                         gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     >  1:780357/1‑68 (MQ=255)
                                         gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     <  1:1957480/68‑1 (MQ=255)
                                         gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                     <  1:271835/68‑1 (MQ=255)
                                          cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCt                    >  1:3131800/1‑68 (MQ=255)
                                          cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCt                    >  1:231634/1‑68 (MQ=255)
                                          cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCt                    >  1:2702082/1‑68 (MQ=255)
                                          cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCt                    >  1:1028346/1‑68 (MQ=255)
                                          cAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCt                    >  1:3005712/1‑68 (MQ=255)
                                           agagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                   >  1:2770644/1‑68 (MQ=255)
                                           agagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                   >  1:1754738/1‑68 (MQ=255)
                                           agagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                   >  1:2338501/1‑68 (MQ=255)
                                            gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttt                  >  1:1515277‑M2/20‑68 (MQ=255)
                                              gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa                >  1:2146733/1‑68 (MQ=255)
                                              gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa                >  1:2220040/1‑68 (MQ=255)
                                              gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa                >  1:875367/1‑68 (MQ=255)
                                              gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa                >  1:2891243/1‑68 (MQ=255)
                                                 ttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCa             <  1:1176465‑M2/54‑1 (MQ=255)
                                                   ctcccaaaccaaCGGTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgg           <  1:2582785‑M2/56‑1 (MQ=255)
                                                   ctcccaaaccaaCCTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgg           >  1:427893‑M2/13‑68 (MQ=255)
                                                    tcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg          >  1:3218361‑M2/12‑68 (MQ=255)
                                                     cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgggg         <  1:328033‑M2/58‑1 (MQ=255)
                                                      ctaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg        <  1:374550‑M2/59‑1 (MQ=255)
                                                       caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:1275736‑M2/60‑1 (MQ=255)
                                                       caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:776469‑M2/60‑1 (MQ=255)
                                                       caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:738824‑M2/60‑1 (MQ=255)
                                                       caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa       <  1:3359645‑M2/60‑1 (MQ=255)
                                                          accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt    >  1:329907‑M2/6‑68 (MQ=255)
                                                            caaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa  >  1:1585870‑M2/4‑68 (MQ=255)
                                                               |                                                                
GAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAA  >  NC_000913/1209556‑1209683

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