Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 1,997,112 Δ26,567 bp [tcyN][fliR] 32 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1997112 2023678 26567 57 [0] [0] 57 [tcyN]–[fliR] [tcyN],tcyL,dcyD,tcyJ,fliZ,fliA,fliC,fliD,fliS,fliT,amyA,yedD,yedE,yedF,yedK,yedL,yedN,intG,fliE,fliF,fliG,fliH,fliI,fliJ,fliK,fliL,fliM,fliN,fliO,fliP,fliQ,[fliR]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 19971110 (0.000)56 (0.590) 40/98 0.4 100% coding (704/753 nt) tcyN cystine ABC transporter ATP binding subunit
?NC_000913 2023679 = 0 (0.000)intergenic (+1/‑289) fliR/rcsA flagellar biosynthesis protein FliR/DNA‑binding transcriptional activator RcsA

AATTTTTAATTTGCTGGCTGATATTATTAGTGAATTGCCATTAATATAATTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAAATAATACC  >  NC_000913/2023630‑2023728
                                                 |                                                 
tattgcagcagaaacttctcgaggaactggcgggtgcgaggctgctcggt                                                   >  1:1874194‑M2/50‑50 (MQ=255)
tattgcagcagaaacttctcgaggaactggcgggggcgcggctgctcggt                                                   >  1:4916196‑M2/50‑50 (MQ=255)
  ttgcagcagaaacttctcgaggaactggcgggtgcgaggctgctcggtt                                                  <  1:3557231‑M2/2‑1 (MQ=255)
    gcagcagaaacttctcgaggaactggcgggtgcgaggctgctcggTTcc                                                <  1:4790050‑M2/4‑1 (MQ=255)
    gcagcagaaacttctcgaggaactggcgggggcgaggctgctcggTTCCg                                               >  1:137852‑M2/46‑50 (MQ=255)
     cagcagaaacttctcgaggaactggcgggtgcgaggctgctcggTTCCg                                               >  1:3473847‑M2/45‑49 (MQ=255)
      agcagaaacttctcgaggaactggcgggtgcgagggtgctcggCTCcgc                                              >  1:6585208‑M2/44‑48 (MQ=255)
      agcagaaacttctcgaggaactggcgggggcgaggctgctcggTTCCGt                                              >  1:6591295‑M2/44‑49 (MQ=255)
       gcagaaacttctcgaggaactggcgggtgcgaggctgctcggTTCCCTaa                                            >  1:5265227‑M2/43‑50 (MQ=255)
        cagaaacttctcgaggaactggcgggggcgaggaggctcggTTCCGTAAc                                           >  1:6357113‑M2/42‑50 (MQ=255)
         agaaacttctcgaggaactggcgagtgcgaggctgctcggTTCCGTAACg                                          <  1:5208741‑M2/10‑1 (MQ=255)
          gaaacttctcgaggaactggcgggtgcgaggctgctcggTTCCGTAACGt                                         <  1:939925‑M2/11‑1 (MQ=255)
          gaaacttctcgaggaactggcgggtgcgaggctgctcggGTCCGTAACGt                                         >  1:617599‑M2/40‑50 (MQ=255)
            aacttctcgaggaactggcgggggcgcggctgctcggTTCCGTAacggt                                        >  1:550879‑M2/38‑47 (MQ=255)
             acttctcgaggaactggcgggtgcgaggctgctcggTTCCGTAACGTTTa                                      >  1:6277667‑M2/37‑50 (MQ=255)
                    gaggaactggcgggtgcgaggctgctcggTTCCGTAACGTTTATCATg                                 <  1:857543‑M2/19‑1 (MQ=255)
                     aggaactggcgggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTa                              <  1:4171181‑M2/22‑1 (MQ=255)
                     aggaactggcgggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTa                              >  1:3068265‑M2/29‑50 (MQ=255)
                        aactggcgggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATcc                           >  1:5678403‑M2/26‑50 (MQ=255)
                          ctggcgggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATCCt                          >  1:1292895‑M2/24‑49 (MQ=255)
                           tggcgggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATCCt                          >  1:3379979‑M2/23‑48 (MQ=255)
                             gcgggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAgg                      <  1:4213618‑M2/30‑1 (MQ=255)
                             gcgggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAg                       >  1:6355874‑M2/21‑49 (MQ=255)
                              cgggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGa                     <  1:4804902‑M2/31‑1 (MQ=255)
                               gggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGAt                    >  1:6624652‑M2/19‑50 (MQ=255)
                                ggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGAtt                   <  1:261406‑M2/33‑1 (MQ=255)
                                ggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGAtt                   <  1:2467606‑M2/33‑1 (MQ=255)
                                ggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGAtt                   <  1:7067245‑M2/33‑1 (MQ=255)
                                 gtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTa                  <  1:1025486‑M2/34‑1 (MQ=255)
                                 gtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTa                  <  1:1627071‑M2/34‑1 (MQ=255)
                                 gtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTa                  >  1:2511531‑M2/17‑50 (MQ=255)
                                   gcgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATc                >  1:3355349‑M2/15‑50 (MQ=255)
                                    cgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATcc               <  1:3897550‑M2/37‑1 (MQ=255)
                                    cgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATcc               <  1:6363056‑M2/37‑1 (MQ=255)
                                    cgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATcc               >  1:2979362‑M2/14‑50 (MQ=255)
                                    cgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATcc               >  1:1247318‑M2/14‑50 (MQ=255)
                                    cgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATc                <  1:1080323‑M2/36‑1 (MQ=255)
                                    cgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATc                <  1:3153901‑M2/36‑1 (MQ=255)
                                     gaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCg              <  1:6736034‑M2/38‑1 (MQ=255)
                                       ggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGaa            <  1:5187497‑M2/40‑1 (MQ=255)
                                       ggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGa             >  1:6254426‑M2/11‑49 (MQ=255)
                                        gctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGaaa           <  1:2482068‑M2/41‑1 (MQ=255)
                                        gctgcgcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGaaa           <  1:4231966‑M2/41‑1 (MQ=255)
                                         ctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGaaaa          <  1:2059690‑M2/42‑1 (MQ=255)
                                           gctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGaaaaa         >  1:4853261‑M2/7‑49 (MQ=255)
                                            ctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAaata       <  1:6271239‑M2/45‑1 (MQ=255)
                                               ggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAaataa      <  1:1863146‑M2/46‑1 (MQ=255)
                                               ggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAaataa      <  1:2694196‑M2/46‑1 (MQ=255)
                                               ggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAaataa      <  1:4528702‑M2/46‑1 (MQ=255)
                                               ggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAaataa      <  1:4640569‑M2/46‑1 (MQ=255)
                                               ggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAAataata    <  1:453798‑M2/48‑1 (MQ=255)
                                               ggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAAataata    <  1:5832316‑M2/48‑1 (MQ=255)
                                               ggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAAataata    >  1:530922‑M2/3‑50 (MQ=255)
                                                gTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAAATAATAc   <  1:2580744‑M2/49‑1 (MQ=255)
                                                gTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAAATAATAc   <  1:2563042‑M2/49‑1 (MQ=255)
                                                gTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAAATAATAc   <  1:5132099‑M2/49‑1 (MQ=255)
                                                 ttCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAAATAATAcc  >  1:2506378/1‑50 (MQ=255)
                                                 |                                                 
AATTTTTAATTTGCTGGCTGATATTATTAGTGAATTGCCATTAATATAATTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAAATAATACC  >  NC_000913/2023630‑2023728

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