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 41 [0] [0] 45 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)41 (0.800) 32/98 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

AAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGG  >  NC_000913/1129349‑1129444
                                               |                                                
aaGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGtt                                                >  1:3577774‑M1/1‑48 (MQ=255)
   aaTTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGttcgc                                             >  1:1741195‑M1/1‑45 (MQ=255)
     tttGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGttcgctt                                           <  1:3332906‑M1/50‑8 (MQ=255)
      ttGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGttcgctt                                           >  1:603034‑M1/1‑42 (MQ=255)
        gCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGttcgcttta                                         <  1:3337357‑M1/49‑10 (MQ=255)
        gCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGttcgcttta                                         <  1:1742398‑M1/49‑10 (MQ=255)
         cccGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGttcgctttaaa                                       <  1:528771‑M1/50‑12 (MQ=255)
         cccGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGttcgctttaaa                                       <  1:231747‑M1/50‑12 (MQ=255)
         cccGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGttcgctttaaa                                       <  1:3578385‑M1/50‑12 (MQ=255)
         cccGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGttcgctttaa                                        <  1:3439092‑M1/49‑11 (MQ=255)
           cgccgGACGGTGTAACAATGCATTCCGGCCTGCAGTGttcgctttaaaac                                     <  1:2083096‑M1/50‑14 (MQ=255)
             ccgGACGGTGTAACAATGCATTCCGGCCTGCAGTGttcgctttaaaacat                                   <  1:2930164‑M1/50‑16 (MQ=255)
              cgGACGGTGTAACAATGCATTCCGGCCTGCAGTGttcgctttaaaacata                                  >  1:3749897‑M1/1‑34 (MQ=255)
               ggACGGTGTAACAATGCATTCCGGCCTGCAGTGttcgctttaaaacatat                                 >  1:3381389‑M1/1‑33 (MQ=255)
               ggACGGTGTAACAATGCATTCCGGCCTGCAGTGttcgctttaaaacatat                                 <  1:3359240‑M1/50‑18 (MQ=255)
               ggACGGTGTAACAATGCATTCCGGCCTGCAGTGttcgctttaaaacatat                                 <  1:646103‑M1/50‑18 (MQ=255)
                   ggTGTAACAATGCATTCCGGCCTGCAGTGttcgctttaaaacatatcatg                             <  1:1786092‑M1/50‑22 (MQ=255)
                   ggTGTAACAATGCATTCCGGCCTGCAGTGttcgctttaaaacatatcatg                             <  1:820607‑M1/50‑22 (MQ=255)
                       tAACAATGCATTCCGGCCTGCAGTGttcgctttaaaacatatcatgaaac                         <  1:879857‑M1/50‑26 (MQ=255)
                       tAACAATGCATTCCGGCCTGCAGTGttcgctttaaaacatatcatgaaa                          <  1:82130‑M1/49‑25 (MQ=255)
                       tAACAATGCATTCCGGCCTGCAGTGttcgctttaaaacatatcatgaaa                          <  1:3159358‑M1/49‑25 (MQ=255)
                         aCAATGCATTCCGGCCTGCAGTGttcgctttaaaacatatcatgaaactg                       <  1:2680759‑M1/50‑28 (MQ=255)
                               cATTCCGGCCTGCAGTGttcgctttaaaacatatcatgaaactgggtat                  <  1:2684902‑M1/49‑33 (MQ=255)
                                aTTCCGGCCTGCAGTGttcgctttaaaacatatcatgaaactgggtatgt                >  1:2862278‑M1/1‑16 (MQ=255)
                                aTTCCGGCCTGCAGTGttcgctttaaaacatatcatgaaactgggtatgt                <  1:2280986‑M1/50‑35 (MQ=255)
                                 ttCCGGCCTGCAGTGttcgctttaaaacatatcatgaaactgggtatgtt               <  1:407068‑M1/50‑36 (MQ=255)
                                   ccGGCCTGCAGTGttcgctttaaaacatatcatgaaactgggtatgtttt             >  1:2827227‑M1/1‑13 (MQ=255)
                                   ccGGCCTGCAGTGttcgctttaaaacatatcatgaaactgggtatgtttt             >  1:188617‑M1/1‑13 (MQ=255)
                                     ggCCTGCAGTGttcgctttaaaacatatcatgaaactgggtatgttttgt           >  1:2760970‑M1/1‑11 (MQ=255)
                                     ggCCTGCAGTGttcgctttaaaacatatcatgaaactgggtatgttttgt           <  1:594949‑M1/50‑40 (MQ=255)
                                      gCCTGCAGTGttcgctttaaaacatatcatgaaactgggtatgttttgtc          <  1:2264615‑M1/50‑41 (MQ=255)
                                      gCCTGCAGTGttcgctttaaaacatatcatgaaactgggtatgttttgtc          <  1:3165250‑M1/50‑41 (MQ=255)
                                      gCCTGCAGTGttcgctttaaaacatatcatgaaactgggtatgttttgtc          >  1:133717‑M1/1‑10 (MQ=255)
                                        cTGCAGTGttcgctttaaaacatatcatgaaactgggtatgttttgtctg        >  1:1833699‑M1/1‑8 (MQ=255)
                                        cTGCAGTGttcgctttaaaacatatcatgaaactgggtatgttttgtctg        <  1:1170041‑M1/50‑43 (MQ=255)
                                         tGCAGTGttcgctttaaaacatatcatgaaactgggtatgttttgtctgc       >  1:1202169‑M1/1‑7 (MQ=255)
                                           cAGTGttcgctttaaaacatatcatgaaactgggtatgttttgtctgcct     >  1:618392‑M1/1‑5 (MQ=255)
                                           cAGTGttcgctttaaaacatatcatgaaactgggtatgttttgtctgcct     >  1:1184069‑M1/1‑5 (MQ=255)
                                           cAGTGttcgctttaaaacatatcatgaaactgggtatgttttgtctgcc      <  1:2119225‑M1/49‑45 (MQ=255)
                                            aGTGttcgctttaaaacatatcatgaaactgggtatgttttgtctgcctg    <  1:1355755‑M1/50‑47 (MQ=255)
                                              tgttcgctttaaaacatatcatgaaactgggtatgttttgtctgcctgct  >  1:2939963‑M1/1‑2 (MQ=255)
                                               |                                                
AAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGG  >  NC_000913/1129349‑1129444

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.