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 166 total aligned reads displayed.

TAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTT  >  NC_000913/1979281‑1979379
                                                 |                                                 
agacatcgagcctttttcactgagttattaaacatactcgcgagcgcgtc                                                   <  1:102272‑M2/1‑1 (MQ=255)
agacatcgagcctttttcactgagttattaaacatactcgcgagcgcgtc                                                   <  1:2602211‑M2/1‑1 (MQ=255)
agacatcgagcctttttcactgagttattaaacatactcgcgagcgcgtc                                                   >  1:5774576‑M2/50‑50 (MQ=255)
agacatcgagcctttttcactgagttattaaacatactcgcgagcgcgtc                                                   <  1:5862494‑M2/1‑1 (MQ=255)
agacatcgagcctttttcactgagttattaaacatactcgcgagcgcgtc                                                   >  1:1139466‑M2/50‑50 (MQ=255)
 gacatcgagcctttttcactgagttattaaacatactcgcgagcgcgtcg                                                  <  1:4266264‑M2/2‑1 (MQ=255)
    atcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgcta                                               >  1:663677‑M2/46‑50 (MQ=255)
    atcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgcta                                               >  1:5055247‑M2/46‑50 (MQ=255)
     tcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgctac                                              <  1:3975568‑M2/6‑1 (MQ=255)
     tcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgctac                                              <  1:4541032‑M2/6‑1 (MQ=255)
       gagcctttttcactgagttattaaacatactcgcgagcgcgtcgctacgc                                            >  1:4721766‑M2/43‑50 (MQ=255)
        agcctttttcactgagttattaaacatactcgcgagcgcgtcgctacgc                                            >  1:4397138‑M2/42‑49 (MQ=255)
         gcctttttcactgagttattaaacatactcgcgagcgcgtcgctacgc                                            >  1:1697998‑M2/41‑48 (MQ=255)
          cctttttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTag                                         <  1:1642546‑M2/11‑1 (MQ=255)
          cctttttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTag                                         >  1:4382185‑M2/40‑50 (MQ=255)
            tttttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagc                                        >  1:1280092‑M2/38‑49 (MQ=255)
            tttttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagc                                        >  1:2430517‑M2/38‑49 (MQ=255)
             ttttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagc                                        >  1:4286237‑M2/37‑48 (MQ=255)
             ttttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagc                                        >  1:4328849‑M2/37‑48 (MQ=255)
              tttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaag                                     >  1:324764‑M2/36‑50 (MQ=255)
              tttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaa                                      <  1:4700499‑M2/14‑1 (MQ=255)
              tttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaa                                      <  1:1554444‑M2/14‑1 (MQ=255)
              tttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagca                                       <  1:5096609‑M2/13‑1 (MQ=255)
               ttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaag                                     >  1:3563239‑M2/35‑49 (MQ=255)
               ttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaag                                     <  1:2335049‑M2/15‑1 (MQ=255)
               ttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaag                                     <  1:2111069‑M2/15‑1 (MQ=255)
               ttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaag                                     >  1:5967356‑M2/35‑49 (MQ=255)
               ttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaag                                     >  1:1792651‑M2/35‑49 (MQ=255)
               ttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaag                                     <  1:5391688‑M2/15‑1 (MQ=255)
                tcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaagc                                    >  1:4989235‑M2/34‑49 (MQ=255)
                 cactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaagca                                   <  1:5899857‑M2/17‑1 (MQ=255)
                  actgagttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCaaa                                 <  1:215371‑M2/19‑1 (MQ=255)
                   ctgagttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCaaa                                 >  1:1740015‑M2/31‑49 (MQ=255)
                   ctgagttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCaaa                                 >  1:5591409‑M2/31‑49 (MQ=255)
                     gagttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAAt                              >  1:4473124‑M2/29‑50 (MQ=255)
                     gagttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAAt                              >  1:1967856‑M2/29‑50 (MQ=255)
                      agttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCaaaaa                               <  1:2779823‑M2/21‑1 (MQ=255)
                      agttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCaaaaa                               <  1:2575629‑M2/21‑1 (MQ=255)
                      agttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAAt                              <  1:1811336‑M2/22‑1 (MQ=255)
                      agttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAAt                              >  1:4756492‑M2/28‑49 (MQ=255)
                      agttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAAt                              <  1:2074675‑M2/22‑1 (MQ=255)
                      agttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAAt                              >  1:3983375‑M2/28‑49 (MQ=255)
                      agttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAAt                              <  1:1020074‑M2/22‑1 (MQ=255)
                      agttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATg                             <  1:5346368‑M2/23‑1 (MQ=255)
                      agttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATg                             <  1:5321763‑M2/23‑1 (MQ=255)
                      agttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATg                             <  1:4957750‑M2/23‑1 (MQ=255)
                      agttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATg                             <  1:979600‑M2/23‑1 (MQ=255)
                       gttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATg                             >  1:1399005‑M2/27‑49 (MQ=255)
                       gttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGa                            <  1:2491464‑M2/24‑1 (MQ=255)
                       gttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGa                            <  1:219300‑M2/24‑1 (MQ=255)
                        ttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGa                            >  1:4029038‑M2/26‑49 (MQ=255)
                        ttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGa                            >  1:4102251‑M2/26‑49 (MQ=255)
                          attaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGaaa                          >  1:2519962‑M2/24‑49 (MQ=255)
                          attaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGaaa                          >  1:947911‑M2/24‑49 (MQ=255)
                          attaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGaaa                          >  1:2859227‑M2/24‑49 (MQ=255)
                          attaaacatacacgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGaaa                          <  1:5664860‑M2/26‑1 (MQ=255)
                           ttaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGaaa                          >  1:905655‑M2/23‑48 (MQ=255)
                           ttaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACa                        <  1:2344610‑M2/28‑1 (MQ=255)
                            taaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAg                       >  1:169835‑M2/22‑50 (MQ=255)
                            taaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAg                       >  1:5694431‑M2/22‑50 (MQ=255)
                             aaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAg                       >  1:1572713‑M2/21‑49 (MQ=255)
                             aaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAg                       >  1:2317535‑M2/21‑49 (MQ=255)
                             aaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAg                       >  1:3850146‑M2/21‑49 (MQ=255)
                              aacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAgg                      >  1:3351674‑M2/20‑49 (MQ=255)
                              aacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAgg                      >  1:700420‑M2/20‑49 (MQ=255)
                              aacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAgg                      <  1:2916908‑M2/30‑1 (MQ=255)
                              aacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGa                     <  1:5198489‑M2/31‑1 (MQ=255)
                               acatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaa                    <  1:1031494‑M2/32‑1 (MQ=255)
                               acatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaa                    >  1:4611799‑M2/19‑50 (MQ=255)
                                catactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaat                   <  1:4059637‑M2/33‑1 (MQ=255)
                                catactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaa                    >  1:1205174‑M2/18‑49 (MQ=255)
                                 atactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaat                   >  1:2601750‑M2/17‑49 (MQ=255)
                                  tactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaataa                 >  1:6019181‑M2/16‑50 (MQ=255)
                                  tactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaataa                 >  1:6000117‑M2/16‑50 (MQ=255)
                                  tactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaataa                 <  1:2713570‑M2/35‑1 (MQ=255)
                                  tactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaataa                 <  1:3133785‑M2/35‑1 (MQ=255)
                                  tactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaataa                 <  1:724328‑M2/35‑1 (MQ=255)
                                    ctcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc               >  1:4513293‑M2/14‑50 (MQ=255)
                                    ctcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc               <  1:2060582‑M2/37‑1 (MQ=255)
                                    ctcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc               <  1:4183624‑M2/37‑1 (MQ=255)
                                    ctcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc               >  1:554653‑M2/14‑50 (MQ=255)
                                      cgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCg              >  1:3293844‑M2/12‑49 (MQ=255)
                                       gcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCg              <  1:2991083‑M2/38‑1 (MQ=255)
                                       gcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGa             >  1:2694748‑M2/11‑49 (MQ=255)
                                        cgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa           <  1:4653705‑M2/41‑1 (MQ=255)
                                         gagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt          >  1:3581854‑M2/9‑50 (MQ=255)
                                         gagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt          >  1:1825920‑M2/9‑50 (MQ=255)
                                         gagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt          >  1:5570371‑M2/9‑50 (MQ=255)
                                         gagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt          >  1:2027268‑M2/9‑50 (MQ=255)
                                         gagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt          >  1:2665676‑M2/9‑50 (MQ=255)
                                         gagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt          <  1:3708041‑M2/42‑1 (MQ=255)
                                         gagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt          >  1:3421579‑M2/9‑50 (MQ=255)
                                          agcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATg         >  1:4356316‑M2/8‑50 (MQ=255)
                                          agcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATg         <  1:31854‑M2/43‑1 (MQ=255)
                                          agcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATg         >  1:865224‑M2/8‑50 (MQ=255)
                                              cgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt     <  1:1266486‑M2/47‑1 (MQ=255)
                                               gtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATg    <  1:4899293‑M2/48‑1 (MQ=255)
                                                tCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGt   >  1:4704713‑M2/2‑50 (MQ=255)
                                                 cgctacgctaGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGtt  >  1:1934881/1‑50 (MQ=255)
                                                 cgctacgctaGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGtt  >  1:393521/1‑50 (MQ=255)
                                                 |                                                 
TAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTT  >  NC_000913/1979281‑1979379

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: