New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi 1207790 =4 (0.070)29 (0.540) 24/266 1.4 87.1% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913_3_bme_pgi 1209619 = 5 (0.090)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

TTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGG  >  NC_000913_3_bme_pgi/1209490‑1209761
                                                                                                                                 |                                                                                                                                              
cggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                                                                                                            >  1:462649‑M2/130‑150 (MQ=255)
cggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                                                                                                            >  1:717494‑M2/130‑150 (MQ=255)
    cAGTTCCGGGTATTCTTCCTCATAAAAATCGGCACCATTGCATTTCAGCCAGCCTGTTTGCGGAGTGGATGAAGGCCACGGAACCGGTACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACa                                                                                                                        <  2:93631/150‑1 (MQ=255)
      tccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGGTTATGAAAATGAAGAAATAACAAg                                                                                                                      <  2:81759‑M2/27‑1 (MQ=255)
       ccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaacccggtgcgtcctaatgtagccgtggctgcattgagttttgctccttcgcccacaccacCGTTTATGAAAATGAAGAAATAACAAGc                                                                                                                     >  2:739494‑M2/123‑150 (MQ=255)
         acccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCaa                                                                                                                   >  1:509483‑M2/121‑150 (MQ=255)
            cccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATg                                                                                                                >  1:761505‑M2/118‑150 (MQ=255)
            cccactgaatgattactgattgctcaataccaccaatcataaccggtatagctaagaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAATTAACAAGCAAATg                                                                                                                >  1:680636‑M2/118‑150 (MQ=255)
                ctgaatgattactgattgcgcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcat                                                                                                            <  2:748096‑M2/37‑1 (MQ=255)
                   aatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatcat                                                                                                         >  1:165531‑M2/111‑150 (MQ=255)
                      gattactgattgctcaataccaccaatcataagcggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                                                                                                      <  2:488230‑M2/43‑1 (MQ=255)
                       gcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGATAATGAAGAAATAACAAGCAAATGGCATCATTCCt                                                                                                     >  2:777539/1‑150 (MQ=255)
                                    caataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                                                                                        <  1:575166‑M2/57‑1 (MQ=255)
                                      ccATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGCACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCATACCACCGTTTATGAAATTGAAGAAATAACAAGCAATTGGCATCATTCCTGCTTTTACCAggcgg                                                                                      >  1:786796/1‑150 (MQ=255)
                                         ccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAtt                                                                                   <  2:462649‑M2/62‑1 (MQ=255)
                                           accaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTa                                                                                 >  1:842892‑M2/87‑150 (MQ=255)
                                               atcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACAt                                                                             <  2:312844‑M2/68‑1 (MQ=255)
                                                         gtatagctatgaaaccggtgcggcctaatggagccgttgctgcattgagtttagctcctgctcccaacccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGc                                                                   <  2:791273‑M2/78‑1 (MQ=255)
                                                               ctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTa                                                             >  1:539786‑M2/67‑150 (MQ=255)
                                                                   gaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAcgcg                                                         <  1:780334‑M2/88‑1 (MQ=255)
                                                                   gaaaccggtgcgtcctaatgtagccgttgcggcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAcgcg                                                         <  1:197335‑M2/88‑1 (MQ=255)
                                                                               tcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTAGGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATg                                             >  2:215129‑M2/51‑150 (MQ=255)
                                                                                    atgtagccgtttctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAg                                        <  1:248398‑M2/105‑1 (MQ=255)
                                                                                          ggACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAAcaca                                  <  1:281132/150‑1 (MQ=255)
                                                                                              acCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATc                              <  2:256862/150‑1 (MQ=255)
                                                                                                 tgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATCTACGCGTATCAACAAATGACCAGAACACAGATCTAc                           >  1:758762‑M2/33‑150 (MQ=255)
                                                                                                  gcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGCTCTACa                          >  1:23458‑M2/32‑150 (MQ=255)
                                                                                                  gcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACa                          <  2:717494‑M2/119‑1 (MQ=255)
                                                                                                  gcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACa                          >  2:811818‑M2/32‑150 (MQ=255)
                                                                                                  gcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACa                          >  2:552174‑M2/32‑150 (MQ=255)
                                                                                                          ttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATg                  <  1:735861‑M2/127‑1 (MQ=255)
                                                                                                          ttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATg                  >  1:453024‑M2/24‑150 (MQ=255)
                                                                                                               ctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAGCAGGGGGATTTAACATGCTGATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTg             >  2:68068‑M2/19‑150 (MQ=255)
                                                                                                               caccttctcccaaaccaaCGTATATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAGTGGCTATCTACGCGTATCAACAAATGACCAGAACACAGACCTACAACGTAATGCGCTg             >  2:924144‑M2/19‑150 (MQ=255)
                                                                                                                tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGa            >  2:68923‑M2/18‑150 (MQ=255)
                                                                                                                    tctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACtg        >  1:368080‑M2/14‑150 (MQ=255)
                                                                                                                      tcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACtgtg      >  1:319346‑M2/12‑150 (MQ=255)
                                                                                                                      tcccaaaccaaCGTTTAGGAAAATGAAGAAATAACCAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGACCACAGATCTACAACGTAATGCGCTGATCtgtg      >  2:240209‑M2/12‑150 (MQ=255)
                                                                                                                          aaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCCGAACACAGATCTACAACGTAATCCGCTGAACTGTGGAgg  >  2:247442‑M2/8‑150 (MQ=255)
                                                                                                                                 |                                                                                                                                              
TTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGG  >  NC_000913_3_bme_pgi/1209490‑1209761

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.