New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 1203246 =19 (0.600)8 (0.320) 6/134 1.4 33.8% coding (290/630 nt) BW25113_RS05990 hypothetical protein
?NZ_CP009273 1205075 = 16 (0.630)coding (40/540 nt) BW25113_RS06005 phage tail protein

ATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCG  >  NZ_CP009273/1204986‑1205162
                                                                                         |                                                                                       
aTTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAAc                                                                                         >  2:127468/1‑90 (MQ=255)
 ccaccaatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCg                                                                                        >  1:352273‑M2/89‑90 (MQ=255)
        tCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGa                                                                                 >  1:363193/1‑90 (MQ=255)
        tCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGa                                                                                 >  1:905156/1‑90 (MQ=255)
                         ggAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACaa                                                                >  2:314911/1‑90 (MQ=255)
                           gaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGc                                                              >  1:134489‑M2/63‑90 (MQ=255)
                                    gcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcat                                                     <  2:655448‑M2/37‑1 (MQ=255)
                                        cACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCAtt                                                 >  2:690879/1‑90 (MQ=255)
                                          cGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTcc                                               <  1:822105/90‑1 (MQ=255)
                                               tagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                                          <  2:544372‑M2/48‑1 (MQ=255)
                                               ccGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                                          <  1:477128/90‑1 (MQ=255)
                                                 gccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtttt                                        >  2:845746‑M2/41‑90 (MQ=255)
                                                 gccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtttt                                        >  2:290964‑M2/41‑90 (MQ=255)
                                                 gccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtttt                                        >  2:86057‑M2/41‑90 (MQ=255)
                                                      acCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATga                                                                           <  1:452231/50‑1 (MQ=255)
                                                      acCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATga                                                                           >  2:452231/1‑50 (MQ=255)
                                                               tAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCa                                    >  2:829987/1‑80 (MQ=255)
                                                               tAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCa                                    <  1:829987/80‑1 (MQ=255)
                                                                 aTGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAAc                        >  2:304932/1‑90 (MQ=255)
                                                                   tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACAt                      <  2:756278‑M2/68‑1 (MQ=255)
                                                                       agCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                                 <  1:453962/75‑1 (MQ=255)
                                                                       agCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                                 >  2:453962/1‑75 (MQ=255)
                                                                       agCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCtt                  <  1:127468/90‑1 (MQ=255)
                                                                       agCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCtt                  >  2:873744/1‑90 (MQ=255)
                                                                         ccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAt                <  1:725263‑M2/74‑1 (MQ=255)
                                                                         ccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAt                <  1:690879‑M2/74‑1 (MQ=255)
                                                                             ctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGc            >  2:268099‑M2/13‑90 (MQ=255)
                                                                             ctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGc            <  1:290964‑M2/78‑1 (MQ=255)
                                                                                  aaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGc                    >  2:359498‑M2/8‑77 (MQ=255)
                                                                                  aaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGc                    <  1:359498‑M2/70‑1 (MQ=255)
                                                                                   aaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTa      >  2:69964‑M2/7‑90 (MQ=255)
                                                                                       aaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAcgcg  >  2:152253‑M2/3‑90 (MQ=255)
                                                                                         |                                                                                       
ATTGCATTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCG  >  NZ_CP009273/1204986‑1205162

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