New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 1203246 =1 (0.080)4 (0.400) 4/134 0.7 74.0% coding (290/630 nt) BW25113_RS05990 hypothetical protein
?NZ_CP009273 1205075 = 2 (0.200)coding (40/540 nt) BW25113_RS06005 phage tail protein

CGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTA  >  NZ_CP009273/1203160‑1203323
                                                                                      |                                                                             
cGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGG                                                                            >  1:218528‑M1/1‑90 (MQ=255)
 gtacatagccaataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGt                                                                           <  1:239610‑M1/5‑1 (MQ=255)
          gcgcTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAA                                                                  <  1:104942‑M1/90‑1 (MQ=255)
          gcgcTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAA                                                                  <  2:151998‑M1/90‑1 (MQ=255)
            ataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGttt                                                                         >  1:239162‑M1/75‑81 (MQ=255)
            ataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGttt                                                                         <  2:239162‑M1/7‑1 (MQ=255)
            ataagcatgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAgg                                                                <  1:46142‑M1/16‑1 (MQ=255)
                  atgttaaatccccctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCaaa                                                          <  1:233385‑M1/22‑1 (MQ=255)
                              cctggtaaaagcaggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATgcagc                                              <  2:182546‑M1/34‑1 (MQ=255)
                                aTCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCaa                                            <  1:232725/90‑1 (MQ=255)
                                      ttGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccgg                                      <  1:320332‑M1/90‑27 (MQ=255)
                                           ggaatgatgccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTa                                 <  1:228379‑M1/47‑1 (MQ=255)
                                                    ccatttgcttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCAcc                        >  2:254887‑M1/35‑90 (MQ=255)
                                                       gATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagcc                     >  1:164862‑M1/1‑47 (MQ=255)
                                                                gtgCAAACGGCTCTCGAAATCCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgcca            <  2:212473‑M1/90‑53 (MQ=255)
                                                                     aaCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacagg       >  2:323652‑M1/1‑33 (MQ=255)
                                                                     aaCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacagg       >  1:163788‑M1/1‑33 (MQ=255)
                                                                      aCGGCTCTCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggc      <  1:227765‑M1/90‑59 (MQ=255)
                                                                          ctctCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggc  >  2:171385‑M1/1‑28 (MQ=255)
                                                                          ctctCGAAAACCTTGGTTTGGGAGAAGGctctgcattacctgttggtgtcccggttccgtggccttcagccactccgccaacaggctggc  >  1:66765‑M1/1‑28 (MQ=255)
                                                                                      |                                                                             
CGCACTTCGCGCGCTGCTGTTAAGCCGCAAGAATCCGTTTGGCGATATCAAATCGGATGGCACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTA  >  NZ_CP009273/1203160‑1203323

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