New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_bme_pgi 1207790 =7 (0.150)22 (0.520) 21/250 1.3 81.2% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913_3_bme_pgi 1209619 = 4 (0.090)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

TTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAG  >  NC_000913_3_bme_pgi/1209497‑1209760
                                                                                                                          |                                                                                                                                             
acacccccatgaaatgattactgattgctcaataccactaatcataacctgtatagctatgaaacccgtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAg                                                                                                                     <  2:416765‑M2/27‑1 (MQ=255)
  acccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccacCGTTTATGAAAATGAAGAAATAACAAGCa                                                                                                                   >  2:624955‑M2/121‑149 (MQ=255)
  acccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCa                                                                                                                   >  2:402988‑M2/121‑149 (MQ=255)
    ccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAACCaaa                                                                                                                 >  2:290507‑M2/119‑149 (MQ=255)
               gattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccagCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc                                                                                                      <  1:142423‑M2/42‑1 (MQ=255)
                   gAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                                                                                                  <  1:798279/149‑1 (MQ=255)
                   actgtttgcccaataccaccaatcataaccggtatagctatgaatccggtgcgtcctaatgtagccgttgctgcattgacttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                                                                                                  <  1:635760‑M2/46‑1 (MQ=255)
                   actgattgctcaataccaccaatcataaccggtatagctatgaaaccagtacgtcctagtgtagccgttgctgcattgagctttgctccatctcccaaaccaaCTTTGATAAAAATGCAGAAATAACAAGCAACAGGCTTAATTccttc                                                                                                  >  2:166097‑M2/104‑147 (MQ=255)
                           ctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGCAGAAATAACAAGCAAATGGCATCATTCCTGCATTTACCa                                                                                          >  1:443714‑M2/96‑149 (MQ=255)
                            tcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCGAATGGCATCATTCCTGCTTTTACCAg                                                                                         <  2:517181‑M2/55‑1 (MQ=255)
                                 accaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa                                                                                    <  2:824005‑M2/60‑1 (MQ=255)
                                   caccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                                                                                                       >  1:607692‑M2/88‑128 (MQ=255)
                                   caccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                                                                                                       <  2:607692‑M2/41‑1 (MQ=255)
                                     ccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTa                                                                                <  1:31313‑M2/64‑1 (MQ=255)
                                       aatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAAc                                                                              <  1:394238‑M2/66‑1 (MQ=255)
                                              cagccTGTTGGCGGGTTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTa                                                                       <  1:449896/149‑1 (MQ=255)
                                               ccggtatagctatgaaaccggttctttcttaggtagcctttgctgcattgagttttgctccttctcccaaaccacCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGTGGGATTTAACATGCTTAt                                                                      <  2:174084‑M2/74‑1 (MQ=255)
                                                         tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTa                                                            >  2:565605‑M2/66‑149 (MQ=255)
                                                          atgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCCATGTAc                                                           >  1:387422‑M2/65‑149 (MQ=255)
                                                            gaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGACATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGCTTTAACATGCTTATTGGCTATGTAcgc                                                         >  1:429186‑M2/63‑149 (MQ=255)
                                                               accggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTa                                                      <  2:292354‑M2/90‑1 (MQ=255)
                                                                 cTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAATGAATTGATGATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATc                                                    <  1:194125/149‑1 (MQ=255)
                                                                        tcctaatgtaggcgttgttgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAGTATGCTTATTGGCTATGTACGCGTATCAACAAAt                                             <  2:615055‑M2/99‑1 (MQ=255)
                                                                                              ttgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATAACCAGAACACAAATCTACAAc                       >  2:369003‑M2/29‑149 (MQ=255)
                                                                                               gTAATGCAGAGCCTTCTCCCAAACCACCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCCTTTACCAGGGGGATTTAACATGCTTATTGACTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACg                      >  2:681277/1‑149 (MQ=255)
                                                                                                         tccttctaccaaaccaaCGTTGATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTg            >  2:311938‑M2/18‑149 (MQ=255)
                                                                                                                cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa                                                                               >  2:412224‑M2/11‑75 (MQ=255)
                                                                                                                cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa                                                                               <  1:412224‑M2/65‑1 (MQ=255)
                                                                                                                   aaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAg  >  1:232379‑M2/8‑149 (MQ=255)
                                                                                                                       caaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAaca                                   <  1:81713‑M2/109‑1 (MQ=255)
                                                                                                                       caaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAaca                                   >  2:81713‑M2/4‑112 (MQ=255)
                                                                                                                          |                                                                                                                                             
TTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAG  >  NC_000913_3_bme_pgi/1209497‑1209760

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.