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 45 [0] [0] 46 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)45 (1.240) 33/94 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

GTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGT  >  NC_000913/1979284‑1979378
                                              |                                                
catcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgc                                                >  1:2007766‑M2/47‑49 (MQ=255)
catcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgc                                                >  1:1574973‑M2/47‑49 (MQ=255)
    gagcctttttcactgagttattaaacatactcgcgagcgcgtcgctacg                                            <  1:1036258‑M2/7‑1 (MQ=255)
     agcctttttcactgagttattaaacatactcgcgagcgcgtcgctacgc                                           >  1:2055072‑M2/42‑49 (MQ=255)
      gcctttttcactgagttattaaacatactcgcgagcgcgtcgctacgct                                          <  1:888524‑M2/9‑1 (MQ=255)
       cctttttcactgagttattaaacatactcgcgagcgcgtcgctacgcta                                         <  1:2557724‑M2/10‑1 (MQ=255)
        ctttttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTag                                        >  1:2569695‑M2/39‑49 (MQ=255)
        ctttttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTag                                        <  1:502236‑M2/11‑1 (MQ=255)
        ctttttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTag                                        >  1:2429256‑M2/39‑49 (MQ=255)
        ctttttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTag                                        >  1:92742‑M2/39‑49 (MQ=255)
        ctttttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTag                                        >  1:1566259‑M2/39‑49 (MQ=255)
           tttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagca                                      <  1:2412754‑M2/13‑1 (MQ=255)
            ttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaa                                     <  1:1854444‑M2/14‑1 (MQ=255)
              cactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaagc                                   <  1:2587294‑M2/16‑1 (MQ=255)
              cactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaagc                                   <  1:1799026‑M2/16‑1 (MQ=255)
               actgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaagca                                  >  1:470200‑M2/32‑48 (MQ=255)
                 tgagttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCaaa                                >  1:657173‑M2/30‑48 (MQ=255)
                   agttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCaaaa                               <  1:1736378‑M2/20‑1 (MQ=255)
                    gttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAAt                             >  1:1702803‑M2/27‑48 (MQ=255)
                    gttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAAt                             >  1:764284‑M2/27‑48 (MQ=255)
                    gttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATg                            <  1:2506210‑M2/23‑1 (MQ=255)
                          aaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAg                      >  1:2013152‑M2/21‑49 (MQ=255)
                          aaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAg                      >  1:829997‑M2/21‑49 (MQ=255)
                           aacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAgg                     <  1:1492177‑M2/30‑1 (MQ=255)
                           aacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAg                      >  1:1186933‑M2/20‑48 (MQ=255)
                           aacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAg                      <  1:65469‑M2/29‑1 (MQ=255)
                            acatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGa                    <  1:109290‑M2/31‑1 (MQ=255)
                              atactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaat                  <  1:2324578‑M2/33‑1 (MQ=255)
                                 ctcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaataat               <  1:2528167‑M2/36‑1 (MQ=255)
                                 ctcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaataa                <  1:2609388‑M2/35‑1 (MQ=255)
                                  tcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc              >  1:1168811‑M2/13‑49 (MQ=255)
                                   cgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc              >  1:753136‑M2/12‑48 (MQ=255)
                                   cgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc              >  1:1068599‑M2/12‑48 (MQ=255)
                                      gagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa          >  1:135696‑M2/9‑49 (MQ=255)
                                        gcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATg        <  1:954979‑M2/43‑1 (MQ=255)
                                         cgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATgg       >  1:2514706‑M2/6‑49 (MQ=255)
                                         cgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATgg       >  1:196036‑M2/6‑49 (MQ=255)
                                         cgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATgg       <  1:1740925‑M2/44‑1 (MQ=255)
                                          gcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATgg       >  1:304310‑M2/5‑48 (MQ=255)
                                          gcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATgg       >  1:191152‑M2/5‑48 (MQ=255)
                                           cgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGa     >  1:2530281‑M2/4‑49 (MQ=255)
                                            gtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt    >  1:1305053‑M2/3‑49 (MQ=255)
                                            gtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt    <  1:1281116‑M2/47‑1 (MQ=255)
                                            gtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt    >  1:1212404‑M2/3‑49 (MQ=255)
                                            gtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt    <  1:1048510‑M2/47‑1 (MQ=255)
                                              cgctacgctaGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGt  >  1:466101/1‑49 (MQ=255)
                                              |                                                
GTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGT  >  NC_000913/1979284‑1979378

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