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 26 [1] [0] 26 [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 = 19971111 (0.030)25 (0.680) 19/96 0.4 98.0% 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

TTTTAATTTGCTGGCTGATATTATTAGTGAATTGCCATTAATATAATTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAAATAATAC  >  NC_000913/2023633‑2023727
                                              |                                                
tgcagcagaaacttctcgaggaactggcgggggcgaggctgctcggTTc                                                >  1:2045954‑M2/47‑49 (MQ=255)
 gcagcagaaacttctcgaggaactggcgggtgcgaggctgctcggTTc                                                <  1:130777‑M2/3‑1 (MQ=255)
 gcagcagaaacttctcgaggaactggcgggtgcgaggctgctcggTTc                                                <  1:100890‑M2/3‑1 (MQ=255)
   agcagaaacttctcgaggaactggcgggtgcgaggctgctcggTTCCGt                                             <  1:2031520‑M2/6‑1 (MQ=255)
    gcagaaacttctcgaggaactggcgggtgcgagggtgctcggTTCCGt                                             >  1:447242‑M2/43‑48 (MQ=255)
    gcagaaacttctcgaggaactggcgggtgcgaggctgctcggTTCCGTa                                            <  1:97744‑M2/7‑1 (MQ=255)
    gcagaaacttctcgaggaactggcgggtgcgaggctgctcggTTCCGTa                                            <  1:2041983‑M2/7‑1 (MQ=255)
      agaaacttctcgaggaactggcgggtgcgaggctgctcggTTCCGTa                                            >  1:834090‑M2/41‑47 (MQ=255)
          acttctcgaggaactggcgggtgcgaggctgctcggTTCCGTAACGttt                                      <  1:1224542‑M2/13‑1 (MQ=255)
              ctcgaggaactggcgggtgcgaggctgctcggTTCCGTAACGTTTATc                                   <  1:2410159‑M2/16‑1 (MQ=255)
                 gaggaactggcgggtgcgaggctgctcggTTCCGTAACGTTTATCATg                                <  1:391405‑M2/19‑1 (MQ=255)
                 gaggaactggcgggtgcgaggctgctcggTTCCGTAACGTTTATCATGt                               <  1:986454‑M2/20‑1 (MQ=255)
                         ggcgggtgcgaggctgctcggTTCCGTAACGTTTATCATGTTATCCt                         >  1:2480730‑M2/22‑47 (MQ=255)
                                 cgaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATc               <  1:2200547‑M2/36‑1 (MQ=255)
                                  gaggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATcc              >  1:262049‑M2/13‑49 (MQ=255)
                                    ggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCg             >  1:1539744‑M2/11‑48 (MQ=255)
                                    ggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCg             >  1:1440504‑M2/11‑48 (MQ=255)
                                    ggctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCg             >  1:1012655‑M2/11‑48 (MQ=255)
                                     gctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGaa           <  1:813767‑M2/40‑1 (MQ=255)
                                      ctgctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGaaa          <  1:2610603‑M2/41‑1 (MQ=255)
                                         ctcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGaaaaa        >  1:2698254‑M2/6‑48 (MQ=255)
                                          tcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAaat       >  1:782529‑M2/5‑48 (MQ=255)
                                          tcggTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAaat       >  1:1330469‑M2/5‑48 (MQ=255)
                                             gTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAAataata   <  1:2204270‑M2/48‑1 (MQ=255)
                                             gTTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAAataata   <  1:863449‑M2/48‑1 (MQ=255)
                                              ttCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAAATAATAc  >  1:373941/1‑49 (MQ=255)
                                              |                                                
TTTTAATTTGCTGGCTGATATTATTAGTGAATTGCCATTAATATAATTCCGTAACGTTTATCATGTTATCCTAAGGATTATCCGAAAAATAATAC  >  NC_000913/2023633‑2023727

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