New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =2 (0.050)31 (0.840) 28/250 0.3 92.8% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 1209619 = 3 (0.080)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

GATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCA  >  NC_000913/1209481‑1209759
                                                                                                                                          |                                                                                                                                            
atgcttttgcggcattccacccccactgaatgattacggattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGttt                                                                                                                                          <  1:664636‑M2/5‑1 (MQ=255)
  tcttgtgctgcattccacccccactgaatgactactgattgttgaataccaccaatcatagccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTAt                                                                                                                                        <  2:205309‑M2/7‑1 (MQ=255)
           gcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaa                                                                                                                               <  1:758273‑M2/16‑1 (MQ=255)
             attccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaacccaaCATTTATgaaaataagga                                                                                                                             >  2:109555‑M2/126‑138 (MQ=255)
               gtgccGGGGACTCTTCAGCAGAAAAAGCGGCACCAGTGCATTGCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGaaa                                                                                                                           <  2:674640/140‑1 (MQ=255)
                ccacccccactgaatgattactgattgctcaagaccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                                                                                                          <  2:166827‑M2/21‑1 (MQ=255)
                  ccccccactggatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATaa                                                                                                                        <  1:334701‑M2/23‑1 (MQ=255)
                     ccccctgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                                                                                                     <  2:670998‑M2/26‑1 (MQ=255)
                            aatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATg                                                                                                              <  1:627175‑M2/33‑1 (MQ=255)
                                         tgctcaataccaccaagcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgcgccctctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGc                                                                                                 <  1:92940‑M2/46‑1 (MQ=255)
                                            tcaatacaaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttt                                                                                              >  2:681005‑M2/95‑143 (MQ=255)
                                              aataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTGATGAAAATGAAGAAATAACAAGCAAATGGAATCATTCTTGCTTTTa                                                                                            >  1:208982‑M2/93‑143 (MQ=255)
                                                  ccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAg                                                                                        <  1:441480‑M2/55‑1 (MQ=255)
                                                   tGCATTTCAGCCAGCCTGTTGGCAGAGTGCCTCAACGCCTCTGTACCTGGACACCAACATGTGATGCAGAGCCTTCTCCCCAACCAACGTGTCTGAAAACGACGAGACAACTAGCAAATGGCATCATTCATGCTTTCACCAgg                                                                                       >  2:269454/1‑143 (MQ=255)
                                                    accaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgcggcattgagttttgctccttctcccaaacaaaCGTTTATGACAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                                                                                      >  1:209586‑M2/87‑143 (MQ=255)
                                                        atcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAt                                                                                  <  1:700610‑M2/61‑1 (MQ=255)
                                                            tagccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcagtgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa                                                                              <  1:18014‑M2/65‑1 (MQ=255)
                                                            taaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctgcttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa                                                                              <  1:382382‑M2/65‑1 (MQ=255)
                                                              accggtatagctatgaaaccggggcgtcctaatgtagccgttgctgcattgagatttgctccttctcaccaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAACTGGCATCATGCCTGCATCTACCAGGGGGATTAAacg                                                                            >  2:629753‑M2/77‑142 (MQ=255)
                                                                              aaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGt                                                            <  1:681005‑M2/83‑1 (MQ=255)
                                                                                ccggtgcgtcctaatgtatccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATCGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAc                                                          <  2:763354‑M2/85‑1 (MQ=255)
                                                                                 tcgtgcgtcttaacgtagccgttgctgcgttgagttttgctcctactcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAcg                                                         <  2:20234‑M2/86‑1 (MQ=255)
                                                                                 cggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAACGAAGAAATAACAAGCAAATGGCATCATTCCTGCTCTTACCAGGGGGATTTAACATGCTTATTGGCTCTGTAcg                                                         >  2:323462‑M2/58‑143 (MQ=255)
                                                                                         cctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTCATTGGCTATGTACGCGTATcaa                                                 >  2:82536‑M2/50‑143 (MQ=255)
                                                                                            aatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATGTAACATGCTTATTGGCTATGTACGCGTAGcaacaa                                              >  2:521828‑M2/47‑143 (MQ=255)
                                                                                            aatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTCTACCAGGAGGCTTTAACATGCTTATTGGCTATGTACGcgtatctccac                                              >  2:659148‑M2/47‑138 (MQ=255)
                                                                                             gAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTCTCATGCTTATCGACTATGTACCCGTATCAACaaa                                             >  1:50487/1‑143 (MQ=255)
                                                                                                     gttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGa                                     >  1:579286‑M2/38‑143 (MQ=255)
                                                                                                      gtgctgcattgagttgtgctccttctcccaaaccacCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGaa                                    <  2:428248‑M2/107‑1 (MQ=255)
                                                                                                           gcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAg                               >  2:390453‑M2/32‑143 (MQ=255)
                                                                                                            tattgagttgtgcgccttcccccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGa                              <  2:708585‑M2/113‑1 (MQ=255)
                                                                                                            cattgaggtttgctccttctcccaaaccaaCGTTTATGAAGATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCATGGAGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGa                              <  2:310386‑M2/113‑1 (MQ=255)
                                                                                                                   ttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACACGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACCCAGCTcaacaa                       >  2:531716‑M2/24‑143 (MQ=255)
                                                                                                                    tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAAc                      <  1:654380‑M2/121‑1 (MQ=255)
                                                                                                                           cttctcccaaaccaaAGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGAGCTACAACGTAGTgc               >  2:551331‑M2/16‑143 (MQ=255)
                                                                                                                            ttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATgcg              >  1:723400‑M2/15‑143 (MQ=255)
                                                                                                                            ttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATACCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATgcg              <  2:455676‑M2/129‑1 (MQ=255)
                                                                                                                               tcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCGACAAATGACCAGGGCACGGATCTACAACGCACTGCGCTg           >  1:43784‑M2/12‑143 (MQ=255)
                                                                                                                                 ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGctgca         >  2:195157‑M2/10‑141 (MQ=255)
                                                                                                                                   aagccaaCGGTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACt       <  2:643476‑M2/136‑1 (MQ=255)
                                                                                                                                        aaCGTCTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGACCTGTGCa  <  1:147089‑M2/141‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                            
GATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCA  >  NC_000913/1209481‑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.