New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 12032610 (0.000)15 (1.030) 12/132 0.2 93.8% coding (305/630 nt) BW25113_RS05990 hypothetical protein
?NZ_CP009273 = 1205058 2 (0.140)coding (57/540 nt) BW25113_RS06005 phage tail protein

GCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCT  >  NZ_CP009273/1204989‑1205145
                                                                     |                                                                                       
gCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggtt                                                                     <  2:365576‑M2/90‑21 (MQ=255)
     tCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcga                                                                <  1:45653‑M2/90‑26 (MQ=255)
        gccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagag                                                             <  1:197531‑M2/90‑29 (MQ=255)
             ccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtt                                                        >  1:467555‑M2/1‑57 (MQ=255)
                  tGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcac                                                   >  2:87451‑M2/1‑52 (MQ=255)
                  tGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcac                                                   >  1:262450‑M2/1‑52 (MQ=255)
                      ggAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtg                                               <  1:177571‑M2/90‑43 (MQ=255)
                      ggAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACaa                                               >  2:282264/1‑90 (MQ=255)
                                    ccACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgat                                 >  1:94971‑M2/1‑34 (MQ=255)
                                                    ccAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagag                                                             >  2:165053‑M2/1‑18 (MQ=255)
                                                    ccAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagag                                                             <  1:165053‑M2/46‑29 (MQ=255)
                                                      aaCAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccggtttgatatcgcaaaacggattctt               <  2:155280‑M2/90‑75 (MQ=255)
                                                         aGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcg            >  1:203332‑M2/1‑13 (MQ=255)
                                                                   gagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagc  <  1:339065‑M2/90‑88 (MQ=255)
                                                                   gagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagc  <  2:196269‑M2/90‑88 (MQ=255)
                                                                   gagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagc  <  2:204059‑M2/90‑88 (MQ=255)
                                                                   gagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCt  >  2:359126/1‑90 (MQ=255)
                                                                     |                                                                                       
GCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCT  >  NZ_CP009273/1204989‑1205145

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