New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi = 12078056 (0.190)21 (0.750) 19/256 0.5 75.9% coding (305/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913_3_bme_pgi = 1209602 8 (0.290)pseudogene (18/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

ATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGG  >  NC_000913_3_bme_pgi/1207659‑1207946
                                                                                                                                                  |                                                                                                                                              
agcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgttgccatttgcttgttatttcttcattttcataaacgTTGGTTTGTGGGAAGGAg                                                                                                                                              <  2:126375‑M1/18‑1 (MQ=255)
   cGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCACACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctc                                                                                                                                           >  1:159963‑M1/1‑144 (MQ=255)
      ttacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccatggtaaaagcaggaatgaagccatttgcttgttatttcttcattttcataaacgTTGGTTTGGAAGCAGGAGCACAAc                                                                                                                                        <  2:221024‑M1/24‑1 (MQ=255)
         cAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattac                                                                                                                                     >  1:298672‑M1/1‑138 (MQ=255)
          gttgtagatccgctttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCaa                                                                                                                                    >  2:175493‑M1/122‑149 (MQ=255)
               ggAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttg                                                                                                                               <  1:192962‑M1/149‑18 (MQ=255)
               ggAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttg                                                                                                                               >  2:352430‑M1/1‑132 (MQ=255)
               agatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgcag                                                                                                                               >  2:280273‑M1/117‑149 (MQ=255)
                     aggctggcacgacCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcc                                                                                                                         <  2:441666‑M1/144‑24 (MQ=255)
                               cacCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtg                                                                                                               <  1:237359‑M1/148‑34 (MQ=255)
                                 tttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcgagcatgatgccatttgcttgttatttcttcattttcataaacgTTAGGTTGGGAGACGGAGCAAAACGCAATGACGCAACGGCCACATgagggc                                                                                                            <  1:242898‑M1/51‑3 (MQ=255)
                                       gatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCCCCg                                                                                                      <  1:149003‑M1/57‑1 (MQ=255)
                                        cttacCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCCAATCTGATGGCAATGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctcttcattacctgttggtgtcctgtttgcgtggccatcacc                                                                                                      >  1:71232‑M1/1‑107 (MQ=255)
                                         ttacCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGt                                                                                                    >  1:362929/1‑149 (MQ=255)
                                           acCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccac                                                                                                   >  2:418637‑M1/1‑104 (MQ=255)
                                           acCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccac                                                                                                   >  1:150301‑M1/1‑104 (MQ=255)
                                             cGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcggccactc                                                                                                 >  2:139907‑M1/1‑102 (MQ=255)
                                             cGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATAGGATGGCACTGTGCAAACGGCTATCGAAACCCTTGGTCTGGGAGAAGGAGACAAACTCAATGCAGCAACGGCTACACTAGGACGACCCGGtccca                                                                                                >  1:112200/1‑145 (MQ=18)
                                                         gctgctGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACATTGGTTTGGGACAAGGctctgcattacctgttggtgtcccggttccgtggccttcagcccctccgccaacaggct                                                                                     >  1:484757‑M1/1‑90 (MQ=255)
                                                             gcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTGCATAGCTATACCGGTTAt                                                                                <  1:141923‑M1/79‑1 (MQ=255)
                                                                tgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGAt                                                                             >  2:382248‑M1/68‑149 (MQ=255)
                                                                    aaatccccctggtaaaagcaggaatgatgccatctgattgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggt                                                                         >  1:238649‑M1/64‑149 (MQ=255)
                                                                       gCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaat                                                                       >  2:363238‑M1/1‑76 (MQ=255)
                                                                           cctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCAACGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTg                                                                  <  1:437651‑M1/93‑1 (MQ=255)
                                                                             aTCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagcccctccgccaacaggctggctgcaatgcaatggtacc                                                                 >  1:207315‑M1/1‑70 (MQ=255)
                                                                              ggtaaaagcaggaatgatgccatttggttgttatttcttcattttcataaacgTTGGTTTGGGAGCAGGCGCAAAACTCAATGCAGCAACGGCTACATTAGGAGGCACCGGTTTCATCACTGTACCGGGTGTGAGTGTTGCTCTTGagt                                                               <  2:444169‑M1/96‑2 (MQ=255)
                                                                                taaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCaa                                                             >  1:517706‑M1/52‑149 (MQ=255)
                                                                                  aaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTAtt                                                                   >  2:399564‑M1/50‑141 (MQ=255)
                                                                                  aaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTAtt                                                                   <  1:399564‑M1/92‑1 (MQ=255)
                                                                                        gATAACAAATCGGACGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgc                                                      <  2:214726‑M1/149‑91 (MQ=255)
                                                                                         atatCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtccctagtccgttgccttcagccactccgccagcaggctggctgaaatgcaatggtgccgctttttctgct                                                     <  1:241118‑M1/149‑92 (MQ=255)
                                                                                            tCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaa                                                  <  1:352430‑M1/149‑95 (MQ=255)
                                                                                               aaTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAg                                              >  1:269420/1‑149 (MQ=255)
                                                                                                gccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGt                                             >  1:327215‑M1/36‑149 (MQ=255)
                                                                                                    gATGGCACTGTGCAAACGGCTCTCGAAAAACATGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactcagccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatcccc                                          >  2:191622‑M1/1‑47 (MQ=255)
                                                                                                        cttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCACTCAGTAACGATTCAGTGGGGGTgg                                     <  1:476769‑M1/122‑1 (MQ=255)
                                                                                                         cACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGGAAAGGCTACATTAGGACGCACAGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTggg                                    >  1:7405/1‑148 (MQ=255)
                                                                                                          aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggacct                                    >  2:295721‑M1/1‑41 (MQ=255)
                                                                                                          aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaact                                    >  2:6652‑M1/1‑41 (MQ=255)
                                                                                                                   ttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAGATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGGGAATGCCGCaaa                          <  1:155594‑M1/133‑1 (MQ=255)
                                                                                                                     ggCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaataccgggaactggcaaaggctt                         >  1:252078‑M1/1‑30 (MQ=255)
                                                                                                                        tttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCAt                     <  2:438181‑M1/138‑1 (MQ=255)
                                                                                                                          tcataaacgTTGGTTTGGGAGAAGGAGCAAACCTCAAGGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCt                   >  2:297897‑M1/10‑149 (MQ=255)
                                                                                                                          tcGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactcccccaacaggctggctgaaatgcaatggtgccgctgtttctgctgaagaatacccggaactggcaaaggcttatcgg                    >  2:19478‑M1/1‑25 (MQ=255)
                                                                                                                            gAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCCCATACTCTAGCATTAGGCAGAACCGTTTTCATACCTTTACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAGAGCATCTGc                 >  2:525035/1‑149 (MQ=255)
                                                                                                                              agacattggtgtGGGAGAAGGatgtgcattacctgttggtgtcccggttccgtggccttcagccactccgcgaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaggggcttatccgacaa                >  2:231225‑M1/6‑21 (MQ=255)
                                                                                                                                     gtgatcGGAGAAGGAGCAAAAGTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCCATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTggggg        >  2:544749‑M1/7‑149 (MQ=255)
                                                                                                                                     gggttGGGAGAAGGTGCGAAACTCAATGCAGCACCGACTACCGTATGTCCGATCTGTGTCATAGCTATAAAAATTAAGCTTGGTATTAGTGATCATTCAGTAATCATTCAGTGGTTGTGGAATGCCGCAAAAGCATCTGCCTCTggggg        <  2:497479/146‑1 (MQ=255)
                                                                                                                                        ttGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcc      <  2:150301‑M1/149‑139 (MQ=255)
                                                                                                                                            gagaAGGctctgcattacttgttggtgttccgggtccgtggccttcagccactccgccaacaggctggctgaaatgcattggtgccgctttttgtgctgaagaatacgcggaactggcaaaggcttattcgagatataagttgcctgat  >  2:64072‑M1/1‑7 (MQ=255)
                                                                                                                                            gagaAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataaattgcctgat  <  2:354696‑M1/149‑143 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
ATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGG  >  NC_000913_3_bme_pgi/1207659‑1207946

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