New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 120780543 (1.210)12 (0.450) 9/98 0.9 33.1% coding (305/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 = 1209602 16 (0.600)pseudogene (18/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

CCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATT  >  NC_000913/1207738‑1207872
                                                                   |                                                                   
ccGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGG                                                                     <  1:670109‑M1/68‑1 (MQ=255)
    ttGGCGATATCATATCGGTTGGCACTGTGCAAACGGCTCTCGAAATCCTTGGTTTGGGAGAAGGctct                                                                 <  1:827496‑M1/68‑5 (MQ=255)
    ttGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctct                                                                 <  1:823559‑M1/68‑5 (MQ=255)
    ttGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctct                                                                 <  2:63690‑M1/68‑5 (MQ=255)
        cGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAAc                                                             <  2:1296100/68‑1 (MQ=255)
        cGATATCAAATCGGATGGCACTGCGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAAc                                                             <  2:219315/68‑1 (MQ=255)
           tatCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattac                                                          <  1:1096757‑M1/68‑12 (MQ=255)
             tCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacct                                                        <  1:320945‑M1/68‑14 (MQ=255)
             tCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAAt                                                        >  2:226665/1‑68 (MQ=255)
             tCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAAt                                                        <  1:49398/68‑1 (MQ=255)
                 aTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgcag                                                    >  2:161589/1‑68 (MQ=255)
                  tctgATGGCACTGTTCAAACGGCTCTCGAACACCTTGGTTTGGGAGAAGGAGCAAAACTCAAtgctgc                                                   <  2:363826/65‑1 (MQ=25)
                  tCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattac                                                          <  1:93718‑M1/61‑12 (MQ=255)
                  tCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattac                                                          >  2:93718‑M1/1‑50 (MQ=255)
                      aTGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACg                                               <  2:7518/68‑1 (MQ=255)
                       tGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                                              <  1:909427/68‑1 (MQ=255)
                        ggCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtccc                                             <  1:920947‑M1/68‑25 (MQ=255)
                         gCACTGTGCAAGCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCt                                            >  2:136255/1‑68 (MQ=255)
                         gCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccg                                            >  2:1044185‑M1/1‑43 (MQ=255)
                           tcTGTGCATACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                          <  2:947860/67‑1 (MQ=255)
                           aCTGTGCAATCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggt                                          <  1:1101007‑M1/68‑28 (MQ=255)
                           aCTGTGCAAACGGCTCTCGAATACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGTCTAc                                          <  2:1010858/68‑1 (MQ=255)
                           aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                          >  1:644662/1‑68 (MQ=255)
                           aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                          >  1:1123969/1‑68 (MQ=255)
                           aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                          >  1:1316281/1‑68 (MQ=255)
                           aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                          <  2:1213140/68‑1 (MQ=255)
                           aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                          <  1:257885/68‑1 (MQ=255)
                           aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                          >  1:129128/1‑68 (MQ=255)
                           aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                          >  2:594153/1‑68 (MQ=255)
                           aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACACAATGCAGCAACGGCTAc                                          >  2:1121085/1‑68 (MQ=255)
                            cTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACa                                         <  2:848000/68‑1 (MQ=255)
                            cTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACa                                         <  1:133101/68‑1 (MQ=255)
                            cTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACa                                         <  1:926983/68‑1 (MQ=255)
                            cTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAAAGCAGCAACTGCTACa                                         <  1:889124/68‑1 (MQ=255)
                                     cGGCTCTCGAAGACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGc                                >  2:658514/1‑68 (MQ=255)
                                     cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttc                                        <  1:658872‑M1/60‑30 (MQ=255)
                                     cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttc                                        >  2:658870‑M1/1‑31 (MQ=255)
                                         tctcGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCg                            <  2:148113/68‑1 (MQ=255)
                                                 acgttGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                    <  1:960247/65‑1 (MQ=255)
                                                 acgttGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                    <  1:924989/65‑1 (MQ=255)
                                                 aCCTTGGTTTGGGAGAATGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                    <  1:466398/68‑1 (MQ=255)
                                                 aCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCTCCGGTTTCATa                    <  2:144724/68‑1 (MQ=255)
                                                 aCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCAtg                    <  1:1176842/68‑2 (MQ=255)
                                                 aCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                    <  1:1216663/68‑1 (MQ=255)
                                                 aCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                    <  1:639689/68‑1 (MQ=255)
                                                 aCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                    <  2:1123968/68‑1 (MQ=255)
                                                 aCCTAGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                    <  1:1222939/68‑1 (MQ=255)
                                                        tttGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggc                                  >  2:330105‑M1/1‑12 (MQ=255)
                                                        tttGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggc                                  <  1:330107‑M1/47‑36 (MQ=255)
                                                         ttGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacagg            <  1:1265571‑M1/68‑58 (MQ=255)
                                                           gggAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCgg          <  1:1061633/68‑1 (MQ=255)
                                                           gggAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCgg          <  2:269458/68‑1 (MQ=255)
                                                           gggAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCgg          <  2:454912/68‑1 (MQ=255)
                                                           gggAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCgg          <  2:754879/68‑1 (MQ=255)
                                                           gggAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCgg          <  1:594155/68‑1 (MQ=255)
                                                            ggAGAAGGAGCAAAACTCAAAGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGt         <  1:804454/68‑1 (MQ=255)
                                                             gagaAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctgg        >  2:522621‑M1/1‑7 (MQ=255)
                                                              agaAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTa       >  2:1185527/1‑68 (MQ=255)
                                                                 aGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctga    <  2:480730‑M1/68‑66 (MQ=255)
                                                                 aGGcgctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctga    <  2:238663‑M1/68‑66 (MQ=255)
                                                                 aGGAGCTATACTCTATGCTGCAATGTCTTCATTAGGTCGCACCGGGTTTATAGCTATACCGGTTATGa    <  2:133949/68‑1 (MQ=255)
                                                                 aGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTGCATAGCTATACCGGTTATGa    <  2:421095/68‑1 (MQ=255)
                                                                   gctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaa  >  1:171310‑M1/1‑1 (MQ=255)
                                                                   |                                                                   
CCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATT  >  NC_000913/1207738‑1207872

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.