Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 1,962,566 Δ16,764 bp flhEinsA‑5 17 genes

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1962566 1979329 16764 75 [0] [0] 73 flhE–insA‑5 flhE,flhA,flhB,cheZ,cheY,cheB,cheR,tap,tar,cheW,cheA,motB,motA,flhC,flhD,insB‑5,insA‑5

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 19625650 (0.000)73 (0.900) 49/96 0.0 100% intergenic (+105/+15) yecT/flhE protein YecT/flagellar protein
?NC_000913 1979330 = 0 (0.000)intergenic (‑115/‑423) insA‑5/uspC IS1 protein InsA/universal stress protein C

TAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATG  >  NC_000913/1979281‑1979377
                                                 |                                               
agacatcgagcctttttcactgagttattaaacatactcgcgagcgcgtc                                                 <  1:4518060‑M2/1‑1 (MQ=255)
agacatcgagcctttttcactgagttattaaacatactcgcgagcgcgtc                                                 <  1:1324216‑M2/1‑1 (MQ=255)
 gacatcgagcctttttcactgagttattaaacatactcgcgagcgcgtcg                                                >  1:737662‑M2/49‑50 (MQ=255)
 gacatcgagcctttttcactgagttattaaacatactcgcgagcgcgtcg                                                >  1:2988010‑M2/49‑50 (MQ=255)
  acatcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgc                                               >  1:259655‑M2/48‑50 (MQ=255)
    atcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgcta                                             >  1:385757‑M2/46‑50 (MQ=255)
     tcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgctac                                            <  1:3443984‑M2/6‑1 (MQ=255)
     tcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgctac                                            <  1:3873592‑M2/6‑1 (MQ=255)
     tcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgctac                                            <  1:5419636‑M2/6‑1 (MQ=255)
     tcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgctac                                            <  1:5631150‑M2/6‑1 (MQ=255)
      cgagcctttttcactgagttattaaacatactcgcgagcgcgtcgctacg                                           >  1:6259723‑M2/44‑50 (MQ=255)
       gagcctttttcactgagttattaaacatactcgcgagcgcgtcgctacgc                                          >  1:3984732‑M2/43‑50 (MQ=255)
       gagcctttttcactgagttattaaacatactcgcgagcgcgtcgctacg                                           <  1:754565‑M2/7‑1 (MQ=255)
       gagcctttttcactgagttattaaacatactcgcgagcgcgtcgctacg                                           <  1:396008‑M2/7‑1 (MQ=255)
          cctttttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTag                                       <  1:4692704‑M2/11‑1 (MQ=255)
          cctttttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTag                                       >  1:5135707‑M2/40‑50 (MQ=255)
           ctttttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTag                                       >  1:5007502‑M2/39‑49 (MQ=255)
           ctttttcactgagttattaaacatactcgcaagcgcgtCGCTACGCTagc                                      <  1:5741273‑M2/12‑1 (MQ=255)
            tttttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagc                                      >  1:3510496‑M2/38‑49 (MQ=255)
            tttttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagc                                      >  1:2451821‑M2/38‑49 (MQ=255)
             ttttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagc                                      >  1:2831598‑M2/37‑48 (MQ=255)
              tttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagca                                     <  1:2524456‑M2/13‑1 (MQ=255)
               ttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaag                                   <  1:478103‑M2/15‑1 (MQ=255)
                tcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaagc                                  <  1:921623‑M2/16‑1 (MQ=255)
                tcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaagc                                  <  1:2102022‑M2/16‑1 (MQ=255)
                tcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaagc                                  >  1:6115486‑M2/34‑49 (MQ=255)
                tcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaagc                                  <  1:831489‑M2/16‑1 (MQ=255)
                 cactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaagca                                 >  1:4180287‑M2/33‑49 (MQ=255)
                  actgagttattaaacatactcgcgagcgcgtCGCTACGCTAgcaagcaa                                <  1:5449706‑M2/18‑1 (MQ=255)
                  actgagttattaaacatactcgcgagcgcgtCGCTACGCTAgcaagcaa                                >  1:5519661‑M2/32‑49 (MQ=255)
                      agttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCaaaaa                             <  1:1047560‑M2/21‑1 (MQ=255)
                      agttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAAt                            <  1:6244089‑M2/22‑1 (MQ=255)
                      agttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAAt                            >  1:2748690‑M2/28‑49 (MQ=255)
                       gttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGa                          <  1:2601846‑M2/24‑1 (MQ=255)
                          attaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGaaa                        >  1:3424442‑M2/24‑49 (MQ=255)
                            taaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAg                     >  1:327757‑M2/22‑50 (MQ=255)
                             aaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAg                     <  1:3368808‑M2/29‑1 (MQ=255)
                             aaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAg                     <  1:1794548‑M2/29‑1 (MQ=255)
                             aaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAg                     >  1:2685401‑M2/21‑49 (MQ=255)
                             aaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAg                     >  1:4504479‑M2/21‑49 (MQ=255)
                              aacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAgg                    >  1:5439760‑M2/20‑49 (MQ=255)
                              aacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAgg                    >  1:1002638‑M2/20‑49 (MQ=255)
                              aacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAgg                    >  1:2875921‑M2/20‑49 (MQ=255)
                                catactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaat                 <  1:713198‑M2/33‑1 (MQ=255)
                                 atactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaata                <  1:1483776‑M2/34‑1 (MQ=255)
                                 atactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaat                 >  1:703491‑M2/17‑49 (MQ=255)
                                  tactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaataa               <  1:4389113‑M2/35‑1 (MQ=255)
                                  tactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaataa               >  1:6244730‑M2/16‑50 (MQ=255)
                                    ctcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc             >  1:1568691‑M2/14‑50 (MQ=255)
                                    ctcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc             <  1:3645775‑M2/37‑1 (MQ=255)
                                    ctcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc             <  1:6136001‑M2/37‑1 (MQ=255)
                                      cgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCg            >  1:5913051‑M2/12‑49 (MQ=255)
                                       gcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCg            <  1:6261571‑M2/38‑1 (MQ=255)
                                        cgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa         <  1:3425098‑M2/41‑1 (MQ=255)
                                        cgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa         <  1:2826162‑M2/41‑1 (MQ=255)
                                         gagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt        >  1:4618023‑M2/9‑50 (MQ=255)
                                         gagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt        >  1:1517994‑M2/9‑50 (MQ=255)
                                         gagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt        >  1:743828‑M2/9‑50 (MQ=255)
                                         gagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt        <  1:650730‑M2/42‑1 (MQ=255)
                                         gagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt        >  1:6129865‑M2/9‑50 (MQ=255)
                                         gagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt        <  1:4169262‑M2/42‑1 (MQ=255)
                                         gagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt        >  1:4832081‑M2/9‑50 (MQ=255)
                                           gcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATgg      <  1:4210496‑M2/44‑1 (MQ=255)
                                           gcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATgg      >  1:3252690‑M2/7‑50 (MQ=255)
                                            cgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATggg     <  1:3022594‑M2/45‑1 (MQ=255)
                                             gcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGa    >  1:4143161‑M2/5‑50 (MQ=255)
                                              cgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt   >  1:1966838‑M2/4‑50 (MQ=255)
                                              cgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt   <  1:2461043‑M2/47‑1 (MQ=255)
                                              cgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt   >  1:2489798‑M2/4‑50 (MQ=255)
                                              cgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt   <  1:4667425‑M2/47‑1 (MQ=255)
                                               gtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATg  >  1:3301030‑M2/3‑50 (MQ=255)
                                               gtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATg  >  1:3283655‑M2/3‑50 (MQ=255)
                                               gtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATg  >  1:4521135‑M2/3‑50 (MQ=255)
                                                 |                                               
TAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATG  >  NC_000913/1979281‑1979377

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