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

CACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGC  >  NC_000913_3_bme_pgi/1209454‑1209748
                                                                                                                                                    |                                                                                                                                                  
cacGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCagag                                                                                                                                                    >  1:295487/1‑149 (MQ=255)
       aTCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctc                                                                                                                                             <  2:354696‑M2/149‑8 (MQ=255)
       aTCAGGCAACTTATATCTCGAATAAGCCTTTGCCAGTTCCGCGTATTCTTCAGCACAAAAAGCGGCACCAATGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGACCCGGAACACCAACAAGTAATGCAGAGccttctc                                                                                                                                             >  2:64072‑M2/1‑142 (MQ=255)
           ggCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaa                                                                                                                                         <  2:150301‑M2/149‑12 (MQ=255)
                     ttGTCGGATAAGCCCCTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTCGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCACATccttctcccacaccaatgtct                                                                                                                               >  2:231225‑M2/1‑128 (MQ=255)
                         ccgatAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAACAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGGGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcga                                                                                                                           >  2:19478‑M2/3‑124 (MQ=255)
                         cGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGa                                                                                                                           <  1:61124/149‑1 (MQ=255)
                              aaGCCTTTGCCAGTTCCCGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagcc                                                                                                                      >  1:252078‑M2/1‑119 (MQ=255)
                                         aGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                                                                                                           >  2:6652‑M2/1‑108 (MQ=255)
                                         aGGTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagt                                                                                                           >  2:295721‑M2/1‑108 (MQ=255)
                                               ggggATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCTGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccatgtttttcgagagccgtttgcacagtgccatc                                                                                                     >  2:191622‑M2/1‑102 (MQ=255)
                                                       ttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttga                                                                                             <  1:352430‑M2/149‑56 (MQ=255)
                                                       ttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                                                                                             >  2:222585/1‑149 (MQ=255)
                                                          agcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAGGGCCACGGAACCGGGAAACCAACAGGTAATGCCGAGCCTTCTCCCAAACCAACGTTTCTGAAAATGAAGAAATAACAAGCAAACGGCGTCATTc                                                                                          >  1:159175/1‑149 (MQ=255)
                                                          agcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGCTGGCGGAGTGGCTGAAGGCAACGGACTAGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatat                                                                                          <  1:241118‑M2/149‑59 (MQ=255)
                                                           gcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccgtccgatttgttatc                                                                                         <  2:214726‑M2/149‑60 (MQ=255)
                                                             agAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg                                                                                       <  1:354696/149‑1 (MQ=255)
                                                                      ggTACCATTGCATTGCAGCCAGCCTGTTGGCGGAGGGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggat                                                                              >  1:207315‑M2/1‑79 (MQ=255)
                                                                            aTTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgc                                                                        >  2:363238‑M2/1‑73 (MQ=255)
                                                                                          agccTGTTGGCGGAGGGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttgtcccaaaccaatgttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagc                                                          >  1:484757‑M2/1‑59 (MQ=255)
                                                                                              tgtgGGCGGAGGAGCGGAAGGACACGGAACCGCGAAACCAACAGGTAAGGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGCAATAACAACCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGc                                                      <  2:164186/149‑1 (MQ=255)
                                                                                                      gAGTGGCTTAAGGCCACGGAACCGGGACACCAACAGGTACTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGACGAAATAACAAGCAAATTGCATCATTCCTGCTTTTACCAGTGCGATTTAACATTCGTATTGGCTATGTAct                                              >  2:436048/1‑148 (MQ=255)
                                                                                                      gAGTGGCCGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcg                                              >  2:139907‑M2/1‑47 (MQ=255)
                                                                                                        gTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggt                                            >  1:150301‑M2/1‑45 (MQ=255)
                                                                                                        gTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggt                                            >  2:418637‑M2/1‑45 (MQ=255)
                                                                                                           ggtgaTGGCCACGCAAACAGGACACCAACAGGTAATGAAGAGccttctcccaaaccaaggttttcgagagccgtttgcacattgccatcagattggatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaag                                         >  1:71232‑M2/3‑42 (MQ=255)
                                                                                                                    cACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtg                                <  1:237359‑M2/149‑117 (MQ=255)
                                                                                                                              ggACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccagcct                      <  2:441666‑M2/149‑127 (MQ=255)
                                                                                                                                    cAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctcc                >  2:352430‑M2/1‑17 (MQ=255)
                                                                                                                                    cAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctcc                <  1:192962‑M2/149‑133 (MQ=255)
                                                                                                                                         ggTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAAc           <  1:222585/149‑1 (MQ=255)
                                                                                                                                          gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctg          >  1:298672‑M2/1‑11 (MQ=255)
                                                                                                                                                gagagccttctcccaaaccaaggttttcgagagccgtgtgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcagcgcgcgaagtgcggtaagtaactggtcgtgccgccccttctccaggctggcaccg    >  1:159963‑M2/2‑5 (MQ=255)
                                                                                                                                                  gagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACCCGTATCAACAAATGACCAGAACACAGATCTACAACGTAACgcgc  >  2:482442/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CACGTAAATCAGGCAATTTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGC  >  NC_000913_3_bme_pgi/1209454‑1209748

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.