New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 120780518 (0.420)13 (0.340) 12/252 2.4 40.6% coding (305/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 = 1209602 22 (0.570)pseudogene (18/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

GGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATG  >  NC_000913/1207663‑1207945
                                                                                                                                              |                                                                                                                                            
cagtactttgtagatctgtgttctggtcatttgttgatacgcgtacctagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAgg                                                                                                                                              >  2:191510‑M1/128‑143 (MQ=255)
 gTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGa                                                                                                                                             <  2:182601/143‑1 (MQ=255)
 ctgcctgccTGGCGACCGCTCGGCACGCCCAGTTACTTACCCCACTTCGCGCGCTGCTGTTAAGCCGCAAGAACCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGa                                                                                                                                             <  2:415808/142‑1 (MQ=17)
       ttgtagatctgtgttctggtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCaaaa                                                                                                                                       >  1:606373‑M1/121‑143 (MQ=255)
                 gtgtgctggtcagttgttgacacgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgcag                                                                                                                             >  2:368687‑M1/111‑143 (MQ=255)
                 ggggCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttg                                                                                                                             <  1:448840‑M1/143‑18 (MQ=255)
                   ggcggcACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggt                                                                                                                           <  1:628768‑M1/143‑20 (MQ=255)
                         gtcatttgttgatacgcgtacatagccaataagcatgttaaatccccctggtaaaagcagcaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGACGGAGCAAAACTCCATGCAGCAACGGCt                                                                                                                     <  2:129995‑M1/41‑1 (MQ=255)
                                gTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGCCACTGTGCAAACGGCTCTCGAACACCCTGGCTTGGGAGCCGGcgccgcatgacctcttcgtgtcccggcgccgt                                                                                                              >  2:748493‑M1/1‑111 (MQ=255)
                                 ttccttgcCGGACTTCGCGGGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAgg                                                                                                             <  2:785853/140‑1 (MQ=255)
                                 ttacttacCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCCCTGTGCAACCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtg                                                                                                             >  1:132761‑M1/1‑110 (MQ=255)
                                  tacttacCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtgg                                                                                                            >  1:599305‑M1/1‑109 (MQ=255)
                                    cttacCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACg                                                                                                          >  1:333702/1‑143 (MQ=255)
                                          gCAGTTCGCGCGCTGGTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCgg                                                                                                    <  2:359715/143‑1 (MQ=255)
                                          gCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTATGACGCACCgg                                                                                                    >  2:287764/1‑143 (MQ=255)
                                              ttCGCTCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGCTATCAAATCGGATGGCACTGTGCAAACGGCACTCGCAACCCATGGTTTGGGAGAAGGcgctgcactacctcttggtgtcccggttccgtggcctgcagccgcc                                                                                                >  2:390204‑M1/1‑97 (MQ=255)
                                              ttCGCGCGCTGCTGTTAAGCCGGAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccact                                                                                                <  1:341058‑M1/143‑47 (MQ=255)
                                                  ccaataagcatgttaactccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAg                                                                                            >  2:20449‑M1/78‑143 (MQ=255)
                                                  ccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAg                                                                                            >  1:417690‑M1/78‑143 (MQ=255)
                                                  ccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGCTTCATAg                                                                                            <  2:23285‑M1/66‑1 (MQ=255)
                                                     gctgctGTTAAGCAGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCta                                                                                         <  2:238107/143‑1 (MQ=255)
                                                     ataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCta                                                                                         >  1:345929‑M1/75‑143 (MQ=255)
                                                       tgctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggagtcccggttccgtggccctcagccaccccgccaaca                                                                                       >  2:510626‑M1/1‑88 (MQ=255)
                                                                              gtcgccgatatCAAGTCGGATGCCACCGTGCGAGCGGCTCTCGAGACCCTTGGGTTGGGAGAAGGctccgcagtacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgc                                                                <  2:417673‑M1/138‑79 (MQ=255)
                                                                              aaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGa                                                                >  1:186020‑M1/50‑143 (MQ=255)
                                                                                   cGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgacatgcaatggtgccgctt                                                           >  1:367930‑M1/1‑60 (MQ=255)
                                                                                    gATTTCAAATCGGATGGTACTGTGCAAACGGCTCTCGAAACCCTTGGTTTGGGAGAAGTAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATc                                                          <  1:455597/143‑1 (MQ=255)
                                                                                    gATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAACACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGACTGGTGCTATTGCGCAATc                                                          >  2:95929/1‑143 (MQ=255)
                                                                                    gATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACGTTAGGACGCACCGGTTTCATAGCTATAACAGTAATGATTAGTGGTATTAAGCAATc                                                          >  1:201614/1‑143 (MQ=38)
                                                                                    gATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATc                                                          >  1:678482/1‑143 (MQ=255)
                                                                                          aaaTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcaataccagttggtgtcccggttccgtggccttcagccactccgccaccaggctggctgacatgcaatggtgccgctttttctgc                                                    >  2:313090‑M1/1‑53 (MQ=255)
                                                                                              catttgcctgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCAtt                                                >  1:235211‑M1/34‑143 (MQ=255)
                                                                                                tttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAAtcatgca                                              <  1:69606‑M1/112‑4 (MQ=255)
                                                                                                   ggCACTGTGCCAACGGCCCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGTACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTg                                           <  1:544948/143‑1 (MQ=255)
                                                                                                    gCAATGCGCAAACGGCACTCGAAAACCTTGGGTTGGGAGAAGGGGCAAAACTCAAGGCAGTAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTgg                                          <  2:241116/143‑1 (MQ=255)
                                                                                                      tgttattgcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTgggg                                        >  1:175981‑M1/26‑143 (MQ=255)
                                                                                                       gttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTCTTGAGCAATCAGCAATCATTCAGTgcggg                                       <  2:383092‑M1/119‑1 (MQ=255)
                                                                                                       cTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgcgaaccggctggctgaaatgcaatggtggcgctgtttctgctgaagaatacccg                                       >  2:539730‑M1/1‑40 (MQ=255)
                                                                                                          atttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGGAATCATTCagtgggggagg                                    <  2:35246‑M1/122‑4 (MQ=255)
                                                                                                              cttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATg                                <  2:443379‑M1/126‑1 (MQ=255)
                                                                                                                     ctcGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggc                         >  1:593432‑M1/1‑26 (MQ=255)
                                                                                                                      tcGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCaaaa                        >  1:173391/1‑143 (MQ=255)
                                                                                                                      ccgcaaaCCTGGGTTTGGGAGAGGGAGCAAAACGCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCaaaa                        <  1:269441/139‑1 (MQ=255)
                                                                                                                         aaaaCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCa                     >  1:533157/1‑143 (MQ=255)
                                                                                                                              cTTGTTCTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGc                <  1:618375/143‑1 (MQ=255)
                                                                                                                               ttGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacacggaactggccaaggcttatccgaca               >  1:183440‑M1/1‑16 (MQ=255)
                                                                                                                                  GTTGGGGAGAAGGAGCAAAACTCAGTGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCctct            >  1:262280‑M1/1‑143 (MQ=255)
                                                                                                                                    ttGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatggtgccgctttttctgctgaagaatacccggaactggcaaaggcttatccgacaaataa          >  2:91721‑M1/1‑11 (MQ=255)
                                                                                                                                         agcaGGTGCAAGACTCAATGCGGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTTTACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTggggggg     <  2:492691/140‑1 (MQ=255)
                                                                                                                                            AGGAGCACAACTCAATGCAGCAACGGCTACAGTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTAGTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATg  >  1:672111‑M1/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                            
GGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATG  >  NC_000913/1207663‑1207945

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