New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =31 (0.780)10 (0.340) 10/96 0.7 24.9% coding (290/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 1209619 = 37 (1.240)pseudogene (37/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

AGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATA  >  NC_000913/1207726‑1207854
                                                                |                                                                
aGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCt                                                                  >  1:2178236/1‑65 (MQ=255)
aGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCt                                                                  >  1:83494/1‑65 (MQ=255)
aGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCt                                                                  >  1:1558447/1‑65 (MQ=255)
 aaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgtt                                                                 <  1:1023854‑M1/2‑1 (MQ=255)
    tccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGt                                                              <  1:3226316‑M1/5‑1 (MQ=255)
     ccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGtt                                                             <  1:1115390‑M1/6‑1 (MQ=255)
         ctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTggg                                                         >  1:3096702‑M1/56‑65 (MQ=255)
             cGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGaa                                                     <  1:1905402/65‑1 (MQ=255)
              aaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAg                                                    >  1:75235‑M1/51‑65 (MQ=255)
              aaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAg                                                    >  1:1176734‑M1/51‑65 (MQ=255)
               aaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAgg                                                   >  1:623066‑M1/50‑65 (MQ=255)
                aagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGa                                                  <  1:1506381‑M1/17‑1 (MQ=255)
                  ggCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGc                                                <  1:2358011/65‑1 (MQ=255)
                  ggCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGc                                                >  1:925428/1‑65 (MQ=255)
                     gATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCaaa                                             <  1:579128/65‑1 (MQ=255)
                     gATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCaaa                                             <  1:799732/65‑1 (MQ=255)
                       aatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAAc                                           <  1:1034928‑M1/24‑1 (MQ=255)
                          cAAATCGGATGGCCCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCa                                        <  1:1804219/65‑1 (MQ=255)
                             gccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATg                                     <  1:532102‑M1/30‑1 (MQ=255)
                              tCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgc                                    >  1:1202130/1‑65 (MQ=255)
                                  aTGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGCGAAGGAGCAAAACTCAATgcagca                                >  1:732690/1‑65 (MQ=255)
                                     gCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACg                             >  1:2202586/1‑65 (MQ=255)
                                     gCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACg                             <  1:296225/65‑1 (MQ=255)
                                     cttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACg                             >  1:1550915‑M1/28‑65 (MQ=255)
                                      ttgttattgcttcattttcataaaggTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                            >  1:1300248‑M1/27‑65 (MQ=255)
                                       tgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGc                           >  1:174142‑M1/26‑65 (MQ=255)
                                       aCTGTGCAAACGGCTCTCGAAAGCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGc                           <  1:297597/65‑1 (MQ=255)
                                       aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGc                           >  1:1533528/1‑65 (MQ=255)
                                       aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGc                           >  1:519427/1‑65 (MQ=255)
                                       aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGc                           >  1:1477806/1‑65 (MQ=255)
                                       aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGc                           >  1:267699/1‑65 (MQ=255)
                                       aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGc                           >  1:63110/1‑65 (MQ=255)
                                       aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGc                           >  1:2169198/1‑65 (MQ=255)
                                           tgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACa                       <  1:803119/65‑1 (MQ=255)
                                           atttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACa                       >  1:2399470‑M1/22‑65 (MQ=255)
                                                ttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAgg                  >  1:2879812‑M1/17‑65 (MQ=255)
                                                 cGGCTCTCGAAAACCTTGGTTTGGGGGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGa                 >  1:1927047/1‑65 (MQ=255)
                                                 cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGa                 >  1:1056963/1‑65 (MQ=255)
                                                 cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGa                 >  1:319133/1‑65 (MQ=255)
                                                 cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGa                 >  1:273673/1‑65 (MQ=255)
                                                 cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGa                 >  1:1642985/1‑65 (MQ=255)
                                                 cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAAATTAGGa                 >  1:1102257/1‑65 (MQ=255)
                                                     tttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCa             >  1:692380‑M1/12‑65 (MQ=255)
                                                     tctcGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCa             <  1:2674732/65‑1 (MQ=255)
                                                      ttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGcag            <  1:1744529‑M1/55‑2 (MQ=255)
                                                        cGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCg          >  1:1920082/1‑65 (MQ=255)
                                                          aaaaCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGt        <  1:1197681/65‑1 (MQ=255)
                                                          aaaaCCTTGGGTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGt        <  1:2647341/65‑1 (MQ=255)
                                                                TTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa  >  1:2763059‑M1/1‑65 (MQ=255)
                                                                TTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa  >  1:2493641‑M1/1‑65 (MQ=255)
                                                                |                                                                
AGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATA  >  NC_000913/1207726‑1207854

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