New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi = 12078054 (0.130)33 (1.240) 29/254 0.1 89.7% coding (305/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913_3_bme_pgi = 1209602 4 (0.150)pseudogene (18/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

ATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACG  >  NC_000913_3_bme_pgi/1209469‑1209740
                                                                                                                                     |                                                                                                                                          
atttatttGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaacca                                                                                                                             <  1:350303‑M2/149‑16 (MQ=255)
   tatttGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaagg                                                                                                                          >  2:301621‑M2/1‑131 (MQ=255)
   gatttgTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGGCAGCGTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaagg                                                                                                                          >  1:73134‑M2/2‑131 (MQ=255)
        gTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttc                                                                                                                     >  1:426628‑M2/1‑126 (MQ=255)
             aTAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagag                                                                                                                >  1:316888‑M2/1‑121 (MQ=255)
                      tGCCAGTTCCGGGTATGCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagcagtttgca                                                                                                       >  1:506750‑M2/1‑112 (MQ=255)
                      ggcCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgca                                                                                                       >  2:453728‑M2/2‑112 (MQ=255)
                         cAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacag                                                                                                    <  2:166010‑M2/149‑41 (MQ=255)
                          aGTGCCGGGCATGCTCCACCAAAAAAAACGGCACCAATGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAACGCAGAGccgtctcccaaaccaatggtttcgagagccgtttgcacagt                                                                                                   >  2:224228‑M2/1‑108 (MQ=255)
                            ttCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagttc                                                                                                 <  2:447510‑M2/149‑44 (MQ=255)
                                gggTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatc                                                                                             >  1:266009‑M2/1‑102 (MQ=255)
                                 ggTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatcc                                                                                            <  1:372409‑M2/149‑49 (MQ=255)
                                    aTTCTTCAGCAGAAAAAGCGGCACCTATGCATTTCAGCCAGCCTGTGGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATTCAGAGGCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGc                                                                                         <  2:135537/149‑1 (MQ=255)
                                           agcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCCAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc                                                                                  >  1:213271/1‑149 (MQ=255)
                                            gcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatc                                                                                 >  1:47265‑M2/1‑90 (MQ=255)
                                                    aGCGGCACCTTTGCGTTGCAGCCAGCCTGTTTGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacg                                                                         >  2:148747‑M2/1‑82 (MQ=255)
                                                    aGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGACACGGAACCGGGACACCAAAAGGTAATGCAGAGcgttctcccaaaccaaggttttcgggagccgtttgcacagtgccgttcggtttgatttcgccaaacg                                                                         <  2:206485‑M2/149‑68 (MQ=255)
                                                      cGGCAACATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCCGCAccttctaccaaaccaaggttttcgacatccgtttgcacagcgccatccgatttgatatcgccaaacgga                                                                       >  2:377754‑M2/1‑80 (MQ=255)
                                                           ccATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattctt                                                                  >  1:475079‑M2/1‑75 (MQ=255)
                                                                gCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggc                                                             >  1:354245‑M2/1‑70 (MQ=255)
                                                                 cATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggct                                                            >  2:83612‑M2/1‑69 (MQ=255)
                                                                      cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATg                                                       <  2:278855/149‑1 (MQ=255)
                                                                           agccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgatagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagc                                                  >  2:417739‑M2/1‑59 (MQ=255)
                                                                           agccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagc                                                  <  2:132467‑M2/149‑91 (MQ=255)
                                                                             ccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgc                                                >  1:132467‑M2/1‑57 (MQ=255)
                                                                                  tGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaa                                           >  2:299314‑M2/1‑52 (MQ=255)
                                                                                  tGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttacacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaa                                           <  2:469076‑M2/149‑98 (MQ=255)
                                                                                    gCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCGGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcgg                                                              >  1:439389‑M2/1‑50 (MQ=255)
                                                                                    gCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCGGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcgg                                                              <  2:439389‑M2/128‑79 (MQ=255)
                                                                                       gAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcg                                      >  1:2685‑M2/1‑47 (MQ=255)
                                                                                             cggaaGGCCACGGAACCGGGACAACAACAGGTAATGCAGAGccttgtcccaaacgagggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagt                                >  2:350303‑M2/3‑41 (MQ=255)
                                                                                             cTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcaccgtgccagctgatgtgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagt                                >  1:365912‑M2/1‑41 (MQ=255)
                                                                                                 aGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaact                            >  1:166010‑M2/1‑37 (MQ=255)
                                                                                                         gAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgc                    <  2:420698‑M2/149‑121 (MQ=255)
                                                                                                             cGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgcc                >  1:420698‑M2/1‑25 (MQ=255)
                                                                                                                     cAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGATTATGAAAATGAAGAAATAACAAGCAAATGACATCATTCCTGCTTTTACCAGGGGGATTTAACATACTTATTGGATATGTACGCGTATCAACAAATGACCAGAACACAGATCt        >  2:465807/1‑149 (MQ=255)
                                                                                                                           gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgtcaaacggattcttgcagcttaacaacagcgcgcgaagtgccgtaagtaactggtcatgccgccccttctacaggctg  <  2:386303‑M2/149‑139 (MQ=255)
                                                                                                                           gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttaatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtacggtaagtaactggtcgtgccgacccttccccaggctg  <  2:386121‑M2/149‑139 (MQ=255)
                                                                                                                                     |                                                                                                                                          
ATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACG  >  NC_000913_3_bme_pgi/1209469‑1209740

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.