New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =4 (0.080)36 (0.860) 30/252 0.4 86.2% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 1209619 = 8 (0.190)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

GGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCA  >  NC_000913/1209480‑1209759
                                                                                                                                           |                                                                                                                                             
cgATAAGCCGCGGCCAGTTCCGTGTCTGCTCCACCAGCAGCAGCGGCACCATGGCATTGCAGCCCGTCTGTTCGCGGAGCGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAAGGtt                                                                                                                                            <  2:641645/142‑1 (MQ=255)
       ccTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGATTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGaaa                                                                                                                                     >  1:183203/1‑143 (MQ=255)
        tgcggcagtccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaaa                                                                                                                                    <  2:817907‑M2/12‑1 (MQ=255)
            gcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaa                                                                                                                               <  1:110596‑M2/16‑1 (MQ=255)
            gcattacaccaccaccgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaa                                                                                                                               <  1:18579‑M2/16‑1 (MQ=255)
                 ccacccccactgaatggttactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaacgaaCGTTTCTGAAAATGAGGAAAt                                                                                                                          >  2:456881‑M2/123‑143 (MQ=255)
                     gggTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACACACGTTTATAAAACTGAAGAAATcaca                                                                                                                      >  2:697084/1‑143 (MQ=255)
                     ccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGCAGAAATAACa                                                                                                                      >  1:414835‑M2/119‑143 (MQ=255)
                      cccactgaatgattactgtttgctccataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                                                                                                     <  2:407301‑M2/26‑1 (MQ=255)
                      cccactgaatgattactcattgctcaataccaccaatcataaccggtatggctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                                                                                                     <  1:756528‑M2/26‑1 (MQ=255)
                                  cagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTTGCTGAAGGCCACGGAACCGGGACACCAACAGGTACTGCAGAGCCTTCTCCCACACCAACGTTTATGAACATGAAGAAATAACAAACAAATGGcctc                                                                                                         >  2:417454/1‑143 (MQ=255)
                                  cagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGCTTGTGAAAATAAGGATAGAACAAGCAAAGGGcatc                                                                                                         >  2:405925/1‑143 (MQ=255)
                                      ggattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgcagccgttgcttcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc                                                                                                     <  2:377676‑M2/42‑1 (MQ=255)
                                        attgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCt                                                                                                   >  1:370365‑M2/100‑143 (MQ=255)
                                                 tcccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAcc                                                                                          <  1:373971‑M2/53‑1 (MQ=255)
                                                 taccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgatgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAAGGGCATCATTCCTGCTTttagc                                                                                          >  1:769478‑M2/91‑141 (MQ=255)
                                                 taccaccaagcagaaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAcc                                                                                          <  1:832205‑M2/53‑1 (MQ=255)
                                                     acccatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                                                                                      <  2:88494‑M2/57‑1 (MQ=255)
                                                        tttCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCACATGGCATCATTCCTGCTTTTACCAGGGGGa                                                                                   >  1:552278/1‑143 (MQ=255)
                                                        aagcacaaccggcagagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa                                                                                   <  1:52290‑M2/60‑1 (MQ=255)
                                                         atcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAt                                                                                  <  2:618836‑M2/61‑1 (MQ=255)
                                                                ccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgccttgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaacgt                                                                           >  1:306929‑M2/76‑141 (MQ=255)
                                                                 cggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcatagagtgttgctcattctcccaaaccaaCGTTTATGAGAATGACGAAATAACACGCAAATGGCATCCGTCCAGCGTTTACAAGGGGGAGTTAACATg                                                                          >  2:58545‑M2/75‑143 (MQ=255)
                                                                      tagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAt                                                                     <  2:707125‑M2/74‑1 (MQ=255)
                                                                        ggCGGAGTGGCTGAAGGCCACGGAACCGTGACCCCAACTTGTAATGCAGAGCCTGCTCCCAAACCAACGTTTATACAAATAAAGAAATAATAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTg                                                                   <  1:396060/143‑1 (MQ=255)
                                                                        gctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAATATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTg                                                                   <  1:457446‑M2/76‑1 (MQ=255)
                                                                          tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGCACTAACAACCAAATGGCATCATTCCGGCTTTTACCAGACGACTGGAACATGCTTAttgcc                                                                 >  2:173478‑M2/66‑141 (MQ=255)
                                                                          tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGc                                                                 >  2:434241‑M2/66‑143 (MQ=255)
                                                                          tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctaccaaacaaaCGATTATGAAAATGAAGAAATAACAAGCAAATGGCACCATTCCTGCTTTTACCAGGGGGATTCACAATGCTTAttgac                                                                 >  2:341992‑M2/66‑141 (MQ=255)
                                                                           atgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATCTAACATGCGTActggat                                                                >  2:462881‑M2/65‑141 (MQ=255)
                                                                                  cggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATACCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAcg                                                         >  1:326597‑M2/58‑143 (MQ=255)
                                                                                     tgcgtcctaaggtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGt                                                      <  1:649585‑M2/89‑1 (MQ=255)
                                                                                       cgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTAt                                                    >  2:191282‑M2/53‑143 (MQ=255)
                                                                                          cctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCGTTTACCAGGGGGATTTAACATGCTTATTGGCTCTGTGCGagtctccg                                                 >  2:451822‑M2/50‑141 (MQ=255)
                                                                                           ctaatgtggccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGCGATTTAACATGCTTATTGGCTATGTACGCGTATcaac                                                <  2:589266‑M2/95‑1 (MQ=255)
                                                                                             aatgtagccgttgctgcattgagttttgctccttctcccaaaccacCGTTTATGACAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACCCGAAtcaacac                                              >  2:632029‑M2/47‑142 (MQ=255)
                                                                                                  agccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGAc                                         <  1:466558‑M2/102‑1 (MQ=255)
                                                                                                    ccgttgctgtattgagttgtgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCa                                       <  2:585441‑M2/104‑1 (MQ=255)
                                                                                                    ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCa                                       >  1:594418‑M2/40‑143 (MQ=255)
                                                                                                    ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAATTAACATGCTTATTGGCTATGTACGCGGATCAAAAAATGACCa                                       >  2:250954‑M2/40‑143 (MQ=255)
                                                                                                         gctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACAAGGGGAATTTAACATGCTTATTGGCTATGTACACGTATCAAGAAATAACCAGAaca                                  >  2:434090‑M2/35‑143 (MQ=255)
                                                                                                           tgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGACcaca                                >  1:147213‑M2/33‑143 (MQ=255)
                                                                                                                   aTGCAGAGCCTTCTCCCAAACCAACGGGAATATAATAATTTAAATAACAAGCAAATGGCATCATTCATGGTTTTAGCTAAAAATGATGACATGCTTATTGGCTATTTTCGCGTTTCGACAAAGGACCATAACACAGATCTACa                        >  2:651497/1‑143 (MQ=255)
                                                                                                                         ctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCTTATCAACAAATGACCAGAACACAGATCTACAACGTaa                  >  1:663563‑M2/19‑143 (MQ=255)
                                                                                                                                tcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATCACCAGAACACAGCTCTACAACGTAATGCGCTg           >  2:346764‑M2/12‑143 (MQ=255)
                                                                                                                                 cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTACCATGCTTATTGACTATGTACGCGTATCCACAAATGCCCAGAACACAGATCTACCACgtaatgctccga          >  1:856523‑M2/11‑138 (MQ=255)
                                                                                                                                      accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACtgt     <  1:179286‑M2/138‑1 (MQ=255)
                                                                                                                                        caaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGAATCAACAAATGACCAGAACACAGATCTACAACGtagtgagctgagctcggc   >  2:615863‑M2/4‑139 (MQ=255)
                                                                                                                                         aaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAGTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCa  <  1:375970‑M2/141‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                             
GGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCA  >  NC_000913/1209480‑1209759

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ 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.