New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi 1207790 =15 (0.270)29 (0.580) 27/246 1.2 76.0% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913_3_bme_pgi 1209619 = 5 (0.100)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

GTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGC  >  NC_000913_3_bme_pgi/1209477‑1209765
                                                                                                                                              |                                                                                                                                                  
gcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatggaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTAt                                                                                                                                              >  1:575001‑M2/143‑149 (MQ=255)
gcagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTAt                                                                                                                                              >  1:698890‑M2/143‑149 (MQ=255)
     tgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaaa                                                                                                                                         <  1:350002‑M2/12‑1 (MQ=255)
        ttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATga                                                                                                                                      >  1:137195‑M2/135‑149 (MQ=255)
        gtttgccgcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATga                                                                                                                                      <  1:47297‑M2/15‑1 (MQ=255)
              ggcattccacccccactgaatgattactgattgctcaataccaccaatcatatccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                                                                                                                <  2:975786‑M2/21‑1 (MQ=255)
               gCCAGATCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTTGCGTAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCATTCTCCCAAACCAACGTTTATGAAAATGGAGAAATa                                                                                                                               >  2:649075/1‑149 (MQ=255)
                       cccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaacgctgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCaa                                                                                                                       >  1:939359‑M2/120‑149 (MQ=255)
                                   gattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc                                                                                                           >  2:390850‑M2/108‑149 (MQ=255)
                                         tgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                                                                                                     <  1:284236‑M2/48‑1 (MQ=255)
                                             tgctcaataccaccaatcataaacggtatagctatgaaaccggtgcgttttaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAc                                                                                                 <  2:698890‑M2/52‑1 (MQ=255)
                                               ggCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCa                                                                                               <  1:653689/149‑1 (MQ=255)
                                                 caataccaccaataataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgg                                                                                             <  2:286687‑M2/56‑1 (MQ=255)
                                                  aataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                                                                                            <  2:740917‑M2/57‑1 (MQ=255)
                                                           aatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                                                                                                         <  2:824120‑M2/44‑1 (MQ=255)
                                                           aatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                                                                                                         >  1:824120‑M2/84‑127 (MQ=255)
                                                                gccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTACTGCTTTTAACAGAGAGATTTAACATGCt                                                                              >  2:410269/1‑149 (MQ=255)
                                                                       tatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAACCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATCCTTCTTGGc                                                                       >  2:486707‑M2/72‑149 (MQ=255)
                                                                                 aaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAcgcg                                                             >  1:2869‑M2/62‑149 (MQ=255)
                                                                                    ccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTAt                                                          <  1:217240‑M2/91‑1 (MQ=255)
                                                                                       gggcgttctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaa                                                       <  2:33668‑M2/94‑1 (MQ=255)
                                                                                                aatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGGTATGTACGCGTATCAACAAATGAcc                                              >  1:635854‑M2/47‑149 (MQ=255)
                                                                                                  tgtagccgttgctgcgttgagttttgctccttctcccaaaccaaCGTTTATGAAGAGGAAGAAATAACAAGGAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTGTTGGCTATGTACGCGTAACAACAAATAGCCAg                                            >  1:298793‑M2/45‑149 (MQ=255)
                                                                                                      gccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAaca                                        >  1:472522‑M2/41‑149 (MQ=255)
                                                                                                             ctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCt                                 >  2:777481‑M2/34‑149 (MQ=255)
                                                                                                             ctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCt                                 <  1:390850‑M2/116‑1 (MQ=255)
                                                                                                              tgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGc                                                                       <  2:333495‑M2/78‑1 (MQ=255)
                                                                                                              tgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGc                                                                       >  1:333495‑M2/33‑110 (MQ=255)
                                                                                                              tgcattgagttttgctccttctcccaaaccaaCGCTTAGGAAAGTGAAGAAATAACCAGCACATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTTGCTATGTACGCGTCTCACCAAATGACGACACCACAGATCTa                                >  2:979709‑M2/33‑149 (MQ=255)
                                                                                                                cattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTAGCAACAAATGACCAGAACACAGATCTACa                              <  2:472522‑M2/119‑1 (MQ=255)
                                                                                                                  ttgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAAc                            >  2:484331‑M2/29‑149 (MQ=255)
                                                                                                                  ggTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAAc                            <  1:344823/149‑1 (MQ=255)
                                                                                                                   gTAATGCAGAGCCTTCTCCCAAACCAACGTCTATGAAAATGAAGAAATAACAAGCAAATGGCATCAATCCTGCTTTTAGCAGGGGGATTTAACATGCTTATTGGATAGGTACGCGTATCAACACATGACCAGAACACAGATCTCCAACg                           >  2:343658/1‑149 (MQ=255)
                                                                                                                          tgctccttctcccaatccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCTTGTACCCCGTTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATgcg                    <  1:486707‑M2/129‑1 (MQ=255)
                                                                                                                                   tcccaaaccaaCGTTTATGAAAATTAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGGTTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACtgt           <  1:922272‑M2/138‑1 (MQ=255)
                                                                                                                                   tcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTCATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGGGCTGGACtgt           >  1:725405‑M2/12‑149 (MQ=255)
                                                                                                                                   acccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACtgt           <  2:137195‑M2/138‑1 (MQ=255)
                                                                                                                                    cccacaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATTTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACtgtg          <  2:298793‑M2/139‑1 (MQ=255)
                                                                                                                                       aaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACTTAATGCTCTGACATGTGCAg       >  2:272719‑M2/8‑149 (MQ=255)
                                                                                                                                         accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGa     >  1:159629‑M2/6‑149 (MQ=255)
                                                                                                                                            aaCTTTTATGAAAATGAAAAAATAACAAGCAAATGGCATCATTCCTGTTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGc  <  2:575001‑M2/147‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
GTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGC  >  NC_000913_3_bme_pgi/1209477‑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.