New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 1203246 =25 (0.730)7 (0.250) 7/138 1.4 31.6% coding (290/630 nt) BW25113_RS05990 hypothetical protein
?NZ_CP009273 1205075 = 10 (0.360)coding (40/540 nt) BW25113_RS06005 phage tail protein

TCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAT  >  NZ_CP009273/1204994‑1205148
                                                                                 |                                                                         
tCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGa                                                                   >  2:36692/1‑90 (MQ=255)
   gccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGaaaa                                                                >  1:384765/1‑90 (MQ=255)
                      ggCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCaaa                                             >  2:58383/1‑90 (MQ=255)
                        cggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAGATg                                           >  2:823228‑M2/58‑90 (MQ=255)
                            gcgccctaatgtagccgttgctgcattgagttttgccccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcat                                       <  1:562785‑M2/37‑1 (MQ=255)
                                 aCGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAAc                                                    <  1:847549/72‑1 (MQ=255)
                                 aCGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAAc                                                    >  2:847549/1‑72 (MQ=255)
                                      gtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCt                             <  1:283610‑M2/47‑1 (MQ=255)
                                       tagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAGATGGCATCATTCCTGCtt                            <  1:823228‑M2/48‑1 (MQ=255)
                                        agccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttt                           >  2:463344‑M2/42‑90 (MQ=255)
                                             ttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCa                      >  2:397372‑M2/37‑90 (MQ=255)
                                                cAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                   <  2:394773/90‑1 (MQ=255)
                                                   cattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGa                <  1:397372‑M2/60‑1 (MQ=255)
                                                       tAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTa            <  1:170616/90‑1 (MQ=255)
                                                           gCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACAt        >  1:467487/1‑90 (MQ=255)
                                                               agCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCtt    >  1:303758/1‑90 (MQ=255)
                                                                gCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGGATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTa   >  2:600369/1‑90 (MQ=255)
                                                                 ccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAt  <  1:81037‑M2/74‑1 (MQ=255)
                                                                          aaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                                      <  2:873469‑M2/38‑1 (MQ=255)
                                                                          aaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                                      >  1:873469‑M2/8‑45 (MQ=255)
                                                                                 |                                                                         
TCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAT  >  NZ_CP009273/1204994‑1205148

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