New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1207790 =43 (0.660)17 (0.310) 16/166 4.0 32.1% coding (290/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 1209619 = 36 (0.660)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.

AGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAAT  >  NC_000913/1209521‑1209717
                                                                                                  |                                                                                                  
aGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACg                                                                                                   <  2:335151/100‑1 (MQ=255)
    gCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTa                                                                                               >  2:830025/1‑100 (MQ=255)
      aCCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAAGGttttcg                                                                                             <  1:1454875/100‑4 (MQ=255)
        cATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAAGGttttcgag                                                                                           <  1:214846/100‑6 (MQ=255)
           tGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAAGGttttcgagagc                                                                                        >  1:429892/1‑92 (MQ=39)
                ttCAGCCGGCCTGTTGGCGGTGTGGCTGAAGGCAACGGAACCGGGACACCAACAGGTAATGCTGTGCCTTCTCCCAATCCAACGTTTATGAAAATgaaga                                                                                   <  2:435004/100‑1 (MQ=255)
                ttCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACTGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaaga                                                                                   <  2:1124172/100‑1 (MQ=255)
                ttCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaaga                                                                                   >  2:525909/1‑100 (MQ=255)
                  cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGaaa                                                                                 <  1:1144439/100‑1 (MQ=255)
                  cagccagcCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGaaa                                                                                 >  2:310339/1‑100 (MQ=255)
                   agccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAAt                                                                                <  2:974283/100‑1 (MQ=255)
                   agccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAAt                                                                                <  2:38330/100‑1 (MQ=255)
                    gccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGCTTATGAAAATGAAGAAATa                                                                               >  1:595941/1‑100 (MQ=255)
                         ccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAg                                                                          <  2:1110135/100‑1 (MQ=255)
                                       ggCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                                                            >  2:368420/1‑100 (MQ=255)
                                        gCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc                                                           >  1:1256375/1‑100 (MQ=255)
                                         cTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                                                          <  1:1062984/100‑1 (MQ=255)
                                           gAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg                                                        <  2:669403/100‑1 (MQ=255)
                                             aGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGAATCATTCCTGCt                                                      >  2:498470/1‑100 (MQ=255)
                                              ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                                                     >  1:39035/1‑100 (MQ=255)
                                              ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                                                     <  1:443458/100‑1 (MQ=255)
                                               gCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttt                                                    >  1:526636/1‑100 (MQ=255)
                                                  aCGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAc                                                 <  2:19580/100‑1 (MQ=255)
                                                   cGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAATAAATAACAAGCAAATGGCATCATTCCTGCTTTTAc                                                 <  1:578644/99‑1 (MQ=255)
                                                        ccGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgggg                                           >  1:1476647/1‑100 (MQ=255)
                                                             acacCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAATAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt                                      <  2:505713/100‑1 (MQ=255)
                                                             acacCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt                                      <  1:310336/100‑1 (MQ=255)
                                                             acacCAACAGGTAATGCAGAGCCTTCTCCCAAACAAACGTTTATGAACATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAttt                                      <  1:432081/100‑1 (MQ=255)
                                                                ccAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAAc                                   <  2:1014308/100‑1 (MQ=255)
                                                                      ggTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCt                              >  2:1467712/1‑99 (MQ=255)
                                                                           tGCAGAGCCTTCTCCCAACCCAACGTTTATGAAAATGAAGACATAACAAGCAAATCGCATCAATCCTGCTATTACAAGGCGGATTTAACATGCTTAtg                          >  1:525897/1‑97 (MQ=255)
                                                                            gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAACAAATGGCATCAATCCTGCTTTTACCAAGGGGAACTAACATGCTTATTGGc                       >  1:936353/1‑100 (MQ=255)
                                                                            gCAGAGCCTTCCCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCGGCTTTTACCAGGGGGATTTAACATGCTTATTGGc                       >  1:842786/1‑100 (MQ=255)
                                                                              agagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTa                     <  1:1191846/100‑1 (MQ=255)
                                                                              agagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTa                     <  2:542497/100‑1 (MQ=255)
                                                                                 gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGt                  >  1:257517/1‑100 (MQ=255)
                                                                                  ccTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCATGCTATTACCAGGGGGATTTAACATGCTTATTGGCTATGTa                 >  1:1339823/1‑100 (MQ=255)
                                                                                   cTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAc                >  1:180351/1‑100 (MQ=255)
                                                                                    ttCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAcg               >  2:637713/1‑100 (MQ=255)
                                                                                     ttacccAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAcg               <  1:1232991/96‑1 (MQ=255)
                                                                                      ctcCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAcgcg             <  2:1443251/100‑1 (MQ=255)
                                                                                      ctcCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAcgcg             <  2:959948/100‑1 (MQ=255)
                                                                                       tcCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGt            >  2:242747/1‑100 (MQ=255)
                                                                                        cccAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTa           >  2:420084/1‑100 (MQ=255)
                                                                                        cccAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTa           >  2:1550685/1‑100 (MQ=255)
                                                                                        cccAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTa           >  1:11605/1‑100 (MQ=255)
                                                                                          cAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATc         >  1:1412998/1‑100 (MQ=255)
                                                                                          cAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATc         >  1:1335269/1‑100 (MQ=255)
                                                                                          cAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATc         <  1:352636/100‑1 (MQ=255)
                                                                                          cAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATc         >  1:923385/1‑100 (MQ=255)
                                                                                          cAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATc         <  2:949079/100‑1 (MQ=255)
                                                                                          cAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTCAACATGCTTATTGGCTATGTACGCGTATc         >  1:963646/1‑100 (MQ=255)
                                                                                             aCCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGTTTATTGGCTATGTACGCGTATcaac      <  2:390607/100‑1 (MQ=255)
                                                                                              ccAACGTTTATGAACATGCAGAAATAACAAGCAAATGGCATCATTCCTGCTATTACGAGGGGGATTTAACATGCTTATTGCCTATGTACGCGTATcaaca     >  1:107737/1‑100 (MQ=255)
                                                                                              ccAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaaca     <  2:87589/100‑1 (MQ=255)
                                                                                              ccAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGAGGGATTTAACATGCTTATTGGCTATGTACGCGTATcaaca     >  1:12008/1‑100 (MQ=255)
                                                                                                 aCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAAt  <  2:768557/100‑1 (MQ=255)
                                                                                                  |                                                                                                  
AGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAAT  >  NC_000913/1209521‑1209717

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.