New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 120780534 (0.360)15 (0.180) 15/252 3.4 32.2% coding (305/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 = 1209602 33 (0.390)pseudogene (18/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related
Rejected: Coverage evenness skew score above cutoff.

CATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATA  >  NC_000913/1207658‑1207951
                                                                                                                                                   |                                                                                                                                                  
cagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAgg                                                                                                                                                    >  3:186655‑M1/133‑148 (MQ=255)
    cGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAggctctg                                                                                                                                               >  3:290773/1‑144 (MQ=255)
         ccAGCCTGGAGAAGGGGCGGCACGACCAGTGACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACt                                                                                                                                           <  8:267416/148‑1 (MQ=255)
         acgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACt                                                                                                                                           >  1:242654‑M1/124‑148 (MQ=255)
          cAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCa                                                                                                                                         <  5:354108/149‑1 (MQ=255)
           gttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattctcataaacgTTGGGTTGGGAGAAGGAGCAAAACTCaa                                                                                                                                        >  2:186963‑M1/122‑149 (MQ=255)
           aGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAACCCGCAAGAATCCGTTTGGCGATATCAAATCCGAATGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCaa                                                                                                                                        <  3:29299/149‑1 (MQ=255)
            gCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAAt                                                                                                                                       >  2:379041/1‑149 (MQ=255)
            gCCTGGAGAAGGGGCGGAACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGATTCCGTTTGGCGATATCAAATCGGATGGCACTGAGAAAACGGCTCTCGAAATCCTTGGGTTGGGAGAAGGAGCAAAACTCAAt                                                                                                                                       >  2:384517/1‑149 (MQ=18)
             tgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATg                                                                                                                                      <  4:29997‑M1/30‑1 (MQ=255)
             ccTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATg                                                                                                                                      <  1:379041/149‑1 (MQ=255)
             ccTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATg                                                                                                                                      <  1:384517/149‑1 (MQ=255)
                 gagaAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgcagc                                                                                                                                  <  4:168967/149‑1 (MQ=255)
                   tctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCaa                                                                                                                                <  4:367627‑M1/36‑1 (MQ=255)
                      ggggCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                                                                                                                             <  1:288195/149‑1 (MQ=255)
                              gtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACAt                                                                                                                       <  6:292793‑M1/45‑1 (MQ=255)
                              aCGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTa                                                                                                                     <  1:59568/149‑1 (MQ=255)
                                 aCCATTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGa                                                                                                                  <  7:91620/149‑1 (MQ=255)
                                     gttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgctttttatttcttcattttcgtaaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCa                                                                                                              <  1:240667‑M1/54‑1 (MQ=255)
                                     gttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCa                                                                                                              <  1:129147‑M1/54‑1 (MQ=255)
                                              cGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCAc                                                                                                             >  8:154579/1‑141 (MQ=255)
                                              cGCACTTCGCGCGCTGCTGTTAAGCAGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCAc                                                                                                             <  7:154579/141‑1 (MQ=255)
                                                      gcgcgcTGCTGTTAAGCCGCAAGAATCCGTTTGGCGACATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCtata                                                                                             <  6:278079/149‑1 (MQ=255)
                                                       ccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATAc                                                                                            <  3:367065‑M1/72‑1 (MQ=255)
                                                                gTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATATCTATAACGGTTATGa                                                                                   >  1:240288/1‑149 (MQ=255)
                                                                atgttaaatccccctggtaaaagcaggaatgatgccagttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGa                                                                                   >  4:202953‑M1/69‑149 (MQ=255)
                                                                   ttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGAtt                                                                                 >  6:96770‑M1/66‑148 (MQ=255)
                                                                   aaGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGAtt                                                                                 <  2:17461/148‑1 (MQ=255)
                                                                     gCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTGTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAACTTAGGACGCACAGCTTTCATAGCTATACAGGTTATGAAtggt                                                                              >  3:86257/1‑149 (MQ=25)
                                                                      aatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggtg                                                                             >  5:292793‑M1/63‑149 (MQ=255)
                                                                        gCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTggtggt                                                                           <  1:9537/149‑1 (MQ=255)
                                                                        gCAAGAATCCGTTTGGCGATATCAAATCCGATGGCACTGTGCAACCGGCTCTCGAAAACCTTGGTGTGGGAGAAGGAGCAAACCTCAATGCAGCATCGGCAACATTAGGAc                                                                                                                 >  1:89048/1‑111 (MQ=21)
                                                                          aaGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTAt                                                                         <  1:183409/149‑1 (MQ=255)
                                                                                         ggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAAt                                                          <  5:130900‑M1/106‑1 (MQ=255)
                                                                                         gATATCAAATCGGATGGAACTGTGCAAACGGCTCTCGAAAAACTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGAACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAAt                                                          >  2:48060/1‑149 (MQ=255)
                                                                                           tatCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCta                                                                                               >  8:401766/1‑110 (MQ=255)
                                                                                           tatCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCta                                                                                               <  7:401766/110‑1 (MQ=255)
                                                                                               aaaTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCAtt                                                      >  7:350526/1‑147 (MQ=255)
                                                                                                        ggCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGc                                                                                                 <  2:214334/95‑1 (MQ=255)
                                                                                                        ggCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGc                                                                                                 >  1:214334/1‑95 (MQ=255)
                                                                                                            gttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAACCTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATAGGTGGGATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAAt                                       <  5:372425‑M1/125‑1 (MQ=255)
                                                                                                               atttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGcc                                    <  7:241127‑M1/128‑1 (MQ=255)
                                                                                                                tttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCg                                   >  3:367627‑M1/21‑149 (MQ=255)
                                                                                                                     cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCaaaa                              <  8:350526/149‑1 (MQ=255)
                                                                                                                      cattttcgtaaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCATCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCaaaa                              >  2:240667‑M1/15‑148 (MQ=255)
                                                                                                                      cattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCaaa                               >  4:367065‑M1/15‑147 (MQ=255)
                                                                                                                            cataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTgggg                                              <  7:25047‑M1/118‑1 (MQ=255)
                                                                                                                            cataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTgggg                                              >  8:25047‑M1/9‑126 (MQ=255)
                                                                                                                            cGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCaaaa                              <  2:236042/142‑1 (MQ=255)
                                                                                                                            cGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCaaaa                              >  1:236042/1‑142 (MQ=255)
                                                                                                                             gAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGc                      >  1:285537/1‑149 (MQ=255)
                                                                                                                             gAAAACCTAGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATAAGGACGCACCGGTATCATAGCTATACCGGTTATGATTGGTTGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCAGCAAAAGCATCtcc                      >  7:276989/1‑147 (MQ=255)
                                                                                                                                   cTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTc                                                     >  4:340419/1‑112 (MQ=255)
                                                                                                                                   cTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTc                                                     <  3:340419/112‑1 (MQ=255)
                                                                                                                                      GGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTggggg             <  1:260703‑M1/149‑1 (MQ=255)
                                                                                                                                       GTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACGTTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTgggggg            >  6:61625‑M1/1‑149 (MQ=255)
                                                                                                                                       GTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTggggg             >  3:29997‑M1/1‑148 (MQ=255)
                                                                                                                                           gggAGAATTAGAAAAACTCAATGCACCAACTTCTACATTAGGCCGCACCGGTTTCTTAGCTATACCGGTTATGATTGGTTGTTTTTATCAATCAGTAATCATTCTGTGGGGGTGGAATGCCGAAAAAGCATCTGCCTCTGGGGGGGAt         <  2:240288/148‑1 (MQ=255)
                                                                                                                                                  GGAGCAAAACTCAATGCAGCAACGGCCACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCACTAATCAATCACTGGGGGTGCAATGCCGAAAAAGCATCTGCCTCTGGGGGGGATGGAAATa  <  1:262518‑M1/148‑1 (MQ=255)
                                                                                                                                                   GAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGGTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGTGGGGATGGAAATa  <  3:256544‑M1/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATA  >  NC_000913/1207658‑1207951

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