New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 19971119 (0.080)59 (0.540) 42/98 1.1 84.9% coding (704/753 nt) tcyN cystine ABC transporter ATP binding subunit
?NC_000913 2023679 = 12 (0.110)intergenic (+1/‑289) fliR/rcsA flagellar biosynthesis protein FliR/DNA‑binding transcriptional activator RcsA

AATTTTTAATTTGCTGGCTGATATTATTAGTGAATTGCCATTAATATAATTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAAATAATAC  >  NC_000913/2023630‑2023727
                                                 |                                                
tattgcagcagaaacttctcgaggaactggcgggggcgaggccgctcggt                                                  >  1:4868775‑M2/50‑50 (MQ=255)
aaTTTTTAATTTGCTGGCTGATATTATTAGTGAATTGCAATTAATATAAt                                                  <  1:4603105/50‑1 (MQ=255)
 attgcagcagaaacttctcgaggaactggcgggggcgaggctgctcggtt                                                 >  1:4378328‑M2/49‑50 (MQ=255)
  tttttAATTTGCTGGCTGATATTATTAGTGAATTGCCATTAATATAATTc                                                <  1:1033311/50‑1 (MQ=255)
    tttAATTTGCTGGCTGATATTATTAGTGAATTGCCATTAATATAATTcc                                               <  1:3086173/49‑1 (MQ=255)
    gcagcagaaacttctcgaggaactggcgggtgcgaggctgctcggTTcc                                               <  1:1120365‑M2/4‑1 (MQ=255)
    gcagcagaaacttctcgaggaactggcgggggcgcggctgctcggTTCCg                                              >  2:3953197‑M2/46‑50 (MQ=255)
     cagcagaaacttctcgaggaactggcgggtgcgaggctgctcggTTCCg                                              <  2:2344993‑M2/5‑1 (MQ=255)
     cagcagaaacttctcgaggaactggcgggtgcgaggctgctcggTTCCGt                                             <  2:1876850‑M2/6‑1 (MQ=255)
       gcagaaacttctcgaggaactggcgggtgcgaggctggtcggTTCCGTa                                            >  2:3277544‑M2/43‑49 (MQ=255)
       aaTTTGCTGGCTGATATTATTAGTGAATTGCCATTAATATAATTCCGTaa                                           >  2:414344/1‑50 (MQ=255)
        cagaaacttctcgaggaactggcgggtgcgaggctgctcggTTCCGTAAc                                          <  1:3562114‑M2/9‑1 (MQ=255)
        cagaaacttctcgaggaactggcgggtgcgaggctgctcggTTCCGGAAc                                          >  2:490030‑M2/42‑50 (MQ=255)
        cagaaacttctcgaggaactggcgggtgcgaggctgctcggGTCCGTAAc                                          >  1:2516190‑M2/42‑50 (MQ=255)
        aTTTGCTGGCTGATATTATTAGTGAATTGCCATTAATATAATTCCGTAAc                                          <  2:1430338/50‑1 (MQ=255)
        aTTTGCTGGCTGATATTATTAGTGAATTGCCATTAATATAATTCCGTAAc                                          <  1:76113/50‑1 (MQ=255)
         agaaacttctcgaggaactggcgggggggaggctgctcggTTCCGTAACg                                         >  2:4988417‑M2/41‑50 (MQ=255)
          gaaacttctcgaggaactggcgggtgcgagggtgctcggTTCCGTAACGt                                        >  2:1735932‑M2/40‑50 (MQ=255)
           aaacttctcgaggaactggcgggtgcgaggctgcccggTTCCGTAACGtt                                       >  2:957217‑M2/39‑50 (MQ=255)
           aaacttctcgaggaactggcgggggggaggctgctcggTTCCGTAACGtt                                       >  2:2356473‑M2/39‑50 (MQ=255)
             acttctcgaggaactggcgggtgcgaggctgctcggTTCCGTAACGTTTa                                     <  1:1165991‑M2/14‑1 (MQ=255)
             acttctcgaggaactggcgggggcgaggctgctcggTTCCGTAACGTTTa                                     >  2:3422014‑M2/37‑50 (MQ=255)
              cttctcgaggaactggcgggggcggggctgctcggTTCCGTAACGTTTAt                                    >  2:481654‑M2/36‑50 (MQ=255)
                tctcgaggaactggcgggtgcgaggctgctcggTTCCGTAACGTTTATCa                                  >  1:4198672‑M2/34‑50 (MQ=255)
                tctcgaggaactggcgggagcgaggctgctcggTTCCCCAACGCTtattt                                  >  2:4714175‑M2/34‑48 (MQ=255)
                  tcgaggaactggcgggttcggggctgctcggGTCCGTAACGCTTATCATg                                >  2:4020820‑M2/32‑50 (MQ=255)
                  tcgaggaactggcgggggcgaggctgctcggTTCCGTAACGTTTATCATg                                >  2:1915767‑M2/32‑50 (MQ=255)
                   cgaggaactggcgggtgcgaggctgctcggGTCCGTAACGTTTATCATGt                               >  2:2329171‑M2/31‑50 (MQ=255)
                   cgaggaactggcgggggcgaggcttctcggTTCCGTAACGTTTTTCATGt                               >  2:3555821‑M2/31‑50 (MQ=255)
                    gaggaactggcgggtgcgaggctgctcggTTCCGTAACGTTTATCATGt                               <  2:543423‑M2/20‑1 (MQ=255)
                    gaggaactggcgggtgcgaggctgctcggTTCCGTAACGTTTATCATGt                               <  2:2677362‑M2/20‑1 (MQ=255)
                    gaggaactggcgggggggaggctgctccgTTCCGTAACGTTTATCATGt                               >  2:4782881‑M2/30‑49 (MQ=255)
                     aggaactggcgggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTa                             <  2:1337446‑M2/22‑1 (MQ=255)
                      ggaactggcgggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTAt                            <  1:2373236‑M2/23‑1 (MQ=255)
                       gaactggcgggggcgaggctgctcggGTCCGTAACGTTTATCATGTTATc                           >  2:3323437‑M2/27‑50 (MQ=255)
                        aactggcgggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATcc                          >  2:2955670‑M2/26‑50 (MQ=255)
                        aactggcgggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATcc                          >  1:4391803‑M2/26‑50 (MQ=255)
                         actggcgggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATCCt                         <  1:1699267‑M2/26‑1 (MQ=255)
                          ttaGTGAATTGCCATTAATATAATTCCGTAACGTTTATCATGTTATccg                         <  2:611187/49‑2 (MQ=255)
                          ttaGTGAATTGCCATTAATATAATTCCGTAACGTTTATCATGTTATCCt                         >  2:232331/1‑49 (MQ=255)
                          ctggcgggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATCCt                         >  1:1372068‑M2/24‑49 (MQ=255)
                           tggcgggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATCCt                         >  1:3927244‑M2/23‑48 (MQ=255)
                           tggcgggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATCCt                         >  2:337826‑M2/23‑48 (MQ=255)
                           tggcgggtgcgaggctgctcggATCCGTAACGTTTATCATGTTATCCt                         >  2:1040042‑M2/23‑48 (MQ=255)
                             gcgggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTaa                       <  2:1369367‑M2/28‑1 (MQ=255)
                             gTGAATTGCCATTAATATAATTCCGTAACGTTTATCATGTTATCCTAAgg                     <  1:3735336/50‑1 (MQ=255)
                               gggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGAt                   <  2:4460054‑M2/32‑1 (MQ=255)
                               gggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGAt                   <  1:3899273‑M2/32‑1 (MQ=255)
                               gggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGAt                   >  2:3510235‑M2/19‑50 (MQ=255)
                               gggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGAt                   <  1:4925827‑M2/32‑1 (MQ=255)
                               gggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGAt                   <  1:1691739‑M2/32‑1 (MQ=255)
                                 gtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTa                 <  2:228455‑M2/34‑1 (MQ=255)
                                    cgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATcc              <  2:1054389‑M2/37‑1 (MQ=255)
                                     gaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATcc              >  1:2747181‑M2/13‑49 (MQ=255)
                                     gaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATcc              >  2:3791152‑M2/13‑49 (MQ=255)
                                     gaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATcc              >  1:3038053‑M2/13‑49 (MQ=255)
                                      cATTAATATAATTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCg             >  1:2804594/1‑49 (MQ=255)
                                        gctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGaaa          <  1:4659796‑M2/41‑1 (MQ=255)
                                         ctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGaaaa         <  1:4761245‑M2/42‑1 (MQ=255)
                                         ctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGaaa          >  1:1030956‑M2/9‑49 (MQ=255)
                                         ctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGaaa          >  1:2410415‑M2/9‑49 (MQ=255)
                                          tgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGaaaaa        >  2:4847669‑M2/8‑50 (MQ=255)
                                             tcggTTCCGTAACGCTTATCATGTTATCCTAAGGATTATCCGAAAaataa     >  2:2649863‑M2/5‑50 (MQ=255)
                                             ataATTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAaataa     >  1:2906063/1‑50 (MQ=255)
                                              cggTTCCGTAACGTTTATCATGTTATCCTAAGGATTCTCCGAAAaataa     >  2:2852622‑M2/4‑49 (MQ=255)
                                               ggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAAataata   >  2:2122041‑M2/3‑50 (MQ=255)
                                               ggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAAataata   >  1:4402284‑M2/3‑50 (MQ=255)
                                                gTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAAATAATAc  <  2:3865821‑M2/49‑1 (MQ=255)
                                                gTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAAATAATAc  >  1:651372‑M2/2‑50 (MQ=255)
                                                gTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAAATAATAc  >  1:4099642‑M2/2‑50 (MQ=255)
                                                aTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAAataata   <  1:2497001/49‑1 (MQ=255)
                                                 |                                                
AATTTTTAATTTGCTGGCTGATATTATTAGTGAATTGCCATTAATATAATTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAAATAATAC  >  NC_000913/2023630‑2023727

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