New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =11 (0.200)40 (0.810) 34/232 NT 83.6% coding (290/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 1209619 = 6 (0.120)pseudogene (37/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

TCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCA  >  NC_000913/1209478‑1209759
                                                                                                                                             |                                                                                                                                            
cagatgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttactaattctcccaaaccaaCg                                                                                                                                             >  1:516298‑M2/142‑143 (MQ=255)
    tgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCgtttg                                                                                                                                         >  1:400922‑M2/138‑142 (MQ=255)
         ctgcggcatcccacccccaccgaatgattactgattgctccataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaa                                                                                                                                    <  2:633059‑M2/11‑1 (MQ=255)
          tgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccagtgcgtcctaatgtagccgttgctgcattgagttttgcgcctgctccccaaccaaCGGatatgcaga                                                                                                                                   >  1:276109‑M2/132‑139 (MQ=255)
           tttGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTTACTGAAGGCCACGGAACCGGCCCAACCGCAGTTAGTGCAGAGCCGAGTCCCACACAAACGTATATGaaact                                                                                                                                  >  2:952235/1‑141 (MQ=255)
              gcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGCAAAtgat                                                                                                                               >  1:215029‑M2/128‑142 (MQ=255)
               cattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaag                                                                                                                              <  1:1106477‑M2/17‑1 (MQ=255)
               cattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccgatgcgtcctaatgtagccgttgctgcattgaggtttgatccttatccgaaacaaaCGTTTATGAAAATggag                                                                                                                              >  1:304878‑M2/127‑143 (MQ=255)
                       ccccactgaatgattactgatttctcaataccaccaaccataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACa                                                                                                                      <  2:928162‑M2/25‑1 (MQ=255)
                        cccactgaatgattactgattgctcaatagcaccagtcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                                                                                                     <  2:688892‑M2/26‑1 (MQ=255)
                         ccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctccccaaccgaCGTTTATGACCATGTAGAAATACCAAg                                                                                                                    >  1:936771‑M2/117‑143 (MQ=255)
                          cactgaatgattactgagtgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccattgctgaattgagttttgctccttctcccaaacacaAGTTGATGAAAAGGAAGAAATAACAAGc                                                                                                                   >  2:367966‑M2/116‑143 (MQ=255)
                               aatgagtactggttgctcaataccaccgatcataccgggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagtttttctccttctcccaaaccatCGTTTATGAAAATGAAGAAATAACAAGCAAATg                                                                                                              <  1:832988‑M2/33‑1 (MQ=255)
                                  gattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtaaccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                                                                                                           >  1:819460‑M2/108‑143 (MQ=255)
                                            tgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCt                                                                                                <  2:507160‑M2/46‑1 (MQ=255)
                                              ctcaataccaccaatcataaccggtatagctatgaaaccggtgtgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                                                                                               >  2:655577‑M2/96‑143 (MQ=255)
                                              ctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAACAACAAGCAAATGGCATCAGTCCTGCtt                                                                                               >  1:1027479‑M2/96‑143 (MQ=255)
                                              ctcaagaccacccagcataaccggtatagctatgaaaccggtgcgtcctaaggtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                                                                                               <  2:807225‑M2/48‑1 (MQ=255)
                                               tcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttt                                                                                              >  1:18581‑M2/95‑143 (MQ=255)
                                                 aataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa                                                                                            >  1:756974‑M2/93‑143 (MQ=255)
                                                       cccaatcatacccgggttagctatgaaaccgttgcgtcctaatgtagccgttgctgcattgagttgtgctccatctcccgtgtcatCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                                                                                      <  2:443520‑M2/57‑1 (MQ=255)
                                                       accaatcatatccggtatagcgatgaaaccggtgcgtcctattgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                                                                                      <  2:868339‑M2/57‑1 (MQ=255)
                                                       accaatcatacccggtatagctatgaaaccggtgcgtcctcatgtagccgttgctgcattgagttttgctccttcttccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                                                                                      <  1:463425‑M2/57‑1 (MQ=255)
                                                        cccatcacaccctgtattgctatgactccggcgcgtcctcatgtagccgttgctgcagtgagtttggctccttctcccaaaccacCGTTTATGAACATGACGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgggg                                                                                     <  2:685022‑M2/58‑1 (MQ=255)
                                                          aatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcgcaaaccaaCGTTTATGGAAATGAAGAAATAACGAGCAAATGGCATCATTCGTGCTTTTACCAGGGGGa                                                                                   >  1:14114‑M2/84‑143 (MQ=255)
                                                               taaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa                                                                              >  1:168355‑M2/79‑143 (MQ=255)
                                                                 accggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACa                                                                            >  1:585476‑M2/77‑143 (MQ=255)
                                                                    ggtatagctatgaaaccggtgcgtcctcatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGc                                                                         <  1:19160‑M2/70‑1 (MQ=255)
                                                                           ctatgaaaccgatgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATAGCATCATTCCTGCTAATACCAGGGGGATTTAACATGCTTATTgg                                                                  >  1:569888‑M2/67‑143 (MQ=255)
                                                                            tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGc                                                                 >  1:135395‑M2/66‑143 (MQ=255)
                                                                            tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGc                                                                 >  1:194273‑M2/66‑143 (MQ=255)
                                                                                   gCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTCTTGGCTATGTAc                                                          >  1:954236/1‑143 (MQ=255)
                                                                                     ggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctccccacccacCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCATGCACCAGGGTGATTTAACATGCTTATTGGCTATGTAcgc                                                        <  2:187717‑M2/87‑1 (MQ=255)
                                                                                       tgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATCCTTATTGGCTATGTACGCGt                                                      >  1:198881‑M2/55‑143 (MQ=255)
                                                                                          gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGAATTTAACATGCATACTGGCTATGTACGCGTATc                                                   >  2:732290/1‑143 (MQ=255)
                                                                                           tcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATca                                                  >  2:449444‑M2/51‑143 (MQ=255)
                                                                                           ccACGGAACCGGGACACCAACAGGTAATGCAGAGCCTGCTCCCAAACCAACGGTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTGATTGGCTATGTACGCGTATca                                                  >  1:1002117/1‑143 (MQ=255)
                                                                                            cctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaa                                                 >  1:1028974‑M2/50‑143 (MQ=255)
                                                                                             ctaatgtagcctttgctgtattgagttttgctccttctcccaaaacaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaac                                                <  2:1106717‑M2/95‑1 (MQ=255)
                                                                                             ctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTatcagc                                                >  2:704213‑M2/49‑141 (MQ=255)
                                                                                             ctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaac                                                <  1:464833‑M2/95‑1 (MQ=255)
                                                                                              taatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaaca                                               <  1:655577‑M2/96‑1 (MQ=255)
                                                                                                    cGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGAc                                         <  1:872448/143‑1 (MQ=255)
                                                                                                     gccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATACACAAATGAcc                                        >  1:661053‑M2/41‑143 (MQ=255)
                                                                                                              aCAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAg                               >  1:900656/1‑143 (MQ=255)
                                                                                                               aattgaattttgctccgcctccaaaaccaaCGTTAATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGa                              <  1:103840‑M2/113‑1 (MQ=255)
                                                                                                                attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACACGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACAAGAACAGAGAt                             >  1:976108‑M2/30‑143 (MQ=255)
                                                                                                                    agttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGCCTATGTACGCGTATCGACAAATGACCAGAACACAGATATAc                         >  2:869897‑M2/26‑143 (MQ=255)
                                                                                                                           ctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTaa                  <  2:18581‑M2/125‑1 (MQ=255)
                                                                                                                              cttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTCCAACGTAATgc               >  1:12629‑M2/16‑143 (MQ=255)
                                                                                                                              cTTCTCCCAAACCAACGTTTAGGAAAATGAAGAAATAACACGCAAATGGCATCATTCATAATTGTACCAGGGGGATTTAACATGCTTATTGCCTATGTACGGGTATCAGCAAATGACAGGAACGCAGATCTACCACGTAATgc               >  2:16848/1‑143 (MQ=255)
                                                                                                                                   cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGa          >  1:1027643‑M2/11‑143 (MQ=255)
                                                                                                                                           atcGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCa  <  1:769218/141‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                            
TCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCA  >  NC_000913/1209478‑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.