New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi 1207790 =5 (0.160)19 (0.680) 19/256 0.5 75.3% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913_3_bme_pgi 1209619 = 8 (0.290)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

GCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCC  >  NC_000913_3_bme_pgi/1207643‑1207932
                                                                                                                                                   |                                                                                                                                               
gCAGCGTGGGTGAGGCATCCGGCGCCAGCCCGGCGAAGGGGCGGCACGACCAGTTACTTCCCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGTCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTT                                                                                                                                                <  2:129519‑M1/149‑1 (MQ=255)
 catcctgcacagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccccggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTg                                                                                                                                               >  1:411718‑M1/147‑149 (MQ=255)
     gtggtgGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGGACAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTT                                                                                                                                           <  2:482431‑M1/149‑1 (MQ=255)
     ctgcacagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcgtgttatttcttcattttcataaacgTTGGttt                                                                                                                                           <  1:232769‑M1/7‑1 (MQ=255)
         acagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccctttgtcgtgtatttcttcattttcataaacgTTGGTTTGGga                                                                                                                                       >  1:540338‑M1/139‑149 (MQ=255)
            gttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGaa                                                                                                                                    >  2:359592‑M1/136‑149 (MQ=255)
                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                                                                                                                 >  2:352430‑M1/1‑132 (MQ=255)
                               ggAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttg                                                                                                                 <  1:192962‑M1/149‑18 (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)
                                                                                                                                                   |                                                                                                                                               
GCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCC  >  NC_000913_3_bme_pgi/1207643‑1207932

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.