New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =8 (0.150)40 (0.850) 31/250 0.5 79.1% coding (290/630 nt) stfP e14 prophage, uncharacterized protein
?NC_000913 1209619 = 14 (0.300)pseudogene (1/501 nt) stfE pseudogene, e14 prophage, side tail fiber protein fragment family,Phage or Prophage Related

GATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCA  >  NC_000913/1209481‑1209759
                                                                                                                                          |                                                                                                                                            
atggttttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGttt                                                                                                                                          <  1:846745‑M2/5‑1 (MQ=255)
      ttgcggcattccacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaa                                                                                                                                    <  1:638910‑M2/11‑1 (MQ=255)
       cTTTGCCGGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGaaaa                                                                                                                                   <  2:777740/143‑1 (MQ=255)
                ccacccccactgaatgatgactgcttgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAAt                                                                                                                          <  2:465820‑M2/21‑1 (MQ=255)
                 cacccccactgaatgattactgattgctcaataccaccaatcataaccggtatagatatgaacccagtgcgtcctaatgtagccgttgctgcattgagttttgctccctctctcaaaccaaCGTTTATGAAAAAGAAGAAATa                                                                                                                         >  1:495249‑M2/122‑143 (MQ=255)
                  acccccactgaatgattcctgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaacccaCGTCTCTGAACGTGAAGCAATaa                                                                                                                        >  2:86719‑M2/121‑143 (MQ=255)
                   cccccactgaatgattactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAAc                                                                                                                       <  2:570331‑M2/24‑1 (MQ=255)
                    gggTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAACTGAAGAAATAACa                                                                                                                      >  1:158784/1‑143 (MQ=255)
                    gggTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAAGGCAGAGCCTTCCCCCACACCAACGCTTACGAACAAGAACAAAAAACa                                                                                                                      >  1:161519/1‑143 (MQ=255)
                     cccactgactgattactgattgctgaatacccccactcataatcggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACaa                                                                                                                     <  2:105908‑M2/26‑1 (MQ=255)
                            ttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTCATGAAAATGAAGAAATAACAAGCAAATg                                                                                                              >  2:640297/1‑143 (MQ=255)
                                gcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcat                                                                                                          >  2:930426/1‑143 (MQ=255)
                                 ttactgattgctcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccatCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                                                                                                         <  2:108920‑M2/38‑1 (MQ=255)
                                           ctccatacccccattcctaaccggtacagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATTACAAGCAAATGGCATCATTCCTGCtt                                                                                               <  1:42735‑M2/48‑1 (MQ=255)
                                            tcaataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttt                                                                                              >  2:176688‑M2/95‑143 (MQ=255)
                                              aataccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttagtgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAACAAATGGCACAATTCCTGcttttg                                                                                            >  2:293774‑M2/93‑142 (MQ=255)
                                                taccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTccc                                                                                          >  1:613544‑M2/91‑143 (MQ=255)
                                                 accaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCa                                                                                         <  2:347760‑M2/54‑1 (MQ=255)
                                                    accaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                                                                                      <  2:392299‑M2/57‑1 (MQ=255)
                                                    accaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                                                                                      >  2:741280‑M2/87‑143 (MQ=255)
                                                    accaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                                                                                      >  1:852874‑M2/87‑143 (MQ=255)
                                                     ccaatcataaccggtatagctatgcaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGCAAATGAAGAAATAACAAGCAAACGGCATCATTCCTGCTTTTAccagcag                                                                                     >  2:501577‑M2/86‑140 (MQ=255)
                                                     cATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgggg                                                                                     <  1:711055/143‑1 (MQ=255)
                                                          cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt                                                                                <  1:514693/143‑1 (MQ=255)
                                                            taaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa                                                                              <  1:609793‑M2/65‑1 (MQ=255)
                                                             aaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctgcgtctcgcaaaccaaCGTTTATGAAAATGAAGACATAACAAGCAGATGGCATCATTCCTGCTTTTACCAGGGGGATTTCAc                                                                             >  1:586687‑M2/78‑143 (MQ=255)
                                                              cagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACaa                                                                                                                     <  1:350053/102‑1 (MQ=255)
                                                              cagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACaa                                                                                                                     >  2:350053/1‑102 (MQ=255)
                                                                   tatagctatgaaaccggtgcgtcctaatgtagccggtgctgcattgagttttgctccttcccccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCtt                                                                       <  2:649470‑M2/72‑1 (MQ=255)
                                                                   tGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGACATAACAAGCAAATGGCATCCTTCCTGCTTTTCCCAGGGCGCTTTAACCTCCtt                                                                       >  1:98561/1‑143 (MQ=255)
                                                                     taactactataccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttcccccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTCTTACCAGGGGGATTTAACATgcttct                                                                     >  2:378676‑M2/70‑141 (MQ=255)
                                                                       gctatgaaaccggtgcgtcctaatgtagccgttgctgcatggagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCGTATTg                                                                   >  1:503024‑M2/68‑143 (MQ=255)
                                                                        ctatgaaaccggtgcgtcctaatgtaggcgttgctgcattgagttttgctccttctcgcaaaccaaCCTTTATGAAAATGAAGAAATAACAAGCAAATGGCGTCATTCCTGCTGTTACCAGGCGGATTGACCAAGCGTATggg                                                                  >  2:58661‑M2/67‑143 (MQ=255)
                                                                         tatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGc                                                                 <  2:431210‑M2/78‑1 (MQ=255)
                                                                          atgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCt                                                                <  1:6662‑M2/79‑1 (MQ=255)
                                                                            ggaggcggtgcgtgctaatgtagccgttgctgcagtgagttttgcgccttctctcagaccatCGGTTAGGAAAATGAAGAAATAACAAGCAAATGGCATCATTGCTGCTTTTACTAGGGGGATTTAACATGCTTATTGGCTAt                                                              <  2:533809‑M2/81‑1 (MQ=255)
                                                                                  tGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTCTACCAGGGGGATTTAAAATGCTTATTGGCTATGTAcgc                                                        <  2:158784/143‑1 (MQ=255)
                                                                                  ggttcgtcctcatgtagcagctgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAcgc                                                        <  2:763684‑M2/87‑1 (MQ=255)
                                                                                  gggggagcgtaaggtagccgttgctgcatcgagttttgctccttctcccaaaccaaCGTGTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTGCTGCGTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAcgc                                                        <  1:851792‑M2/87‑1 (MQ=255)
                                                                                  ggggcgtcctaatgtagccgttgctgcattgagttttgctccttgtcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAcgc                                                        <  2:741777‑M2/87‑1 (MQ=255)
                                                                                  agtcccttctcatgtagccgttgtttcattgagttttgctccttcgcccaaaccatTGTTTATGAAAATTAAGTAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGAAcgc                                                        <  2:109980‑M2/87‑1 (MQ=255)
                                                                                     gcggcctaatgtagccgttgctgccttgagttttgctccttctccccaaccatCGTTTATGAAAACGAAGAAATACCAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTa                                                     <  1:531480‑M2/90‑1 (MQ=255)
                                                                                        tcctagtgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATca                                                  >  1:450424‑M2/51‑143 (MQ=255)
                                                                                                 agccgttgctgcattgagttttgctcgttctcccagttgagCTTTCATGAGAATGAAGAAAGGACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGAc                                         <  2:905387‑M2/102‑1 (MQ=255)
                                                                                                 agccgttgctgcattgagttttgctccgtctcccaaaccaaACTCTCTTGATATTAAGAACTAAAAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGAc                                         <  1:635474‑M2/102‑1 (MQ=255)
                                                                                                   acgatgtcacaaggagtttttctccttctcccaatccatCGTTTATGAAAATGAAGAAATTACAATCTAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCa                                       <  2:548958‑M2/104‑1 (MQ=255)
                                                                                                      cacCAACAGGTAAGGCAGAGCCTTCTCCCAAACCAACGTTTAGGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGaa                                    <  1:72816/143‑1 (MQ=255)
                                                                                                       tgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTgg                                                                  >  1:184539‑M2/36‑112 (MQ=255)
                                                                                                       tgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTgg                                                                  <  2:184539‑M2/77‑1 (MQ=255)
                                                                                                         cAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGCacac                                 >  1:910734/1‑143 (MQ=255)
                                                                                                             atttagttttgcgccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGAt                             <  2:10561‑M2/114‑1 (MQ=255)
                                                                                                              ttgagttttgctccttctcccaaaccaaCGTTTATGAACATGAAGAAATAACCCGCACATGGCCTCCTTCCTGCTTTTACCAGGGGGCTTTAACATGCTTATTGGCTATGTACGCGTATCAACACATAACCAGCACAGAAATc                            >  1:139907‑M2/29‑143 (MQ=255)
                                                                                                               gTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCt                           >  1:393219/1‑143 (MQ=255)
                                                                                                                         tccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAAt                 >  1:79733‑M2/18‑143 (MQ=255)
                                                                                                                               tcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATAAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTg           >  2:728244‑M2/12‑143 (MQ=255)
                                                                                                                               tcCCAAACCAACGTTAATGAACATGAAGAAATAACAAGCAAATAGCATCATTCCTACTTTTACCAGGGGGATTTAACATGCTTACTGGCTATGTACGCGTATAAAGAAATGACAAGAACACCGACCTCCAACATAATGCGCta           >  2:147508/1‑142 (MQ=255)
                                                                                                                                 ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGaa         <  1:741280‑M2/134‑1 (MQ=255)
                                                                                                                                 ccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTACTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCACAACACACATATACAACGCCATGCActgcc         >  2:746606‑M2/10‑141 (MQ=255)
                                                                                                                                      ccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTGTTGAAAGGGGGAGTTAATATGTGCATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGGACtgtg    >  2:111539‑M2/5‑143 (MQ=255)
                                                                                                                                        aaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTGATTGGCTATGTACGCGCATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCa  <  1:930426‑M2/141‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                            
GATAAGCCTTTGCCAGTTCCGGGTATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCA  >  NC_000913/1209481‑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.