New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 1203246 =10 (0.600)12 (0.900) 12/134 0.2 63.0% coding (290/630 nt) BW25113_RS05990 hypothetical protein
?NZ_CP009273 1205075 = 6 (0.450)coding (40/540 nt) BW25113_RS06005 phage tail protein

GCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGT  >  NZ_CP009273/1203168‑1203331
                                                                              |                                                                                     
gccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGgag                                                                            <  1:490306‑M1/12‑1 (MQ=255)
  gcgcTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAA                                                                          <  2:490497‑M1/90‑1 (MQ=255)
         tgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCaa                                                                   <  1:12883/90‑1 (MQ=255)
              aGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTc                                                              <  1:472641/90‑1 (MQ=255)
                 cGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAAt                                                           >  2:103517/1‑90 (MQ=255)
                    aaGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgtt                                                        >  2:182641‑M1/1‑74 (MQ=255)
                     ccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgcag                                                       <  1:151865‑M1/33‑1 (MQ=255)
                         ggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgcagca                                                     <  2:37237‑M1/35‑1 (MQ=255)
                         ggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgcagca                                                     >  1:37237‑M1/54‑88 (MQ=255)
                             tttGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCt                                               <  1:344385/90‑1 (MQ=255)
                              aagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTa                                              >  2:479854‑M1/49‑90 (MQ=255)
                                  aggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACAtt                                          >  2:490306‑M1/45‑90 (MQ=255)
                                    atatCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgt                                        >  1:210988‑M1/1‑58 (MQ=255)
                                      atCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtgg                                      >  1:458070‑M1/1‑56 (MQ=255)
                                       tgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGAc                                     >  1:133861‑M1/40‑90 (MQ=255)
                                            ccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCAcc                                >  1:58536‑M1/35‑90 (MQ=255)
                                               gATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagcc                             >  2:310091‑M1/1‑47 (MQ=255)
                                                aTGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagcca                            <  1:271790‑M1/90‑45 (MQ=255)
                                                  gcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTc                          >  1:130654‑M1/29‑90 (MQ=255)
                                                        gtgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTc                                                              >  1:328549/1‑48 (MQ=255)
                                                        gtgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTc                                                              <  2:328549/48‑1 (MQ=255)
                                                        atgCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAACAAAACTCAATGCAGCAACGGCTACATTAGGACACACCGGTTTCCTAGCt                    <  2:281041/89‑1 (MQ=38)
                                                            tcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATAc                >  1:173046‑M1/19‑90 (MQ=255)
                                                             aaCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATAcc               >  2:190171/1‑90 (MQ=255)
                                                              aCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggc              <  2:210988‑M1/90‑59 (MQ=255)
                                                               cGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggct             <  2:17734‑M1/90‑60 (MQ=255)
                                                                ggCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctg            >  1:234647‑M1/1‑30 (MQ=255)
                                                                 gCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGtt           <  1:479543/90‑1 (MQ=255)
                                                                     tcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGa       >  1:231533‑M1/10‑90 (MQ=255)
                                                                       ataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGAtt     <  1:431412‑M1/83‑1 (MQ=255)
                                                                          aacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggt  >  1:218060‑M1/5‑90 (MQ=255)
                                                                          aaCCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggt  <  1:103517/90‑1 (MQ=255)
                                                                              |                                                                                     
GCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGT  >  NZ_CP009273/1203168‑1203331

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