New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 12032616 (0.250)17 (0.870) 11/136 0.5 83.2% coding (305/630 nt) BW25113_RS05990 hypothetical protein
?NZ_CP009273 = 1205058 2 (0.100)coding (57/540 nt) BW25113_RS06005 phage tail protein

AGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCT  >  NZ_CP009273/1204977‑1205145
                                                                                 |                                                                                       
aGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcc                                                                                 >  1:666793‑M2/1‑82 (MQ=255)
        cATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaa                                                                         >  1:480923‑M2/1‑74 (MQ=255)
                 tCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGa                                                                >  1:104461/1‑90 (MQ=255)
                    gccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagag                                                             >  2:229379‑M2/1‑62 (MQ=255)
                         ccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtt                                                        >  1:352390‑M2/1‑57 (MQ=255)
                         ccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtt                                                        >  1:642173‑M2/1‑57 (MQ=255)
                          cTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATgaaga                                                       <  2:66136/90‑1 (MQ=255)
                              tGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcac                                                   <  2:154769‑M2/90‑39 (MQ=255)
                                        gCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatcc                                         >  2:265168‑M2/1‑42 (MQ=255)
                                             aGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgattt                                    <  1:23241‑M2/90‑54 (MQ=255)
                                                        ccGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaa                         <  2:517987‑M2/90‑65 (MQ=255)
                                                                 cAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattct                <  1:694139‑M2/90‑74 (MQ=255)
                                                                       gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggc          >  2:503524‑M2/1‑11 (MQ=255)
                                                                          aTGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggc          <  1:287880‑M2/87‑80 (MQ=255)
                                                                          aTGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggc          >  2:287880‑M2/1‑8 (MQ=255)
                                                                          aTGCAGAGccttctcccaaaccaaggttttcgagag                                                             <  2:693395‑M2/36‑29 (MQ=255)
                                                                          aTGCAGAGccttctcccaaaccaaggttttcgagag                                                             >  1:693395‑M2/1‑8 (MQ=255)
                                                                              agagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacag   <  2:326571‑M2/90‑87 (MQ=255)
                                                                              agagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacag   <  1:223650‑M2/90‑87 (MQ=255)
                                                                               gagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagc  <  2:526752‑M2/90‑88 (MQ=255)
                                                                               gagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagc  <  1:593091‑M2/90‑88 (MQ=255)
                                                                                 |                                                                                       
AGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCT  >  NZ_CP009273/1204977‑1205145

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