New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =24 (0.190)46 (0.380) 41/548 3.6 71.0% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 1209619 = 15 (0.120)pseudogene (1/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related
Rejected: Coverage evenness skew score above cutoff.

AGTAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCA‑AAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTTTATTACGGTGCGC  >  NC_000913/1207493‑1208081
                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                    
aGTAATCCAAGACACGTAAGGACACGAAAACCGGGCAGAAGGAGAAGTTCGGCCCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGGCATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTgg                                                                                                                                                                                                                                                                                                   <  1:983125/301‑1 (MQ=255)
   aatccatgagtaggataaccccGAAAACCTCGCAGAAGGATAGGTTCGGTGCGGGTAACACCGGTTTTACCCGTGGTAACCTCCAGACCGAGACGACTGCCACCGGTCTGGATGATGGCTACTTTGGCATGTTGCAGGAGGAACTTTGCAGCGTGGTTGAGGCATCCTGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAGGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGttt                                                                                                                                                                                                                                                                                                <  1:356971/280‑1 (MQ=255)
      gcaagccgcggataGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGGCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTgg                                                                                                                                                                                                                                                                                              <  1:110929/289‑1 (MQ=255)
      ccAGGCAGCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTggg                                                                                                                                                                                                                                                                                             <  2:132270/301‑1 (MQ=255)
           gagagtaacgacaacaaaaCCGCGAAGAAGGATAAGTTCGGCGCGGGGCAGAACGGTTTTACCCGTGGTAACCCCCAAACCGGCACGCCTTCCACCGATCTGGATGATAACTACTGTGACATGTTGCAGGGGTAACTTTGCAGCGTTGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGAACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAg                                                                                                                                                                                                                                                                                        <  1:387900/286‑1 (MQ=17)
           cATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCCCTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTGTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCGCTCGAAAACCTTGGTTTGGgaca                                                                                                                                                                                                                                                                                          >  2:273360/1‑297 (MQ=255)
                tataGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCCAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGtagatcggagagaac                                                                                                                                                                                                                                                                                    >  2:959123/1‑286 (MQ=255)
                tataGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGt                                                                                                                                                                                                                                                                                                  <  1:958871/286‑1 (MQ=255)
                 ccggtctcGAAAACCGGGCAGAAGGATAAGTGGGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCa                                                                                                                                                                                                                                                                                   <  1:235927/293‑1 (MQ=255)
                      aacGAAAACCGGGCAGAAGGAAAAGTGGGCGGCGGGTAAGAAGGGGTTTACCCGTGGTAACCCCCAGACCGGCACGCATGCCACCGAGCTGGATGATGACTACTTTGAAATGTTGCAGGAGGAACTTTGGAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAACCCGTATGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGAA‑AAACTc                                                                                                                                                                                                                                                                            <  2:862981/300‑1 (MQ=255)
                       acGAAAACCGCGCAGAATGATAAGTTCGGCGCGGGTAAGAACGGTGTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAggctctgcattac                                                                                                                                                                                                                                                                            <  2:67097/301‑12 (MQ=255)
                        cGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCCCGAAAACCTTGGTTTGGGAGAAGGAGCa                                                                                                                                                                                                                                                                                   >  2:235975/1‑293 (MQ=255)
                         gAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGACggctcttcattacct                                                                                                                                                                                                                                                                          >  2:789204/1‑288 (MQ=255)
                                gcgcAGAAGGATGAGTTCGGCGCGGGTAAGAACGGTGTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTAGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAggctctgcattacctgttggtg                                                                                                                                                                                                                                                                   <  2:29935/301‑21 (MQ=255)
                                  gcAGAAGGATAAGTTCGGCCAGGGTAAAAACGGGTTTGCCCGGGGTCAACCCCCGACCGGCACGCCTGCCACCGAGCTGGATGATGACTACTTGGACAGGTTGCAGGAGGAACGTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTGACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAggctctgcattacctgttggtgtc                                                                                                                                                                                                                                                                 <  1:116011/301‑23 (MQ=17)
                                    aGAAGGCTAAGTTCGGGACGGGTAAGAGCGGTTGTGCCCGGGGTAACCCCCACACCGGCACGCCTGCCACCGATGTGGATGATGACTACTTTGAAAGGTGGCAGGAGGAACGTTGCAGCGGGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCGGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAggctctgcattacctgttggtgtccc                                                                                                                                                                                                                                                               <  1:273304/301‑25 (MQ=17)
                                                  ggCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCA‑AAACTCAATGCAGCAACGGCTACATTAGGAc                                                                                                                                                                                                                                                   <  2:673353/298‑1 (MQ=255)
                                                  ggCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCCCTCGAAAACCTTGGTTTGGGAGAAGGAGCA‑AAACTCAATGCAGCAAAGGCTACATTAGGaa                                                                                                                                                                                                                                                   >  1:673160/1‑297 (MQ=255)
                                                                                     ccAGACCCACACTCCGGCCACCCATCGGGATGATGACTACTGTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCA‑AAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTggtggt                                                                                                                                                                                                             <  2:720499/301‑1 (MQ=255)
                                                                                                                              ttcaCATGTTGCAGCAGGAACTTTGCAGCGGGGTGGAGGCCGCCGGTGCCAGCCTGGAGAAGGGTCTGCACGACCAGTTACTTACCGCACTTCGCGCTCTGCTGTTAAGCCGAAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAGACCTTGGTTTGGGAGAAGGAGCA‑AAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCg                                                                                                                                                                     <  2:1142114/297‑1 (MQ=255)
                                                                                                                                                                     cATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCA‑AAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATTGGCCTCTGGGGGGGATGGAAATCCAGgtgt                                                                                                                             >  1:485295/1‑301 (MQ=255)
                                                                                                                                                                               cAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGCAGGAGCA‑AAACTCAATGCAACAAAGGCTACATTAGGACGCACCGGTATCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGGAATCATTCAGTGGGGGGGGAATGCCGCGAAAGTATCGGCCCCTGGGGGGGAGGGAAAGACAGTTGTGTTCCCggcg                                                                                                                   >  1:599094/1‑299 (MQ=37)
                                                                                                                                                                                                                                                            gggATAACAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCA‑AAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACggg                                      <  2:485425/299‑1 (MQ=255)
                                                                                                                                                                                                                                                                                    atgtcttcattttcataaaCGTTGGTTTGGGAGAAGGAGCA‑AAACTCAATGCAGCAACGGCGACAGTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACa                                                                        <  2:291475/227‑1 (MQ=255)
                                                                                                                                                                                                                                                                                          tcattttcataaaCGTTGGTTTGGGAGAAGGAGCA‑AAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCGGTGTTTCCGTGGTTGCAAATTATGAACATGTCAGCGCACCTATCATTGCAGTGGCAACGGGGGGATATAGAACCAACCAGTTGTTATTAcc        >  1:42341/11‑300 (MQ=255)
                                                                                                                                                                                                                                                                                               tgcataaaGCTTGGTTTGGG‑GAAAGAGGACAAATTCACAGGAGGACCGGATACATGAGTACGCACGGGTTTCAAAGCTTTACCAGGTATGAGTGGTGTGATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCAGCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTTTATTACGGTgc    <  1:1057899/295‑1 (MQ=255)
                                                                                                                                                                                                                                                                                               ctcGAAAACCTTGGTTTGGGAGAAGGAGCA‑AAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCAc                                                            >  2:212084/1‑244 (MQ=255)
                                                                                                                                                                                                                                                                                                  ataaaCGTTGGTTTGGGAGAAGGAGCA‑AAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTAGTTTTTATTACGGTacgc  >  2:19996/3‑299 (MQ=255)
                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                                    
AGTAATCCATGCATCGTATAGACACGAAAACCGCGCAGAAGGATAAGTTCGGCGCGGGTAAGAACGGTTTTACCCGTGGTAACCCCCAGACCGGCACGCCTGCCACCGATCTGGATGATGACTACTTTGACATGTTGCAGGAGGAACTTTGCAGCGTGGTGGAGGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCA‑AAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGTGGCAACGGGGGGATATACAACCACTTCGTTTTTATTACGGTGCGC  >  NC_000913/1207493‑1208081

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