New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =1 (0.040)16 (0.950) 12/100 0.2 84.8% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 1209619 = 5 (0.300)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

AGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATA  >  NC_000913/1207726‑1207854
                                                                |                                                                
aGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGG                                                               >  2:199853‑M1/1‑68 (MQ=255)
aGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGG                                                               >  1:392250‑M1/1‑68 (MQ=255)
 gCCGCAAGAATCCGTTTGGCGATATCAAACCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGT                                                              <  1:156812‑M1/68‑1 (MQ=255)
 aaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGt                                                              >  2:410647‑M1/64‑68 (MQ=255)
 aaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGt                                                              <  1:535775‑M1/5‑1 (MQ=255)
 aaatccccctggtaaaagcaggaatgatgccatttgcttgttattacttcattttcataaacgTTGGt                                                              <  2:708366‑M1/5‑1 (MQ=255)
   atccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGttt                                                            <  2:481275‑M1/7‑1 (MQ=255)
    tccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTg                                                           >  1:162090‑M1/61‑68 (MQ=255)
             taaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGa                                                  >  1:143550‑M1/52‑68 (MQ=255)
             cGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGc                                                  <  1:396815‑M1/68‑2 (MQ=255)
               aaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGc                                                >  1:483419‑M1/50‑68 (MQ=255)
                  gcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCaaa                                             <  2:727137‑M1/22‑1 (MQ=255)
                   gCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgca                                            <  2:135711‑M1/68‑8 (MQ=255)
                       aatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCa                                        >  2:203379‑M1/42‑68 (MQ=255)
                         tCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGT                                                              <  2:761408‑M1/44‑1 (MQ=255)
                         tCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGT                                                              >  1:761407‑M1/1‑44 (MQ=255)
                            cctcgGACGGCACTGTGCAAACGGCTCTCGACTCCCTTGGTTTGGGAGAAGGctctgcactacctgtc                                   <  1:568255‑M1/66‑17 (MQ=255)
                             gccatttgcttgttatatcttcattatcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgcag                                  <  1:545463‑M1/33‑1 (MQ=255)
                              ccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgcagc                                 <  1:394680‑M1/34‑1 (MQ=255)
                               cGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggt                                >  2:791814‑M1/1‑49 (MQ=255)
                                     cttgttatttcttcattttcataaacgTTGGTTTGGGCGAAGGAGCAAAACTCAATGCAGCAACGGCt                          <  1:239166‑M1/41‑1 (MQ=255)
                                       tgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGgcttc                        <  2:370289‑M1/43‑3 (MQ=255)
                                       tgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                        >  1:19312‑M1/26‑68 (MQ=255)
                                       tgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                        >  1:241086‑M1/26‑68 (MQ=255)
                                       aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggt                        >  2:450877‑M1/1‑41 (MQ=255)
                                        gttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACa                       >  1:774806‑M1/25‑68 (MQ=255)
                                         ttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCa                                        >  2:285508‑M1/24‑50 (MQ=255)
                                         ttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCa                                        <  1:285507‑M1/27‑1 (MQ=255)
                                               aaCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggc                <  1:243501‑M1/68‑36 (MQ=255)
                                                 cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggcct              >  2:652586‑M1/1‑31 (MQ=255)
                                                 cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggcct              >  1:663225‑M1/1‑31 (MQ=255)
                                                 cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggcct              >  2:189925‑M1/1‑31 (MQ=255)
                                                 cGGCTCTCGAAAAACTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGc              >  2:312166/1‑68 (MQ=255)
                                                       tcataaacgTTGGCTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGt        >  1:204811‑M1/10‑68 (MQ=255)
                                                        cGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccc       >  2:76651‑M1/1‑24 (MQ=255)
                                                             acgttGGTTTTGGAGAATGATCAATACTCATTGCATCAACGGCTACATTTGGACGCACCGTTTTCTTa  <  2:109263/65‑1 (MQ=255)
                                                              cgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGt        >  2:294395‑M1/3‑61 (MQ=255)
                                                              cgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGt        <  1:294394‑M1/59‑1 (MQ=255)
                                                                |                                                                
AGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATA  >  NC_000913/1207726‑1207854

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.