New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 120780542 (0.640)18 (0.330) 17/166 3.7 34.6% coding (305/630 nt) stfP e14 prophage; uncharacterized protein
?NC_000913 = 1209602 33 (0.600)pseudogene (18/501 nt) stfE pseudogene, e14 prophage; side tail fiber protein fragment family;Phage or Prophage Related
Rejected: Coverage evenness skew score above cutoff.

ATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGT  >  NC_000913/1209505‑1209701
                                                                                                 |                                                                                                   
aTTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGcc                                                                                                   >  1:949071/1‑100 (MQ=255)
 ttcttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGcc                                                                                                   <  2:984297/99‑1 (MQ=255)
   cttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTtc                                                                                                >  1:912005/1‑100 (MQ=255)
   cttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTtc                                                                                                >  1:87588/1‑100 (MQ=255)
   cttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTtc                                                                                                >  1:660311/1‑100 (MQ=255)
   cttcAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACAGGGACACCAACAGGTAATGCAGAGCCTtc                                                                                                >  1:959939/1‑100 (MQ=255)
      cagcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTcc                                                                                             <  2:1540036/100‑1 (MQ=255)
       agcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTccc                                                                                            <  2:912012/100‑1 (MQ=255)
       agcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTccc                                                                                            <  2:660316/100‑1 (MQ=255)
       agcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTccc                                                                                            >  2:463372/1‑100 (MQ=255)
       agcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTccc                                                                                            >  2:1365299/1‑100 (MQ=255)
        gcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCa                                                                                           >  2:809142/1‑100 (MQ=255)
          agAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCaaa                                                                                         >  2:352639/1‑100 (MQ=255)
          agAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCaaa                                                                                         <  1:878404/100‑1 (MQ=255)
           gAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACAGCAACAGGTAATGCAGAGCCTTCTCCCAAAc                                                                                        >  1:249283/1‑100 (MQ=255)
                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     <  2:542497/100‑1 (MQ=255)
                                                                                              agagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTa     <  1:1191846/100‑1 (MQ=255)
                                                                                                 gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGt  >  1:257517/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
ATTCTTCAGCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGT  >  NC_000913/1209505‑1209701

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.