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 91 [0] [0] 93 [dgcJ] [dgcJ]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 18720290 (0.000)88 (1.080) 57/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

AGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTGCCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCG  >  NC_000913/1871980‑1872077
                                                 |                                                
aGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAG                                                  <  1:5510205‑M1/50‑1 (MQ=255)
  gCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGat                                                >  1:689415‑M1/1‑48 (MQ=255)
  gCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGat                                                <  1:3275130‑M1/50‑3 (MQ=255)
   ccGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatg                                               <  1:1241595‑M1/50‑4 (MQ=255)
   ccGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatg                                               >  1:6043998‑M1/1‑47 (MQ=255)
     gCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcg                                             >  1:2514144‑M1/1‑45 (MQ=255)
      cAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgt                                            <  1:4171847‑M1/50‑7 (MQ=255)
      cAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcg                                             >  1:4513614‑M1/1‑44 (MQ=255)
       aaTCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgta                                           >  1:6093451‑M1/1‑43 (MQ=255)
       aaTCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgta                                           >  1:4264372‑M1/1‑43 (MQ=255)
        aTCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtat                                          <  1:2354020‑M1/50‑9 (MQ=255)
        aTCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgta                                           >  1:584928‑M1/1‑42 (MQ=255)
          cGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtata                                         >  1:4177399‑M1/1‑40 (MQ=255)
           gAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaa                                       <  1:6386538‑M1/50‑12 (MQ=255)
           gAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaa                                       <  1:6066400‑M1/50‑12 (MQ=255)
           gAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaa                                       <  1:4136538‑M1/50‑12 (MQ=255)
           gAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaa                                       <  1:1781170‑M1/50‑12 (MQ=255)
           gAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaa                                       <  1:5694726‑M1/50‑12 (MQ=255)
           gAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaa                                       <  1:1056101‑M1/50‑12 (MQ=255)
           gAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtata                                         >  1:3580804‑M1/1‑39 (MQ=255)
           gAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtata                                         >  1:2699411‑M1/1‑39 (MQ=255)
             aCTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagc                                     <  1:5813098‑M1/50‑14 (MQ=255)
             aCTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagc                                     <  1:4762334‑M1/50‑14 (MQ=255)
             aCTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagc                                     <  1:2462285‑M1/50‑14 (MQ=255)
              cTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagct                                    <  1:3687179‑M1/50‑15 (MQ=255)
               tttATCTGGTTTTCTCGTTTCACTAACCGCAGGAGatgcgtataaagctt                                   >  1:3720672‑M1/1‑35 (MQ=255)
               tttATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagctt                                   >  1:1989325‑M1/1‑35 (MQ=255)
               tttATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagctt                                   >  1:6226600‑M1/1‑35 (MQ=255)
               tttATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagctt                                   >  1:4705668‑M1/1‑35 (MQ=255)
                  aTCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttcc                                 <  1:2203392‑M1/49‑18 (MQ=255)
                   tCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccga                               >  1:5891986‑M1/1‑31 (MQ=255)
                   tCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccga                               >  1:269014‑M1/1‑31 (MQ=255)
                      ggTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgat                              >  1:6334966‑M1/1‑28 (MQ=255)
                       gTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgag                           <  1:3045022‑M1/50‑24 (MQ=255)
                       gTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgag                           <  1:5179618‑M1/50‑24 (MQ=255)
                       gTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatga                            <  1:3064168‑M1/49‑23 (MQ=255)
                       gTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatga                            <  1:6217784‑M1/49‑23 (MQ=255)
                       gTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgat                              <  1:6382647‑M1/47‑21 (MQ=255)
                       gTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccga                               <  1:5400380‑M1/46‑20 (MQ=255)
                        ttttCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagc                          <  1:4747574‑M1/50‑25 (MQ=255)
                        ttttCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagc                          <  1:466976‑M1/50‑25 (MQ=255)
                         tttCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcg                         >  1:6331455‑M1/1‑25 (MQ=255)
                         tttCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcg                         >  1:5333650‑M1/1‑25 (MQ=255)
                          ttCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgt                        >  1:1477168‑M1/1‑24 (MQ=255)
                          ttCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcg                         <  1:3434643‑M1/49‑26 (MQ=255)
                          ttCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcg                         <  1:3025297‑M1/49‑26 (MQ=255)
                           tctcGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgt                        <  1:4301554‑M1/49‑27 (MQ=255)
                            ctcGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtt                       <  1:1679030‑M1/49‑28 (MQ=255)
                            ctcGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtt                       >  1:4775433‑M1/1‑22 (MQ=255)
                            ctcGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtt                       <  1:35734‑M1/49‑28 (MQ=255)
                            ctcGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtt                       >  1:3511847‑M1/1‑22 (MQ=255)
                             tcGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgttta                     >  1:589611‑M1/1‑21 (MQ=255)
                               gTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttata                   >  1:543079‑M1/1‑19 (MQ=255)
                               gTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttata                   <  1:5523941‑M1/50‑32 (MQ=255)
                               gTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttata                   >  1:1317684‑M1/1‑19 (MQ=255)
                               gTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttata                   >  1:2592691‑M1/1‑19 (MQ=255)
                               gTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttat                    <  1:5416123‑M1/49‑31 (MQ=255)
                               gTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgttta                     <  1:654479‑M1/48‑30 (MQ=255)
                                tttCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatat                  >  1:2585313‑M1/1‑18 (MQ=255)
                                 ttCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatg                 >  1:1494721‑M1/1‑17 (MQ=255)
                                 ttCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatg                 >  1:2432311‑M1/1‑17 (MQ=255)
                                  tCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgt                >  1:2330345‑M1/1‑16 (MQ=255)
                                   cACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgt                >  1:4571345‑M1/1‑15 (MQ=255)
                                   cACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgt                >  1:1825289‑M1/1‑15 (MQ=255)
                                   cACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgt                >  1:5244943‑M1/1‑15 (MQ=255)
                                    aCTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtca              <  1:3625999‑M1/50‑37 (MQ=255)
                                    aCTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtc               <  1:451589‑M1/49‑36 (MQ=255)
                                    aCTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtc               >  1:2215334‑M1/1‑14 (MQ=255)
                                    aCTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtc               >  1:3229175‑M1/1‑14 (MQ=255)
                                    aCTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtc               >  1:281762‑M1/1‑14 (MQ=255)
                                     cTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaa             >  1:5757062‑M1/1‑13 (MQ=255)
                                      tAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaa             >  1:1137765‑M1/1‑12 (MQ=255)
                                       aaCCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaat            >  1:1837659‑M1/1‑11 (MQ=255)
                                        aCCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataa          <  1:4219715‑M1/50‑41 (MQ=255)
                                        aCCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataa          <  1:1253157‑M1/50‑41 (MQ=255)
                                        aCCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataa          <  1:2092681‑M1/50‑41 (MQ=255)
                                        aCCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaat            >  1:4198497‑M1/1‑10 (MQ=255)
                                        aCCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaat            <  1:6262125‑M1/48‑39 (MQ=255)
                                         ccGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataag         <  1:820062‑M1/50‑42 (MQ=255)
                                         ccGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataag         >  1:6048522‑M1/1‑9 (MQ=255)
                                          cGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataagc        <  1:590051‑M1/50‑43 (MQ=255)
                                          cGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataag         >  1:2399301‑M1/1‑8 (MQ=255)
                                           gAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataag         >  1:4420459‑M1/1‑7 (MQ=255)
                                             aGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataagcaga     <  1:629101‑M1/50‑46 (MQ=255)
                                             aGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataagcag      >  1:3668994‑M1/1‑5 (MQ=255)
                                             aGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataagcag      >  1:1702372‑M1/1‑5 (MQ=255)
                                              ggAGatgcgtataaagcttccgatgagcgtttatatgtcaataagcagaa    <  1:1964321‑M1/50‑47 (MQ=255)
                                              ggAGatgcgtataaagcttccgatgagcgtttatatgtcaataagcagaa    >  1:40941‑M1/1‑4 (MQ=255)
                                              ggAGatgcgtataaagcttccgatgagcgtttatatgtcaataagcagaa    <  1:4586360‑M1/50‑47 (MQ=255)
                                                 gatgcgtataaagcttccgatgagcgtttatatgtcaataagcagaaca  >  1:2774127‑M1/1‑1 (MQ=255)
                                                 gatgcgtataaagcttccgatgagcgtttatatgtcaataagcagaaca  >  1:4113247‑M1/1‑1 (MQ=255)
                                                 |                                                
AGGCCGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGTGCCATTTATCATGAAATTGCACCATAGAATGCTCCGGCATTTTATCG  >  NC_000913/1871980‑1872077

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