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

CATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTAATTTTTTTGTCCTTCAGCGATAATTCACAATCACTTCATTACGTTG  >  NC_000913/1962520‑1962612
                                             |                                               
cATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgc                                              >  1:2007766‑M1/1‑46 (MQ=255)
cATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgc                                              >  1:1574973‑M1/1‑46 (MQ=255)
    gAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacg                                          <  1:1036258‑M1/49‑8 (MQ=255)
     aGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgc                                         >  1:2055072‑M1/1‑41 (MQ=255)
      gCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgct                                        <  1:888524‑M1/49‑10 (MQ=255)
       ccTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgcta                                       <  1:2557724‑M1/49‑11 (MQ=255)
        cTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctag                                      >  1:2569695‑M1/1‑38 (MQ=255)
        cTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctag                                      <  1:502236‑M1/49‑12 (MQ=255)
        cTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctag                                      >  1:2429256‑M1/1‑38 (MQ=255)
        cTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctag                                      >  1:92742‑M1/1‑38 (MQ=255)
        cTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctag                                      >  1:1566259‑M1/1‑38 (MQ=255)
           tttCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagca                                    <  1:2412754‑M1/48‑14 (MQ=255)
            ttCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaa                                   <  1:1854444‑M1/48‑15 (MQ=255)
              cACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagc                                 <  1:2587294‑M1/48‑17 (MQ=255)
              cACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagc                                 <  1:1799026‑M1/48‑17 (MQ=255)
               aCTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagca                                >  1:470200‑M1/1‑31 (MQ=255)
                 tGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaa                              >  1:657173‑M1/1‑29 (MQ=255)
                   aGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaa                             <  1:1736378‑M1/47‑21 (MQ=255)
                    gTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatg                          <  1:2506210‑M1/49‑24 (MQ=255)
                    gTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaat                           >  1:1702803‑M1/1‑26 (MQ=255)
                    gTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaat                           >  1:764284‑M1/1‑26 (MQ=255)
                          aaaCATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacag                    >  1:2013152‑M1/1‑20 (MQ=255)
                          aaaCATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacag                    >  1:829997‑M1/1‑20 (MQ=255)
                           aaCATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacagg                   <  1:1492177‑M1/49‑31 (MQ=255)
                           aaCATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacag                    >  1:1186933‑M1/1‑19 (MQ=255)
                           aaCATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacag                    <  1:65469‑M1/48‑30 (MQ=255)
                            aCATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacagga                  <  1:109290‑M1/49‑32 (MQ=255)
                              aTACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaat                <  1:2324578‑M1/49‑34 (MQ=255)
                                 cTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataat             <  1:2528167‑M1/49‑37 (MQ=255)
                                 cTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataa              <  1:2609388‑M1/48‑36 (MQ=255)
                                  tCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatc            >  1:1168811‑M1/1‑12 (MQ=255)
                                   cgcgAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatc            >  1:1068599‑M1/1‑11 (MQ=255)
                                   cgcgAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatc            >  1:753136‑M1/1‑11 (MQ=255)
                                      gAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaa        >  1:135696‑M1/1‑8 (MQ=255)
                                        gcgcGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatg      <  1:954979‑M1/49‑44 (MQ=255)
                                         cgcgTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgg     >  1:2514706‑M1/1‑5 (MQ=255)
                                         cgcgTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgg     >  1:196036‑M1/1‑5 (MQ=255)
                                         cgcgTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgg     <  1:1740925‑M1/49‑45 (MQ=255)
                                          gcgTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgg     >  1:304310‑M1/1‑4 (MQ=255)
                                          gcgTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgg     >  1:191152‑M1/1‑4 (MQ=255)
                                           cgTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggga   >  1:2530281‑M1/1‑3 (MQ=255)
                                            gTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggat  >  1:1305053‑M1/1‑2 (MQ=255)
                                            gTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggat  <  1:1281116‑M1/49‑48 (MQ=255)
                                            gTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggat  >  1:1212404‑M1/1‑2 (MQ=255)
                                            gTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggat  <  1:1048510‑M1/49‑48 (MQ=255)
                                             |                                               
CATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTAATTTTTTTGTCCTTCAGCGATAATTCACAATCACTTCATTACGTTG  >  NC_000913/1962520‑1962612

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.