New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =6 (0.120)41 (0.910) 33/230 NT 85.2% coding (290/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 1209619 = 9 (0.200)pseudogene (37/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

GGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCA  >  NC_000913/1209480‑1209759
                                                                                                                                           |                                                                                                                                            
ggATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACAAACGtt                                                                                                                                           >  1:1013603/1‑143 (MQ=255)
 gATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGttt                                                                                                                                          >  2:366666/1‑143 (MQ=255)
     ttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagtcttgctccgtctcccaaaccacCGTTTAAGa                                                                                                                                      >  2:800703‑M2/135‑143 (MQ=255)
        cTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCACACCAACGTTGATGaaaa                                                                                                                                   >  1:840779/1‑143 (MQ=255)
          cggcattccacccccactgaatgattactgattgctcaataccacaaatcataacagggatagcaatgaaaccggtgcttcctaatatcgccgttgatgcattgagttttgcgcagtctgccaaaccaaCGTTTatgaaagta                                                                                                                                 >  2:369621‑M2/130‑140 (MQ=255)
            ccatgccacccccactcaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaa                                                                                                                               <  2:422537‑M2/16‑1 (MQ=255)
                  cacccccactgagcgatgactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                                                                                                         <  2:656259‑M2/22‑1 (MQ=255)
                  cacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgggttttgctcctgttgccaaaccaaCGGTTATgagattgaatagatg                                                                                                                         >  1:701107‑M2/122‑137 (MQ=255)
                            gaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaacaggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcacaaacgaaCGGTTATGAAAATGACGAACTAAGGAgcaata                                                                                                               >  2:750874‑M2/112‑141 (MQ=255)
                            cttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAAt                                                                                                               >  1:60817/1‑143 (MQ=255)
                                  ttactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttcccccaaaccaaCGTTTATGAAAGTGAAGAAATAACAAGCAAACGGcatc                                                                                                         >  2:871238‑M2/106‑143 (MQ=255)
                                  ttactgattgctcaataccaccaatcataaccggtatagctatgaaaacgctgcgtcctaatgtagccgttgctgcattgaggtctgctccctctcccaaaccaaCCTTTATGAAAACAAAGAAATCACAAGCAAATAGcatc                                                                                                         >  2:172087‑M2/106‑143 (MQ=255)
                                   tactgagtgctcactaccaccaatcatacccgctatagctagccaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                                                                                                        <  2:941032‑M2/39‑1 (MQ=255)
                                     ctgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                                                                                                      >  1:1010097‑M2/103‑143 (MQ=255)
                                         ttgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg                                                                                                  >  1:917889‑M2/99‑143 (MQ=255)
                                            ctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATCACAAGCAAATGGCATCATTCCTGCtt                                                                                               >  1:571998‑M2/96‑143 (MQ=255)
                                            ctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                                                                                               >  1:461360‑M2/96‑143 (MQ=255)
                                               aataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgccttgagttttgctccttctcccaaacgatCGAGCATAAAAATGAAGACATAACAAGCAAATGGCATCATTCCTGATTTTa                                                                                            >  2:719116‑M2/93‑143 (MQ=255)
                                               aataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa                                                                                            >  1:920987‑M2/93‑143 (MQ=255)
                                                  accaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgggttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTtacaa                                                                                         >  1:710687‑M2/90‑141 (MQ=255)
                                                         atcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCAACATTCCTGCTTTTACCAGGGGGAt                                                                                  >  1:784885‑M2/83‑143 (MQ=255)
                                                          tcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAtt                                                                                 <  1:41976‑M2/62‑1 (MQ=255)
                                                               accggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACa                                                                            <  1:290979‑M2/67‑1 (MQ=255)
                                                                    tGTTGGCGGAGTGGCTGAAGGCCAGGCAACCGCAACACCAACTGTTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCtt                                                                       <  1:956380/143‑1 (MQ=255)
                                                                       agctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAtt                                                                    <  1:630850‑M2/75‑1 (MQ=255)
                                                                            tgcaaacggtgcgccgtagtgtagccgttgctgcattaagttttgccccttctcccaaaccaaCGTTTACGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTa                                                               <  2:266418‑M2/80‑1 (MQ=255)
                                                                               aaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGt                                                            <  1:891845‑M2/83‑1 (MQ=255)
                                                                                   tGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAcgc                                                        <  2:1013603/143‑1 (MQ=255)
                                                                                       cgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTAt                                                    <  1:885089‑M2/91‑1 (MQ=255)
                                                                                       cgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTAt                                                    <  2:917889‑M2/91‑1 (MQ=255)
                                                                                          cctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTCCGCGTATcaa                                                 >  1:879403‑M2/50‑143 (MQ=255)
                                                                                                 tagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTCGCTATGTACGCGTATCAACAAATGa                                          <  2:998127‑M2/101‑1 (MQ=255)
                                                                                                   gggACACCAACAGATAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCATAACAACAAATGaac                                        >  2:621043/1‑141 (MQ=255)
                                                                                                    ggACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGATATGTAAGCATATCAACAAATGACCa                                       >  1:629530/1‑143 (MQ=255)
                                                                                                    ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCa                                       <  1:693589‑M2/104‑1 (MQ=255)
                                                                                                    ccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCGATGTACGCGTATCAACAAATGACCa                                       >  1:367269‑M2/40‑143 (MQ=255)
                                                                                                     cgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAg                                      >  1:833279‑M2/39‑143 (MQ=255)
                                                                                                          ctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACACATGACCAGAacac                                 >  1:179328‑M2/34‑143 (MQ=255)
                                                                                                          ctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAacac                                 >  1:155407‑M2/34‑143 (MQ=255)
                                                                                                          cAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTCTTGGCTATGTACGCGCCGCAACAAATGACCAGAacac                                 >  2:694853/1‑143 (MQ=255)
                                                                                                           tgcgttgaggtttgctcctggtcccaaaccaaCGTTTGTGAAAATGAAGAAATAACAAGCAAAGGGCATCATTCGTGCCTTTACCAGGGGGATTTAAAATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAAcaca                                >  2:819379‑M2/33‑143 (MQ=255)
                                                                                                           tgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTCGCTATGTACACGTATCAACAAATGACCAGACcaca                                >  1:570689‑M2/33‑143 (MQ=255)
                                                                                                           tgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACAAGAGGGATTTAACATGCTTATTGGCTATGTACGCGCATCAACAAATGACCAGAAcaca                                >  2:592730‑M2/33‑143 (MQ=255)
                                                                                                            gcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAg                               >  2:331869‑M2/32‑143 (MQ=255)
                                                                                                             cattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAAAACAGa                              >  2:844683‑M2/31‑143 (MQ=255)
                                                                                                              attgagttttgctccttctcacaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGAt                             <  2:851260‑M2/114‑1 (MQ=255)
                                                                                                              aGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGAt                             <  1:653027/143‑1 (MQ=255)
                                                                                                               ttgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTCTTGGCTATGTACGCGTATCAACAAATGACCACAACACAGATc                            >  1:944538‑M2/29‑143 (MQ=255)
                                                                                                                tgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCt                           >  1:885120‑M2/28‑143 (MQ=255)
                                                                                                                tgagttttgctcattctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAGTGGCTATGTACGCGTAGCAACAAATGACAAGAACACAGATCt                           >  2:79845‑M2/28‑143 (MQ=255)
                                                                                                                 tAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTa                          >  2:851276/1‑143 (MQ=255)
                                                                                                                    ttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACaa                       <  2:463940‑M2/120‑1 (MQ=255)
                                                                                                                        gctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAGATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTa                   >  1:964787‑M2/20‑143 (MQ=255)
                                                                                                                         ctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTGAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTaa                  >  1:512681‑M2/19‑143 (MQ=255)
                                                                                                                            cttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATgc               >  2:337867‑M2/16‑143 (MQ=255)
                                                                                                                              tctcccaacccagCGTTTATGAGAATGAAGAAAGAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATgcgc             <  2:461360‑M2/130‑1 (MQ=255)
                                                                                                                              tctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATgcgc             <  2:571998‑M2/130‑1 (MQ=255)
                                                                                                                              tctcCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAAGGGGCTTTAACATGCTTATTAGCCCAGCACGCATATCCACACAAGACCAGAACACCCATCTCCAACGTAGTgcac             >  2:836348/1‑141 (MQ=255)
                                                                                                                                 cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGa          >  2:774140‑M2/11‑143 (MQ=255)
                                                                                                                                        caaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGc   <  1:366666‑M2/140‑1 (MQ=255)
                                                                                                                                        caaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGc   >  1:652404‑M2/4‑143 (MQ=255)
                                                                                                                                        caaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAATGAACATACTTCTTGGCTATGTACGAGTATCAAAAAATGACCAGAAAAGAGATCTACAACGTAATGTGAAGCACTGTGc   >  2:168650‑M2/4‑143 (MQ=255)
                                                                                                                                         accGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCa  <  2:840779/141‑1 (MQ=255)
                                                                                                                                         aaCGTTTAGGGAGATGAAGCAAGAACAAGCACATGGCAGCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGTCCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCa  <  2:632341‑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 < 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.