New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 1203246 =5 (0.280)5 (0.350) 4/136 1.0 62.4% coding (290/630 nt) BW25113_RS05990 hypothetical protein
?NZ_CP009273 1205075 = 2 (0.140)coding (40/540 nt) BW25113_RS06005 phage tail protein

GCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTA  >  NZ_CP009273/1203161‑1203323
                                                                                     |                                                                             
gtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGt                                                                           <  2:103742‑M1/5‑1 (MQ=255)
gCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGT                                                                           >  1:84664‑M1/1‑90 (MQ=255)
     tCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGG                                                                      >  2:268977‑M1/1‑90 (MQ=255)
             tgctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGct                                                              >  1:149129‑M1/1‑88 (MQ=255)
                 atgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCaaa                                                          <  1:119097‑M1/22‑1 (MQ=255)
                 atgttaaatccccctggtaaaagcaggaatgataccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCaaa                                                          <  2:37568‑M1/22‑1 (MQ=255)
                        cGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctc                                                             <  1:247370‑M1/80‑4 (MQ=255)
                        cGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctc                                                             >  2:247370‑M1/1‑77 (MQ=255)
                              aaTCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggt                                             <  2:214425‑M1/90‑20 (MQ=255)
                               tggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCaa                                            >  1:367575‑M1/55‑90 (MQ=255)
                                      tGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGc                                                             >  1:149073/1‑66 (MQ=255)
                                      tGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGc                                                             <  2:149073/66‑1 (MQ=255)
                                      tGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTAc                                     >  1:164307/1‑90 (MQ=255)
                                          ggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTa                                 <  2:354994‑M1/47‑1 (MQ=255)
                                            tatCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtg                               <  1:225780‑M1/90‑34 (MQ=255)
                                                   tCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCAcc                        <  2:164307/90‑1 (MQ=255)
                                                         ggCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTc                  <  2:173728/90‑1 (MQ=255)
                                                         gcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTc                  >  2:356135‑M1/29‑90 (MQ=255)
                                                          gCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGG                                                                >  2:230142‑M1/1‑43 (MQ=255)
                                                          gCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGG                                                                <  1:230142‑M1/43‑1 (MQ=255)
                                                               gtgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgcca            >  2:398333‑M1/1‑38 (MQ=255)
                                                                         ctctCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggc  >  2:53792‑M1/1‑28 (MQ=255)
                                                                                     |                                                                             
GCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTA  >  NZ_CP009273/1203161‑1203323

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