New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =4 (0.140)13 (0.520) 10/232 NT 85.2% coding (290/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 1209619 = 1 (0.040)pseudogene (37/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

GTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTG  >  NC_000913/1207651‑1207930
                                                                                                                                           |                                                                                                                                            
cacagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttaattctcatgaacgTTgg                                                                                                                                           <  1:502913‑M1/4‑1 (MQ=255)
cacagttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttggtatttcttcattttcataaacgTTgg                                                                                                                                           <  1:207192‑M1/4‑1 (MQ=255)
 acagttcagctcatttcgttgtaaatctatgttcttgtcctttgttattattcgtacatagccaataatcatgttttttccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGt                                                                                                                                          >  2:197017‑M1/139‑143 (MQ=255)
    aGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTG                                                                                                                                       <  1:336900‑M1/143‑1 (MQ=255)
     ttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTgg                                                                                                                                      >  1:301866‑M1/135‑143 (MQ=255)
     ttcagcgcattacgttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTgg                                                                                                                                      >  2:288779‑M1/135‑143 (MQ=255)
                      tagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTc                                                                                                                     >  2:349462‑M1/118‑143 (MQ=255)
                       agatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCa                                                                                                                    >  1:230679‑M1/117‑143 (MQ=255)
                         agccgtgttctggtcctttgttgatacgcgtacatagcgaataagcatgttaaattcccctggtaaaaccaggaatgatgccatttggttgttatttcttcattttcatacacgTTGGTTTGGGAGAAGGAGCAAAACTCAAt                                                                                                                  >  2:236032‑M1/115‑143 (MQ=255)
                           gagtattttggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTCGGTTTGGGGGCAGGAGCAAAACTCAATgc                                                                                                                <  2:232791‑M1/31‑1 (MQ=255)
                             ggggCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttg                                                                                                              <  2:341826‑M1/143‑18 (MQ=255)
                                           tgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTa                                                                                                >  1:454786‑M1/97‑143 (MQ=255)
                                           tggtgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcagtttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTa                                                                                                >  1:174330‑M1/97‑143 (MQ=255)
                                                    gcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCAcc                                                                                       <  1:269856‑M1/56‑1 (MQ=255)
                                                      gCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAGAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgttccggttccgtgttcttctgc                                                                                     >  1:413132‑M1/1‑101 (MQ=255)
                                                                    tccatgttcactccccctggtcaaagcgggaatgatgccatttgcttgttatttcttcattttcataaacgTTGTTTTGGGAGCAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATAc                                                                       >  1:69007‑M1/72‑143 (MQ=255)
                                                                      cctgttaaagccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCg                                                                     >  1:376839‑M1/70‑143 (MQ=255)
                                                                       gTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGA                                                                                                                                    <  2:530230‑M1/79‑1 (MQ=255)
                                                                       gTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGA                                                                                                                                    >  1:530230‑M1/1‑79 (MQ=255)
                                                                                   tgacccGTGTGGCGAGATCAACTCGGAGGGCACTGTGCAAACGGCTCGCGACAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATGAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggtg                                                        <  2:374458/139‑1 (MQ=38)
                                                                                    ctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgGTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTGATGattggtag                                                       <  1:479527‑M1/88‑3 (MQ=255)
                                                                                    aaTCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaa                                                       >  1:537645‑M1/1‑71 (MQ=255)
                                                                                        cGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggt                                                   <  2:390868‑M1/143‑77 (MQ=255)
                                                                                         gTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGGTATGATTGGTGGTATTg                                                  >  1:416072/1‑143 (MQ=255)
                                                                                           ttGGCGATGTCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgcc                                                <  2:366298‑M1/143‑80 (MQ=255)
                                                                                             ggCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCa                                              <  1:508637/143‑1 (MQ=255)
                                                                                                ggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAAATCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGGAATc                                           <  2:38225‑M1/100‑1 (MQ=255)
                                                                                                  aatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAg                                         >  2:110565‑M1/42‑143 (MQ=255)
                                                                                                     cAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctg                                      <  1:559893‑M1/143‑90 (MQ=255)
                                                                                                          cGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagcaactccgccaacaggctggctgaaatgcaaaagtgccgcttcttctgctgaa                                 >  2:288596‑M1/1‑49 (MQ=255)
                                                                                                              tGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgcttgttcttctgaagaat                             >  1:75917‑M1/1‑45 (MQ=255)
                                                                                                               ggCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATGAGGACGCACCGGTTTCATAACTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATGCAGTg                            >  2:110839/1‑143 (MQ=255)
                                                                                                                      tgCAAACGGCTGTCGAAACCCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaa                     <  2:332026‑M1/143‑107 (MQ=255)
                                                                                                                          aaCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactgg                 <  1:198854‑M1/143‑111 (MQ=255)
                                                                                                                          aaCGGCTCTCGAAAACCTTGGTGTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactgc                 >  2:124746‑M1/1‑33 (MQ=255)
                                                                                                                                  tcacaaactTTGGTTTGGGAGAAGGAGCAGAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCaaaa         >  2:334824‑M1/10‑143 (MQ=255)
                                                                                                                                  tcGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggcactggcataggct         >  1:479758‑M1/1‑25 (MQ=255)
                                                                                                                                  tcGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgccttttctgctgaagcatacccggaactggccaaggct         >  1:187866‑M1/1‑25 (MQ=255)
                                                                                                                                    gAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccacgccgccaacaggctggctgcagtgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggctta       >  1:420960‑M1/1‑23 (MQ=255)
                                                                                                                                         ccTGGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccga  <  2:500943‑M1/143‑126 (MQ=255)
                                                                                                                                           |                                                                                                                                            
GTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTG  >  NC_000913/1207651‑1207930

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.