New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 1203246 =39 (0.650)6 (0.120) 6/132 2.5 17.3% coding (290/630 nt) BW25113_RS05990 hypothetical protein
?NZ_CP009273 1205075 = 26 (0.540)coding (40/540 nt) BW25113_RS06005 phage tail protein
Rejected: Frequency below/above cutoff threshold.

GCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGT  >  NZ_CP009273/1204997‑1205163
                                                                              |                                                                                        
taaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGaaaa                                                                               >  2:365480‑M2/79‑90 (MQ=255)
gccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGaaaa                                                                               >  2:410570/1‑90 (MQ=255)
   agccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATga                                                                            <  2:673724/90‑1 (MQ=255)
               gAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAg                                                                <  2:1703788/90‑1 (MQ=255)
                aGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGc                                                               >  1:62320/1‑90 (MQ=255)
                  tGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCaa                                                             >  2:1548771/1‑90 (MQ=255)
                    gCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAAt                                                           <  2:1623155/90‑1 (MQ=255)
                          ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                                                     <  1:1404477/90‑1 (MQ=255)
                          ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                                                     <  2:571724/90‑1 (MQ=255)
                             cACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                                                  >  1:1172263/1‑90 (MQ=255)
                              aCGGAACCGGTACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTc                                                 >  1:1016824/1‑90 (MQ=255)
                               taatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                                                <  2:1556963‑M2/43‑1 (MQ=255)
                               cGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                                                <  1:889972/90‑1 (MQ=255)
                                    ccGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcat                                                      <  1:1109723/79‑1 (MQ=255)
                                    ccGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcat                                                      >  2:1109723/1‑79 (MQ=255)
                                    ccGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                                           <  2:34060/90‑1 (MQ=255)
                                     agccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaaga                                                                         >  1:790092‑M2/42‑59 (MQ=255)
                                     agccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATgaaga                                                                         <  2:790092‑M2/18‑1 (MQ=255)
                                         gttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTAcc                                      >  1:910705‑M2/38‑90 (MQ=255)
                                          ttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCa                                     >  1:1209443‑M2/37‑90 (MQ=255)
                                            acAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgg                                   <  2:100123/89‑1 (MQ=255)
                                             cAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACa                                                                  <  1:1779188/58‑1 (MQ=255)
                                             cAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACa                                                                  >  2:1779188/1‑58 (MQ=255)
                                             cAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                                  >  2:1310700/1‑90 (MQ=255)
                                                 aGGTAATGCAGAGCCTTATCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAt                              >  1:646643/1‑90 (MQ=255)
                                                 aGGTAATGCAGAGCCTTATCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAt                              >  1:953120/1‑90 (MQ=255)
                                                       tgaagagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACa                        <  1:1670483/87‑1 (MQ=255)
                                                        gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACAt                       <  1:71359/90‑1 (MQ=255)
                                                           gagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCt                    >  1:1585091/1‑90 (MQ=255)
                                                            agCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                                  <  2:96165/75‑1 (MQ=255)
                                                            agCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                                  >  1:96165/1‑75 (MQ=255)
                                                            agCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCtt                   >  1:974857/1‑90 (MQ=255)
                                                              ccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAt                 <  2:113480‑M2/74‑1 (MQ=255)
                                                              ccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAt                 <  1:1310700‑M2/74‑1 (MQ=255)
                                                                ttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTg               >  1:478079‑M2/15‑90 (MQ=255)
                                                                  ctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGc             <  1:410570‑M2/78‑1 (MQ=255)
                                                                  ctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGc             <  2:669361‑M2/78‑1 (MQ=255)
                                                                  ctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGc             >  1:1488133‑M2/13‑90 (MQ=255)
                                                                    cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTa           <  2:1714892‑M2/80‑1 (MQ=255)
                                                                       aaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGt        >  2:1461351‑M2/8‑90 (MQ=255)
                                                                         accaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAc      <  2:478079‑M2/85‑1 (MQ=255)
                                                                           caaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAcgc    <  1:1548771‑M2/87‑1 (MQ=255)
                                                                             aCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGt  >  1:882921‑M2/2‑90 (MQ=255)
                                                                             aCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGt  >  1:357971‑M2/2‑90 (MQ=255)
                                                                              |                                                                                        
GCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCGT  >  NZ_CP009273/1204997‑1205163

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