New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 12078058 (0.310)21 (1.080) 15/100 0.2 84.0% coding (305/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 = 1209602 2 (0.100)pseudogene (18/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

CGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGA  >  NC_000913/1207739‑1207870
                                                                  |                                                                 
cGTTTGGCGATATCAAATCGGATGGCACTGTTCATACGGCGCTCGAAAACCTTGGTTTGTGAGAAGGc                                                                  <  1:761751‑M1/68‑2 (MQ=255)
cGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGc                                                                  <  1:39045‑M1/68‑2 (MQ=255)
   ttGGCGATATCATATCGGATGGCACTGTGCAAACGGCTCTCGAATACCTTGGTTTGGGAGAAGGctct                                                               <  1:905777‑M1/68‑5 (MQ=255)
    agcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCaa                                                              <  1:814122‑M1/21‑1 (MQ=255)
     ggCGAAATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCaaa                                                             <  2:68575/68‑1 (MQ=255)
     gcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCaaa                                                             <  2:417835‑M1/22‑1 (MQ=255)
     gcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCaaa                                                             <  2:501315‑M1/22‑1 (MQ=255)
      gCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGcgctgca                                                            <  2:326856‑M1/68‑8 (MQ=255)
       aggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAAc                                                           >  2:846001‑M1/45‑68 (MQ=255)
          aatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCa                                                        >  1:533247‑M1/42‑68 (MQ=255)
            tgatgccatttgcttgttatttcttcattttcataatcgTTGGTTTGGGAGAATGAGCAAAACTCAAt                                                      >  1:932599‑M1/40‑68 (MQ=255)
            tCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacct                                                      <  1:96534‑M1/68‑14 (MQ=255)
            tCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGG                                                                   <  2:178271‑M1/55‑1 (MQ=255)
            tCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGG                                                                   >  1:178271‑M1/1‑55 (MQ=255)
                    gATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgt                                              <  1:783218‑M1/68‑22 (MQ=255)
                     aTGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtc                                             <  2:655128‑M1/68‑23 (MQ=255)
                     aTGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACg                                             <  2:556345/68‑1 (MQ=255)
                      tGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                                            <  2:706501/68‑1 (MQ=255)
                       gcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGc                                           <  2:408316‑M1/40‑1 (MQ=255)
                          tgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                        >  1:660697‑M1/26‑68 (MQ=255)
                          aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggt                                        >  1:346235‑M1/1‑41 (MQ=255)
                          aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggcgtcccggc                                        <  1:948777‑M1/68‑28 (MQ=255)
                          aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                        >  2:639207/1‑68 (MQ=255)
                          aCTGTGCAAACGGCTCTCGAAAACCTTGGATTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                        >  1:100615/1‑68 (MQ=255)
                           gttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACa                                       >  1:773043‑M1/25‑68 (MQ=255)
                           gttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACa                                       >  2:54150‑M1/25‑68 (MQ=255)
                           cTGTGCAAACGGCTCTCGAAAACCTTGTTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACa                                       <  1:639862/68‑1 (MQ=255)
                            tgtgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttc                                      <  1:415919‑M1/68‑30 (MQ=255)
                            tgtgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttc                                      <  1:568213‑M1/68‑30 (MQ=255)
                              atttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTa                                    <  1:601911‑M1/47‑1 (MQ=255)
                              attacttcattttcataaacgTAGGTTTGGGAGAAGGAGCAAAACTCAATGCATCAACGGCTACattt                                    <  2:738098‑M1/47‑2 (MQ=255)
                               gCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgt                                   >  2:712908‑M1/1‑36 (MQ=255)
                                  aaCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggc                                <  2:914998‑M1/68‑36 (MQ=255)
                                    cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggcct                              >  2:556088‑M1/1‑31 (MQ=255)
                                    cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgca                                                            <  2:521623‑M1/38‑8 (MQ=255)
                                    cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgca                                                            >  1:521620‑M1/1‑31 (MQ=255)
                                    cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgc                                                             <  1:143646‑M1/37‑7 (MQ=255)
                                    cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgc                                                             >  2:143646‑M1/1‑31 (MQ=255)
                                    cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGc                              >  1:762840/1‑68 (MQ=255)
                                           cataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGtt                       >  1:201851‑M1/9‑68 (MQ=255)
                                             aaaaCCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccact                     >  1:210080‑M1/1‑22 (MQ=255)
                                              aaaCCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactc                    <  1:281986‑M1/68‑48 (MQ=255)
                                                aCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                  <  2:762844/68‑1 (MQ=255)
                                                  cTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgcc                <  1:470606‑M1/68‑52 (MQ=255)
                                                     GGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCtat             >  2:703857‑M1/1‑68 (MQ=255)
                                                          GGGAGAAGGAGCAAAACTCAATGCAGCAACGTCTACATTAGGTCGCACCGGTTTCATAGCTTTACcg       >  1:305644‑M1/1‑67 (MQ=255)
                                                          GGGAGAAGGAGCAAAACTCAATGCAGCAACGTCTACATTAGGACGCACCGGTTTCATAGCTATACCgg        >  2:254429‑M1/1‑68 (MQ=255)
                                                          GGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCgg        >  2:100615‑M1/1‑68 (MQ=255)
                                                          GGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCgg        >  1:316165‑M1/1‑68 (MQ=255)
                                                           GGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGt       >  1:262036‑M1/1‑68 (MQ=255)
                                                            gagaAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctgg      <  2:918567‑M1/68‑62 (MQ=255)
                                                              gaAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggct    <  1:231370‑M1/68‑64 (MQ=255)
                                                              gaAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggct    <  1:847008‑M1/68‑64 (MQ=255)
                                                                aGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctga  <  1:786727‑M1/68‑66 (MQ=255)
                                                                aGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctga  <  2:210080‑M1/68‑66 (MQ=255)
                                                                  gctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggc     <  2:61711‑M1/63‑63 (MQ=255)
                                                                  gctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggc     >  1:61711‑M1/1‑1 (MQ=255)
                                                                  |                                                                 
CGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGA  >  NC_000913/1207739‑1207870

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