Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 1,129,397 Δ11,591 bp flgNflgL flgN, flgM, flgA, flgB, flgC, flgD, flgE, flgF, flgG, flgH, flgI, flgJ, flgK, flgL

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1129397 1140987 11591 62 [0] [0] 67 flgN–flgL flgN,flgM,flgA,flgB,flgC,flgD,flgE,flgF,flgG,flgH,flgI,flgJ,flgK,flgL

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 11293960 (0.000)61 (0.840) 35/68 0.1 100% intergenic (+22/+18) murJ/flgN putative lipid II flippase MurJ/flagellar biosynthesis protein FlgN
?NC_000913 1140988 = 0 (0.000)intergenic (+2/+194) flgL/rne flagellar hook‑filament junction protein 2/ribonuclease E

AGGGATTGTCGCTCTTCCAGCTCAGCAAATAATTTCGCTTTAAAACATATCATGAAACTGGGTATGTT  >  NC_000913/1140955‑1141022
                                 |                                  
ggacggtgtaacaatgcattccggcctgcagtgtt                                   <  1:7791102‑M2/2‑1 (MQ=255)
ggacggtgtaacaatgcattccggcctgcagtgtt                                   <  1:5061211‑M2/2‑1 (MQ=255)
 gacggtgtaacaatgcattccggcctgcagtgTTc                                  <  1:7743311‑M2/3‑1 (MQ=255)
  acggtgtaacaatgcattccggcctgcagtgTTCg                                 <  1:1598681‑M2/4‑1 (MQ=255)
  acggtgtaacaatgcattccggcctgcagtgTTCg                                 <  1:2052799‑M2/4‑1 (MQ=255)
   cggtgtaacaatgcattccggcctgcagtgTTCGc                                <  1:7552983‑M2/5‑1 (MQ=255)
   cggtgtaacaatgcattccggcctgcagtgTTCGc                                <  1:6086569‑M2/5‑1 (MQ=255)
    ggtgtaacaatgcattccggcctgcagtgTTCGCt                               <  1:5320001‑M2/6‑1 (MQ=255)
    ggtgtaacaatgcattccggcctgcagtgTTCGCt                               <  1:2723992‑M2/6‑1 (MQ=255)
     gtgtaacaatgcattccggcctgcagtgTTCGCtt                              >  1:1019237‑M2/29‑35 (MQ=255)
     gtgtaacaatgcattccggcctgcagtgTTCGCtt                              >  1:5613085‑M2/29‑35 (MQ=255)
     gtgtaacaatgcattccggcctgcagtgTTCGCtt                              >  1:2250341‑M2/29‑35 (MQ=255)
     gtgtaacaatgcattccggcctgcagtgTTCGCtt                              <  1:2588676‑M2/7‑1 (MQ=255)
      tgtaacaatgcattccggcctgcagtgTTCGCttt                             >  1:5824471‑M2/28‑35 (MQ=255)
      tgtaacaatgcattccggcctgcagtgTTCGCttt                             <  1:6687848‑M2/8‑1 (MQ=255)
       gtaacaatgcattccggcctgcagtgTTCGCTTTa                            <  1:5702684‑M2/9‑1 (MQ=255)
       gtaacaatgcattccggcctgcagtgTTCGCTTTa                            >  1:2142607‑M2/27‑35 (MQ=255)
       gtaacaatgcattccggcctgcagtgTTCGCTTTa                            <  1:5376991‑M2/9‑1 (MQ=255)
        taacaatgcattccggcctgcagtgTTCGCTTTaa                           >  1:5699387‑M2/26‑35 (MQ=255)
        taacaatgcattccggcctgcagtgTTCGCTTTaa                           >  1:6957876‑M2/26‑35 (MQ=255)
         aacaatgcattccggcctgcagtgTTCGCTTTaaa                          <  1:5403484‑M2/11‑1 (MQ=255)
          acaatgcattccggcctgcagtgTTCGCTTTaaaa                         <  1:1535590‑M2/12‑1 (MQ=255)
          acaatgcattccggcctgcagtgTTCGCTTTaaaa                         >  1:4695271‑M2/24‑35 (MQ=255)
          acaatgcattccggcctgcagtgTTCGCTTTaaaa                         <  1:4871983‑M2/12‑1 (MQ=255)
            aatgcattccggcctgcagtgTTCGCTTTAAAACa                       <  1:3746741‑M2/14‑1 (MQ=255)
            aatgcattccggcctgcagtgTTCGCTTTAAAACa                       <  1:2205293‑M2/14‑1 (MQ=255)
            aatgcattccggcctgcagtgTTCGCTTTAAAACa                       >  1:5276355‑M2/22‑35 (MQ=255)
            aatgcattccggcctgcagtgTTCGCTTTAAAACa                       <  1:4812090‑M2/14‑1 (MQ=255)
             atgcattccggcctgcagtgTTCGCTTTAAAACat                      >  1:612619‑M2/21‑35 (MQ=255)
             atgcattccggcctgcagtgTTCGCTTTAAAACat                      <  1:734778‑M2/15‑1 (MQ=255)
             atgcattccggcccgcagtgTTCGCTTTAAAACat                      >  1:5436419‑M2/21‑35 (MQ=255)
              tgcattccggcctgcagtgTTCGCTTTAAAACata                     <  1:7174007‑M2/16‑1 (MQ=255)
               gcattccggcctgcagtgTTCGCTTTAAAACatat                    <  1:3744842‑M2/17‑1 (MQ=255)
               gcattccggcctgcagtgTTCGCTTTAAAACatat                    <  1:3013696‑M2/17‑1 (MQ=255)
               gcattccggcctgcagtgTTCGCTTTAAAACatat                    <  1:5722151‑M2/17‑1 (MQ=255)
                cattccggcctgcagtgTTCGCTTTAAAACATATc                   <  1:689185‑M2/18‑1 (MQ=255)
                cattccggcctgcagtgTTCGCTTTAAAACATATc                   <  1:2960151‑M2/18‑1 (MQ=255)
                cattccggcctgcagtgTTCGCTTTAAAACATATc                   >  1:2838053‑M2/18‑35 (MQ=255)
                 attccggcctgcagtgTTCGCTTTAAAACATATCa                  >  1:5027586‑M2/17‑35 (MQ=255)
                  ttccggcctgcagtgTTCGCTTTAAAACATATCAt                 <  1:1762601‑M2/20‑1 (MQ=255)
                  ttccggcctgcagtgTTCGCTTTAAAACATATCAt                 <  1:2803789‑M2/20‑1 (MQ=255)
                  ttccggcctgcagtgTTCGCTTTAAAACATATCAt                 <  1:5564375‑M2/20‑1 (MQ=255)
                  ttccggcctgcagtgTTCGCTTTAAAACATATCAt                 <  1:7507863‑M2/20‑1 (MQ=255)
                  ttccggcctgcagtgTTCGCTTTAAAACATATCAt                 <  1:1574455‑M2/20‑1 (MQ=255)
                    ccggcctgcagtgTTCGCTTTAAAACATATCATGa               <  1:1659058‑M2/22‑1 (MQ=255)
                    ccggcctgcagtgTTCGCTTTAAAACATATCATGa               <  1:3714726‑M2/22‑1 (MQ=255)
                    ccggcctgcagtgTTCGCTTTAAAACATATCATGa               <  1:6402101‑M2/22‑1 (MQ=255)
                    ccggcctgcagtgTTCGCTTTAAAACATATCATGa               <  1:6860191‑M2/22‑1 (MQ=255)
                       gcctgcagtgTTCGCTTTAAAACATATCATGAAAc            >  1:5523673‑M2/11‑35 (MQ=255)
                           gcagtgTTCGCTTTAAAACATATCATGAAACTggg        <  1:1324833‑M2/29‑1 (MQ=255)
                           gcagtgTTCGCTTTAAAACATATCATGAAACTggg        <  1:621254‑M2/29‑1 (MQ=255)
                            cagtgTTCGCTTTAAAACATATCATGAAACTGGGt       >  1:7083294‑M2/6‑35 (MQ=255)
                             agtgTTCGCTTTAAAACATATCATGAAACTGGGTa      <  1:6513174‑M2/31‑1 (MQ=255)
                             agtgTTCGCTTTAAAACATATCATGAAACTGGGTa      >  1:4191625‑M2/5‑35 (MQ=255)
                             agtgTTCGCTTTAAAACATATCATGAAACTGGGTa      <  1:7576940‑M2/31‑1 (MQ=255)
                              gtgTTCGCTTTAAAACATATCATGAAACTGGGTAt     <  1:6975440‑M2/32‑1 (MQ=255)
                              gtgTTCGCTTTAAAACATATCATGAAACTGGGTAt     >  1:4223830‑M2/4‑35 (MQ=255)
                               tgTTCGCTTTAAAACATATCATGAAACTGGGTATg    <  1:6088663‑M2/33‑1 (MQ=255)
                               tgTTCGCTTTAAAACATATCATGAAACTGGGTATg    <  1:7258682‑M2/33‑1 (MQ=255)
                               tgTTCGCTTTAAAACATATCATGAAACTGGGTATg    <  1:2478567‑M2/33‑1 (MQ=255)
                                gTTCGCTTTAAAACATATCATGAAACTGGGTATGt   <  1:6812949‑M2/34‑1 (MQ=255)
                                 ttCGCTTTAAAACATATCATGAAACTGGGTATGtt  >  1:6015666/1‑35 (MQ=255)
                                 ttCGCTTTAAAACATATCATGAAACTGGGTATGtt  >  1:374674/1‑35 (MQ=255)
                                 ttCGCTTTAAAACATATCATGAAACTGGGTATGtt  >  1:7731018/1‑35 (MQ=255)
                                 ttCGCTTTAAAACATATCATGAAACTGGGTATGtt  <  1:3199254/35‑1 (MQ=255)
                                 ttCGCTTTAAAACATATCATGAAACTGGGTATGtt  >  1:6817851/1‑35 (MQ=255)
                                 ttCGCTTTAAAACATATCATGAAACTGCGTATGtt  >  1:695538/1‑35 (MQ=255)
                                 |                                  
AGGGATTGTCGCTCTTCCAGCTCAGCAAATAATTTCGCTTTAAAACATATCATGAAACTGGGTATGTT  >  NC_000913/1140955‑1141022

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