New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 1203246 =0 (0.000)19 (0.830) 12/132 0.4 92.7% coding (290/630 nt) BW25113_RS05990 hypothetical protein
?NZ_CP009273 1205075 = 3 (0.130)coding (40/540 nt) BW25113_RS06005 phage tail protein

GCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGA  >  NZ_CP009273/1203161‑1203326
                                                                                     |                                                                                
gCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGT                                                                              >  1:42620‑M1/1‑90 (MQ=255)
     tagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTggg                                                                         <  2:550708‑M1/10‑1 (MQ=255)
       gccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGgag                                                                       <  2:531928‑M1/12‑1 (MQ=255)
       gccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGgag                                                                       <  2:224389‑M1/12‑1 (MQ=255)
         caataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGaa                                                                     <  1:154194‑M1/14‑1 (MQ=255)
             tgctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGct                                                                 <  1:8296‑M1/90‑3 (MQ=255)
              gctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctc                                                                <  2:42620‑M1/90‑4 (MQ=255)
                     taaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTc                                                         >  1:550708‑M1/65‑90 (MQ=255)
                     taaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTc                                                         >  2:398646‑M1/65‑90 (MQ=255)
                      gCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattac                                                        >  1:823242‑M1/1‑79 (MQ=255)
                            ccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgcag                                                  <  2:94457‑M1/33‑1 (MQ=255)
                            ccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgcag                                                  <  2:447103‑M1/33‑1 (MQ=255)
                            ccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgcag                                                  <  2:316549‑M1/33‑1 (MQ=255)
                             gAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttgg                                                 <  2:345075‑M1/90‑19 (MQ=255)
                              ctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgcagca                                                >  2:784554‑M1/56‑90 (MQ=255)
                                    aaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCt                                          >  1:54994‑M1/50‑90 (MQ=255)
                                          ggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACAtt                                     >  1:679300‑M1/44‑89 (MQ=255)
                                          ggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACAtt                                     <  2:679300‑M1/46‑1 (MQ=255)
                                          ggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTa                                    <  2:728848‑M1/47‑1 (MQ=255)
                                          ggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTa                                    <  1:401858‑M1/47‑1 (MQ=255)
                                            tatCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtg                                  <  1:108648‑M1/90‑34 (MQ=255)
                                            tatCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtg                                  <  2:823242‑M1/90‑34 (MQ=255)
                                               gatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACg                               <  2:576383‑M1/52‑1 (MQ=255)
                                               gatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACg                               <  2:107624‑M1/52‑1 (MQ=255)
                                                 aaTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtccc                                           <  1:555974‑M1/76‑25 (MQ=255)
                                                 aaTCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtccc                                           >  2:555974‑M1/1‑52 (MQ=255)
                                                         gcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTc                     >  1:107624‑M1/29‑90 (MQ=255)
                                                           ttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCAt                   >  2:401858‑M1/27‑90 (MQ=255)
                                                           ttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACAGCTACATTAGGACGCACCGGTTTCAt                   >  1:531928‑M1/27‑90 (MQ=255)
                                                               gtgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgcca               >  1:727197‑M1/1‑38 (MQ=255)
                                                                      tcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCtat             <  2:726277‑M1/70‑1 (MQ=255)
                                                                      tcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCtat             >  1:726277‑M1/16‑85 (MQ=255)
                                                                         ctctCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggc     >  2:540642‑M1/1‑28 (MQ=255)
                                                                            tcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGa  >  2:66712‑M1/10‑90 (MQ=255)
                                                                                     |                                                                                
GCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGA  >  NZ_CP009273/1203161‑1203326

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