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 44 [0] [0] 44 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)44 (1.480) 32/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

AAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATG  >  NC_000913/1979282‑1979377
                                                |                                               
gacatcgagcctttttcactgagttattaaacatactcgcgagcgcgtcg                                                <  1:1186659‑M2/2‑1 (MQ=255)
   atcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgcta                                             <  1:1118940‑M2/5‑1 (MQ=255)
    tcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgctac                                            <  1:463101‑M2/6‑1 (MQ=255)
    tcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgctac                                            <  1:2613259‑M2/6‑1 (MQ=255)
    tcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgctac                                            <  1:2593871‑M2/6‑1 (MQ=255)
    tcgagcctttttcactgagttattaaacatactcgcgagcgcgtcgctac                                            <  1:1723816‑M2/6‑1 (MQ=255)
     cgagcctttttcactgagttattaaacatactcgcgagcgcgtcgctacg                                           >  1:891163‑M2/44‑50 (MQ=255)
         cctttttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTag                                       <  1:1023175‑M2/11‑1 (MQ=255)
         cctttttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTag                                       <  1:507581‑M2/11‑1 (MQ=255)
          ctttttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagc                                      <  1:2276801‑M2/12‑1 (MQ=255)
          ctttttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTag                                       >  1:1961909‑M2/39‑49 (MQ=255)
              ttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaag                                   >  1:565541‑M2/35‑49 (MQ=255)
              ttcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaag                                   >  1:670047‑M2/35‑49 (MQ=255)
               tcactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaagc                                  <  1:2505840‑M2/16‑1 (MQ=255)
                cactgagttattaaacatactcgcgagcgcgtCGCTACGCTagcaagca                                 >  1:220484‑M2/33‑49 (MQ=255)
                 actgagttattaaacatactcgcgagcgcgtCGCTACGCTAgcaagcaa                                <  1:827444‑M2/18‑1 (MQ=255)
                   tgagttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCaaaaa                             >  1:466251‑M2/30‑50 (MQ=255)
                     agttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAAt                            <  1:2452560‑M2/22‑1 (MQ=255)
                     agttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATg                           <  1:2222362‑M2/23‑1 (MQ=255)
                      gttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATg                           >  1:399422‑M2/27‑49 (MQ=255)
                       ttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGa                          >  1:1437363‑M2/26‑49 (MQ=255)
                       ttattaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGa                          >  1:627496‑M2/26‑49 (MQ=255)
                          ttaaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACa                      <  1:948295‑M2/28‑1 (MQ=255)
                            aaacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAg                     >  1:1533281‑M2/21‑49 (MQ=255)
                             aacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAgg                    >  1:41809‑M2/20‑49 (MQ=255)
                             aacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAgg                    >  1:2679344‑M2/20‑49 (MQ=255)
                             aacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAgg                    >  1:432989‑M2/20‑49 (MQ=255)
                             aacatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGa                   <  1:686332‑M2/31‑1 (MQ=255)
                              acatactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaa                  <  1:2153649‑M2/32‑1 (MQ=255)
                               catactcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaat                 <  1:1527865‑M2/33‑1 (MQ=255)
                                  actcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaataat              >  1:1420286‑M2/15‑50 (MQ=255)
                                  actcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaataat              >  1:2115614‑M2/15‑50 (MQ=255)
                                   ctcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaataat              <  1:73133‑M2/36‑1 (MQ=255)
                                   ctcgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGaataat              <  1:106‑M2/36‑1 (MQ=255)
                                     cgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCg            >  1:2189809‑M2/12‑49 (MQ=255)
                                     cgcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCg            >  1:120991‑M2/12‑49 (MQ=255)
                                      gcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCg            <  1:496061‑M2/38‑1 (MQ=255)
                                      gcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGaa          <  1:2519725‑M2/40‑1 (MQ=255)
                                      gcgagcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGa           >  1:2391689‑M2/11‑49 (MQ=255)
                                         agcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATg       <  1:354665‑M2/43‑1 (MQ=255)
                                         agcgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATg       <  1:164682‑M2/43‑1 (MQ=255)
                                           cgcgtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATggg     <  1:314876‑M2/45‑1 (MQ=255)
                                              gtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATg  <  1:67185‑M2/48‑1 (MQ=255)
                                              gtCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATg  >  1:244776‑M2/3‑50 (MQ=255)
                                                |                                               
AAGTCATACTTTTGTTTTGGGTGTATTTCAATCTGTTAAAAAGTTTTTCGCTACGCTAGCAAGCAAAAATGAAACAGGAATAATCGAAATGGGATG  >  NC_000913/1979282‑1979377

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: