New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =2 (0.080)17 (0.870) 11/100 0.4 90.6% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 1209619 = 2 (0.100)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

AGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGC  >  NC_000913/1207726‑1207856
                                                                |                                                                  
aGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGG                                                                 <  1:9670‑M1/68‑1 (MQ=255)
 gCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGT                                                                <  1:755579‑M1/68‑1 (MQ=255)
 aaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGt                                                                <  1:888536‑M1/5‑1 (MQ=255)
 aaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGt                                                                >  2:857200‑M1/64‑68 (MQ=255)
    tccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTg                                                             <  2:444747‑M1/8‑1 (MQ=255)
    tccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTg                                                             <  2:522512‑M1/8‑1 (MQ=255)
             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:740446‑M1/1‑31 (MQ=255)
                                                 cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggcct                >  1:596900‑M1/1‑31 (MQ=255)
                                                 cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggcct                >  1:961471‑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:959556‑M1/10‑68 (MQ=255)
                                                       tcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGt          >  2:382300‑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)
                                                                |                                                                  
AGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGC  >  NC_000913/1207726‑1207856

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.