New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =1 (0.030)33 (1.130) 28/232 NT 95.8% coding (290/630 nt) ycfK e14 prophage; protein StfP
?NC_000913 1209619 = 2 (0.070)pseudogene (37/537 nt) stfE e14 prophage; putative side tail fiber protein fragment

ATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCA  >  NC_000913/1209482‑1209759
                                                                                                                                         |                                                                                                                                            
tgtctttgcggcgtgccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTa                                                                                                                                         <  2:207964‑M2/6‑1 (MQ=255)
        cggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATg                                                                                                                                 >  1:597948‑M2/130‑143 (MQ=255)
          gcattccaccccccctgaatgattactgctcgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgtggctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaa                                                                                                                               <  2:173499‑M2/16‑1 (MQ=255)
          gcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaa                                                                                                                               <  2:363743‑M2/16‑1 (MQ=255)
               ccacccccactgaatgattactgattgctcaataccaccaagcataaccggtatagctatgaaaacggtgcgtcctaatttaaccgttgctgcattaagtgttgctccctctccgaaaccgaCGTGTACGAAGATAAAgaaag                                                                                                                          >  2:533527‑M2/123‑142 (MQ=255)
                    cccactgaatgatcactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                                                                                                     <  1:358091‑M2/26‑1 (MQ=255)
                         tgaattattacgtattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCaaa                                                                                                                <  2:134747‑M2/31‑1 (MQ=255)
                          gaaggagtactgagtgctcagtaccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAAt                                                                                                               <  1:613068‑M2/32‑1 (MQ=255)
                              gattactgattgctcactaccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                                                                                                           <  2:246629‑M2/36‑1 (MQ=255)
                                ttactgattgctcaataccaccaatcataaccggtataactatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATgccacc                                                                                                         >  2:360346‑M2/106‑141 (MQ=255)
                                    tgattgctcaatcccaccaatcataaccggtatagcgatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc                                                                                                     <  2:44874‑M2/42‑1 (MQ=255)
                                     gattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatctagccgttgctgcaatgagttttgcgccttctcccaaaccaaCGCTTATGAAAATGAAGAAATAACACGCAAATGGCATcattgc                                                                                                    >  2:5230‑M2/101‑141 (MQ=255)
                                        tgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccagCGTTTATGACAATGAAGAAATAACAAGCAAATGGCCTCATTCCTGc                                                                                                 >  1:549105‑M2/98‑143 (MQ=255)
                                          ctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttctgccccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                                                                                               >  1:90555‑M2/96‑143 (MQ=255)
                                                  caccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgg                                                                                       >  2:398925‑M2/88‑143 (MQ=255)
                                                   accaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctgcttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttttacacggg                                                                                      >  2:285650‑M2/87‑138 (MQ=255)
                                                     cactcatgaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAGAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg                                                                                    <  2:653335‑M2/59‑1 (MQ=255)
                                                     caagcataacagctactgctatgaaactggtgcgccctaatgtagccgttgctgcattgagttttgctccttcactcaaaccaaCGTTTACGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg                                                                                    <  2:454993‑M2/59‑1 (MQ=255)
                                                          ataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgatgcattgagttttgctccttctcccaacCCAACGTTTCTGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTa                                                                               >  1:292304‑M2/76‑143 (MQ=255)
                                                             accggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACa                                                                            <  2:125452‑M2/67‑1 (MQ=255)
                                                                 gtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACAGGCt                                                                        >  1:174288‑M2/73‑143 (MQ=255)
                                                                         gtgagaccggtgcgtcctagcgtgcccttggctgcattgagttttgctccttctcccaaaccaaCGTTGATCAAAATGAAGACATAGCAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCt                                                                <  2:131621‑M2/79‑1 (MQ=255)
                                                                          gAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTa                                                               <  2:90699/143‑1 (MQ=255)
                                                                                 ggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccacacccaCTTCCATGAAAATGAAGAAATACCAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAcgc                                                        <  1:149843‑M2/87‑1 (MQ=255)
                                                                                        cctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaa                                                 <  2:297916‑M2/94‑1 (MQ=255)
                                                                                        cctaaggcagccgtcgctgcattgagttttgctcctgctcctaatccaaCGTTTATGAAAATGATGAAATCACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaa                                                 <  2:529501‑M2/94‑1 (MQ=255)
                                                                                         ctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaac                                                <  1:256303‑M2/95‑1 (MQ=255)
                                                                                           agtgtagccgttgctgtattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaacaa                                              <  2:90555‑M2/97‑1 (MQ=255)
                                                                                            atgtagccgttgctgcattgagttttgctccttctcccaacCCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTATCAGGCTTATTGGCTATGTACGCGTATCAACaaa                                             <  1:467092‑M2/102‑1 (MQ=255)
                                                                                            atgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACaaa                                             <  1:313958‑M2/98‑1 (MQ=255)
                                                                                             tgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCACATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAAt                                            <  1:4022‑M2/99‑1 (MQ=255)
                                                                                                     cacCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGaa                                    >  1:614031/1‑143 (MQ=255)
                                                                                                         tgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATAACCAGcacaca                                >  2:542352‑M2/33‑143 (MQ=255)
                                                                                                            attgagttatgcaccttctcgcaacCCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGAt                             <  2:597948‑M2/118‑1 (MQ=255)
                                                                                                             ttgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTCCGCGTATAAACAAATGACCAGAACACagacc                            >  1:54921‑M2/29‑141 (MQ=255)
                                                                                                             ttgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCGGAGCGCagacc                            >  1:67041‑M2/29‑141 (MQ=255)
                                                                                                             ttgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATc                            >  1:549610‑M2/29‑143 (MQ=255)
                                                                                                                   tttgcgccttctcccaaaccaaCGTTTATGAAAATGAACAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAAc                      <  1:463094‑M2/121‑1 (MQ=255)
                                                                                                                           ttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACCCAGATCTACAACGTAATccg              >  2:171456‑M2/15‑143 (MQ=255)
                                                                                                                              tcccaaaccaaCGTTTATGAAAATGAAGAAATAACGAGCAAATGGCATCATTCCTGCTTTTGCCGGGGGGATTCAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAgtgcgcgg           >  2:396446‑M2/12‑141 (MQ=255)
                                                                                                                                    ctcAACGTTTATGCAAAGCAAGAAAGAGCAAGCATGTGGCTTCATTCCTGCTTTTACCGGGGGGATTTAACATGCCTATTCGCTATGTACCCGTATCAACAAATGACCAGACCACAGATCTACAGCGTAATGCGCTGAACtgt     <  2:101613/141‑1 (MQ=255)
                                                                                                                                      caaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCCCTGAACTGTGc   >  2:330932‑M2/4‑143 (MQ=255)
                                                                                                                                       aaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCa  <  1:455875‑M2/141‑1 (MQ=255)
                                                                                                                                       aaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATATACGCGTATCCACAAATGACCAGAACACAGATCTACAACGTAATGCACTGAACAGTGCa  >  2:578652‑M2/3‑143 (MQ=255)
                                                                                                                                         |                                                                                                                                            
ATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCA  >  NC_000913/1209482‑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.