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 59 [0] [0] 60 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)59 (1.020) 45/98 0.0 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

AAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGG  >  NC_000913/1129348‑1129445
                                                |                                                 
aaaGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGt                                                  >  1:4218480‑M1/1‑49 (MQ=255)
 aaGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGtt                                                 >  1:554535‑M1/1‑48 (MQ=255)
  aGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGttc                                                >  1:3856883‑M1/1‑47 (MQ=255)
    aaTTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGttcgc                                              <  1:2549721‑M1/50‑6 (MQ=255)
      tttGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGttcgctt                                            <  1:2171193‑M1/50‑8 (MQ=255)
          cccGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGttcgctttaaa                                        <  1:2500822‑M1/50‑12 (MQ=255)
          cccGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGttcgctttaa                                         >  1:3237527‑M1/1‑39 (MQ=255)
          cccGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGttcgctttaa                                         <  1:3044494‑M1/49‑11 (MQ=255)
          cccGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGttcgctttaa                                         <  1:291403‑M1/49‑11 (MQ=255)
          cccGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGttcgctttaa                                         <  1:2901911‑M1/49‑11 (MQ=255)
          cccGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGttcgctttaa                                         >  1:211120‑M1/1‑39 (MQ=255)
           ccgccgGACGGTGTAACAATGCATTCCGGCCTGCAGTGttcgctttaaaa                                       <  1:3804433‑M1/50‑13 (MQ=255)
            cgccgGACGGTGTAACAATGCATTCCGGCCTGCAGTGttcgctttaaaa                                       >  1:1630196‑M1/1‑37 (MQ=255)
            cgccgGACGGTGTAACAATGCATTCCGGCCTGCAGTGttcgctttaaaa                                       >  1:3104290‑M1/1‑37 (MQ=255)
             gccgGACGGTGTAACAATGCATTCCGGCCTGCAGTGttcgctttaaaaca                                     >  1:4063823‑M1/1‑36 (MQ=255)
              ccgGACGGTGTAACAATGCATTCCGGCCTGCAGTGttcgctttaaaacat                                    <  1:2021425‑M1/50‑16 (MQ=255)
              ccgGACGGTGTAACAATGCATTCCGGCCTGCAGTGttcgctttaaaacat                                    <  1:3542432‑M1/50‑16 (MQ=255)
              ccgGACGGTGTAACAATGCATTCCGGCCTGCAGTGttcgctttaaaacat                                    <  1:4314759‑M1/50‑16 (MQ=255)
              ccgGACGGTGTAACAATGCATTCCGGCCTGCAGTGttcgctttaaaaca                                     >  1:319674‑M1/1‑35 (MQ=255)
               cgGACGGTGTAACAATGCATTCCGGCCTGCAGTGttcgctttaaaacata                                   <  1:3924127‑M1/50‑17 (MQ=255)
               cgGACGGTGTAACAATGCATTCCGGCCTGCAGTGttcgctttaaaacata                                   <  1:966403‑M1/50‑17 (MQ=255)
               cgGACGGTGTAACAATGCATTCCGGCCTGCAGTGttcgctttaaaacata                                   >  1:2257090‑M1/1‑34 (MQ=255)
                ggACGGTGTAACAATGCATTCCGGCCTGCAGTGttcgctttaaaacatat                                  <  1:1346318‑M1/50‑18 (MQ=255)
                 gACGGTGTAACAATGCATTCCGGCCTGCAGTGttcgctttaaaacatat                                  >  1:793959‑M1/1‑32 (MQ=255)
                  aCGGTGTAACAATGCATTCCGGCCTGCAGTGttcgctttaaaacatatca                                >  1:1242902‑M1/1‑31 (MQ=255)
                    ggTGTAACAATGCATTCCGGCCTGCAGTGttcgctttaaaacatatcat                               >  1:2876875‑M1/1‑29 (MQ=255)
                     gtgtAACAATGCATTCCGGCCTGCAGTGttcgctttaaaacatatcatga                             <  1:1737416‑M1/50‑23 (MQ=255)
                     gtgtAACAATGCATTCCGGCCTGCAGTGttcgctttaaaacatatcatg                              >  1:3722045‑M1/1‑28 (MQ=255)
                       gtAACAATGCATTCCGGCCTGCAGTGttcgctttaaaacatatcatgaaa                           >  1:26481‑M1/1‑26 (MQ=255)
                        tAACAATGCATTCCGGCCTGCAGTGttcgctttaaaacatatcatgaaac                          <  1:1016923‑M1/50‑26 (MQ=255)
                          aCAATGCATTCCGGCCTGCAGTGttcgctttaaaacatatcatgaacctg                        >  1:2812049‑M1/1‑23 (MQ=255)
                          aCAATGCATTCCGGCCTGCAGTGttcgctttaaaacatatcatgaaactg                        <  1:385004‑M1/50‑28 (MQ=255)
                          aCAATGCATTCCGGCCTGCAGTGttcgctttaaaacatatcatgaaact                         <  1:3503788‑M1/49‑27 (MQ=255)
                          aCAATGCATTCCGGCCTGCAGTGttcgctttaaaacatatcatgaaact                         <  1:3118430‑M1/49‑27 (MQ=255)
                           cAATGCATTCCGGCCTGCAGTGttcgctttaaaacatatcatgaaactgg                       <  1:323640‑M1/50‑29 (MQ=255)
                            aaTGCATTCCGGCCTGCAGTGttcgctttaaaacatatcatgaaactggg                      >  1:838036‑M1/1‑21 (MQ=255)
                            aaTGCATTCCGGCCTGCAGTGttcgctttaaaacatatcatgaaactggg                      <  1:3158204‑M1/50‑30 (MQ=255)
                            aaTGCATTCCGGCCTGCAGTGttcgctttaaaacatatcatgaaactggg                      <  1:2393039‑M1/50‑30 (MQ=255)
                             aTGCATTCCGGCCTGCAGTGttcgctttaaaacatatcatgaaactggg                      >  1:957723‑M1/1‑20 (MQ=255)
                               gCATTCCGGCCTGCAGTGttcgctttaaaacatatcatgaaactgggtat                   >  1:1667206‑M1/1‑18 (MQ=255)
                                cATTCCGGCCTGCAGTGttcgctttaaaacatatcatgaaactgggtatg                  <  1:3563635‑M1/50‑34 (MQ=255)
                                 aTTCCGGCCTGCTGTGttcgctttaaaacatatcatgaaactgggtatgc                 >  1:1054820‑M1/1‑16 (MQ=255)
                                    ccGGCCTGCAGTGttcgctttaaaacatatcatgaaactgggtatgtttt              >  1:3335821‑M1/1‑13 (MQ=255)
                                    ccGGCCTGCAGTGttcgctttaaaacatatcatgaaactgggtatgtttt              >  1:3412202‑M1/1‑13 (MQ=255)
                                    ccGGCCTGCAGTGttcgctttaaaacatatcatgaaactgggtatgttt               <  1:3982638‑M1/49‑37 (MQ=255)
                                      ggCCTGCAGTGttcgctttaaaacatatcatgaaactgggtatgttttgt            >  1:2572081‑M1/1‑11 (MQ=255)
                                      ggCCTGCAGTGttcgctttaaaacatatcatgaaactgggtatgttttgt            >  1:1880141‑M1/1‑11 (MQ=255)
                                      ggCCTGCAGTGttcgctttaaaacatatcatgaaactgggtatgttttgt            <  1:1710711‑M1/50‑40 (MQ=255)
                                        ccTGCAGTGttcgctttaaaacatatcatgaaactgggtatgttttgtct          >  1:1528229‑M1/1‑9 (MQ=255)
                                        ccTGCAGTGttcgctttaaaacatatcatgaaactgggtatgttttgtct          <  1:942411‑M1/50‑42 (MQ=255)
                                           gCAGTGttcgctttaaaacatatcatgaaactgggtatgttttgtctgcc       >  1:4067980‑M1/1‑6 (MQ=255)
                                           gCAGTGttcgctttaaaacatatcatgaaactgggtatgttttgtctgcc       >  1:787557‑M1/1‑6 (MQ=255)
                                           gCAGTGttcgctttaaaacatatcatgaaactgggtatgttttgtctgcc       >  1:3490749‑M1/1‑6 (MQ=255)
                                            cAGTGttcgctttaaaacatatcatgaaactgggtatgttttgtctgcct      <  1:3334679‑M1/50‑46 (MQ=255)
                                            cAGTGttcgctttaaaacatatcatgaaactgggtatgttttgtctgcc       <  1:1982035‑M1/49‑45 (MQ=255)
                                            cAGTGttcgctttaaaacatatcatgaaactgggtatgttttgtctgc        <  1:2026799‑M1/48‑44 (MQ=255)
                                              gtgttcgctttaaaacatatcatggaactgggtatgttttgtctgcctgc    <  1:474071‑M1/50‑48 (MQ=255)
                                               tgttcgctttaaaacatatcatgaaactgggtatgttttgtctgcctgct   >  1:2942966‑M1/1‑2 (MQ=255)
                                                gttcgctttaaaacatatcatgaaactgggtatgttttgtctgcctgctc  >  1:2815928‑M1/1‑1 (MQ=255)
                                                |                                                 
AAAGAATTTGCCCGCCGGACGGTGTAACAATGCATTCCGGCCTGCAGTGCAGGCCGGAGATAATCTTCAGATCGAAATCTTTTTACCGCCGCGATGGG  >  NC_000913/1129348‑1129445

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ 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.