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 58 [1] [0] 55 [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.010)54 (0.740) 42/98 0.1 99.1% 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

TTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAACGCTTTTGCGGCCCCCTGCTCGACTATCCGCCCCT  >  NC_000913/1997062‑1997159
                                                 |                                                
ttattGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTTGTCgg                                                  >  1:2019564/1‑50 (MQ=255)
ttattGCAGCAGAAACTTCTCGAGGAACTGGCGGGGGCGAGGCTGCTCgg                                                  >  1:581546/1‑50 (MQ=255)
 tattGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCggg                                                 >  1:4863904/1‑50 (MQ=255)
 tattGCAGCAGAAACTTCTCGAGGAACTGGCGGGCGCCGGGCGGCTCggc                                                 >  1:955581/1‑49 (MQ=255)
 tATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGt                                                 >  1:1249116‑M1/1‑49 (MQ=255)
     gcagcaGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGttcc                                              <  1:4496105‑M1/49‑5 (MQ=255)
       agcagAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGttccgta                                           <  1:4868304‑M1/50‑8 (MQ=255)
       agcagAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGttccgta                                           <  1:4005632‑M1/50‑8 (MQ=255)
        gcagAAACTTCTCGAGGAACTGGCGGGTGGGAGGCTGCTCGGttccgtaa                                          <  1:232391‑M1/50‑9 (MQ=255)
        gcagAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGttccgtaa                                          <  1:2906211‑M1/50‑9 (MQ=255)
         cagAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGttccgtaac                                         <  1:3595952‑M1/50‑10 (MQ=255)
          agAAACTTCTCGAGGAACTGGCGGGTGCGCGGCTGCTCGGgtccgtgacg                                        >  1:722178‑M1/1‑40 (MQ=255)
          agAAACTTCTCGAGGAACTGGCGGGTGCGAGGGTGCTCGGttccgtaacg                                        >  1:1557972‑M1/1‑40 (MQ=255)
          agAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGttccgtaacg                                        <  1:1357350‑M1/50‑11 (MQ=255)
          agAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGgtccgtaaag                                        >  1:3998804‑M1/1‑40 (MQ=255)
          agAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCCGttccgtaacg                                        >  1:1187171‑M1/1‑40 (MQ=255)
           gAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGgtccctaacgt                                       >  1:4178455‑M1/1‑39 (MQ=255)
            aaaCTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGgtccgtaacgtt                                      >  1:3873636‑M1/1‑38 (MQ=255)
            aaaCTTCTCGAGGAACTGGCGGGGGCGAGGCGGCTCGGgtccgtaaagtt                                      >  1:410421‑M1/1‑38 (MQ=255)
             aaCTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGttccgtaacgttt                                     <  1:2679507‑M1/50‑14 (MQ=255)
               cTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGttccgtaacgtttat                                   <  1:5115461‑M1/50‑16 (MQ=255)
                  ctcGAGGAACTGGCGGGTGCGAGGCTGCTCGGttccgtaacgtttatca                                 <  1:3328250‑M1/49‑18 (MQ=255)
                    cGAGGAACTGGCGGGTGCGAGGCTGCTCGGttccgtaacgtttatcatgt                              >  1:23126‑M1/1‑30 (MQ=255)
                     gAGGAACTGGCGGGTGCGAGGCTGCTCGGttccgtaacgtttatcatgt                              >  1:2949172‑M1/1‑29 (MQ=255)
                     gAGGAACTGGCGGGTGCGAGGCTGCTCGGttccgtaacgtttatcatg                               <  1:4422267‑M1/48‑20 (MQ=255)
                     gAGGAACTGGCGGGTGCGAGGCTGCTCGGttccgtaacgtttatcatg                               <  1:5339828‑M1/48‑20 (MQ=255)
                     gAGGAACTGGCGGGTGCGAGGCTGCTCGGttccgtaacgtttatcatg                               <  1:5521473‑M1/48‑20 (MQ=255)
                     gAGGAACTGGCGGGTGCGAGGCTGCTCGGttccgtaacgtttatcat                                <  1:1484753‑M1/47‑19 (MQ=255)
                       ggAACTGGCGGGTGCGAGGCTGCTCGGttccgtaacgtttatcatgttat                           <  1:5142854‑M1/50‑24 (MQ=255)
                         aaCTGGCGGGTGCGAGGCTGCTCGGttccgtaacgtttatcatgttatcc                         >  1:2512364‑M1/1‑25 (MQ=255)
                           cTGGCGGGTGCGAGGCTGCTCGGttccgtaacgtttatcatgttatccta                       <  1:2813595‑M1/50‑28 (MQ=255)
                           cTGGCGGGTGCGAGGCTGCTCGGttccgtaacgtttatcatgttatcct                        >  1:2133347‑M1/1‑23 (MQ=255)
                            tGGCGGGTGCGAGGCTGCTCGGttccgtaacgtttatcatgttatcct                        >  1:4071887‑M1/1‑22 (MQ=255)
                              gCGGGTGCGAGGCTGCTCGGttccgtaacgtttatcatgttatcctaagg                    <  1:4247962‑M1/50‑31 (MQ=255)
                              gCGGGTGCGAGGCTGCTCGGttccgtaacgtttatcatgttatcctaag                     <  1:3940138‑M1/49‑30 (MQ=255)
                               cGGGTGCGAGGCTGCTCGGttccgtaacgtttatcatgttatcctaagga                   <  1:4475866‑M1/50‑32 (MQ=255)
                                gggTGCGAGGCTGCTCGGttccgtaacgtttatcatgttatcctaaggat                  >  1:1087888‑M1/1‑18 (MQ=255)
                                gggTGCGAGGCTGCTCGGttccgtaacgtttatcatgttatcctaaggat                  <  1:2779504‑M1/50‑33 (MQ=255)
                                  gTGCGAGGCTGCTCGGttccgtaacgtttatcatgttatcctaaggatta                <  1:3199117‑M1/50‑35 (MQ=255)
                                   tGCGAGGCTGCTCGGttccgtaacgtttatcatgttatcctaaggattat               <  1:2474042‑M1/50‑36 (MQ=255)
                                     cGAGGCTGCTCGGttccgtaacgtttatcatgttatcctaaggattatcc             <  1:1076298‑M1/50‑38 (MQ=255)
                                     cGAGGCTGCTCGGttccgtaacgtttatcatgttatcctaaggattatc              <  1:4037348‑M1/49‑37 (MQ=255)
                                      gAGGCTGCTCGGttccgtaacgtttatcatgttatcctaaggattatccg            <  1:4915654‑M1/50‑39 (MQ=255)
                                      gAGGCTGCTCGGttccgtaacgtttatcatgttatcctaaggattatcc             >  1:4745156‑M1/1‑12 (MQ=255)
                                      gAGGCTGCTCGGttccgtaacgtttatcatgttatcctaaggattatcc             >  1:1505759‑M1/1‑12 (MQ=255)
                                       aGGCTGCTCGGttccgtaacgtttatcatgttatcctaaggattatccga           <  1:2753517‑M1/50‑40 (MQ=255)
                                       aGGCTGCTCGGttccgtaacgtttatcatgttatcctaaggattatccga           <  1:197094‑M1/50‑40 (MQ=255)
                                       aGGCTGCTCGGttccgtaacgtttatcatgttatcctaaggattatccg            >  1:2801963‑M1/1‑11 (MQ=255)
                                         gctgctCGGttccgtaacgtttatcatgttatcctaaggattatccgaaa         <  1:3005669‑M1/50‑42 (MQ=255)
                                         gctgctCGGttccgtaacgtttatcatgttatcctaaggattatccgaa          >  1:1938313‑M1/1‑9 (MQ=255)
                                           tgctCGGttccgtaacgtttatcatgttatcctaaggattatccgaaaaa       >  1:3209030‑M1/1‑7 (MQ=255)
                                             ctCGGttccgtaacgtttatcatgttatcctaaggattatccgaaaaata     <  1:4739306‑M1/50‑46 (MQ=255)
                                             ctCGGttccgtaacgtttatcatgttatcctaaggattatccgaaaaat      <  1:1388619‑M1/49‑45 (MQ=255)
                                             ctCGGttccgtaacgtttatcatgttatcctaaggattatccgaaaaat      <  1:594915‑M1/49‑45 (MQ=255)
                                             ctCGGttccgtaacgtttatcatgttatcctaaggattatccgaaaaa       <  1:5365141‑M1/48‑44 (MQ=255)
                                                ggttccgtaacgtttatcatgttatcctaaggattatccgaaaaataata  <  1:1158961‑M1/50‑49 (MQ=255)
                                                ggttccgtaacgtttatcatgttatcctaaggattatccgaaaaataa    <  1:2102798‑M1/48‑47 (MQ=255)
                                                ggttccgtaacgtttatcatgttatcctaaggattatccgaaaaataa    <  1:2695876‑M1/48‑47 (MQ=255)
                                                 |                                                
TTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAACGCTTTTGCGGCCCCCTGCTCGACTATCCGCCCCT  >  NC_000913/1997062‑1997159

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