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 102 [0] [0] 100 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)97 (1.360) 58/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

TAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTT  >  NC_000913/1979281‑1979379
                                                 |                                                 
agacatcgagcctttttcactgagttattaaacatactcgcgagcgcgtc                                                   >  1:4282517‑M2/50‑50 (MQ=255)
agacatcgagcctttttcactgagttattaaacatactcgcgagcgcgtc                                                   >  1:3022717‑M2/50‑50 (MQ=255)
  acatcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgc                                                 <  1:3448500‑M2/3‑1 (MQ=255)
    atcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgcta                                               <  1:751214‑M2/5‑1 (MQ=255)
     tcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgctac                                              >  1:66840‑M2/45‑50 (MQ=255)
     tcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgctac                                              >  1:1257114‑M2/45‑50 (MQ=255)
     tcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgctac                                              <  1:4627450‑M2/6‑1 (MQ=255)
     tcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgctac                                              <  1:3851241‑M2/6‑1 (MQ=255)
     tcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgctac                                              <  1:2985554‑M2/6‑1 (MQ=255)
       gagcctttttcactgagttattaaacatactcgcgagcgcgtcgctacgc                                            >  1:1027649‑M2/43‑50 (MQ=255)
       gagcctttttcactgagttattaaacatactcgcgagcgcgtcgctacgc                                            >  1:319708‑M2/43‑50 (MQ=255)
       gagcctttttcactgagttattaaacatactcgcgagcgcgtcgctacgc                                            >  1:2907927‑M2/43‑50 (MQ=255)
        ggcctttttcactgagttattaaacatactcgcgagcgcgtcgctacgct                                           <  1:2842396‑M2/9‑1 (MQ=255)
        agcctttttcactgagttattaaacatactcgcgagcgcgtcgctacgct                                           <  1:5183172‑M2/9‑1 (MQ=255)
        agcccttttcactgagttattaaacatactcgcgagcgcgtcgctacgc                                            >  1:4799778‑M2/42‑49 (MQ=255)
         gcctttttcactgagttattaaacatactcgcgagcgcgtcgctacgcta                                          <  1:996410‑M2/10‑1 (MQ=255)
         gcctttttcactgagttattaaacatactcgcgagcgcgtcgctacgc                                            >  1:936934‑M2/41‑48 (MQ=255)
         gcctttttcactgagttattaaacatactcgcgagcgcgtcgctacgc                                            >  1:3588405‑M2/41‑48 (MQ=255)
          cctttttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTag                                         <  1:5322100‑M2/11‑1 (MQ=255)
           ctttttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTag                                         >  1:2507694‑M2/39‑49 (MQ=255)
            tttttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagc                                        >  1:4649301‑M2/38‑49 (MQ=255)
              tttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaag                                     >  1:4579089‑M2/36‑50 (MQ=255)
               ttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaag                                     <  1:2401908‑M2/15‑1 (MQ=255)
               ttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaag                                     >  1:2073908‑M2/35‑49 (MQ=255)
                tcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaagc                                    <  1:2029043‑M2/16‑1 (MQ=255)
                tcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaagc                                    >  1:339603‑M2/34‑49 (MQ=255)
                tcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaagc                                    >  1:3231227‑M2/34‑49 (MQ=255)
                tcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaagc                                    >  1:4000469‑M2/34‑49 (MQ=255)
                tcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaagc                                    >  1:5230907‑M2/34‑49 (MQ=255)
                 cactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaagca                                   >  1:1761798‑M2/33‑49 (MQ=255)
                 cactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaagca                                   >  1:5020385‑M2/33‑49 (MQ=255)
                  actgagttattaaacatactcgcgagcgcgtCGCTACGCTAgcaagcaa                                  >  1:227098‑M2/32‑49 (MQ=255)
                  actgagttattaaacatactcgcgagcgcgtCGCTACGCTAgcaagcaa                                  >  1:4833557‑M2/32‑49 (MQ=255)
                   ctgagttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCaaaa                                <  1:1281580‑M2/20‑1 (MQ=255)
                    tgagttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCaaaaa                               >  1:1738011‑M2/30‑50 (MQ=255)
                     gagttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAAt                              >  1:743175‑M2/29‑50 (MQ=255)
                      agttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAAt                              <  1:4057378‑M2/22‑1 (MQ=255)
                      agttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATg                             <  1:2642766‑M2/23‑1 (MQ=255)
                       gttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATg                             >  1:3820711‑M2/27‑49 (MQ=255)
                         tattaaacatactcgcgagcgcgtCGCTACGCTAGTAAGCAAAAATGaaa                          >  1:2400900‑M2/25‑50 (MQ=255)
                         tattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGaaa                          >  1:2532063‑M2/25‑50 (MQ=255)
                         tattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGaaa                          >  1:2466507‑M2/25‑50 (MQ=255)
                          attaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGaaa                          >  1:4919159‑M2/24‑49 (MQ=255)
                          attaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGaaa                          >  1:2030484‑M2/24‑49 (MQ=255)
                           ttaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGaaa                          >  1:3539824‑M2/23‑48 (MQ=255)
                            taaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAg                       >  1:4076578‑M2/22‑50 (MQ=255)
                            taaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAg                       >  1:4694332‑M2/22‑50 (MQ=255)
                              aacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAgg                      <  1:4853047‑M2/30‑1 (MQ=255)
                              aacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAgg                      <  1:3366646‑M2/30‑1 (MQ=255)
                              aacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGa                     <  1:4280251‑M2/31‑1 (MQ=255)
                              aacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGa                     <  1:1126312‑M2/31‑1 (MQ=255)
                              aacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGa                     <  1:2874808‑M2/31‑1 (MQ=255)
                               acatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaa                    <  1:1621017‑M2/32‑1 (MQ=255)
                               acatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaa                    >  1:2718324‑M2/19‑50 (MQ=255)
                               acatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaa                    <  1:4869006‑M2/32‑1 (MQ=255)
                                catactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaat                   <  1:1179006‑M2/33‑1 (MQ=255)
                                catactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaa                    >  1:3506515‑M2/18‑49 (MQ=255)
                                catactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaa                    >  1:3776452‑M2/18‑49 (MQ=255)
                                 atactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaat                   >  1:1400853‑M2/17‑49 (MQ=255)
                                 atactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaat                   >  1:5523231‑M2/17‑49 (MQ=255)
                                 atactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaat                   >  1:3075831‑M2/17‑49 (MQ=255)
                                  tactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaataa                 >  1:1585962‑M2/16‑50 (MQ=255)
                                  tactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaataa                 <  1:5334585‑M2/35‑1 (MQ=255)
                                  tactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaataa                 <  1:1977190‑M2/35‑1 (MQ=255)
                                  tactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaataa                 <  1:1958341‑M2/35‑1 (MQ=255)
                                  tactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaataa                 <  1:2818726‑M2/35‑1 (MQ=255)
                                   actcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaataat                >  1:3582972‑M2/15‑50 (MQ=255)
                                   actcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaataat                >  1:1969898‑M2/15‑50 (MQ=255)
                                    ctcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaataat                <  1:3288580‑M2/36‑1 (MQ=255)
                                    ctcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc               <  1:4186021‑M2/37‑1 (MQ=255)
                                    ctcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc               >  1:1898306‑M2/14‑50 (MQ=255)
                                       gcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCg              <  1:2149482‑M2/38‑1 (MQ=255)
                                       gcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaa            <  1:2574995‑M2/40‑1 (MQ=255)
                                       gcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGa             >  1:1111999‑M2/11‑49 (MQ=255)
                                       gcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGa             <  1:2542839‑M2/39‑1 (MQ=255)
                                       gcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGa             >  1:853438‑M2/11‑49 (MQ=255)
                                        cgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa           <  1:1989908‑M2/41‑1 (MQ=255)
                                        cgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa           >  1:3475195‑M2/10‑50 (MQ=255)
                                        cgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa           <  1:3516017‑M2/41‑1 (MQ=255)
                                         gagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt          >  1:2261170‑M2/9‑50 (MQ=255)
                                         gagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt          >  1:863860‑M2/9‑50 (MQ=255)
                                         gagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt          >  1:903203‑M2/9‑50 (MQ=255)
                                         gagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt          >  1:5044600‑M2/9‑50 (MQ=255)
                                         gagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt          >  1:4106596‑M2/9‑50 (MQ=255)
                                         gagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt          >  1:1873239‑M2/9‑50 (MQ=255)
                                         gagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt          >  1:4702381‑M2/9‑50 (MQ=255)
                                         gagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt          >  1:3388057‑M2/9‑50 (MQ=255)
                                           gcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATgg        >  1:4533595‑M2/7‑50 (MQ=255)
                                           gcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATgg        <  1:5432872‑M2/44‑1 (MQ=255)
                                            cgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATgg        >  1:5272051‑M2/6‑49 (MQ=255)
                                             gcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGa      <  1:1291629‑M2/46‑1 (MQ=255)
                                              cgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt     >  1:5091066‑M2/4‑50 (MQ=255)
                                              cgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt     <  1:3419890‑M2/47‑1 (MQ=255)
                                              cgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt     <  1:4517792‑M2/47‑1 (MQ=255)
                                              cgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATATCGAAATGGGATg    <  1:3456088‑M2/47‑1 (MQ=255)
                                               gtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATg    >  1:5264793‑M2/3‑50 (MQ=255)
                                               gtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATg    <  1:3500038‑M2/48‑1 (MQ=255)
                                                tCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAACTGGGATGt   >  1:3433369‑M2/2‑50 (MQ=255)
                                                tCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGt   <  1:1886869‑M2/49‑1 (MQ=255)
                                                 cgctacgctaGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGtt  >  1:3306676/1‑50 (MQ=255)
                                                 |                                                 
TAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTT  >  NC_000913/1979281‑1979379

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: