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 171 [0] [0] 166 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)163 (2.100) 72/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 171 total aligned reads displayed.

AAGACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTAATTTTTTTGTCCTTCAGCGATAATTCACAATCACTTCATTACGTTGTA  >  NC_000913/1962516‑1962614
                                                 |                                                 
aaGACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGT                                                   <  1:478817‑M1/50‑1 (MQ=255)
aaGACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGT                                                   <  1:4101034‑M1/50‑1 (MQ=255)
aaGACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGT                                                   >  1:3295934‑M1/1‑50 (MQ=255)
aaGACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGT                                                   >  1:1675764‑M1/1‑50 (MQ=255)
aaGACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGT                                                   >  1:5447885‑M1/1‑50 (MQ=255)
aaGACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGT                                                   >  1:2716799‑M1/1‑50 (MQ=255)
 aGACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTc                                                  <  1:102272‑M1/50‑2 (MQ=255)
 aGACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTc                                                  <  1:278080‑M1/50‑2 (MQ=255)
  gACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcg                                                 <  1:3937566‑M1/50‑3 (MQ=255)
  gACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcg                                                 <  1:4266264‑M1/50‑3 (MQ=255)
     aTCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgcta                                              <  1:3549568‑M1/50‑6 (MQ=255)
     aTCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgcta                                              <  1:6000058‑M1/50‑6 (MQ=255)
     aTCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgcta                                              <  1:2682426‑M1/50‑6 (MQ=255)
     aTCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgcta                                              >  1:663677‑M1/1‑45 (MQ=255)
      tCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctac                                             <  1:3975568‑M1/50‑7 (MQ=255)
      tCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctac                                             <  1:4541032‑M1/50‑7 (MQ=255)
      tCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctac                                             <  1:2121000‑M1/50‑7 (MQ=255)
         aGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgc                                           >  1:4397138‑M1/1‑41 (MQ=255)
             tttttCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagc                                       >  1:2430517‑M1/1‑37 (MQ=255)
              ttttCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagc                                       >  1:4286237‑M1/1‑36 (MQ=255)
              ttttCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagc                                       >  1:4328849‑M1/1‑36 (MQ=255)
              ttttCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagc                                       >  1:4848616‑M1/1‑36 (MQ=255)
              ttttCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagc                                       >  1:1349882‑M1/1‑36 (MQ=255)
               tttCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaag                                    >  1:324764‑M1/1‑35 (MQ=255)
               tttCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaa                                     <  1:4700499‑M1/49‑15 (MQ=255)
               tttCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaa                                     <  1:1554444‑M1/49‑15 (MQ=255)
               tttCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagca                                      <  1:2054576‑M1/48‑14 (MQ=255)
               tttCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagca                                      <  1:5096609‑M1/48‑14 (MQ=255)
                ttCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaag                                    <  1:5391688‑M1/49‑16 (MQ=255)
                 tCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagc                                   >  1:4989235‑M1/1‑33 (MQ=255)
                 tCACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagc                                   >  1:1187351‑M1/1‑33 (MQ=255)
                  cACTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagca                                  <  1:5899857‑M1/49‑18 (MQ=255)
                   aCTGAGTTATTAAACATGCTCGCGAGCGCGTcgctacgctagcaagcaaa                                <  1:5769323‑M1/50‑20 (MQ=255)
                   aCTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaa                                 >  1:5856331‑M1/1‑31 (MQ=255)
                    cTGAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaa                                >  1:1740015‑M1/1‑30 (MQ=255)
                      gAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaat                             >  1:4473124‑M1/1‑28 (MQ=255)
                      gAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaat                             >  1:1967856‑M1/1‑28 (MQ=255)
                      gAGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaat                             >  1:2591388‑M1/1‑28 (MQ=255)
                       aGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatg                            <  1:5346368‑M1/50‑24 (MQ=255)
                       aGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatg                            <  1:979600‑M1/50‑24 (MQ=255)
                       aGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatg                            <  1:5321763‑M1/50‑24 (MQ=255)
                       aGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatg                            <  1:4957750‑M1/50‑24 (MQ=255)
                       aGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaat                             <  1:1020074‑M1/49‑23 (MQ=255)
                       aGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaat                             <  1:2074675‑M1/49‑23 (MQ=255)
                       aGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaat                             <  1:1811336‑M1/49‑23 (MQ=255)
                       aGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaat                             >  1:1597551‑M1/1‑27 (MQ=255)
                       aGTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaa                              <  1:2779823‑M1/48‑22 (MQ=255)
                        gTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatga                           <  1:219300‑M1/50‑25 (MQ=255)
                        gTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatga                           <  1:1653670‑M1/50‑25 (MQ=255)
                        gTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatga                           <  1:2491464‑M1/50‑25 (MQ=255)
                        gTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatg                            >  1:1399005‑M1/1‑26 (MQ=255)
                        gTTATTAAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatg                            >  1:4692021‑M1/1‑26 (MQ=255)
                         ttattaAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatga                           >  1:405291‑M1/1‑25 (MQ=255)
                         ttattaAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatga                           >  1:4102251‑M1/1‑25 (MQ=255)
                           attaAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaa                         >  1:3836332‑M1/1‑23 (MQ=255)
                           attaAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaa                         >  1:947911‑M1/1‑23 (MQ=255)
                            ttaAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaaca                       <  1:2344610‑M1/50‑29 (MQ=255)
                             taAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacag                      >  1:169835‑M1/1‑21 (MQ=255)
                             taAACATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacag                      >  1:5694431‑M1/1‑21 (MQ=255)
                              aaaCATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacag                      >  1:398408‑M1/1‑20 (MQ=255)
                              aaaCATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacag                      >  1:2317535‑M1/1‑20 (MQ=255)
                               aaCATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacagga                    <  1:5198489‑M1/50‑32 (MQ=255)
                               aaCATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacagg                     <  1:2916908‑M1/49‑31 (MQ=255)
                               aaCATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacagg                     >  1:2617622‑M1/1‑19 (MQ=255)
                                aCATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaa                   >  1:4611799‑M1/1‑18 (MQ=255)
                                aCATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaa                   <  1:1031494‑M1/50‑33 (MQ=255)
                                 cATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaa                   >  1:1205174‑M1/1‑17 (MQ=255)
                                 cATACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaa                   >  1:2090879‑M1/1‑17 (MQ=255)
                                  ttaCTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaat                  >  1:2321139‑M1/2‑16 (MQ=255)
                                  aTACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaat                  >  1:2574161‑M1/1‑16 (MQ=255)
                                  aTACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaat                  >  1:2601750‑M1/1‑16 (MQ=255)
                                  aTACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaat                  >  1:176267‑M1/1‑16 (MQ=255)
                                   tACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataa                <  1:724328‑M1/50‑36 (MQ=255)
                                   tACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataa                >  1:6019181‑M1/1‑15 (MQ=255)
                                   tACTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataa                <  1:3133785‑M1/50‑36 (MQ=255)
                                     cTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatc              >  1:554653‑M1/1‑13 (MQ=255)
                                     cTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatc              <  1:4183624‑M1/50‑38 (MQ=255)
                                     cTCGCGAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataat               <  1:5575739‑M1/49‑37 (MQ=255)
                                       cgcgAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcg             >  1:3293844‑M1/1‑11 (MQ=255)
                                        gcgAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaa           <  1:5598365‑M1/50‑41 (MQ=255)
                                        gcgAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcga            >  1:4520783‑M1/1‑10 (MQ=255)
                                        gcgAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcga            >  1:2412823‑M1/1‑10 (MQ=255)
                                        gcgAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcga            >  1:2694748‑M1/1‑10 (MQ=255)
                                        gcgAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcg             <  1:2991083‑M1/48‑39 (MQ=255)
                                          gAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaat         >  1:5052083‑M1/1‑8 (MQ=255)
                                          gAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaat         >  1:896098‑M1/1‑8 (MQ=255)
                                          gAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaat         <  1:3390669‑M1/50‑43 (MQ=255)
                                          gAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaat         >  1:3581854‑M1/1‑8 (MQ=255)
                                          gAGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaat         >  1:1825920‑M1/1‑8 (MQ=255)
                                           aGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatg        >  1:4356316‑M1/1‑7 (MQ=255)
                                           aGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatg        >  1:675265‑M1/1‑7 (MQ=255)
                                           aGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatg        >  1:865224‑M1/1‑7 (MQ=255)
                                           aGCGCGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatg        <  1:31854‑M1/50‑44 (MQ=255)
                                            gcgcGTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgg       <  1:5931861‑M1/50‑45 (MQ=255)
                                              gcgTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggga     <  1:3598045‑M1/50‑47 (MQ=255)
                                              gcgTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggga     >  1:4725986‑M1/1‑4 (MQ=255)
                                              gcgTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatggga     >  1:1987768‑M1/1‑4 (MQ=255)
                                               cgTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggat    <  1:1266486‑M1/50‑48 (MQ=255)
                                               cgTcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggat    >  1:4388336‑M1/1‑3 (MQ=255)
                                                 tcgctacgctagcaagcaaaaatgaaacaggaataatcgaaatgggatgt  >  1:4704713‑M1/1‑1 (MQ=255)
                                                 |                                                 
AAGACATCGAGCCTTTTTCACTGAGTTATTAAACATACTCGCGAGCGCGTAATTTTTTTGTCCTTCAGCGATAATTCACAATCACTTCATTACGTTGTA  >  NC_000913/1962516‑1962614

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