New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 12078052 (0.080)12 (0.620) 8/100 0.7 87.2% 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

CGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGAT  >  NC_000913/1207739‑1207871
                                                                  |                                                                  
cGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGc                                                                   <  1:870156‑M1/68‑2 (MQ=255)
   ttGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctct                                                                <  1:308699‑M1/68‑5 (MQ=255)
   ccggcgATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctct                                                                <  2:952959‑M1/66‑5 (MQ=255)
     gcaggaatgatgccatttgcttgttatttcttcattttcatatacgTTGGTTTGGGAGAAGGAGCaaa                                                              <  1:413091‑M1/22‑1 (MQ=255)
     gcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCaaa                                                              >  1:791000‑M1/47‑68 (MQ=255)
      gCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgca                                                             <  1:286905‑M1/68‑8 (MQ=255)
       aggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAAc                                                            >  2:779485‑M1/45‑68 (MQ=255)
                gccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgcag                                                   <  2:774351‑M1/33‑1 (MQ=255)
                 tCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttgg                                                  >  1:927169‑M1/1‑50 (MQ=255)
                     tagcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACg                                              >  1:760305‑M1/31‑68 (MQ=255)
                     aTGGCACTGTGCAAACGGCTCTCGAAAACCTTGATTTGGGAGAAGGctctgcattacctgtaggtgtc                                              >  2:914534‑M1/1‑46 (MQ=255)
                          tgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                         >  2:888533‑M1/26‑68 (MQ=255)
                          tgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                         >  2:364059‑M1/26‑68 (MQ=255)
                           gttatttcttcattttcataaacgTTGGCTTGGGAGAAGGAGCAACACTCAATGCAGCAACGGCTACa                                        >  2:363405‑M1/25‑68 (MQ=255)
                               gCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgt                                    >  1:220365‑M1/1‑36 (MQ=255)
                                    cGGCTCTCGAAGACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggcct                               >  1:890623‑M1/1‑31 (MQ=255)
                                    cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggcct                               >  1:961471‑M1/1‑31 (MQ=255)
                                    cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggcct                               >  1:740446‑M1/1‑31 (MQ=255)
                                    cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggcct                               >  1:596900‑M1/1‑31 (MQ=255)
                                    cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGc                               >  2:447649/1‑68 (MQ=255)
                                      gCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttc                             >  1:703164‑M1/1‑29 (MQ=255)
                                          tcatatacgCTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACACTAGGACGCACCGGt                         >  2:271449‑M1/10‑68 (MQ=255)
                                          tcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGt                         >  2:382300‑M1/10‑68 (MQ=255)
                                          tcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGt                         >  2:959556‑M1/10‑68 (MQ=255)
                                           cataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGtt                        >  1:559532‑M1/9‑68 (MQ=255)
                                           cataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGtt                        >  1:522512‑M1/9‑68 (MQ=255)
                                                acgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                   >  1:716228‑M1/4‑68 (MQ=255)
                                                 ccTTTGTTTGGTAGAAGGATCAATACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAg                  <  1:119647/68‑1 (MQ=255)
                                                  gTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGc                 >  1:89645‑M1/2‑68 (MQ=255)
                                                  gTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGc                 <  1:79323‑M1/67‑1 (MQ=255)
                                                  gTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGc                 <  1:166608‑M1/67‑1 (MQ=255)
                                                          GGGAGAAGGAGCAAAACTCATTGCTGCAACGTCTACTTTAGGACTCACCGGTTTCATAGCTATACCgg         >  2:79323‑M1/1‑68 (MQ=255)
                                                            gagaAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctgg       <  1:749196‑M1/68‑62 (MQ=255)
                                                            gagaAGGccctgcattacctgttggcgtcccggttccgtggccttcagccactccgccaacaggctgg       <  2:808329‑M1/68‑62 (MQ=255)
                                                             agaAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggcaggc      >  2:234010‑M1/1‑6 (MQ=255)
                                                             AGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTa      <  1:897692‑M1/68‑1 (MQ=255)
                                                                aGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctga   <  2:200506‑M1/68‑66 (MQ=255)
                                                                 GGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGAt  <  1:567068‑M1/68‑1 (MQ=255)
                                                                  |                                                                  
CGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGAT  >  NC_000913/1207739‑1207871

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