New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 1203246 =4 (0.080)22 (0.580) 17/136 0.6 89.4% coding (290/630 nt) BW25113_RS05990 hypothetical protein
?NZ_CP009273 1205075 = 2 (0.050)coding (40/540 nt) BW25113_RS06005 phage tail protein

TTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCG  >  NZ_CP009273/1204992‑1205162
                                                                                   |                                                                                       
aatcataaccggtatagctatgaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTAt                                                                                   <  1:986273‑M2/7‑1 (MQ=255)
  tCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGa                                                                                 >  2:239165/1‑90 (MQ=255)
                     gaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGc                                                              >  2:49655‑M2/63‑90 (MQ=255)
                     gaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGc                                                              >  2:125783‑M2/63‑90 (MQ=255)
                     gaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGc                                                              >  2:118995‑M2/63‑90 (MQ=255)
                     gaaaccggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaAGTTTATGAAAATGAAGAAATAACAAGc                                                              >  1:133160‑M2/63‑90 (MQ=255)
                           ggtgcgtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATgg                                                        >  2:408221‑M2/57‑90 (MQ=255)
                                gtcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatca                                                   >  2:71226‑M2/52‑90 (MQ=255)
                                 tcctaatgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGcatcat                                                  >  2:295640‑M2/51‑90 (MQ=255)
                                      atgtagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTg                                             <  1:1215404‑M2/45‑1 (MQ=255)
                                         tagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                                          <  2:597428‑M2/48‑1 (MQ=255)
                                         tagccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtt                                          <  1:1091931‑M2/48‑1 (MQ=255)
                                          agccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCttt                                         >  2:1359090‑M2/42‑90 (MQ=255)
                                           gccgttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCtttt                                        >  1:1330558‑M2/41‑90 (MQ=255)
                                               ttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCa                                    >  2:676350‑M2/37‑90 (MQ=255)
                                               ttgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCa                                    >  2:661338‑M2/37‑90 (MQ=255)
                                                tgctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAg                                   >  1:549337‑M2/36‑90 (MQ=255)
                                                 gctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAgg                                  <  1:1341797‑M2/56‑1 (MQ=255)
                                                  ctgcattgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAggg                                 <  1:49655‑M2/57‑1 (MQ=255)
                                                      attgagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGAt                             >  2:260071‑M2/30‑90 (MQ=255)
                                                         gagttttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTa                          >  2:994196‑M2/27‑90 (MQ=255)
                                                             tttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACAt                      <  1:676350‑M2/68‑1 (MQ=255)
                                                              ttgctccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATg                     >  1:896460‑M2/22‑90 (MQ=255)
                                                               agagCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGc                    >  2:1232940/1‑90 (MQ=255)
                                                                   ccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAt                <  1:260071‑M2/74‑1 (MQ=255)
                                                                   ccttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAt                <  2:520934‑M2/74‑1 (MQ=255)
                                                                    cttctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTAtt               <  2:415920‑M2/75‑1 (MQ=255)
                                                                       ctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGc            <  1:71226‑M2/78‑1 (MQ=255)
                                                                       ctcccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATAGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGc            <  2:414614‑M2/78‑1 (MQ=255)
                                                                       ctcccaaaccaaAGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGc            <  2:133160‑M2/78‑1 (MQ=255)
                                                                         cccaaaccaaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTa          <  2:343450‑M2/80‑1 (MQ=255)
                                                                                 aaCGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTAcgcg  >  2:150053‑M2/3‑90 (MQ=255)
                                                                                   |                                                                                       
TTTCAGCCAGCCTGTTGGCGGAGTGGCTGAAGGCCACGGAACCGGGACACCAACAGGTAATGCAGAGCCTTCTCCCAAACCAACGTTTATGAAAATGAAGAAATAACAAGCAAATGGCATCATTCCTGCTTTTACCAGGGGGATTTAACATGCTTATTGGCTATGTACGCG  >  NZ_CP009273/1204992‑1205162

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