New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi 1207790 =17 (0.410)31 (0.840) 29/258 0.3 76.4% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913_3_bme_pgi 1209619 = 4 (0.110)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

TTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGA  >  NC_000913_3_bme_pgi/1209471‑1209767
                                                                                                                                                    |                                                                                                                                                    
ccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaac                                                                                                                                                      >  1:448230‑M2/149‑149 (MQ=255)
    aggcagatgcttttgcggcattccatccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGttt                                                                                                                                                  <  1:258168‑M2/5‑1 (MQ=255)
      gTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCGTGTCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTAt                                                                                                                                                >  1:391228/1‑149 (MQ=255)
              ttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGCTTATGAAAATga                                                                                                                                        >  1:567374‑M2/135‑149 (MQ=255)
               gCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCAGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGACTTCTCCCCAACCAACGTTTATGAAAATgaa                                                                                                                                       >  2:304543/1‑149 (MQ=255)
                    gtgattccacccccactgcatgattactgattgctcaataccaacaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                                                                                                                  <  2:467117‑M2/21‑1 (MQ=255)
                    ggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgcagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTTTGAAAATAAAGAAAt                                                                                                                                  >  2:651355‑M2/129‑149 (MQ=255)
                      ccAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATCAAAAAATac                                                                                                                                >  2:713939/1‑148 (MQ=255)
                            acccccactgaatgattactgattgctcaatacttcagatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgccttgagttttgatccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCa                                                                                                                          >  1:358831‑M2/121‑149 (MQ=255)
                            acccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCa                                                                                                                          >  2:426003‑M2/121‑149 (MQ=255)
                                ccactgaatgtttactgattgcttaatactatcaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgaggtttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATg                                                                                                                      <  1:258646‑M2/33‑1 (MQ=255)
                                    tgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGt                                                                                                                                                    <  1:118025‑M2/3‑1 (MQ=255)
                                    tgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGt                                                                                                                                                    >  2:118025‑M2/113‑115 (MQ=255)
                                       atgattactgattgctcaacaccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatcat                                                                                                               >  1:223592‑M2/110‑149 (MQ=255)
                                        tgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                                                                                                              >  1:188807‑M2/109‑149 (MQ=255)
                                             actgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                                                                                                         >  2:503832‑M2/104‑149 (MQ=255)
                                               tgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtacccgttgctgcattgagttttgctcttgtgatccaaccaaCGTTTATGAAAATGAAGAAATACCAAGCAACTGGCATCATTCCTGCtt                                                                                                       >  1:352027‑M2/102‑149 (MQ=255)
                                                aaaGCGGCACCATGGCATTTCAGCCAGCCTGTTGGCGGAGGGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttt                                                                                                      <  2:482352/149‑1 (MQ=255)
                                                 attgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtttt                                                                                                     >  1:433816‑M2/100‑149 (MQ=255)
                                                             caccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATAGCATCATTCCTGCTTTTACCAGAGGCAtt                                                                                         >  2:40272‑M2/88‑149 (MQ=255)
                                                              accaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt                                                                                        >  1:271126‑M2/87‑149 (MQ=255)
                                                              accaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt                                                                                        >  2:128354‑M2/87‑149 (MQ=255)
                                                                        accggtatagctatgaaacctgtgcgtcctaatgtagccgttgctccattgagttttgctccttctccccaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCCTa                                                                              >  2:16446‑M2/77‑149 (MQ=255)
                                                                            gtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTgg                                                                          <  1:241955‑M2/77‑1 (MQ=255)
                                                                                     tgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGAATTTCACATACTTATTGTCTATGTAcg                                                                 >  2:441928‑M2/64‑149 (MQ=255)
                                                                                          ctggtgcgtcctaatgtagccgttgctgcattgagttttactccttctcccaaaccaaCGTTTATGAAAATGAATAAATAACAAGCAAATGGCATCAGTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTAt                                                            <  1:374708‑M2/91‑1 (MQ=255)
                                                                                             gtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATACTTATTGGCTATATACGCGTATcaa                                                         >  1:183958‑M2/56‑149 (MQ=255)
                                                                                             gcgcgtcctaatgcagccgttgctgcattgaggtttgctcctgctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaa                                                         <  1:470735‑M2/94‑1 (MQ=255)
                                                                                                       atgtagccgttgctgcattgagtttggcgccttctcccaaaccaaCGTGTATGAAAATGAAGAAATAACAAGCAAAGGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCa                                               <  1:341637‑M2/104‑1 (MQ=255)
                                                                                                         gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGa                                             >  1:615996‑M2/44‑149 (MQ=255)
                                                                                                          tcgccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGaa                                            <  1:58410‑M2/107‑1 (MQ=255)
                                                                                                            gccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACGAGGGGGATTTAACATGCTAATTGGCTATGTACGCGTATCAACAAATGACAagaaaa                                          >  1:314788‑M2/41‑147 (MQ=255)
                                                                                                                   ctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAAAACACAGATCt                                   >  2:51175‑M2/34‑149 (MQ=255)
                                                                                                                        ttgagttttgctccttctcccaaaccaaCTGTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGTATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAAc                              <  2:646132‑M2/121‑1 (MQ=255)
                                                                                                                                tgttccttctcccaaaccaaCGATTATGAAAATGAAGAACTAACGAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATgcg                      <  1:686246‑M2/129‑1 (MQ=255)
                                                                                                                                     cttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACACCGTAATGCGCTGaa                 >  2:463547‑M2/16‑149 (MQ=255)
                                                                                                                                          cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGGAAACAGATTTACAACGTAATGCGCTGAAAtgtg            >  2:310140‑M2/11‑149 (MQ=255)
                                                                                                                                           ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGc           <  1:206947‑M2/140‑1 (MQ=255)
                                                                                                                                            caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCa          >  1:109777‑M2/9‑149 (MQ=255)
                                                                                                                                              aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCCACAAAAGACCAGAACCCAAATCTACAACGTAATGCGCTGAACTCTGCAgg        >  1:80525‑M2/7‑149 (MQ=255)
                                                                                                                                               accaaCGTTTATGAAAATGAAGAAATAACAAGCAAGTGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGa       >  1:290032‑M2/6‑149 (MQ=255)
                                                                                                                                               accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAAt                                                    >  2:53431‑M2/6‑104 (MQ=255)
                                                                                                                                               accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAAt                                                    <  1:53431‑M2/99‑1 (MQ=255)
                                                                                                                                                    cGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAAGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGa  <  1:461205/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGA  >  NC_000913_3_bme_pgi/1209471‑1209767

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.