New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =30 (1.360)8 (0.480) 8/100 0.5 29.8% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 1209619 = 15 (0.900)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

AGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCT  >  NC_000913/1207726‑1207857
                                                                |                                                                   
taaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTgg                                                                  >  2:278480‑M1/65‑68 (MQ=255)
aGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTgg                                                                  >  1:28979/1‑68 (MQ=255)
 gCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGt                                                                 <  2:274906/68‑1 (MQ=255)
    tccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTg                                                              <  1:663208‑M1/8‑1 (MQ=255)
    tccacctggtaaaagaaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTg                                                              <  1:240999‑M1/8‑1 (MQ=255)
            gtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAgg                                                      >  2:730016‑M1/53‑68 (MQ=255)
             taatagcatgaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGgagaatgc                                                     >  2:434252‑M1/52‑65 (MQ=255)
             taaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGa                                                     >  1:6181‑M1/52‑68 (MQ=255)
             cGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGa                                                     <  2:305021/68‑1 (MQ=255)
               tttGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGc                                                   <  2:666586/68‑1 (MQ=255)
                    cGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAAc                                              <  2:687729/68‑1 (MQ=255)
                        atgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCaa                                          <  1:200425‑M1/28‑1 (MQ=255)
                         tgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAAt                                         >  1:785737‑M1/40‑68 (MQ=255)
                                  ctGGCACTGTGCAATCGGCTCTCGAATACCTTGGTTTGGGAGAAGGTGCAAAACTCAATGCAGCAACg                                <  1:550931/67‑1 (MQ=25)
                                   tgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                               >  1:254315‑M1/30‑68 (MQ=255)
                                   tGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                               <  1:128333/68‑1 (MQ=255)
                                       aCTGTGCAAACGGCTCTCGAACACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGAGACGGCAAc                           >  2:521317/1‑68 (MQ=25)
                                       aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAACCTCAATGCAGCACCGGCTAc                           >  1:782518/1‑68 (MQ=255)
                                       aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                           <  1:422196/68‑1 (MQ=255)
                                       aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                           >  1:656194/1‑68 (MQ=255)
                                       aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                           >  1:353687/1‑68 (MQ=255)
                                       aCTGTGCAAACGGCTCTCGAAAACCTTGGCTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                           <  1:129370/68‑1 (MQ=255)
                                        cTGTTCTAACGTCTCTCGAAAACCTAGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACa                          <  1:388935/68‑1 (MQ=37)
                                        cTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACa                          <  2:700481/68‑1 (MQ=255)
                                         tgtgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACAt                         <  1:96823/68‑1 (MQ=255)
                                           tgCAAAAGGCTCTCGAAAACCTTGGTATGGGAGAAGGAGCAAAACTCAATGCCGCAACGGCTACATTa                       >  1:374707/1‑68 (MQ=25)
                                                 tcattttcatatacgTTGGTTTGGGAGAAGGAGCAATACTCAATGCAGCAACGGCTACATTAGGACGc                 >  2:373757‑M1/16‑68 (MQ=255)
                                                 cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCACTGCAGCAACGGCTACATTAGGACGc                 >  2:25656/1‑68 (MQ=255)
                                                 cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGc                 >  1:186177/1‑68 (MQ=255)
                                                 cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGc                 >  2:332064/1‑68 (MQ=255)
                                                 cGGCACTCGAAAAACTTGGTATGGGAGAAGGAGCAAAACTCAATGCAGCAACGTCTACATTAGGACGc                 >  2:449710/1‑68 (MQ=255)
                                                         gAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGttt         <  2:274909/68‑1 (MQ=255)
                                                         gAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGttt         <  2:186177/68‑1 (MQ=255)
                                                          aaaaCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                           >  1:617688/1‑49 (MQ=255)
                                                          aaaaCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                           <  2:617688/49‑1 (MQ=255)
                                                          aaaaCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTc        >  2:223231/1‑68 (MQ=255)
                                                             acgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                           >  1:46240‑M1/4‑46 (MQ=255)
                                                             acgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                           <  2:46240‑M1/43‑1 (MQ=255)
                                                             aCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACTCACCGGTTTCATa     <  2:277983/68‑1 (MQ=255)
                                                             aCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa     <  2:14524/68‑1 (MQ=255)
                                                              cgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAg    >  2:200425‑M1/3‑68 (MQ=255)
                                                               cTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGc   <  1:389669/68‑1 (MQ=255)
                                                               cTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCGACATTAGGACGCACCGGTTTCATAGc   <  1:620020/68‑1 (MQ=255)
                                                                TTGGTGTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCt  >  2:565421‑M1/1‑68 (MQ=255)
                                                                |                                                                   
AGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCT  >  NC_000913/1207726‑1207857

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