New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =8 (0.310)13 (0.670) 10/100 0.5 76.4% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 1209619 = 2 (0.100)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

TAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGC  >  NC_000913/1207724‑1207856
                                                                  |                                                                  
tAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTT                                                                   >  1:168995‑M1/1‑68 (MQ=255)
  aGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGG                                                                 >  2:613113‑M1/1‑68 (MQ=255)
   gCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGT                                                                >  1:556342‑M1/1‑68 (MQ=255)
   aaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGt                                                                >  1:244769‑M1/64‑68 (MQ=255)
     ctcaccctggtaaaagcaggaatgatgccatttgcttgttatatcttcattatcataaacgTAGGttt                                                              <  1:860111‑M1/7‑1 (MQ=255)
      tccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTg                                                             <  1:915325‑M1/8‑1 (MQ=255)
      tccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTg                                                             >  2:113523‑M1/61‑68 (MQ=255)
       ccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTgg                                                            <  2:201851‑M1/9‑1 (MQ=255)
               cGTTTGGCGATATCAAATCGGATGGCACTGTTCATACGGCGCTCGAAAACCTTGGTTTGTGAGAAGGc                                                    <  1:761751‑M1/68‑2 (MQ=255)
               cGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGc                                                    <  1:39045‑M1/68‑2 (MQ=255)
                  ttGGCGATATCATATCGGATGGCACTGTGCAAACGGCTCTCGAATACCTTGGTTTGGGAGAAGGctct                                                 <  1:905777‑M1/68‑5 (MQ=255)
                   agcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCaa                                                <  1:814122‑M1/21‑1 (MQ=255)
                    ggCGAAATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCaaa                                               <  2:68575/68‑1 (MQ=255)
                    gcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCaaa                                               <  2:417835‑M1/22‑1 (MQ=255)
                    gcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCaaa                                               <  2:501315‑M1/22‑1 (MQ=255)
                     gCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGcgctgca                                              <  2:326856‑M1/68‑8 (MQ=255)
                      aggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAAc                                             >  2:846001‑M1/45‑68 (MQ=255)
                         aatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCa                                          >  1:533247‑M1/42‑68 (MQ=255)
                           tgatgccatttgcttgttatttcttcattttcataatcgTTGGTTTGGGAGAATGAGCAAAACTCAAt                                        >  1:932599‑M1/40‑68 (MQ=255)
                           tCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacct                                        <  1:96534‑M1/68‑14 (MQ=255)
                           tCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGG                                                     <  2:178271‑M1/55‑1 (MQ=255)
                           tCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGG                                                     >  1:178271‑M1/1‑55 (MQ=255)
                                tCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGA                                                        <  1:635679‑M1/47‑1 (MQ=255)
                                tCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGA                                                        >  2:635682‑M1/1‑47 (MQ=255)
                                   gATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgt                                <  1:783218‑M1/68‑22 (MQ=255)
                                    aTGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtc                               <  2:655128‑M1/68‑23 (MQ=255)
                                    aTGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACg                               <  2:556345/68‑1 (MQ=255)
                                     tGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                              <  2:706501/68‑1 (MQ=255)
                                      gcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGc                             <  2:408316‑M1/40‑1 (MQ=255)
                                         tgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                          >  1:660697‑M1/26‑68 (MQ=255)
                                         aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggt                          >  1:346235‑M1/1‑41 (MQ=255)
                                         aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggcgtcccggc                          <  1:948777‑M1/68‑28 (MQ=255)
                                         aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                          >  2:639207/1‑68 (MQ=255)
                                         aCTGTGCAAACGGCTCTCGAAAACCTTGGATTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                          >  1:100615/1‑68 (MQ=255)
                                          gttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACa                         >  1:773043‑M1/25‑68 (MQ=255)
                                          gttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACa                         >  2:54150‑M1/25‑68 (MQ=255)
                                          cTGTGCAAACGGCTCTCGAAAACCTTGTTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACa                         <  1:639862/68‑1 (MQ=255)
                                           tgtgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttc                        <  1:568213‑M1/68‑30 (MQ=255)
                                           tgtgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttc                        <  1:415919‑M1/68‑30 (MQ=255)
                                             atttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTa                      <  1:601911‑M1/47‑1 (MQ=255)
                                             attacttcattttcataaacgTAGGTTTGGGAGAAGGAGCAAAACTCAATGCATCAACGGCTACattt                      <  2:738098‑M1/47‑2 (MQ=255)
                                              gCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgt                     >  2:712908‑M1/1‑36 (MQ=255)
                                                 aaCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggc                  <  2:914998‑M1/68‑36 (MQ=255)
                                                   cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggcct                >  2:556088‑M1/1‑31 (MQ=255)
                                                   cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgca                                              <  2:521623‑M1/38‑8 (MQ=255)
                                                   cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgca                                              >  1:521620‑M1/1‑31 (MQ=255)
                                                   cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgc                                               >  2:143646‑M1/1‑31 (MQ=255)
                                                   cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgc                                               <  1:143646‑M1/37‑7 (MQ=255)
                                                   cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGc                >  1:762840/1‑68 (MQ=255)
                                                          cataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGtt         >  1:201851‑M1/9‑68 (MQ=255)
                                                            aaaaCCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccact       >  1:210080‑M1/1‑22 (MQ=255)
                                                             aaaCCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactc      <  1:281986‑M1/68‑48 (MQ=255)
                                                               aCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa    <  2:762844/68‑1 (MQ=255)
                                                                 cTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgcc  <  1:470606‑M1/68‑52 (MQ=255)
                                                                  |                                                                  
TAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGC  >  NC_000913/1207724‑1207856

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.