New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =68 (0.970)15 (0.260) 15/166 4.7 21.9% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 1209619 = 50 (0.850)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.

GCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATG  >  NC_000913/1209522‑1209718
                                                                                                 |                                                                                                   
gCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGt                                                                                                   >  2:724250/1‑100 (MQ=255)
 cGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAAGGt                                                                                                   >  1:1211275/1‑99 (MQ=255)
 cGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGtt                                                                                                  <  1:208910/100‑1 (MQ=255)
   gCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAAGGtttt                                                                                                >  2:1698316/1‑99 (MQ=255)
    cACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAAGGttttc                                                                                               <  2:1283928/100‑3 (MQ=255)
    cACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTa                                                                                                >  1:395371/1‑99 (MQ=255)
    cACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTAt                                                                                               <  2:59983/100‑1 (MQ=255)
     aCCATTGCAGTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTAt                                                                                               <  1:1219643/99‑1 (MQ=255)
        aTTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGaaa                                                                                           >  2:1555412/1‑100 (MQ=255)
        aTTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGaaa                                                                                           <  1:1555405/100‑1 (MQ=255)
            cATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAAGGttttcgagagccg                                                                                       <  2:1019297/100‑11 (MQ=38)
              tttCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaag                                                                                     <  2:195395/100‑1 (MQ=255)
               ttCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaaga                                                                                    >  2:185621/1‑100 (MQ=255)
                    ccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACAGCAAAAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATaa                                                                               >  1:867772/1‑100 (MQ=255)
                       gccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACaa                                                                            >  2:615984/1‑100 (MQ=255)
                        ccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAg                                                                           <  2:64819/100‑1 (MQ=255)
                        ccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAg                                                                           <  2:1292891/100‑1 (MQ=255)
                         cTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGc                                                                          >  2:1043563/1‑100 (MQ=255)
                          tGTTGGCGGAGTGGCTGAAAGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCa                                                                         >  1:1064667/1‑100 (MQ=255)
                           gTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCaa                                                                        >  1:104577/1‑100 (MQ=255)
                                 ggAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGca                                                                  >  1:536697/1‑100 (MQ=255)
                                  gAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcat                                                                 >  2:764321/1‑100 (MQ=255)
                                   aGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGc                                                                   <  1:1696787/97‑1 (MQ=255)
                                      ggCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAAGCAACGTTTATGAGAATGAAGAAATAACAAGCAAATGGcatcat                                                              >  2:1352719/1‑99 (MQ=255)
                                      ggCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                                                             >  1:38400/1‑100 (MQ=255)
                                      ggCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                                                             >  2:1359213/1‑100 (MQ=255)
                                      ggCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                                                             >  1:710449/1‑100 (MQ=255)
                                       gCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc                                                            >  1:623689/1‑100 (MQ=255)
                                        cTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCACCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                                                           >  1:869005/1‑100 (MQ=255)
                                        cTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                                                           <  1:1662309/100‑1 (MQ=255)
                                        cTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                                                           >  2:1030577/1‑100 (MQ=255)
                                        cTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATAATTc                                                            >  1:1209169/1‑99 (MQ=255)
                                               ccACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtttt                                                    <  1:170521/100‑1 (MQ=255)
                                                cACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTa                                                   <  2:215103/100‑1 (MQ=255)
                                                   ggAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCa                                                >  1:762420/1‑100 (MQ=255)
                                                        cGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggggg                                           >  2:734671/1‑100 (MQ=255)
                                                         gggACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGATATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa                                          <  1:932763/100‑1 (MQ=255)
                                                         gggACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa                                          <  1:1091765/100‑1 (MQ=255)
                                                          ggACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAt                                         <  1:1359206/100‑1 (MQ=255)
                                                          ggACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAt                                         >  1:1453879/1‑100 (MQ=255)
                                                              acCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTaa                                     <  1:1533815/100‑1 (MQ=255)
                                                               ccAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAAc                                    >  2:1167937/1‑100 (MQ=255)
                                                                    aGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCt                               <  1:1122845/100‑1 (MQ=255)
                                                                     ggTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCtc                              >  2:1037053/1‑99 (MQ=255)
                                                                     ggTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATAGCATCATTCCTGCTTTTACCAGGGGGATTTAACCTGCtt                              >  1:429838/1‑100 (MQ=255)
                                                                       gagttttgctccTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAt                            >  1:1340056/11‑100 (MQ=255)
                                                                        aaTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAtt                           >  1:1083224/1‑100 (MQ=255)
                                                                         aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTCATTg                          <  1:1037049/100‑1 (MQ=255)
                                                                           gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGc                        <  1:295119/100‑1 (MQ=255)
                                                                              gagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTAt                     <  1:823981/100‑1 (MQ=255)
                                                                              gagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAGTTAACATGCTTATTGGCTAt                     <  1:1043559/100‑1 (MQ=255)
                                                                                   ttCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAcg                <  2:516835/100‑1 (MQ=255)
                                                                                   ttCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAcg                >  2:162054/1‑100 (MQ=255)
                                                                                     ctcCCAAACCAACGTTTATGAAAATGAAGAAATAACAATCAAATGGCATCATTCCTGCCTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAcgcg              <  1:928643/100‑1 (MQ=255)
                                                                                     ctcCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAcgcg              <  1:1064981/100‑1 (MQ=255)
                                                                                      tcCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGt             >  2:1630841/1‑100 (MQ=255)
                                                                                       cccAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGt             >  2:150143/1‑99 (MQ=255)
                                                                                       cccAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTa            >  1:451705/1‑100 (MQ=255)
                                                                                         cAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATc          <  1:969051/100‑1 (MQ=255)
                                                                                               aaCGTTTATGAAAATGAAGAAATAACAAGCAAACGGCATCATTCCTGCTTTTACCAGGGGGAATTAACATGCTTATTGGCTATGTACGCGTATCAACaaa    >  1:1122522/1‑100 (MQ=255)
                                                                                                 cGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACAAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAGCAAATg  >  1:713810/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
GCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATG  >  NC_000913/1209522‑1209718

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ 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.