New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi 1207790 =10 (0.210)38 (0.900) 37/252 0.2 84.6% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913_3_bme_pgi 1209619 = 5 (0.120)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

TTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGC  >  NC_000913_3_bme_pgi/1209471‑1209765
                                                                                                                                                    |                                                                                                                                                  
ccagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaac                                                                                                                                                    <  1:748001‑M2/1‑1 (MQ=255)
 cagaggcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCg                                                                                                                                                   >  2:331710‑M2/148‑149 (MQ=255)
          atgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaa                                                                                                                                          <  2:546740‑M2/11‑1 (MQ=255)
            gcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAAt                                                                                                                                        <  1:331710‑M2/13‑1 (MQ=255)
               gctgcggcagcccccacccactgaatgtttgctgattgcgcaataccagcaatcataaccggtatagctatgaaagcggtgcgtcctaatgtagccgttgctgcatggagttttggtccttctcccaaaccaaCGTTTATGAAAATgaa                                                                                                                                     <  2:124889‑M2/16‑1 (MQ=255)
                  gcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaagaa                                                                                                                                  <  1:330596‑M2/19‑1 (MQ=255)
                  gcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTGTATGAAAATgaagaa                                                                                                                                  >  1:209987‑M2/131‑149 (MQ=255)
                    ggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaag                                                                                                                                    >  1:383691‑M2/129‑145 (MQ=255)
                    ggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaag                                                                                                                                    <  2:383691‑M2/17‑1 (MQ=255)
                        ttccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtacgtcctaatgtaggcgccgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACa                                                                                                                            >  1:194893‑M2/125‑149 (MQ=255)
                         gTTCCGCGTAGTCTTCAGCAGAAAAAGGGGCACCATTGCATTTCAGCCAGGCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACaa                                                                                                                           <  2:108880/149‑1 (MQ=255)
                            acccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttactgcattgagttttgctccttctcccaaaccacCGTTTATGAAATTAAAGAAATACAAAGCa                                                                                                                        >  2:527496‑M2/121‑149 (MQ=255)
                                ccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGCAAATGAAGAAATAAAAAGCAAATg                                                                                                                    >  1:325626‑M2/117‑149 (MQ=255)
                                         gattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTCATGAAAATGAAGAAATAACAAGCAAATGGCACCATTc                                                                                                           >  2:256546‑M2/108‑149 (MQ=255)
                                              ctgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctattgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCAgcg                                                                                                      >  1:139214‑M2/103‑148 (MQ=255)
                                              ctgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATATCAAGCAAATGGCATCATTCCTGCt                                                                                                      >  1:139166‑M2/103‑149 (MQ=255)
                                              ctgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgccacattgagttctgaactcactcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCt                                                                                                      <  1:823441‑M2/47‑1 (MQ=255)
                                               tgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCCTGAatccttcgt                                                                                                     >  1:570316‑M2/102‑145 (MQ=255)
                                                 attgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAACAAATAACAAGCAAAAGGCATCCTTACTGCtttt                                                                                                   >  2:717549‑M2/100‑149 (MQ=255)
                                                   tgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAc                                                                                                 <  1:22504‑M2/52‑1 (MQ=255)
                                                   tgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTGTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCGGCTTTTAc                                                                                                 >  1:344701‑M2/98‑149 (MQ=255)
                                                     ggCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCa                                                                                               <  1:734738/149‑1 (MQ=255)
                                                     ctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCa                                                                                               >  1:536394‑M2/96‑149 (MQ=255)
                                                      tcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAg                                                                                              >  1:18831‑M2/95‑149 (MQ=255)
                                                                caatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa                                                                                    <  1:475115‑M2/65‑1 (MQ=255)
                                                                 aatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATTGCATCATTCCTGCTTTTACCAGGGGGATTTAAc                                                                                   <  1:496704‑M2/66‑1 (MQ=255)
                                                                  atcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATACCAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACa                                                                                  >  1:359154‑M2/83‑149 (MQ=255)
                                                                  atcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACa                                                                                  <  1:567925‑M2/67‑1 (MQ=255)
                                                                    cataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATg                                                                                <  1:423488‑M2/69‑1 (MQ=255)
                                                                           ggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTg                                                                         >  2:133851‑M2/74‑149 (MQ=255)
                                                                              atagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAGTGGCATCATTCCTGCGTTTACCAGGGGGATTTAACATGCTTTTTGGCt                                                                      >  1:21287‑M2/71‑149 (MQ=255)
                                                                                          gCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAGTTAACATGCTTATTGGCTATGTACGCGTAt                                                          >  1:499008/1‑149 (MQ=255)
                                                                                              aaGGCCAGGGGACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTacatag                                                      <  2:831023/149‑6 (MQ=255)
                                                                                                     taatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGtt                                                                                                                                                 >  1:326257‑M2/48‑51 (MQ=255)
                                                                                                     taatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGtt                                                                                                                                                 <  2:326257‑M2/4‑1 (MQ=255)
                                                                                                       gAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGCCGGGATAACGAGGAATAAGCGGTTATCCTGCTTTTACGAGGGGCATTTGACACGCGTAGTGGCGATGTACGCGTAGCAACAGATGACCa                                             >  1:197803/1‑149 (MQ=255)
                                                                                                        tgtagccgttgctgcattgagttttgctccttctcccaagccgaCTTTGATGAAAATGCAGAAATCCCAATCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATGAACAAATGACCAg                                            >  1:785471‑M2/45‑149 (MQ=255)
                                                                                                          tagccgttgctacattgagtttttctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGaa                                          <  1:300769‑M2/107‑1 (MQ=255)
                                                                                                            gccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATTTAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAaca                                        >  2:765686‑M2/41‑149 (MQ=255)
                                                                                                               gttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCACGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAg                                     <  1:28926‑M2/112‑1 (MQ=255)
                                                                                                                   ctgcattgagttttgcaccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCt                                 >  1:731195‑M2/34‑149 (MQ=255)
                                                                                                                             ttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGCCTATGTACGCGTATCAACAAATGACCAGAACACGGATCTACAACgtagt                       >  2:330955‑M2/24‑147 (MQ=255)
                                                                                                                                  ctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCt                  <  1:720312‑M2/131‑1 (MQ=255)
                                                                                                                                     cttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTCCGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTTaa               <  1:245489‑M2/134‑1 (MQ=255)
                                                                                                                                        ctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACtg            >  1:536809‑M2/13‑149 (MQ=255)
                                                                                                                                          cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCTGAACACAGATCTACAACGTAATGCGCTGGActgtt          <  2:731195‑M2/139‑2 (MQ=255)
                                                                                                                                           ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCCCTGAACTGTGc         >  2:50901‑M2/10‑149 (MQ=255)
                                                                                                                                            caaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCa        <  2:536809‑M2/141‑1 (MQ=255)
                                                                                                                                                  aaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGc  <  1:717549‑M2/147‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TTATTTGTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGC  >  NC_000913_3_bme_pgi/1209471‑1209765

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.