New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =46 (0.910)5 (0.110) 5/250 5.9 9.1% coding (290/630 nt) stfP e14 prophage, uncharacterized protein
?NC_000913 1209619 = 59 (1.310)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.
Rejected: Frequency below/above cutoff threshold.

GGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGC  >  NC_000913/1207656‑1207925
                                                                                                                                      |                                                                                                                                       
ggCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGCTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTgg                                                                                                                                 <  2:871431/143‑1 (MQ=255)
     ccGGTGCCTGCCTGGGGAAGGGTCGGCCCGCGCCGTTACTTATCGCACTGCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGaa                                                                                                                            <  1:424289/143‑1 (MQ=17)
      cgctGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAg                                                                                                                           <  2:709298/140‑1 (MQ=255)
        gggCCAACCTCCAGAAGGGGCCGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGa                                                                                                                         <  1:858014/141‑1 (MQ=255)
           ccAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATTTCAAATCGGATGGCACGGTGCAAACGGCTCTCGAAAACCTTGGTGTGGGAGGAGGAGCa                                                                                                                      >  2:554524/1‑143 (MQ=255)
            cAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTTGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCaa                                                                                                                     <  2:86513/143‑1 (MQ=255)
             aGCCTGGAGAAGGGGCGTCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAggctctgc                                                                                                                    <  2:783486/143‑7 (MQ=255)
             aGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCGGCTGTTAAGCCGCAAGAGTCCGTTTGGCGATATCAAAGCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCaaa                                                                                                                    <  1:592742/143‑1 (MQ=255)
                     gaAGGGGCGGCACCACCAGTTACTTACCGCACTTCGCGCGCTGCTCTTAAGCCGCAAGAACCCGTTTGGCGATATCAAATCGGCTGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAACCTCAATg                                                                                                            <  2:612421/143‑1 (MQ=18)
                      aaGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGGACGAGGAAATCTCAATgc                                                                                                           >  2:549594/1‑143 (MQ=255)
                       aGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATgca                                                                                                          <  2:63676/143‑1 (MQ=255)
                         gggCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTGAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTATGGGAGAAGGAGCAAAAGTCAATgcagc                                                                                                        >  1:429312/1‑143 (MQ=255)
                           gcggcACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCaa                                                                                                      >  2:775118/1‑143 (MQ=255)
                              gcACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTACATGCAGCCACgg                                                                                                   >  1:783900/1‑143 (MQ=255)
                              gcACGACCAGTTACTTACCGCACTTCGCGCGCTGCGGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGACAACCCTGGGTTGGGCGCAGCGGCAACACTCACTacagcaccac                                                                                                   >  1:621700/1‑139 (MQ=17)
                              gcACGACCAGTTACTTACCGCACTGCGGGCGCTGCTGTTAAGGCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACgg                                                                                                   <  1:379098/143‑1 (MQ=255)
                                  tttttttttcttacCGAACGGCGCGCGCAGCTGGTAAGCCGCAAGAAGCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                                                                               <  2:10668/134‑1 (MQ=21)
                                                  gaggtCGCGCGTGGCGGGTAAGCCGCAAGACTCCGTGTGGCGATATCAAATCGGATGGCACTGTGCTAACGGCTCTCGAAAACCTTGGTTTGAGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCTGt                                                                               <  1:850728/139‑1 (MQ=25)
                                                        gcgcggTGCTGGGAAGGCGCAAGAATCCGTTTGGCGAGATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATa                                                                         <  1:557289/143‑1 (MQ=25)
                                                          gcgcTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCACATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGCAGGAGCAAAACTCAATGCAGCAAAGGCGACATTAGGACGGAACGGGTTCATacc                                                                       >  1:579437/1‑141 (MQ=21)
                                                                  gTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCgg                                                               >  2:510985/1‑143 (MQ=255)
                                                                       gCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTAAATGCAGCAACGGCTAAATTAGGACGAACCGGTTTCATAGCTATACAGGTCATg                                                          >  2:293830/1‑143 (MQ=25)
                                                                                   cGTTTGGCGATATCCAAGCGGATGGCACGGTGCAAACGGCGCTCGAAAGCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTAtt                                              <  2:464816/143‑1 (MQ=255)
                                                                                    gtgtGGCGATATCAAAGCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAGTGCAGCAACGGCGACATTAGGACGCACCGGTTTCCTAGCTATACCGGTTATGATTGGTGGTATTg                                             <  1:188058/140‑1 (MQ=38)
                                                                                     tgtGGCGATATGAGATCGGATGGCACTGAGCAAACGGCTCTCGAAAACCTTGGTGTGGGAGAAGGAGCAAAACGCAATGCAGCAGCGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGa                                            <  2:179917/141‑1 (MQ=255)
                                                                                      tgggCGATATCAGATCGGCGGGCACGGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCCACACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAg                                           <  1:145816/141‑1 (MQ=38)
                                                                                        ggCGATATCGAAGCGGATGGCACTGGGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCa                                         <  2:22625/143‑1 (MQ=255)
                                                                                          ggATATCAACTGGGGTGGCACGGGGCGAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCGACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAAt                                       <  1:799090/142‑1 (MQ=38)
                                                                                                cAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTCTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTaa                                 <  2:136718/143‑1 (MQ=255)
                                                                                                   accgtttGGCACTGTGCAAACGGCTCTCGAAAACCTTGGGATAAAAGAAGGAGTAAAGCTAATTTTCTCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCa                              <  2:605145/137‑1 (MQ=255)
                                                                                                   aTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCa                              >  1:814503/1‑143 (MQ=255)
                                                                                                      ggATGGCAGTGTGCAAACGGCTCGCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTc                           <  2:579437/143‑1 (MQ=255)
                                                                                                       gATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGGAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTcc                          >  1:519185/1‑142 (MQ=255)
                                                                                                          ggCACTGTGCAAGCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTg                       <  1:358473/143‑1 (MQ=255)
                                                                                                            cACTGCGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCGATACCGGTTATGACTGGTGGTATTGAGCAATCAGTAATCACTCAGTggg                     >  1:738446/1‑143 (MQ=255)
                                                                                                             aCTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGCgggg                    >  1:427297/1‑143 (MQ=255)
                                                                                                                   cAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGGTTGGTGGTATTGTGCAATCAGTAATCATTCCGTGTGGGTGGaa              >  2:429488/1‑143 (MQ=255)
                                                                                                                       cagCTCTCGACAACCGTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGcc          <  1:375841/141‑1 (MQ=255)
                                                                                                                         cctctGGAGAACCTTCGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGc        <  1:114674/142‑1 (MQ=255)
                                                                                                                             tcGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCTGTGGGGGTGGAATGCCGCaaaa    >  2:522293/1‑143 (MQ=255)
                                                                                                                             tcGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCaaaa    >  2:525112/1‑143 (MQ=255)
                                                                                                                             tcGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTACTGAGCAATCAGTAATCATTCAGTGGGGGGGGAATGCCGCaaaa    >  1:593479/1‑143 (MQ=255)
                                                                                                                             gcGAACACCGTGGTCTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGCACGCACCGGTTTCATAGCGATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCaaaa    <  2:806665/142‑1 (MQ=255)
                                                                                                                               gAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGc  >  2:799924/1‑143 (MQ=255)
                                                                                                                               gAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGc  <  2:864423/143‑1 (MQ=255)
                                                                                                                               gAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGATATACCGGTTATGATTGGTGGTAGTGAGCAATCAATAATCACTCAGCGGGGGTGGAAGCCAGCAAAAGc  >  1:214593/1‑143 (MQ=255)
                                                                                                                                      |                                                                                                                                       
GGCATCCGGTGCCAGCCTGGAGAAGGGGCGGCACGACCAGTTACTTACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGC  >  NC_000913/1207656‑1207925

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.