New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 1203246 =3 (0.140)14 (0.820) 11/134 0.4 83.8% coding (290/630 nt) BW25113_RS05990 hypothetical protein
?NZ_CP009273 1205075 = 3 (0.180)coding (40/540 nt) BW25113_RS06005 phage tail protein

TACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTA  >  NZ_CP009273/1203157‑1203335
                                                                                         |                                                                                         
tACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCT                                                                                           >  1:340131‑M1/1‑90 (MQ=255)
 tccGCCCCTCGCGCGCCGCTGTTAAGCCGCAAGGATCCGTTTGGCGATAGCAAAACGGATGGCACTTTGCAAACGGCACTCGGAATCCtt                                                                                          <  1:282924/89‑1 (MQ=17)
  ccGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTG                                                                                         >  2:485094‑M1/1‑90 (MQ=255)
    gCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGT                                                                                       >  1:113105‑M1/1‑90 (MQ=255)
         tagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTggg                                                                                  <  1:511382‑M1/10‑1 (MQ=255)
             caataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGaa                                                                              >  2:546375‑M1/77‑90 (MQ=255)
               ataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAgg                                                                            <  2:318708‑M1/16‑1 (MQ=255)
                 tgctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGct                                                                          >  1:174877‑M1/1‑88 (MQ=255)
                 tgctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGct                                                                          <  1:485094‑M1/90‑3 (MQ=255)
                 tgctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGct                                                                          <  2:189175‑M1/90‑3 (MQ=255)
                   ctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctct                                                                        <  1:471224‑M1/90‑5 (MQ=255)
                    catgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCaa                                                                       >  2:541521‑M1/70‑90 (MQ=255)
                     atgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCaaa                                                                      <  1:342348‑M1/22‑1 (MQ=255)
                     atgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCaaa                                                                      <  1:575287‑M1/22‑1 (MQ=255)
                       tAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAAc                                                                    <  2:183491/90‑1 (MQ=255)
                         aGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcatta                                                                  <  2:211125‑M1/90‑11 (MQ=255)
                            atccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAAt                                                               >  2:575287‑M1/62‑90 (MQ=255)
                                ccctggtaaaagcaggagtgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgcag                                                           <  1:401933‑M1/33‑1 (MQ=255)
                                ccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgcag                                                           <  2:157152‑M1/33‑1 (MQ=255)
                                ccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgcag                                                           <  2:35159‑M1/33‑1 (MQ=255)
                                  aaTCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggt                                                         <  1:552322‑M1/90‑20 (MQ=255)
                                    tCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgt                                                       >  2:110487‑M1/1‑69 (MQ=255)
                                    ggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAAc                                                       >  1:157152‑M1/54‑90 (MQ=255)
                                     gtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACg                                                      >  2:610720‑M1/53‑90 (MQ=255)
                                         aagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTa                                                  <  2:435941‑M1/42‑1 (MQ=255)
                                                tatCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctc                                                                         <  1:425244‑M1/60‑4 (MQ=255)
                                                tatCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctc                                                                         >  2:425244‑M1/1‑57 (MQ=255)
                                                 atCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGa                                          >  1:236611/1‑90 (MQ=255)
                                                       ccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCAcc                                    >  1:35159‑M1/35‑90 (MQ=255)
                                                        catttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCg                                   <  1:236025‑M1/57‑1 (MQ=255)
                                                          gATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagcc                                 >  1:483941‑M1/1‑47 (MQ=255)
                                                                      ttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCtata                     <  1:71757‑M1/71‑1 (MQ=255)
                                                                      ttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCtata                     >  2:401933‑M1/20‑90 (MQ=255)
                                                                          cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccaccaacaggct                 <  1:522495‑M1/90‑60 (MQ=255)
                                                                             ctctCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggc              >  1:109972‑M1/1‑28 (MQ=255)
                                                                                tcGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctga           <  2:113105‑M1/90‑66 (MQ=255)
                                                                                tcGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctga           <  2:483941‑M1/90‑66 (MQ=255)
                                                                                tcGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctga           <  1:110487‑M1/90‑66 (MQ=255)
                                                                                         ttGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggctgaaatgcaatg  >  2:565303‑M1/1‑16 (MQ=255)
                                                                                         |                                                                                         
TACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTA  >  NZ_CP009273/1203157‑1203335

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