New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 1203246 =10 (0.500)11 (0.690) 7/134 0.7 66.6% coding (290/630 nt) BW25113_RS05990 hypothetical protein
?NZ_CP009273 1205075 = 3 (0.190)coding (40/540 nt) BW25113_RS06005 phage tail protein

TACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGT  >  NZ_CP009273/1203157‑1203331
                                                                                         |                                                                                     
tacCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCt                                                                                       >  1:570949/1‑90 (MQ=255)
 cgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgtt                                                                                      >  1:276453‑M1/89‑90 (MQ=255)
 cgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgtt                                                                                      <  1:303279‑M1/2‑1 (MQ=255)
 cgcgtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgtt                                                                                      >  2:495181‑M1/89‑90 (MQ=255)
    gtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGt                                                                                   <  2:343261‑M1/5‑1 (MQ=255)
    gCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGt                                                                                   >  1:340548/1‑90 (MQ=255)
           gccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGgag                                                                            <  2:89994‑M1/12‑1 (MQ=255)
            cgcgcTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGgaga                                                                           >  1:176552/1‑90 (MQ=255)
            ccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGgaga                                                                           >  1:302728‑M1/78‑90 (MQ=255)
               gctgctGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAgg                                                                        <  2:340881/90‑1 (MQ=255)
               ataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAgg                                                                        >  2:257086‑M1/75‑90 (MQ=255)
               ataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAgg                                                                        <  2:370545‑M1/16‑1 (MQ=255)
               ataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAgg                                                                        <  2:54535‑M1/16‑1 (MQ=255)
                 tgctgTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAggct                                                                      <  1:429637/90‑3 (MQ=255)
                    catgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCaa                                                                   >  1:152335‑M1/70‑90 (MQ=255)
                     atgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCaaa                                                                  <  2:561775‑M1/22‑1 (MQ=255)
                                       aaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGc                                                >  2:380818‑M1/51‑90 (MQ=255)
                                       aaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGc                                                >  1:188939‑M1/51‑90 (MQ=255)
                                       aaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGc                                                >  1:161800‑M1/51‑90 (MQ=255)
                                                        catttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCg                               >  2:303279‑M1/34‑90 (MQ=255)
                                                         ggATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAggc                                                                       >  1:15436/1‑48 (MQ=255)
                                                         ggATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAggc                                                                       <  2:15436/49‑2 (MQ=255)
                                                                       tcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATAc                <  2:502758‑M1/72‑1 (MQ=255)
                                                                       aaaCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCCCCGGTTTCATAGCTATAc                <  2:340548/90‑1 (MQ=255)
                                                                             ttttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTa          >  1:30637‑M1/13‑90 (MQ=255)
                                                                             ctctCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTa          >  1:583284/1‑90 (MQ=255)
                                                                             ctctCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTa          >  1:576849/1‑90 (MQ=255)
                                                                                     aacgTTGGTTTGGGAGAAGGAGCGAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggt  >  1:89994‑M1/5‑90 (MQ=255)
                                                                                     aacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggt  >  1:561775‑M1/5‑90 (MQ=255)
                                                                                     aacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATtggt  >  1:502758‑M1/5‑90 (MQ=255)
                                                                                         |                                                                                     
TACCGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGT  >  NZ_CP009273/1203157‑1203331

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