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 98 [0] [0] 96 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)94 (0.980) 57/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

TAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGT  >  NC_000913/1979281‑1979378
                                                 |                                                
agacatcgagcctttttcactgagttattaaacatactcgcgggcgcgtc                                                  <  1:373052‑M2/1‑1 (MQ=255)
agacatcgagcctttttcactgagttattaaacatactcgcgagcgcgtc                                                  >  1:1374074‑M2/50‑50 (MQ=255)
agacatcgagcctttttcactgagttattaaacatactcgcgagcgcgtc                                                  <  1:3366773‑M2/1‑1 (MQ=255)
agacatcgagcctttttcactgagttattaaacatactcgcgagcgcgtc                                                  <  1:1624705‑M2/1‑1 (MQ=255)
agacatcgagcctttttcactgagttattaaacatactcgcgagcgcgtc                                                  <  1:2336379‑M2/1‑1 (MQ=255)
agacatcgagcctttttcactgagttattaaacatactcgcgagcgcgtc                                                  <  1:6401832‑M2/1‑1 (MQ=255)
 gacatcgagcctttttcactgagttattaaacatactcgcgagcgcgtcg                                                 >  1:6346766‑M2/49‑50 (MQ=255)
  acatcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgc                                                >  1:7060043‑M2/48‑50 (MQ=255)
    atcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgcta                                              <  1:2161147‑M2/5‑1 (MQ=255)
    atcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgcta                                              >  1:2060841‑M2/46‑50 (MQ=255)
    atcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgcta                                              <  1:1845624‑M2/5‑1 (MQ=255)
     tcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgctac                                             <  1:1868635‑M2/6‑1 (MQ=255)
     tcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgctac                                             >  1:6779467‑M2/45‑50 (MQ=255)
     tcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgctac                                             <  1:2370265‑M2/6‑1 (MQ=255)
       gagcctttttcactgagttattaaacatactcgcgagcgcgtcgctacgc                                           >  1:3132635‑M2/43‑50 (MQ=255)
       gagcctttttcactgagttattaaacatactcgcgagcgcgtcgctacgc                                           >  1:5587448‑M2/43‑50 (MQ=255)
       gagcctttttcactgagttattaaacatactcgcgagcgcgtcgctacgc                                           >  1:1668972‑M2/43‑50 (MQ=255)
       gagcctttttcactgagttattaaacatactcgcgagcgcgtcgctacg                                            <  1:2920289‑M2/7‑1 (MQ=255)
        agcctttttcactgagttattaaacatactcgcgagcgcgtcgctacgct                                          <  1:4931795‑M2/9‑1 (MQ=255)
        agcctttttcactgagttattaaacatactcgcgagcgcgtcgctacgct                                          <  1:690294‑M2/9‑1 (MQ=255)
         gcctttttcactgagttattaaacatactcgcgagcgcgtcgctacgcta                                         <  1:2442059‑M2/10‑1 (MQ=255)
          cctttttcactgagttattaaacatactcgcgagcgtgtCGCTACGCTag                                        >  1:1219817‑M2/40‑50 (MQ=255)
           ctttttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTag                                        >  1:4096608‑M2/39‑49 (MQ=255)
              tttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaag                                    >  1:3206470‑M2/36‑50 (MQ=255)
              tttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaa                                     <  1:5607085‑M2/14‑1 (MQ=255)
              tttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagca                                      <  1:5876067‑M2/13‑1 (MQ=255)
              tttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagca                                      <  1:2061757‑M2/13‑1 (MQ=255)
               ttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaag                                    >  1:5168482‑M2/35‑49 (MQ=255)
                 cactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaagca                                  >  1:301715‑M2/33‑49 (MQ=255)
                 cactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaagca                                  >  1:5289358‑M2/33‑49 (MQ=255)
                  actgagttattaaacatactcgcgagcgcgtCGCTACGCTAgcaagcaa                                 <  1:6014828‑M2/18‑1 (MQ=255)
                  actgagttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCaaa                                <  1:1932164‑M2/19‑1 (MQ=255)
                  actgagttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCaaa                                <  1:1418252‑M2/19‑1 (MQ=255)
                   ctgagttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCaaaa                               <  1:3991067‑M2/20‑1 (MQ=255)
                   ctgagttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCaaa                                >  1:2645715‑M2/31‑49 (MQ=255)
                   ctgagttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCaaa                                >  1:3915401‑M2/31‑49 (MQ=255)
                    tgagttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCaaaaa                              >  1:7110918‑M2/30‑50 (MQ=255)
                    tgagttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCaaaaa                              >  1:2124606‑M2/30‑50 (MQ=255)
                      agttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCaaaaa                              <  1:6827739‑M2/21‑1 (MQ=255)
                      agttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAAt                             <  1:2478515‑M2/22‑1 (MQ=255)
                      agttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAAt                             <  1:7010381‑M2/22‑1 (MQ=255)
                      agttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAAt                             >  1:2068738‑M2/28‑49 (MQ=255)
                      agttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAAt                             <  1:7361746‑M2/22‑1 (MQ=255)
                      agttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAAt                             <  1:6953916‑M2/22‑1 (MQ=255)
                       gttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGa                           <  1:4127694‑M2/24‑1 (MQ=255)
                        ttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGa                           >  1:4644293‑M2/26‑49 (MQ=255)
                        ttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGa                           >  1:807908‑M2/26‑49 (MQ=255)
                          attaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGaaa                         >  1:4326664‑M2/24‑49 (MQ=255)
                           ttaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAAc                        <  1:6436586‑M2/27‑1 (MQ=255)
                           ttaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAAc                        <  1:2771360‑M2/27‑1 (MQ=255)
                           ttaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAAc                        <  1:860136‑M2/27‑1 (MQ=255)
                           ttaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACa                       <  1:1198694‑M2/28‑1 (MQ=255)
                            taaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAg                      >  1:7131187‑M2/22‑50 (MQ=255)
                             aaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAg                      >  1:906882‑M2/21‑49 (MQ=255)
                             aaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAg                      >  1:474886‑M2/21‑49 (MQ=255)
                             aaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAg                      >  1:4742572‑M2/21‑49 (MQ=255)
                              aacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAgg                     <  1:4496526‑M2/30‑1 (MQ=255)
                               acatactcgcgagcgcgtCGCTACGTTAGCAAGCAAAAACGAAACAGGaa                   >  1:202755‑M2/19‑50 (MQ=255)
                               acatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaa                   <  1:741463‑M2/32‑1 (MQ=255)
                                catactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaa                   >  1:2122683‑M2/18‑49 (MQ=255)
                                catactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaa                   >  1:3887470‑M2/18‑49 (MQ=255)
                                 atactcgggagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaata                 <  1:3041844‑M2/34‑1 (MQ=255)
                                 atactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaata                 <  1:6557748‑M2/34‑1 (MQ=255)
                                 atactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaata                 <  1:6624847‑M2/34‑1 (MQ=255)
                                 atactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaata                 <  1:3901242‑M2/34‑1 (MQ=255)
                                 atactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaata                 <  1:6067002‑M2/34‑1 (MQ=255)
                                  tactcgcgagcgcgtCGCTACGCTAGCCAGCAAAAATGAAACAGGaataa                >  1:139861‑M2/16‑50 (MQ=255)
                                  tactcgcgagcgcgtCGCTACGCTAGCAAGCAGAAATGAAACAGGaataa                <  1:5019063‑M2/35‑1 (MQ=255)
                                  tactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaataa                <  1:449405‑M2/35‑1 (MQ=255)
                                  tactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaataa                >  1:6304605‑M2/16‑50 (MQ=255)
                                  tactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaataa                <  1:6203657‑M2/35‑1 (MQ=255)
                                  tactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaataa                >  1:2189573‑M2/16‑50 (MQ=255)
                                  tactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaataa                <  1:3810437‑M2/35‑1 (MQ=255)
                                   actcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaataat               >  1:1802376‑M2/15‑50 (MQ=255)
                                    ctcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc              <  1:3211691‑M2/37‑1 (MQ=255)
                                    ctcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc              >  1:4565759‑M2/14‑50 (MQ=255)
                                      cgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCg             >  1:4143649‑M2/12‑49 (MQ=255)
                                       gcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGa            >  1:968932‑M2/11‑49 (MQ=255)
                                       gcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGa            >  1:4550225‑M2/11‑49 (MQ=255)
                                        cgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa          >  1:6192094‑M2/10‑50 (MQ=255)
                                        cgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa          <  1:7462140‑M2/41‑1 (MQ=255)
                                         gagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt         >  1:2139115‑M2/9‑50 (MQ=255)
                                         gagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt         >  1:6939619‑M2/9‑50 (MQ=255)
                                          agcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATg        <  1:3505668‑M2/43‑1 (MQ=255)
                                            cgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATggg      <  1:3969505‑M2/45‑1 (MQ=255)
                                            cgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATgg       >  1:60948‑M2/6‑49 (MQ=255)
                                             gcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGa     >  1:4142222‑M2/5‑50 (MQ=255)
                                             gcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGa     <  1:871998‑M2/46‑1 (MQ=255)
                                              cgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt    >  1:2445119‑M2/4‑50 (MQ=255)
                                              cgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt    >  1:564022‑M2/4‑50 (MQ=255)
                                              cgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt    <  1:5683610‑M2/47‑1 (MQ=255)
                                               gtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATg   <  1:4941205‑M2/48‑1 (MQ=255)
                                               gtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATg   <  1:2904828‑M2/48‑1 (MQ=255)
                                               gtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATg   <  1:5654956‑M2/48‑1 (MQ=255)
                                                tCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGt  <  1:6571059‑M2/49‑1 (MQ=255)
                                                tCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGt  <  1:516669‑M2/49‑1 (MQ=255)
                                                 |                                                
TAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGT  >  NC_000913/1979281‑1979378

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