New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 357534827 (1.080)7 (0.280) 3/168 2.4 20.6% coding (126/1179 nt) hcp Hcp1 family type VI secretion system effector
?NZ_CP009273 3576788 = NA (NA)coding (391/417 nt) BW25113_RS17825 YrhA family protein

TGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGT  >  NZ_CP009273/3576709‑3576866
                                                                               |                                                                              
tgatTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGAtgctgc                                                                      >  2:324062/1‑90 (MQ=255)
                gATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAg                                                      <  1:89189/90‑1 (MQ=255)
                     aCAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCa                                                                    <  2:437472/71‑1 (MQ=255)
                     aCAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCa                                                                    >  1:437472/1‑71 (MQ=255)
                        aGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGAttt                                                        <  1:275934/80‑1 (MQ=255)
                        aGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGAttt                                                        >  2:275934/1‑80 (MQ=255)
                        aGAACGGTTAGTGATCGGGGATTATAGCATTTCAATACTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTatgat                                              >  2:499496/1‑90 (MQ=255)
                          aaCGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATgg                                            >  1:634576/1‑90 (MQ=255)
                           aCGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGgt                                           <  2:634576/90‑1 (MQ=255)
                                tAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGtttt                                      >  2:471611/1‑90 (MQ=255)
                                       cGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTg                               >  1:536448/1‑90 (MQ=255)
                                        ggggATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTAc                                                              >  1:599228/1‑58 (MQ=255)
                                        ggggATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTAc                                                              <  2:599228/58‑1 (MQ=255)
                                            aTTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCa                                                                    <  2:635041/48‑1 (MQ=255)
                                            aTTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCa                                                                    >  1:635041/1‑48 (MQ=255)
                                                taGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTgg                      >  2:469755/1‑90 (MQ=255)
                                                  gCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCt                    >  1:704288/1‑90 (MQ=255)
                                                    aTTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCtt                   <  1:494630/89‑1 (MQ=255)
                                                    aTTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCtt                   >  2:494630/1‑89 (MQ=255)
                                                       tCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGt               >  1:563333/1‑90 (MQ=255)
                                                               tACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTa          >  1:400003/1‑87 (MQ=255)
                                                               tACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTa          <  2:400003/87‑1 (MQ=255)
                                                                    aTGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTatgat                                              >  2:119023/1‑46 (MQ=255)
                                                                    aTGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTatgat                                              <  1:119023/46‑1 (MQ=255)
                                                                    aTGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGt  >  1:321137/1‑90 (MQ=255)
                                                                    aTGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGt  >  1:511559/1‑90 (MQ=255)
                                                                               |                                                                              
TGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGT  >  NZ_CP009273/3576709‑3576866

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