New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =17 (0.380)42 (1.060) 38/252 0.1 78.5% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 1209619 = 8 (0.200)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related

GTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAG  >  NC_000913/1209477‑1209760
                                                                                                                                              |                                                                                                                                             
gcagatgcttttgcggcattccaaccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaac                                                                                                                                               <  1:392438‑M2/1‑1 (MQ=255)
    atgcttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgcaaccggtgcgtcctaatgtagccgttgctgcattgagtcttgccccttctcccaacccaaCGttt                                                                                                                                           >  1:202519‑M2/139‑143 (MQ=255)
     tgcttttgcggcattccacccccactgaatgattactgattgctcaataccacccatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccggtgctgcattgagttttgctccttctcccaaaccaaCGTTTa                                                                                                                                          <  2:101772‑M2/6‑1 (MQ=255)
          ttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaa                                                                                                                                     <  1:592116‑M2/11‑1 (MQ=255)
             ttGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATg                                                                                                                                  >  1:482249/1‑143 (MQ=255)
              tGCCAGTGCCGGGGATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCACCGTTTATGAAAATga                                                                                                                                 <  2:289080/143‑1 (MQ=255)
              ggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccagtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccacagcaaCGTCtatgaaactaa                                                                                                                                 >  2:503582‑M2/129‑139 (MQ=255)
               gcattccccccccactgaatgattactgatcgctcaacaccaccaatcataaccggtatagctatgaacccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaa                                                                                                                                <  1:252968‑M2/16‑1 (MQ=255)
               gcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaacccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaa                                                                                                                                <  1:530296‑M2/16‑1 (MQ=255)
               gCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCCTTTCAGCCAGCCTGTGGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaa                                                                                                                                <  2:393988/143‑1 (MQ=255)
                     cacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                                                                                                          >  2:145333‑M2/122‑143 (MQ=255)
                     cacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATa                                                                                                                          <  1:573092‑M2/22‑1 (MQ=255)
                        ccccactgaatgattactgattgctcaataccaccactcataaccggtatagctatgagaccgctgcatcctaatgtagccattgctgcattgagttttgctccttctcccaaacccaCGTCTAGGAAAATGAAAAAATcaca                                                                                                                       >  2:630563‑M2/119‑143 (MQ=255)
                                     ttactggttgctcaagaccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                                                                                                          <  1:792102‑M2/38‑1 (MQ=255)
                                     ttactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGATTATGAAAATGAAGAAATAACACGCAAATGGcatc                                                                                                          >  1:105325‑M2/106‑143 (MQ=255)
                                          gattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcagtgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                                                                                                     >  2:106727‑M2/101‑143 (MQ=255)
                                              gctcagtaccgccaatcataaccggtatagctatgaaaccggtgcgccctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCt                                                                                                 <  2:525166‑M2/47‑1 (MQ=255)
                                               ctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                                                                                                >  2:596171‑M2/96‑143 (MQ=255)
                                                tcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCATATGGCATCCATCCTGCttt                                                                                               >  1:712838‑M2/95‑143 (MQ=255)
                                                tcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttt                                                                                               >  2:159532‑M2/95‑143 (MQ=255)
                                                tcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAACCAACAGCCGGCATTCCTGCttt                                                                                               >  1:294024‑M2/95‑143 (MQ=255)
                                                  aataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctctcaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa                                                                                             <  1:156851‑M2/51‑1 (MQ=255)
                                                       caccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgg                                                                                        >  1:465441‑M2/88‑143 (MQ=255)
                                                        accaatcataaccggtatagctatgaaacaggtacgtcctaatgtagccgttgctgcattgagttttgctacctctcccaaaccaaCGTTGATGAAAACGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                                                                                       >  2:15768‑M2/87‑143 (MQ=255)
                                                         ccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgggg                                                                                      <  1:625708‑M2/58‑1 (MQ=255)
                                                                taaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa                                                                               <  1:170498‑M2/65‑1 (MQ=255)
                                                                 aaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAAc                                                                              >  2:66465‑M2/78‑143 (MQ=255)
                                                                   ccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACAt                                                                            <  1:186844‑M2/68‑1 (MQ=255)
                                                                         tagttatgaaaccggtgcgtcctaatgtagccgttgctgtattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAt                                                                      <  2:492983‑M2/74‑1 (MQ=255)
                                                                                gaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTAt                                                               >  1:660501‑M2/63‑143 (MQ=255)
                                                                                  aaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGt                                                             >  2:330648‑M2/61‑143 (MQ=255)
                                                                                   accggtgcttcctaccgtagccgttgctgcattcagttttgctccttctcccaaaccaaCGTTTATGAGTATGAAGCAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTa                                                            <  2:132025‑M2/84‑1 (MQ=255)
                                                                                      ggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAcgc                                                         <  1:596171‑M2/87‑1 (MQ=255)
                                                                                         gcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTa                                                      >  1:242191‑M2/54‑143 (MQ=255)
                                                                                            tcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTCTGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCATTTACCAGGGGGATTTAACATGCTTATTGGTTATGTTCGCATATca                                                   >  2:592051‑M2/51‑143 (MQ=255)
                                                                                            tcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTACCTACGCGTATca                                                   >  2:106725‑M2/51‑143 (MQ=255)
                                                                                             cctactgtagccgttgctgcattgatttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaa                                                  <  2:35657‑M2/94‑1 (MQ=255)
                                                                                              aCGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATccac                                                 >  2:265763/1‑143 (MQ=255)
                                                                                                 gAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACaaa                                              >  2:409367/1‑143 (MQ=255)
                                                                                                 atgtagccgttgctgcatggagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACaaa                                              <  2:24798‑M2/98‑1 (MQ=255)
                                                                                                   gtagccgttgctgcattgatttttgctccttctcccaaaccaaCTTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATg                                            <  2:204995‑M2/100‑1 (MQ=255)
                                                                                                         gttgcagcattgagttctgctccttctcccaaaccagCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATGCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGa                                      <  1:9873‑M2/106‑1 (MQ=255)
                                                                                                           acCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAac                                    >  1:601344/1‑143 (MQ=255)
                                                                                                             cAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAAGAAATGACCAGCacac                                  >  1:531444/1‑143 (MQ=255)
                                                                                                             cAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAacac                                  >  1:180549/1‑143 (MQ=255)
                                                                                                              tgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCCGAAcaca                                 >  2:142342‑M2/33‑143 (MQ=255)
                                                                                                                  ttgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACCACCAAATGGCATCATTCCTCCTTTTACCAGGGGGCTCTACCATGCTTATTGCCTCTGTACGCGTATCACCAAATGCCCAGACCACAGATc                             >  2:334977‑M2/29‑143 (MQ=255)
                                                                                                                   tgattcttgctcctcctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAACTGACCAGAACACAGATCt                            <  1:11001‑M2/116‑1 (MQ=255)
                                                                                                                     agttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTAc                          <  1:693217‑M2/118‑1 (MQ=255)
                                                                                                                      gttttggtgcttctcccaaaccaaCGTTTATGAATATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACGAATGACCAGAACACAGATCTACa                         <  2:350064‑M2/119‑1 (MQ=255)
                                                                                                                              ccttctcccaagccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATg                 <  2:40335‑M2/127‑1 (MQ=255)
                                                                                                                                 tctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATgcgc              <  2:437724‑M2/130‑1 (MQ=255)
                                                                                                                                  ctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCt             <  1:68812‑M2/131‑1 (MQ=255)
                                                                                                                                     ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTgac          >  1:703541‑M2/10‑142 (MQ=255)
                                                                                                                                     ccAAACCAACGTTTATGAAAATGAAGAACTAACAAGCAAACGGCATCATTACTGGTAGCACCAGGGGGATTTAAAATGCTTATTGGCTATTTACGCGTATGAACAAATGACCAGAACACAGATCTAACAGCTAATGCGCTGaa          >  2:355851/1‑143 (MQ=255)
                                                                                                                                           caaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACCGTGc    >  1:441025‑M2/4‑143 (MQ=255)
                                                                                                                                            aaCGTTTATGAAAATGAAGGAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCa   <  2:506568‑M2/141‑1 (MQ=255)
                                                                                                                                             aCGTTTATGAAAATGAAGACATTCTCCACGATTTACATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAg  <  1:716811‑M2/142‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                             
GTCGGATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAG  >  NC_000913/1209477‑1209760

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.