New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 12032611 (0.030)22 (0.880) 15/136 0.4 90.1% coding (305/630 nt) BW25113_RS05990 hypothetical protein
?NZ_CP009273 = 1205058 4 (0.160)coding (57/540 nt) BW25113_RS06005 phage tail protein

GCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTA  >  NZ_CP009273/1204969‑1205147
                                                                                         |                                                                                         
gcagAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCagag                                                                                           >  2:263989/1‑90 (MQ=255)
     aaaaGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGGAATGCAGAGccttc                                                                                      <  2:863467‑M2/90‑6 (MQ=255)
           ggCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaa                                                                                <  2:593267‑M2/90‑12 (MQ=255)
                   tGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggt                                                                        >  2:385017‑M2/1‑71 (MQ=255)
                   tGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggt                                                                        >  2:317163‑M2/1‑71 (MQ=255)
                   tGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggt                                                                        >  1:74671‑M2/1‑71 (MQ=255)
                         tCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcga                                                                  <  2:252562‑M2/90‑26 (MQ=255)
                            gccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagag                                                               >  2:249008‑M2/1‑62 (MQ=255)
                            gccagccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagag                                                               >  2:736393‑M2/1‑62 (MQ=255)
                                gccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgt                                                           <  2:383313‑M2/90‑33 (MQ=255)
                                 ccTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtt                                                          >  2:110696‑M2/1‑57 (MQ=255)
                                      tGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcac                                                     >  2:108656‑M2/1‑52 (MQ=255)
                                      tGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcac                                                     >  1:593267‑M2/1‑52 (MQ=255)
                                          ggAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtg                                                 <  1:745789‑M2/90‑43 (MQ=255)
                                          ggAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACaa                                                 >  1:824198/1‑90 (MQ=255)
                                              tGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCaa                                             >  2:212589/1‑90 (MQ=255)
                                                       gCCACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttga                                    <  1:492156‑M2/90‑56 (MQ=255)
                                                         cACGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgata                                  <  1:512853‑M2/90‑58 (MQ=255)
                                                          aCGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagag                                                               <  1:853542‑M2/60‑29 (MQ=255)
                                                          aCGGAACCGGGACACCAACAGGTAATGCAGAGccttctcccaaaccaaggttttcgagag                                                               >  2:853542‑M2/1‑32 (MQ=255)
                                                                             aGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttt                          <  1:167003/78‑1 (MQ=255)
                                                                             aGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttt                          >  2:167003/1‑78 (MQ=255)
                                                                               gTAATGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggc            >  2:467348‑M2/1‑11 (MQ=255)
                                                                                  aTGCAGAGccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcaccaaacggattcttgcggctta         <  1:789265‑M2/90‑83 (MQ=255)
                                                                                      agagccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacag     <  1:744590‑M2/90‑87 (MQ=255)
                                                                                        agccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagca   <  2:325105‑M2/90‑89 (MQ=255)
                                                                                        agccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagca   <  2:668798‑M2/90‑89 (MQ=255)
                                                                                        agccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagca   <  2:791388‑M2/90‑89 (MQ=255)
                                                                                         gccttctcccaaaccaaggttttcgagagccgtttgcacagtgccatccgatttgatatcgccaaacggattcttgcggcttaacagcag  <  2:650316‑M2/90‑90 (MQ=255)
                                                                                         |                                                                                         
GCAGAAAAAGCGGCACCATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTA  >  NZ_CP009273/1204969‑1205147

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