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 153 [0] [0] 151 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)148 (1.390) 70/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 151 total aligned reads displayed.

TAAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGTT  >  NC_000913/1979281‑1979379
                                                 |                                                 
agacatcgagcctttttcactgagttattaaacatactcgcgagcgcgtc                                                   <  1:88484‑M2/1‑1 (MQ=255)
 gacatcgagcctttttcactgagttattaaacatactcgcgagcgcgtcg                                                  >  1:5305533‑M2/49‑50 (MQ=255)
 gacatcgagcctttttcactgagttattaaacatactcgcgagcgcgtcg                                                  >  1:2292238‑M2/49‑50 (MQ=255)
  acatcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgc                                                 >  1:1322355‑M2/48‑50 (MQ=255)
  acatcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgc                                                 >  1:4939861‑M2/48‑50 (MQ=255)
  acatcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgc                                                 >  1:1575601‑M2/48‑50 (MQ=255)
   catcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgct                                                <  1:1225429‑M2/4‑1 (MQ=255)
   catcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgc                                                 >  1:5252681‑M2/47‑49 (MQ=255)
    atcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgcta                                               >  1:7438059‑M2/46‑50 (MQ=255)
    atcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgcta                                               >  1:6427787‑M2/46‑50 (MQ=255)
     tcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgctac                                              >  1:4041804‑M2/45‑50 (MQ=255)
     tcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgctac                                              <  1:5507314‑M2/6‑1 (MQ=255)
      cgagcctttttcactgagttattaaacatactcgcgagcgcgtcgctacg                                             >  1:4449036‑M2/44‑50 (MQ=255)
      cgagcctttttcactgagttattaaacatactcgcgagcgcgtcgctacg                                             >  1:3367400‑M2/44‑50 (MQ=255)
      cgagcctttttcactgagttattaaacatactcgcgagcgcgtcgctacg                                             >  1:8148150‑M2/44‑50 (MQ=255)
       gagcctttttcactgagttattaaacatactcgcgagcgcgtcgctacgc                                            >  1:3371743‑M2/43‑50 (MQ=255)
       gagcctttttcactgagttattaaacatactcgcgagcgcgtcgctacgc                                            >  1:2357158‑M2/43‑50 (MQ=255)
        agcctttttcactgagttattaaacatactcgcgagcgcgtcgctacgct                                           <  1:1877252‑M2/9‑1 (MQ=255)
        agcctttttcactgagttattaaacatactcgcgagcgcgtcgctacgct                                           <  1:4670424‑M2/9‑1 (MQ=255)
        agcctttttcactgagttattaaacatactcgcgagcgcgtcgctacgct                                           <  1:7697739‑M2/9‑1 (MQ=255)
         gcctttttcactgagttattaaacatactcgcgagcgcgtcgctacgcta                                          <  1:6427749‑M2/10‑1 (MQ=255)
         gcctttttcactgagttattaaacatactcgcgagcgcgtcgctacgc                                            >  1:1760296‑M2/41‑48 (MQ=255)
           ctttttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagc                                        <  1:4791249‑M2/12‑1 (MQ=255)
            tttttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagc                                        >  1:4965143‑M2/38‑49 (MQ=255)
             ttttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagc                                        >  1:839486‑M2/37‑48 (MQ=255)
              tttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaag                                     >  1:733822‑M2/36‑50 (MQ=255)
              tttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaa                                      <  1:5522189‑M2/14‑1 (MQ=255)
              tttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagca                                       <  1:6929762‑M2/13‑1 (MQ=255)
              tttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagca                                       <  1:4539387‑M2/13‑1 (MQ=255)
              tttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagca                                       <  1:3579494‑M2/13‑1 (MQ=255)
                tcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaagc                                    >  1:4700870‑M2/34‑49 (MQ=255)
                 cactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaagca                                   >  1:205501‑M2/33‑49 (MQ=255)
                  actgagttattaaacatactcgcgagcgcgtCGCTACGCTAgtaagcaa                                  <  1:1914154‑M2/18‑1 (MQ=255)
                  actgagttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCaaa                                 <  1:4149310‑M2/19‑1 (MQ=255)
                    tgagttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCaaaaa                               >  1:4333627‑M2/30‑50 (MQ=255)
                    tgagttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCaaaaa                               >  1:958156‑M2/30‑50 (MQ=255)
                      agttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCaaaaa                               <  1:2776246‑M2/21‑1 (MQ=255)
                      agttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCaaaaa                               <  1:3725131‑M2/21‑1 (MQ=255)
                      agttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAAt                              >  1:4887557‑M2/28‑49 (MQ=255)
                      agttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAAt                              >  1:2614253‑M2/28‑49 (MQ=255)
                      agttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAAt                              <  1:7800654‑M2/22‑1 (MQ=255)
                      agttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATg                             <  1:5736407‑M2/23‑1 (MQ=255)
                       gttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATg                             >  1:3523790‑M2/27‑49 (MQ=255)
                       gttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGa                            <  1:7990440‑M2/24‑1 (MQ=255)
                       gttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGa                            <  1:8217031‑M2/24‑1 (MQ=255)
                       gttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGa                            <  1:2200894‑M2/24‑1 (MQ=255)
                         tattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGaaa                          >  1:7720567‑M2/25‑50 (MQ=255)
                         tattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGaaa                          >  1:2975707‑M2/25‑50 (MQ=255)
                          attaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGaaa                          <  1:2231659‑M2/26‑1 (MQ=255)
                          attaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGaaa                          >  1:3123486‑M2/24‑49 (MQ=255)
                           ttaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGaaa                          >  1:4318594‑M2/23‑48 (MQ=255)
                           ttaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACa                        <  1:437740‑M2/28‑1 (MQ=255)
                            taaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAg                       >  1:1562760‑M2/22‑50 (MQ=255)
                             aaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAg                       >  1:1821968‑M2/21‑49 (MQ=255)
                             aaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAg                       >  1:6586208‑M2/21‑49 (MQ=255)
                             aaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAg                       >  1:6208839‑M2/21‑49 (MQ=255)
                              aacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAgg                      <  1:291808‑M2/30‑1 (MQ=255)
                               acatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaa                    <  1:3236805‑M2/32‑1 (MQ=255)
                               acatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaa                    <  1:2905329‑M2/32‑1 (MQ=255)
                               acatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaa                    <  1:3619021‑M2/32‑1 (MQ=255)
                                 atactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaata                  <  1:6819915‑M2/34‑1 (MQ=255)
                                 atactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaata                  <  1:5042128‑M2/34‑1 (MQ=255)
                                 atactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaat                   >  1:5445323‑M2/17‑49 (MQ=255)
                                 atactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaat                   >  1:3769844‑M2/17‑49 (MQ=255)
                                 atactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaat                   >  1:2864444‑M2/17‑49 (MQ=255)
                                  tactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaataa                 >  1:4808318‑M2/16‑50 (MQ=255)
                                  tactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaataa                 <  1:3980565‑M2/35‑1 (MQ=255)
                                    ctcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaataat                <  1:7156728‑M2/36‑1 (MQ=255)
                                    ctcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaataat                <  1:3487256‑M2/36‑1 (MQ=255)
                                    ctcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaataat                <  1:3848856‑M2/36‑1 (MQ=255)
                                    ctcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc               >  1:1046024‑M2/14‑50 (MQ=255)
                                    ctcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATc               <  1:3872862‑M2/37‑1 (MQ=255)
                                      cgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCg              >  1:4742814‑M2/12‑49 (MQ=255)
                                       gcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGa             >  1:1307383‑M2/11‑49 (MQ=255)
                                       gcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGa             >  1:7981426‑M2/11‑49 (MQ=255)
                                        cgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa           >  1:8168793‑M2/10‑50 (MQ=255)
                                        cgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa           <  1:5072135‑M2/41‑1 (MQ=255)
                                        cgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaaa           >  1:659237‑M2/10‑50 (MQ=255)
                                         gagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGTAATAATCCAAAt          >  1:4519840‑M2/9‑50 (MQ=255)
                                         gagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt          >  1:6163313‑M2/9‑50 (MQ=255)
                                         gagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAAt          <  1:6229728‑M2/42‑1 (MQ=255)
                                           gcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATgg        <  1:1661533‑M2/44‑1 (MQ=255)
                                            cgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATggg       <  1:6123889‑M2/45‑1 (MQ=255)
                                            cgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATggg       <  1:4556443‑M2/45‑1 (MQ=255)
                                            cgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATgg        >  1:8218394‑M2/6‑49 (MQ=255)
                                             gcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGa      <  1:5255256‑M2/46‑1 (MQ=255)
                                              cgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt     <  1:2953631‑M2/47‑1 (MQ=255)
                                              cgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt     <  1:60217‑M2/47‑1 (MQ=255)
                                              cgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt     <  1:3884146‑M2/47‑1 (MQ=255)
                                              cgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt     <  1:4224711‑M2/47‑1 (MQ=255)
                                              cgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt     <  1:3143435‑M2/47‑1 (MQ=255)
                                              cgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGAt     <  1:6753229‑M2/47‑1 (MQ=255)
                                               gtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATg    <  1:4551534‑M2/48‑1 (MQ=255)
                                               gtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATg    <  1:1744820‑M2/48‑1 (MQ=255)
                                               gtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATg    <  1:7674442‑M2/48‑1 (MQ=255)
                                               gtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATg    >  1:7604551‑M2/3‑50 (MQ=255)
                                               gtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATg    >  1:6948001‑M2/3‑50 (MQ=255)
                                                tCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGt   <  1:3065209‑M2/49‑1 (MQ=255)
                                                tCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGt   >  1:4602235‑M2/2‑50 (MQ=255)
                                                 cgctacgctaGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATGtt  >  1:3941854/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 < 33 ≤ 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.