Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 1,872,030 Δ1,437 bp [dgcJ] [dgcJ]

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 1872030 1873466 1437 30 [0] [0] 30 [dgcJ] [dgcJ]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 18720290 (0.000)30 (1.010) 27/96 0.0 100% intergenic (+169/‑12) cdgI/dgcJ putative c‑di‑GMP binding protein CdgI/putative diguanylate cyclase DgcJ
?NC_000913 1873467 = 0 (0.000)coding (1427/1491 nt) dgcJ putative diguanylate cyclase DgcJ

GCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTGCCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCG  >  NC_000913/1871982‑1872077
                                               |                                                
gCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGat                                                <  1:2622826‑M1/50‑3 (MQ=255)
 ccGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatg                                               >  1:27985‑M1/1‑47 (MQ=255)
   gCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcg                                             >  1:1587829‑M1/1‑45 (MQ=255)
     aaTCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgta                                           >  1:1390587‑M1/1‑43 (MQ=255)
      aTCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtat                                          <  1:1567560‑M1/50‑9 (MQ=255)
      aTCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgta                                           >  1:470261‑M1/1‑42 (MQ=255)
      aTCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgta                                           >  1:132712‑M1/1‑42 (MQ=255)
        cGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtata                                         >  1:930425‑M1/1‑40 (MQ=255)
          aaCTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaag                                      <  1:1938772‑M1/50‑13 (MQ=255)
           aCTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagc                                     <  1:2549356‑M1/50‑14 (MQ=255)
               tATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttcc                                 >  1:349626‑M1/1‑33 (MQ=255)
                aTCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttcc                                 <  1:1817914‑M1/49‑18 (MQ=255)
                  cTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgat                              >  1:884021‑M1/1‑30 (MQ=255)
                   tGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgat                              >  1:2014547‑M1/1‑29 (MQ=255)
                     gTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatga                            <  1:1988871‑M1/49‑23 (MQ=255)
                       tttCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcg                         >  1:2076982‑M1/1‑25 (MQ=255)
                       tttCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcg                         >  1:1123868‑M1/1‑25 (MQ=255)
                        ttCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgt                        >  1:2349208‑M1/1‑24 (MQ=255)
                        ttCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcg                         <  1:1080644‑M1/49‑26 (MQ=255)
                             gTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgttta                     <  1:1407335‑M1/48‑30 (MQ=255)
                              tttCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatat                  >  1:754477‑M1/1‑18 (MQ=255)
                               ttCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatg                 >  1:2617605‑M1/1‑17 (MQ=255)
                                 cACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgt                >  1:1250293‑M1/1‑15 (MQ=255)
                                  aCTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtc               >  1:2091203‑M1/1‑14 (MQ=255)
                                  aCTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtc               >  1:1557987‑M1/1‑14 (MQ=255)
                                    tAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaa             >  1:1773661‑M1/1‑12 (MQ=255)
                                     aaCCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaat            >  1:1007298‑M1/1‑11 (MQ=255)
                                           aGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataagcag      >  1:1069394‑M1/1‑5 (MQ=255)
                                             gagatgcgtataaagcttccgatgagcgtttatatgtcaataagcagaa    >  1:891890‑M1/1‑3 (MQ=255)
                                              agatgcgtataaagcttccgatgagcgtttatatgtcaataagcagaaca  <  1:1318055‑M1/50‑49 (MQ=255)
                                               |                                                
GCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTGCCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCG  >  NC_000913/1871982‑1872077

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: 
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.