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 107 [0] [0] 107 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)105 (0.940) 59/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

Only 100 of 107 total aligned reads displayed.

AGACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTAATTTTTTTGTCCTTCAGCGATAATTCACAATCACTTCATTACGTTGTA  >  NC_000913/1962517‑1962614
                                                |                                                 
aGACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTc                                                  <  1:486818‑M1/50‑2 (MQ=255)
 gACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcg                                                 >  1:442850‑M1/1‑48 (MQ=255)
 gACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcg                                                 >  1:4611927‑M1/1‑48 (MQ=255)
   cATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgc                                                >  1:619802‑M1/1‑46 (MQ=255)
    aTCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgcta                                              >  1:3371420‑M1/1‑45 (MQ=255)
    aTCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgcta                                              >  1:4684186‑M1/1‑45 (MQ=255)
    aTCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgcta                                              >  1:2496971‑M1/1‑45 (MQ=255)
     tCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctac                                             >  1:4290944‑M1/1‑44 (MQ=255)
     tCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctac                                             <  1:4504060‑M1/50‑7 (MQ=255)
     tCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctac                                             <  1:456543‑M1/50‑7 (MQ=255)
      cGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacg                                            >  1:1532197‑M1/1‑43 (MQ=255)
       gAGCCTTTTTCTCTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgc                                           >  1:6932911‑M1/1‑42 (MQ=255)
       gAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgc                                           >  1:5162332‑M1/1‑42 (MQ=255)
        aGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgc                                           >  1:3386374‑M1/1‑41 (MQ=255)
         gCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgcta                                         <  1:2452517‑M1/50‑11 (MQ=255)
         gCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgcta                                         <  1:4794793‑M1/50‑11 (MQ=255)
         gCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgcta                                         <  1:7897977‑M1/50‑11 (MQ=255)
          ccTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctag                                        <  1:1528660‑M1/50‑12 (MQ=255)
          ccTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctag                                        <  1:6628341‑M1/50‑12 (MQ=255)
           cTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagc                                       <  1:2793851‑M1/50‑13 (MQ=255)
           cTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagc                                       <  1:2825455‑M1/50‑13 (MQ=255)
           cTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctag                                        >  1:6517205‑M1/1‑38 (MQ=255)
           cTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctag                                        >  1:4512251‑M1/1‑38 (MQ=255)
            tttttCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagc                                       >  1:3628109‑M1/1‑37 (MQ=255)
              tttCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaag                                    >  1:4054301‑M1/1‑35 (MQ=255)
              tttCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaag                                    >  1:8558200‑M1/1‑35 (MQ=255)
              tttCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaag                                    >  1:7996879‑M1/1‑35 (MQ=255)
              tttCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagca                                      <  1:3978758‑M1/48‑14 (MQ=255)
               ttCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaag                                    <  1:7538366‑M1/49‑16 (MQ=255)
                tCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagc                                   >  1:4040980‑M1/1‑33 (MQ=255)
                tCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagc                                   <  1:8121106‑M1/49‑17 (MQ=255)
                tCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagc                                   >  1:6704081‑M1/1‑33 (MQ=255)
                tCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaaga                                   <  1:2674432‑M1/49‑17 (MQ=255)
                 cACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagca                                  <  1:696035‑M1/49‑18 (MQ=255)
                 cACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagca                                  <  1:4232038‑M1/49‑18 (MQ=255)
                  aCTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaa                                 <  1:3993607‑M1/49‑19 (MQ=255)
                   cTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaa                               <  1:5667460‑M1/50‑21 (MQ=255)
                   cTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaa                               <  1:5631567‑M1/50‑21 (MQ=255)
                   cTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaa                               <  1:4154761‑M1/50‑21 (MQ=255)
                   cTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaa                               <  1:3841228‑M1/50‑21 (MQ=255)
                   cTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaa                                >  1:8401814‑M1/1‑30 (MQ=255)
                    tGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaa                              >  1:3605713‑M1/1‑29 (MQ=255)
                     gAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaat                             >  1:126821‑M1/1‑28 (MQ=255)
                     gAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaat                             >  1:6136724‑M1/1‑28 (MQ=255)
                       gTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatg                            >  1:3501904‑M1/1‑26 (MQ=255)
                        ttattaAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaa                          <  1:8235528‑M1/50‑26 (MQ=255)
                        ttattaAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaa                          <  1:1715050‑M1/50‑26 (MQ=255)
                          attaAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaa                         >  1:4737648‑M1/1‑23 (MQ=255)
                          attaAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaa                         >  1:1652885‑M1/1‑23 (MQ=255)
                           ttaAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaac                        <  1:7564654‑M1/49‑28 (MQ=255)
                             aaaCATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacag                      <  1:1058190‑M1/49‑30 (MQ=255)
                             aaaCATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacag                      <  1:4953275‑M1/49‑30 (MQ=255)
                              aaCATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacagg                     <  1:2158378‑M1/49‑31 (MQ=255)
                              aaCATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacagg                     <  1:6670173‑M1/49‑31 (MQ=255)
                               aCATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaa                   <  1:6660535‑M1/50‑33 (MQ=255)
                               aCATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaa                   <  1:1749628‑M1/50‑33 (MQ=255)
                               aCATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaa                   >  1:7015534‑M1/1‑18 (MQ=255)
                               aCATACTCGCGAGCGCGTcgctacgctagcaagcaaaaataaaacaggaa                   <  1:6704995‑M1/50‑33 (MQ=255)
                                cATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaat                  <  1:5001615‑M1/50‑34 (MQ=255)
                                cATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaat                  <  1:5779888‑M1/50‑34 (MQ=255)
                                cATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaat                  <  1:8505907‑M1/50‑34 (MQ=255)
                                 aTACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaata                 <  1:3826883‑M1/50‑35 (MQ=255)
                                 aTACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaat                  >  1:7582996‑M1/1‑16 (MQ=255)
                                 aTACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaat                  >  1:8227442‑M1/1‑16 (MQ=255)
                                 aTACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaat                  >  1:7884641‑M1/1‑16 (MQ=255)
                                  tACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataa                <  1:8310133‑M1/50‑36 (MQ=255)
                                  tACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataa                >  1:3082518‑M1/1‑15 (MQ=255)
                                  tACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataa                <  1:3145629‑M1/50‑36 (MQ=255)
                                  tACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataa                <  1:4244657‑M1/50‑36 (MQ=255)
                                   aCTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataat               >  1:4831501‑M1/1‑14 (MQ=255)
                                    cTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatc              <  1:6069341‑M1/50‑38 (MQ=255)
                                    cTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataat               <  1:5935094‑M1/49‑37 (MQ=255)
                                     tCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatc              >  1:8330278‑M1/1‑12 (MQ=255)
                                       gcgAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcg             <  1:4128878‑M1/48‑39 (MQ=255)
                                        cgAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaa          >  1:1157857‑M1/1‑9 (MQ=255)
                                        cgAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaa          <  1:3521715‑M1/50‑42 (MQ=255)
                                        cgAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaa          >  1:4027740‑M1/1‑9 (MQ=255)
                                         gAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaat         >  1:7783511‑M1/1‑8 (MQ=255)
                                         gAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaat         >  1:1234780‑M1/1‑8 (MQ=255)
                                         gAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaat         >  1:2760625‑M1/1‑8 (MQ=255)
                                         gAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaat         <  1:4563528‑M1/50‑43 (MQ=255)
                                          aGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatg        >  1:6239502‑M1/1‑7 (MQ=255)
                                          aGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatg        >  1:6235893‑M1/1‑7 (MQ=255)
                                          aGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatg        <  1:8334034‑M1/50‑44 (MQ=255)
                                          aGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatg        >  1:6456302‑M1/1‑7 (MQ=255)
                                          aGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatg        >  1:5675967‑M1/1‑7 (MQ=255)
                                           gcgcGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgg       <  1:5082027‑M1/50‑45 (MQ=255)
                                           gcgcGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgg       <  1:683904‑M1/50‑45 (MQ=255)
                                           gcgcGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgg       <  1:5244835‑M1/50‑45 (MQ=255)
                                            cgcgTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggg      <  1:2658239‑M1/50‑46 (MQ=255)
                                             gcgTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggga     <  1:3370102‑M1/50‑47 (MQ=255)
                                             gcgTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggga     >  1:7413445‑M1/1‑4 (MQ=255)
                                              cgTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggat    <  1:576318‑M1/50‑48 (MQ=255)
                                              cgTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggat    >  1:8356586‑M1/1‑3 (MQ=255)
                                              cgTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggat    >  1:1674267‑M1/1‑3 (MQ=255)
                                              cgTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggat    <  1:4816409‑M1/50‑48 (MQ=255)
                                               gTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatg   <  1:5521887‑M1/50‑49 (MQ=255)
                                               gTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatg   >  1:5167119‑M1/1‑2 (MQ=255)
                                                tcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgt  <  1:2776928‑M1/50‑50 (MQ=255)
                                                tcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgt  <  1:8359193‑M1/50‑50 (MQ=255)
                                                |                                                 
AGACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTAATTTTTTTGTCCTTCAGCGATAATTCACAATCACTTCATTACGTTGTA  >  NC_000913/1962517‑1962614

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