New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 120780515 (0.770)8 (0.550) 8/100 0.4 42.9% coding (305/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 = 1209602 10 (0.680)pseudogene (18/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

CCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTA  >  NC_000913/1207738‑1207867
                                                                   |                                                              
gtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAgg                                                                >  1:204948‑M1/53‑68 (MQ=255)
 cGTTTGGCGATATCAAACCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGc                                                               <  2:315333‑M1/68‑2 (MQ=255)
     tGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctg                                                           <  2:409828‑M1/68‑6 (MQ=255)
       gCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgca                                                         <  1:78134‑M1/68‑8 (MQ=255)
             tgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAAt                                                   >  1:300473‑M1/40‑68 (MQ=255)
                aaTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgtt                                                <  1:156621‑M1/68‑17 (MQ=255)
                 aTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgcag                                               >  2:418427/1‑68 (MQ=255)
                   catttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgcagca                                             <  1:638021‑M1/35‑1 (MQ=255)
                      aTGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACg                                          <  2:360017/68‑1 (MQ=255)
                       tgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                                         >  2:221943‑M1/30‑68 (MQ=255)
                       agcttgttatttcttcattttcataaacgTTGGTTTGTGAGAAGGAGCAAAACTCAATTCAGCAAcgt                                         >  2:4955‑M1/30‑67 (MQ=255)
                         gCACTGTGCAAACGGCTCTCGAACACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCt                                       <  2:348197/68‑1 (MQ=255)
                         gCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccg                                       >  2:173326‑M1/1‑43 (MQ=255)
                         gCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCt                                       <  2:199459/68‑1 (MQ=255)
                         gCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCt                                       <  1:650065/68‑1 (MQ=255)
                           tgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                     >  1:133329‑M1/26‑68 (MQ=255)
                           aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggt                                     >  2:424795‑M1/1‑41 (MQ=255)
                           aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAACCTCAATGCCGCAGCGGCTAc                                     >  2:445049/1‑68 (MQ=25)
                           aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                     >  1:155139/1‑68 (MQ=255)
                           aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                     >  1:165525/1‑68 (MQ=255)
                                   aaCGGCTCTCGAAAACCTTGGTTTGGGAGAATGctctgcagtacctgtttgtgtcccggttccgtggc                             <  1:294045‑M1/68‑36 (MQ=255)
                                     cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggcct                           >  1:666940‑M1/1‑31 (MQ=255)
                                     cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGc                           >  2:183301/1‑68 (MQ=255)
                                         tctcGATAACCTTGGTTTGGTAGAAGGAGCAATACTCAATGCAGCAACGTCTACATTATGACGCACCg                       <  2:171183/68‑1 (MQ=37)
                                            cGATAACCTTGGGTTGGGAGAAGGAGCAAAACTCATTGCAGCAACGGCTACTTTAGGACGCACCGGtt                    <  1:261158/68‑1 (MQ=255)
                                                 acgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa               >  1:167732‑M1/4‑68 (MQ=255)
                                                 aCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa               <  1:159233/68‑1 (MQ=255)
                                                 aCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa               <  1:222209/68‑1 (MQ=255)
                                                 aCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa               <  1:186683/68‑1 (MQ=255)
                                                   cTTGGTTTGGGAGAAGGAGCAATACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGc             <  2:627689/68‑1 (MQ=255)
                                                    TTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCt            >  1:256918‑M1/1‑68 (MQ=255)
                                                    TTGGTTGGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCTCCGGTTTCATAGCt            >  2:638022‑M1/1‑68 (MQ=255)
                                                        tttGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacag        >  2:151619‑M1/1‑12 (MQ=255)
                                                          TGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCg      >  2:165525‑M1/1‑68 (MQ=255)
                                                            GGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGt    >  1:311387‑M1/1‑68 (MQ=255)
                                                             gagaAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctgg   >  1:419362‑M1/1‑7 (MQ=255)
                                                              agaAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggc  >  2:12973‑M1/1‑6 (MQ=255)
                                                                   |                                                              
CCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTA  >  NC_000913/1207738‑1207867

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.