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

TTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAACGCTTTTGCGGCCCCCTGCTCGACTATCCGCCCCTG  >  NC_000913/1997062‑1997160
                                                 |                                                 
ttattGCAGCAGAAACTTCTCGAGGAACTGGCGGGGGCGAGGCTGCTCgg                                                   >  1:1375257/1‑50 (MQ=255)
 tATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGt                                                  >  1:1874194‑M1/1‑49 (MQ=255)
 tATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGGGCGCGGCTGCTCGGt                                                  >  1:4916196‑M1/1‑49 (MQ=255)
   ttGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGtt                                                 <  1:3557231‑M1/49‑3 (MQ=255)
     gcagcaGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGttcc                                               <  1:4790050‑M1/49‑5 (MQ=255)
     gcagcaGAAACTTCTCGAGGAACTGGCGGGGGCGAGGCTGCTCGGttccg                                              >  1:137852‑M1/1‑45 (MQ=255)
      cagcagAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGttccg                                              >  1:3473847‑M1/1‑44 (MQ=255)
       agcagAAACTTCTCGAGGAACTGGCGGGTGCGAGGGTGCTCGGctccgc                                             >  1:6585208‑M1/1‑43 (MQ=255)
       agcagAAACTTCTCGAGGAACTGGCGGGGGCGAGGCTGCTCGGttccgt                                             >  1:6591295‑M1/1‑43 (MQ=255)
        gcagAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGttccctaa                                           >  1:5265227‑M1/1‑42 (MQ=255)
         cagAAACTTCTCGAGGAACTGGCGGGGGCGAGGAGGCTCGGttccgtaac                                          >  1:6357113‑M1/1‑41 (MQ=255)
          agAAACTTCTCGAGGAACTGGCGAGTGCGAGGCTGCTCGGttccgtaacg                                         <  1:5208741‑M1/50‑11 (MQ=255)
           gAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGttccgtaacgt                                        <  1:939925‑M1/50‑12 (MQ=255)
           gAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGgtccgtaacgt                                        >  1:617599‑M1/1‑39 (MQ=255)
             aaCTTCTCGAGGAACTGGCGGGGGCGCGGCTGCTCGGttccgtaacggt                                       >  1:550879‑M1/1‑37 (MQ=255)
              aCTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGttccgtaacgttta                                     >  1:6277667‑M1/1‑36 (MQ=255)
                     gAGGAACTGGCGGGTGCGAGGCTGCTCGGttccgtaacgtttatcatg                                <  1:857543‑M1/48‑20 (MQ=255)
                      aGGAACTGGCGGGTGCGAGGCTGCTCGGttccgtaacgtttatcatgtta                             >  1:3068265‑M1/1‑28 (MQ=255)
                      aGGAACTGGCGGGTGCGAGGCTGCTCGGttccgtaacgtttatcatgtta                             <  1:4171181‑M1/50‑23 (MQ=255)
                         aaCTGGCGGGTGCGAGGCTGCTCGGttccgtaacgtttatcatgttatcc                          >  1:5678403‑M1/1‑25 (MQ=255)
                           cTGGCGGGTGCGAGGCTGCTCGGttccgtaacgtttatcatgttatcct                         >  1:1292895‑M1/1‑23 (MQ=255)
                            tGGCGGGTGCGAGGCTGCTCGGttccgtaacgtttatcatgttatcct                         >  1:3379979‑M1/1‑22 (MQ=255)
                              gCGGGTGCGAGGCTGCTCGGttccgtaacgtttatcatgttatcctaagg                     <  1:4213618‑M1/50‑31 (MQ=255)
                              gCGGGTGCGAGGCTGCTCGGttccgtaacgtttatcatgttatcctaag                      >  1:6355874‑M1/1‑20 (MQ=255)
                               cGGGTGCGAGGCTGCTCGGttccgtaacgtttatcatgttatcctaagga                    <  1:4804902‑M1/50‑32 (MQ=255)
                                gggTGCGAGGCTGCTCGGttccgtaacgtttatcatgttatcctaaggat                   >  1:6624652‑M1/1‑18 (MQ=255)
                                 ggTGCGAGGCTGCTCGGttccgtaacgtttatcatgttatcctaaggatt                  <  1:261406‑M1/50‑34 (MQ=255)
                                 ggTGCGAGGCTGCTCGGttccgtaacgtttatcatgttatcctaaggatt                  <  1:7067245‑M1/50‑34 (MQ=255)
                                 ggTGCGAGGCTGCTCGGttccgtaacgtttatcatgttatcctaaggatt                  <  1:2467606‑M1/50‑34 (MQ=255)
                                  gTGCGAGGCTGCTCGGttccgtaacgtttatcatgttatcctaaggatta                 <  1:1025486‑M1/50‑35 (MQ=255)
                                  gTGCGAGGCTGCTCGGttccgtaacgtttatcatgttatcctaaggatta                 >  1:2511531‑M1/1‑16 (MQ=255)
                                  gTGCGAGGCTGCTCGGttccgtaacgtttatcatgttatcctaaggatta                 <  1:1627071‑M1/50‑35 (MQ=255)
                                    gCGAGGCTGCTCGGttccgtaacgtttatcatgttatcctaaggattatc               >  1:3355349‑M1/1‑14 (MQ=255)
                                     cGAGGCTGCTCGGttccgtaacgtttatcatgttatcctaaggattatcc              <  1:3897550‑M1/50‑38 (MQ=255)
                                     cGAGGCTGCTCGGttccgtaacgtttatcatgttatcctaaggattatcc              <  1:6363056‑M1/50‑38 (MQ=255)
                                     cGAGGCTGCTCGGttccgtaacgtttatcatgttatcctaaggattatcc              >  1:2979362‑M1/1‑13 (MQ=255)
                                     cGAGGCTGCTCGGttccgtaacgtttatcatgttatcctaaggattatcc              >  1:1247318‑M1/1‑13 (MQ=255)
                                     cGAGGCTGCTCGGttccgtaacgtttatcatgttatcctaaggattatc               <  1:1080323‑M1/49‑37 (MQ=255)
                                     cGAGGCTGCTCGGttccgtaacgtttatcatgttatcctaaggattatc               <  1:3153901‑M1/49‑37 (MQ=255)
                                      gAGGCTGCTCGGttccgtaacgtttatcatgttatcctaaggattatccg             <  1:6736034‑M1/50‑39 (MQ=255)
                                        ggCTGCTCGGttccgtaacgtttatcatgttatcctaaggattatccgaa           <  1:5187497‑M1/50‑41 (MQ=255)
                                        ggCTGCTCGGttccgtaacgtttatcatgttatcctaaggattatccga            >  1:6254426‑M1/1‑10 (MQ=255)
                                         gctgctCGGttccgtaacgtttatcatgttatcctaaggattatccgaaa          <  1:2482068‑M1/50‑42 (MQ=255)
                                         gctgcgCGGttccgtaacgtttatcatgttatcctaaggattatccgaaa          <  1:4231966‑M1/50‑42 (MQ=255)
                                          ctgctCGGttccgtaacgtttatcatgttatcctaaggattatccgaaaa         <  1:2059690‑M1/50‑43 (MQ=255)
                                            gctCGGttccgtaacgtttatcatgttatcctaaggattatccgaaaaa        >  1:4853261‑M1/1‑6 (MQ=255)
                                             ctCGGttccgtaacgtttatcatgttatcctaaggattatccgaaaaata      <  1:6271239‑M1/50‑46 (MQ=255)
                                                ggttccgtaacgtttatcatgttatcctaaggattatccgaaaaataata   >  1:530922‑M1/1‑2 (MQ=255)
                                                ggttccgtaacgtttatcatgttatcctaaggattatccgaaaaataata   <  1:5832316‑M1/50‑49 (MQ=255)
                                                ggttccgtaacgtttatcatgttatcctaaggattatccgaaaaataata   <  1:453798‑M1/50‑49 (MQ=255)
                                                ggttccgtaacgtttatcatgttatcctaaggattatccgaaaaataa     <  1:1863146‑M1/48‑47 (MQ=255)
                                                ggttccgtaacgtttatcatgttatcctaaggattatccgaaaaataa     <  1:2694196‑M1/48‑47 (MQ=255)
                                                ggttccgtaacgtttatcatgttatcctaaggattatccgaaaaataa     <  1:4528702‑M1/48‑47 (MQ=255)
                                                ggttccgtaacgtttatcatgttatcctaaggattatccgaaaaataa     <  1:4640569‑M1/48‑47 (MQ=255)
                                                 gttccgtaacgtttatcatgttatcctaaggattatccgaaaaataatac  <  1:2580744‑M1/50‑50 (MQ=255)
                                                 gttccgtaacgtttatcatgttatcctaaggattatccgaaaaataatac  <  1:2563042‑M1/50‑50 (MQ=255)
                                                 gttccgtaacgtttatcatgttatcctaaggattatccgaaaaataatac  <  1:5132099‑M1/50‑50 (MQ=255)
                                                 |                                                 
TTATTGCAGCAGAAACTTCTCGAGGAACTGGCGGGTGCGAGGCTGCTCGGGGTCGGCAAATAACGCTTTTGCGGCCCCCTGCTCGACTATCCGCCCCTG  >  NC_000913/1997062‑1997160

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.