New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 12032612 (0.090)15 (0.880) 12/134 0.3 86.7% coding (305/630 nt) BW25113_RS05990 hypothetical protein
?NZ_CP009273 = 1205058 3 (0.180)coding (57/540 nt) BW25113_RS06005 phage tail protein

CATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTT  >  NZ_CP009273/1204985‑1205146
                                                                         |                                                                                        
cATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaa                                                                          >  2:565303‑M2/1‑74 (MQ=255)
         tCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcga                                                                 <  1:110487‑M2/90‑26 (MQ=255)
         tCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcga                                                                 <  2:483941‑M2/90‑26 (MQ=255)
         tCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcga                                                                 <  2:113105‑M2/90‑26 (MQ=255)
            gccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagag                                                              >  1:109972‑M2/1‑62 (MQ=255)
               agccTGTTGGTGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccg                                                           <  1:522495‑M2/90‑32 (MQ=255)
                               ggCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatc                                           >  1:483941‑M2/1‑43 (MQ=255)
                                      ggCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatc                                    <  2:198385/90‑1 (MQ=255)
                                           cGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                               <  1:266261/90‑1 (MQ=255)
                                                     acacCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacgga                     >  2:110487‑M2/1‑21 (MQ=255)
                                                       acCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggatt                   <  1:552322‑M2/90‑72 (MQ=255)
                                                         cAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                 <  1:8118/90‑1 (MQ=255)
                                                                tAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggct          <  2:211125‑M2/90‑81 (MQ=255)
                                                                      agagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacag    <  1:471224‑M2/90‑87 (MQ=255)
                                                                       gagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgata                                 <  1:425244‑M2/60‑58 (MQ=255)
                                                                       gagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgata                                 >  2:425244‑M2/1‑3 (MQ=255)
                                                                        agccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagca  <  1:485094‑M2/90‑89 (MQ=255)
                                                                        agccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagca  <  2:189175‑M2/90‑89 (MQ=255)
                                                                        agccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagca  >  1:174877‑M2/1‑2 (MQ=255)
                                                                         |                                                                                        
CATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTT  >  NZ_CP009273/1204985‑1205146

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