New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 120780525 (0.780)8 (0.330) 8/100 1.0 32.8% coding (305/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 = 1209602 14 (0.580)pseudogene (18/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

CGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATT  >  NC_000913/1207739‑1207872
                                                                  |                                                                   
taaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGa                                                                    <  1:1181027‑M1/17‑1 (MQ=255)
taaaagcaggaatgatgcaatttgcttgttatttcatcattttcataaacgTTGGTTTGGGAGAaggg                                                                    <  2:54956‑M1/17‑2 (MQ=255)
cGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGc                                                                    <  1:226293‑M1/68‑2 (MQ=255)
     ggCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAGCCTTGGTTTGGGAGAAGGctctgc                                                               >  2:130682‑M1/1‑62 (MQ=255)
     ggCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgc                                                               <  2:252371‑M1/68‑7 (MQ=255)
       cGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAAc                                                             <  2:1095458/68‑1 (MQ=255)
       cGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAAc                                                             <  2:1096846/68‑1 (MQ=255)
       cGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAAc                                                             <  1:13356/68‑1 (MQ=255)
               tgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgca                                                     <  1:1032170‑M1/32‑1 (MQ=255)
               aaTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgtt                                                     <  2:79671‑M1/68‑17 (MQ=255)
                gccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgcag                                                    <  1:1049366‑M1/33‑1 (MQ=255)
                     aTGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACg                                               <  2:1103664/68‑1 (MQ=255)
                     aTGGCACTGTGCAAACGGCGCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACg                                               <  2:1135403/68‑1 (MQ=255)
                      tGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                                              <  1:687775/68‑1 (MQ=255)
                        tgaCTTTGCAATCGGCTCTCGAAAACCTTTTTTTGGTAGAAGGAGCATTACTCAATGCTTCAACGGCt                                            <  2:751864/66‑1 (MQ=16)
                        gCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCt                                            <  2:51847/68‑1 (MQ=255)
                          tgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                          >  2:1181025‑M1/26‑68 (MQ=255)
                          tgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                          >  2:944131‑M1/26‑68 (MQ=255)
                          tgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                          >  1:526511‑M1/26‑68 (MQ=255)
                          aCTGTGCAATCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                          <  2:752595/68‑1 (MQ=255)
                          aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggt                                          >  1:444574‑M1/1‑41 (MQ=255)
                          aCTGTGCAAACGGCTCTCGAAAACCTAGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                          <  2:490699/68‑1 (MQ=255)
                          aCTGTGAAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                          >  1:1040334/1‑68 (MQ=255)
                           gttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACa                                         >  2:1049365‑M1/25‑68 (MQ=255)
                           gttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACa                                         >  2:1032169‑M1/25‑68 (MQ=255)
                           ctttgCAAACGGCGCTCGAAAACCTTTGTTAGGGAGATGGAGCAAAACTCAATGCAGCAACGGCTACa                                         <  2:374362/65‑1 (MQ=37)
                           cTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACa                                         <  1:515854/68‑1 (MQ=255)
                           cTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACa                                         <  1:439090/68‑1 (MQ=255)
                           cTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACa                                         <  1:1178312/68‑1 (MQ=255)
                           cTGTGCAAACGGCTCTCGAAAACCTTGGCTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACa                                         <  2:338183/68‑1 (MQ=255)
                               gCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgt                                     <  1:804416‑M1/68‑33 (MQ=255)
                                    ctgCTCTCGAAAACCTTGGTTAGGGAGAAGGAGCAAAACTCAAAGAAGCAACGGCTACAATAGGACGc                                >  1:1092429/3‑68 (MQ=25)
                                    cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggcct                                >  2:187729‑M1/1‑31 (MQ=255)
                                    cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgcagc                                                   <  1:748369/49‑1 (MQ=255)
                                    cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgcagc                                                   >  2:748369/1‑49 (MQ=255)
                                    cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGc                                >  1:53103/1‑68 (MQ=255)
                                     cattttcataaacgTTGGTTTGTGAGAAGGAGCAAAACTCAATGCAGCATCGGCTACATTAGGACGCa                               >  2:136758‑M1/15‑68 (MQ=255)
                                           cGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCTATTCAGCAACGGCTACATTAGGACGCACCGGtt                         <  1:767438/68‑1 (MQ=255)
                                             aaaaCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGttcc                       <  2:713504/68‑3 (MQ=255)
                                                acgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                    >  1:981976‑M1/4‑68 (MQ=255)
                                                aCCTTGTTTTGGGATAATGAGCAAAACTCAATGCATCAACGGCTACATTATGACGCACCGGTTTCAtt                    <  1:232435/68‑2 (MQ=255)
                                                 cgTTGGTTTGGGAGAAGGAGCAAAACTCAATgcagc                                                   >  1:626024‑M1/3‑36 (MQ=255)
                                                 cgTTGGTTTGGGAGAAGGAGCAAAACTCAATgcagc                                                   <  2:626024‑M1/34‑1 (MQ=255)
                                                  cTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgcc                  <  1:734800‑M1/68‑52 (MQ=255)
                                                  cTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGc                  <  2:923868/68‑1 (MQ=255)
                                                  cTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGc                  >  1:1021279/1‑68 (MQ=255)
                                                    tGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaa                <  1:112880‑M1/68‑54 (MQ=255)
                                                         TGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCg           >  1:627272‑M1/1‑68 (MQ=255)
                                                         TGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCg           >  1:625669‑M1/1‑68 (MQ=255)
                                                          gtGAGAAGGAGCAAAACTCAATGCAGCTACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCgg          >  2:1042648‑M1/3‑68 (MQ=255)
                                                          GGGAGAAGGATCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCgg          >  2:1070380‑M1/1‑68 (MQ=255)
                                                           GGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGt         >  1:93656‑M1/1‑68 (MQ=255)
                                                             AGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCGCCGGTTTCATAGCTATACCGGTTa       <  2:340656‑M1/68‑1 (MQ=255)
                                                              gaAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggct      <  1:113564‑M1/68‑64 (MQ=255)
                                                                 GGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGAt   <  2:239552‑M1/68‑1 (MQ=255)
                                                                 GGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGAt   >  1:477250‑M1/1‑68 (MQ=255)
                                                                  gctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaa  >  1:1117229‑M1/1‑1 (MQ=255)
                                                                  gctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaa  >  1:929837‑M1/1‑1 (MQ=255)
                                                                  |                                                                   
CGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATT  >  NC_000913/1207739‑1207872

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