New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 12078054 (0.100)18 (0.620) 12/98 0.7 87.8% coding (305/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 = 1209602 2 (0.070)pseudogene (18/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

CGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTA  >  NC_000913/1207739‑1207867
                                                                  |                                                              
taaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGa                                                               <  2:1214231‑M1/17‑1 (MQ=255)
cGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGc                                                               <  1:1269041‑M1/68‑2 (MQ=255)
  aaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGc                                                             >  1:127689‑M1/50‑68 (MQ=255)
   ttGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctct                                                            <  1:662122‑M1/68‑5 (MQ=255)
   ttGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctct                                                            <  1:1408297‑M1/68‑5 (MQ=255)
    tGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctg                                                           <  1:147322‑M1/68‑6 (MQ=255)
     ggCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgc                                                          <  1:711392‑M1/68‑7 (MQ=255)
     gcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCaaa                                                          <  1:167204‑M1/22‑1 (MQ=255)
      gCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgca                                                         >  1:443825‑M1/1‑61 (MQ=255)
      gCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCaaaa                                                         >  2:854695/1‑68 (MQ=255)
          tatCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattac                                                     <  2:477570‑M1/68‑12 (MQ=255)
          tatCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattac                                                     <  1:978142‑M1/68‑12 (MQ=255)
          tatCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattac                                                     <  1:441976‑M1/68‑12 (MQ=255)
          aatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCa                                                     <  2:766294‑M1/27‑1 (MQ=255)
               aaTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgtt                                                <  1:1288779‑M1/68‑17 (MQ=255)
                 tCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGcactgcattaccagttgg                                              >  2:1134241‑M1/1‑50 (MQ=255)
                 ccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgcagc                                              <  1:243424‑M1/34‑1 (MQ=255)
                    tttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAAc                                           >  1:766298‑M1/32‑68 (MQ=255)
                      tgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                                         >  2:696217‑M1/30‑68 (MQ=255)
                        cttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCt                                       >  2:288472‑M1/28‑68 (MQ=255)
                          tgttttttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                     >  2:950221‑M1/26‑68 (MQ=255)
                          tgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                     >  1:1214235‑M1/26‑68 (MQ=255)
                          tgttatttcttcattttcataaacgTCGGTTTGGGAGAAGGAGCAATACTCAATGCAGCAACGGCTAc                                     >  2:935483‑M1/26‑68 (MQ=255)
                          aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggt                                     >  1:1162324‑M1/1‑41 (MQ=255)
                          aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                     <  1:1290013/68‑1 (MQ=255)
                           gttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACa                                    >  2:48418‑M1/25‑68 (MQ=255)
                           gttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACa                                    >  2:1093492‑M1/25‑68 (MQ=255)
                            tgtgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttc                                   <  1:664588‑M1/68‑30 (MQ=255)
                                    cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggcct                           >  1:464227‑M1/1‑31 (MQ=255)
                                    cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggcct                           >  2:439736‑M1/1‑31 (MQ=255)
                                    cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggcct                           >  2:11972‑M1/1‑31 (MQ=255)
                                    cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattaactgttggtgtcccggttccgtggcct                           >  1:305789‑M1/1‑31 (MQ=255)
                                    cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGc                           >  1:905238/1‑68 (MQ=255)
                                    cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGc                           >  2:329156/1‑68 (MQ=255)
                                             aaaaCCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccact                  >  2:143002‑M1/1‑22 (MQ=255)
                                                gcgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCatg               <  2:1152083‑M1/65‑2 (MQ=255)
                                                acgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa               >  1:636617‑M1/4‑68 (MQ=255)
                                                acgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa               >  1:1304056‑M1/4‑68 (MQ=255)
                                                acgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa               >  1:829014‑M1/4‑68 (MQ=255)
                                                acgTTGGTGTGGGAGAAGGAGCTAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa               >  1:1470770‑M1/4‑68 (MQ=255)
                                                 cgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAg              <  1:40235‑M1/66‑1 (MQ=255)
                                                   ttGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgcca            <  1:242886‑M1/68‑53 (MQ=255)
                                                      gTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaaca         <  1:814010‑M1/68‑56 (MQ=255)
                                                       TTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATAc        <  1:1383694‑M1/68‑1 (MQ=255)
                                                          GGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCgg     >  1:666550‑M1/1‑68 (MQ=255)
                                                           GGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                     <  1:1370855‑M1/35‑1 (MQ=255)
                                                           GGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                     >  2:1370850‑M1/1‑35 (MQ=255)
                                                            GAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGGCGCACCGGTTTCATAGCTATACCGGtt   <  2:1396041‑M1/68‑1 (MQ=255)
                                                            GAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGtt   <  1:405935‑M1/68‑1 (MQ=255)
                                                             AGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTa  <  2:220923‑M1/68‑1 (MQ=255)
                                                                  |                                                              
CGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTA  >  NC_000913/1207739‑1207867

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.