New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 120780530 (1.060)9 (0.420) 6/100 1.1 34.9% coding (305/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 = 1209602 11 (0.510)pseudogene (18/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

CGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGAT  >  NC_000913/1207739‑1207871
                                                                  |                                                                  
cGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGa                                                                   <  2:990232/68‑1 (MQ=255)
     gcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCaaa                                                              <  1:271007‑M1/22‑1 (MQ=255)
       cGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAAc                                                            <  1:655776/68‑1 (MQ=255)
          tatCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattac                                                         <  2:234877‑M1/68‑12 (MQ=255)
          tatCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattac                                                         <  1:855102‑M1/68‑12 (MQ=255)
          tatCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCa                                                         <  1:740472/68‑1 (MQ=255)
               aaTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgtt                                                    <  2:176605‑M1/68‑17 (MQ=255)
                gccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgcag                                                   <  2:290212‑M1/33‑1 (MQ=255)
                 tCGGATGGCACTGTGCAAACGGCTCTCGAAATCCTTGGTTTGGGAGAAGGAGCAAAACTCAATgcagc                                                  <  1:162084/68‑1 (MQ=255)
                  cGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgcagca                                                 >  2:322696/1‑68 (MQ=255)
                  cGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgcagca                                                 >  2:356539/1‑68 (MQ=255)
                    gATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAAc                                               >  2:490462/1‑68 (MQ=255)
                     ttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACg                                              >  2:899982‑M1/31‑68 (MQ=255)
                      tgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                                             >  2:313637‑M1/30‑68 (MQ=255)
                      tGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                                             <  2:309049/68‑1 (MQ=255)
                        gCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCt                                           <  2:423969/68‑1 (MQ=255)
                        gCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCt                                           <  2:423970/68‑1 (MQ=255)
                        cttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCt                                           >  2:11637‑M1/28‑68 (MQ=255)
                         cACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccgg                                          >  2:757243‑M1/1‑42 (MQ=255)
                          tgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgca                                                    <  2:473440‑M1/32‑1 (MQ=255)
                          tgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgca                                                    >  1:473440‑M1/26‑57 (MQ=255)
                          aCTGTGCAAGCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                         >  1:663021/1‑68 (MQ=255)
                          aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggt                                         >  2:243145‑M1/1‑41 (MQ=255)
                          aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggt                                         >  1:245801‑M1/1‑41 (MQ=255)
                          aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtc                                              >  2:241931‑M1/1‑41 (MQ=255)
                          aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtc                                              <  1:241931‑M1/63‑23 (MQ=255)
                          aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                         >  2:600839/1‑68 (MQ=255)
                          aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                         >  2:11778/1‑68 (MQ=255)
                           gttatttcttcattttcataaacgTTGGTGTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACa                                        >  1:62502‑M1/25‑68 (MQ=255)
                           cTGTGCAAACGGCTCTCGAAATCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACa                                        <  2:241296/68‑1 (MQ=255)
                           cTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACa                                        <  1:822299/68‑1 (MQ=255)
                           cTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACa                                        <  1:405505/68‑1 (MQ=255)
                                 tcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGAAACGGCTACATTAGGa                                  <  2:665763‑M1/50‑1 (MQ=255)
                                  aaCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGAc                                 >  1:720978/1‑68 (MQ=255)
                                    cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGc                               >  2:158224/1‑68 (MQ=255)
                                    cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGc                               >  2:218522/1‑68 (MQ=255)
                                     ggCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCa                              <  1:922779/68‑1 (MQ=255)
                                       ttttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCAcc                            <  1:242841‑M1/56‑1 (MQ=255)
                                          tcGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGt                         <  1:158224/68‑1 (MQ=255)
                                          tcGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGt                         <  2:51069/68‑1 (MQ=255)
                                           cGAATACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCATCGGCTACATTAGGACGCACCGGtt                        <  1:483207/68‑1 (MQ=255)
                                           cGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagcca                        <  1:97562‑M1/68‑45 (MQ=255)
                                            gAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGttt                       <  2:805170/68‑1 (MQ=255)
                                             aaaaCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTc                      >  1:904632/1‑68 (MQ=255)
                                                cccTTGGTATGGGAGAAGGAGCAAAACTCTATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                   <  1:832245/67‑1 (MQ=255)
                                                acgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                   >  1:364878‑M1/4‑68 (MQ=255)
                                                acgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                   >  2:391783‑M1/4‑68 (MQ=255)
                                                acgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                   >  1:665763‑M1/4‑68 (MQ=255)
                                                aCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                   <  2:335937/68‑1 (MQ=255)
                                                aCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                   <  1:279706/68‑1 (MQ=255)
                                                aCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                   <  2:727468/68‑1 (MQ=255)
                                                aCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                   >  2:703090/1‑68 (MQ=255)
                                                 cgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAg                  <  2:566166‑M1/66‑1 (MQ=255)
                                                  gTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                   >  1:312586‑M1/2‑66 (MQ=255)
                                                  gTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                   <  2:312586‑M1/65‑1 (MQ=255)
                                                     GGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCtat              >  1:326444‑M1/1‑68 (MQ=255)
                                                      gcTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCtata             >  1:655402‑M1/3‑68 (MQ=255)
                                                        ttGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacagg           <  1:988695‑M1/68‑58 (MQ=255)
                                                          GGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCgg         >  2:242841‑M1/1‑68 (MQ=255)
                                                            GAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                   >  1:73803‑M1/1‑56 (MQ=255)
                                                            GAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                   >  2:365154‑M1/1‑56 (MQ=255)
                                                            GAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                   <  2:73803‑M1/56‑1 (MQ=255)
                                                            GAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                   <  1:365154‑M1/56‑1 (MQ=255)
                                                                aGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctga   <  1:854668‑M1/68‑66 (MQ=255)
                                                                 GGAGCAAATCTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGAt  >  2:319571‑M1/1‑68 (MQ=255)
                                                                  |                                                                  
CGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGAT  >  NC_000913/1207739‑1207871

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