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 169 [0] [0] 168 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)166 (1.120) 73/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 168 total aligned reads displayed.

TAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTT  >  NC_000913/1979281‑1979379
                                                 |                                                 
agacatcgagcctttttcactgagttattaaacatactcgcgagcgcgtc                                                   <  1:5144822‑M2/1‑1 (MQ=255)
agacatcgagcctttttcactgagttattaaacatactcgcgagcgcgtc                                                   <  1:694916‑M2/1‑1 (MQ=255)
agacatcgagcctttttcactgagttattaaacatactcgcgagcgcgtc                                                   <  1:4071657‑M2/1‑1 (MQ=255)
   catcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgct                                                <  2:334975‑M2/4‑1 (MQ=255)
   catcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgc                                                 >  2:676168‑M2/47‑49 (MQ=255)
   catcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgc                                                 >  1:5766194‑M2/47‑49 (MQ=255)
   catcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgc                                                 >  2:5107067‑M2/47‑49 (MQ=255)
    atcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgcta                                               <  2:4203022‑M2/5‑1 (MQ=255)
    atcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgcta                                               >  1:3522312‑M2/46‑50 (MQ=255)
     tcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgctac                                              <  1:3751094‑M2/6‑1 (MQ=255)
     tcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgctac                                              <  2:4103536‑M2/6‑1 (MQ=255)
     tcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgctac                                              <  1:4212586‑M2/6‑1 (MQ=255)
     tcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgctac                                              <  2:1198138‑M2/6‑1 (MQ=255)
     tcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgctac                                              >  2:1033594‑M2/45‑50 (MQ=255)
      cgagcctttttcactgagttattaaacatactcgcgagcgcgtcgctacg                                             >  2:3947555‑M2/44‑50 (MQ=255)
       gagcctttttcactgagttattaaccatactcgcgagcgcgtcgctacg                                             <  2:753874‑M2/7‑1 (MQ=255)
       gagcctttttcactgagttattaaacatactcgcgagcgcgtcgctacgc                                            <  1:3619715‑M2/8‑1 (MQ=255)
       gagcctttttcactgagttattaaacatactcgcgagcgcgtcgctacgc                                            >  1:3601771‑M2/43‑50 (MQ=255)
       gagcctttttcactgagttattaaacatactcgcgagcgcgtcgctacgc                                            >  2:3195043‑M2/43‑50 (MQ=255)
        agcctttttcactgagttattaaacatactcgcgagcgcgtcgctacgct                                           <  2:165060‑M2/9‑1 (MQ=255)
        agcctttttcactgagttattaaacatactcgcgagcgcgtcgctacgct                                           <  1:762851‑M2/9‑1 (MQ=255)
        agcctttttcactgagttattaaacatactcgcgagcgcgtcgctacgct                                           <  2:3630820‑M2/9‑1 (MQ=255)
        agcctttttcactgagttattaaacatactcgcgagcgcgtcgctacgc                                            >  2:5430991‑M2/42‑49 (MQ=255)
         gcctttttcactgagttattaaacatactcgcgagcgcgtcgctacgcta                                          <  1:1210582‑M2/10‑1 (MQ=255)
         gcctttttcactgagttattaaacatactcgcgagcgcgtcgctacgcta                                          <  1:843030‑M2/10‑1 (MQ=255)
         gcctttttcactgagttattaaacatactcgcgagcgcgtcgctacgcta                                          <  2:660258‑M2/10‑1 (MQ=255)
          cctttttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTag                                         <  1:3322411‑M2/11‑1 (MQ=255)
          cctttttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTag                                         <  1:4591909‑M2/11‑1 (MQ=255)
           ctttttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTag                                         >  2:1362059‑M2/39‑49 (MQ=255)
            tttttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagc                                        >  1:3670602‑M2/38‑49 (MQ=255)
            tttttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagc                                        >  1:3840151‑M2/38‑49 (MQ=255)
              tttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagca                                       <  1:2993026‑M2/13‑1 (MQ=255)
               ttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaag                                     <  1:5349421‑M2/15‑1 (MQ=255)
               ttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaag                                     <  2:1467746‑M2/15‑1 (MQ=255)
               ttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaag                                     <  2:705789‑M2/15‑1 (MQ=255)
                tcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaagc                                    >  2:2270225‑M2/34‑49 (MQ=255)
                 cactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaagca                                   >  1:2390382‑M2/33‑49 (MQ=255)
                  actgagttattaaacatactcgcgagcgcgtCGCTACGCTAgcaagcaa                                  <  2:5798925‑M2/18‑1 (MQ=255)
                  actgagttattaaacatactcgcgagcgcgtCGCTACGCTAgcaagcaa                                  <  1:5077041‑M2/18‑1 (MQ=255)
                   ctgagttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCaaaa                                <  2:92622‑M2/20‑1 (MQ=255)
                   ctgagttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCaaaa                                <  2:4817824‑M2/20‑1 (MQ=255)
                    tgagttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCaaaaa                               >  1:2205634‑M2/30‑50 (MQ=255)
                     gagttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAAt                              >  2:143120‑M2/29‑50 (MQ=255)
                      agttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCaaaaa                               <  1:3148540‑M2/21‑1 (MQ=255)
                      agttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCaaaaa                               <  2:4872533‑M2/21‑1 (MQ=255)
                      agttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCaaaaa                               <  1:5416618‑M2/21‑1 (MQ=255)
                      agttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAAt                              <  1:5747522‑M2/22‑1 (MQ=255)
                      agttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAAt                              >  1:5661263‑M2/28‑49 (MQ=255)
                      agttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAAt                              <  1:5059736‑M2/22‑1 (MQ=255)
                      agttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATg                             <  1:4094493‑M2/23‑1 (MQ=255)
                       gttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGa                            <  2:3047324‑M2/24‑1 (MQ=255)
                        ttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGaa                           <  1:2247803‑M2/25‑1 (MQ=255)
                        ttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGaa                           <  2:3728730‑M2/25‑1 (MQ=255)
                         tattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGaaa                          >  1:5838734‑M2/25‑50 (MQ=255)
                          attaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGaaa                          <  2:5309318‑M2/26‑1 (MQ=255)
                          attaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGaaa                          <  2:4179805‑M2/26‑1 (MQ=255)
                           ttaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGaaa                          >  1:5518444‑M2/23‑48 (MQ=255)
                           ttaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGaaa                          >  1:2630911‑M2/23‑48 (MQ=255)
                             aaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAg                       >  1:4124334‑M2/21‑49 (MQ=255)
                              aacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAgg                      <  1:4782882‑M2/30‑1 (MQ=255)
                              aacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAgg                      >  2:4122714‑M2/20‑49 (MQ=255)
                              aacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGa                     <  1:4428831‑M2/31‑1 (MQ=255)
                               acatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaa                    <  1:1334395‑M2/32‑1 (MQ=255)
                               acatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaa                    <  2:5420216‑M2/32‑1 (MQ=255)
                               acatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaa                    <  2:1100500‑M2/32‑1 (MQ=255)
                                catactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaat                   <  1:4771093‑M2/33‑1 (MQ=255)
                                 atactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaata                  <  1:5600584‑M2/34‑1 (MQ=255)
                                 atactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaata                  <  1:4332639‑M2/34‑1 (MQ=255)
                                 atactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaat                   >  1:3899665‑M2/17‑49 (MQ=255)
                                  tactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaataa                 <  1:1773982‑M2/35‑1 (MQ=255)
                                   actcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaataat                >  2:320551‑M2/15‑50 (MQ=255)
                                    ctcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc               <  1:3217780‑M2/37‑1 (MQ=255)
                                    ctcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc               >  1:4535422‑M2/14‑50 (MQ=255)
                                     tcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc               >  2:5178442‑M2/13‑49 (MQ=255)
                                       gcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCg              <  1:519953‑M2/38‑1 (MQ=255)
                                       gcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCg              <  2:1355504‑M2/38‑1 (MQ=255)
                                       gcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCg              <  2:586990‑M2/38‑1 (MQ=255)
                                       gcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGa             <  1:1571337‑M2/39‑1 (MQ=255)
                                       gcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGa             >  2:2611874‑M2/11‑49 (MQ=255)
                                        cgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa           <  1:4554814‑M2/41‑1 (MQ=255)
                                        cgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa           <  1:2294089‑M2/41‑1 (MQ=255)
                                        cgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa           >  1:4653989‑M2/10‑50 (MQ=255)
                                         gagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt          >  1:2088644‑M2/9‑50 (MQ=255)
                                         gagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt          >  2:4941717‑M2/9‑50 (MQ=255)
                                         gagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt          >  2:3616644‑M2/9‑50 (MQ=255)
                                          agcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATg         >  2:2735401‑M2/8‑50 (MQ=255)
                                          agcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATg         <  2:1483478‑M2/43‑1 (MQ=255)
                                          agcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATg         <  1:257097‑M2/43‑1 (MQ=255)
                                           gcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATgg        <  1:3492543‑M2/44‑1 (MQ=255)
                                           gcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATgg        <  1:5315845‑M2/44‑1 (MQ=255)
                                           gcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATgg        <  2:5653834‑M2/44‑1 (MQ=255)
                                            cgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATggg       <  1:3562496‑M2/45‑1 (MQ=255)
                                             gcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGa      <  2:3334517‑M2/46‑1 (MQ=255)
                                             gcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGa      <  1:5514496‑M2/46‑1 (MQ=255)
                                             gcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGa      <  1:3249953‑M2/46‑1 (MQ=255)
                                              cgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt     <  1:2838194‑M2/47‑1 (MQ=255)
                                              cgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt     <  2:1931346‑M2/47‑1 (MQ=255)
                                               gtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATg    <  2:4334221‑M2/48‑1 (MQ=255)
                                                tCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGt   <  1:161459‑M2/49‑1 (MQ=255)
                                                 cgctacgctaGCAAGCAAAACTGAAACAGGAATAATCGAAATGGGCTGtt  >  2:1971421/1‑50 (MQ=255)
                                                 |                                                 
TAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTT  >  NC_000913/1979281‑1979379

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