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 100 [0] [0] 101 [dgcJ] [dgcJ]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 18720290 (0.000)98 (1.020) 64/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:5121473‑M1/50‑1 (MQ=255)
   ccGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatg                                               <  1:3455111‑M1/50‑4 (MQ=255)
   ccGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatg                                               >  1:6962254‑M1/1‑47 (MQ=255)
    cGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgc                                              <  1:3825912‑M1/50‑5 (MQ=255)
    cGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgc                                              >  1:2107444‑M1/1‑46 (MQ=255)
      cAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgt                                            <  1:6915215‑M1/50‑7 (MQ=255)
      cAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcg                                             >  1:3072998‑M1/1‑44 (MQ=255)
       aaTCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgta                                           <  1:4720623‑M1/50‑8 (MQ=255)
        aTCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtat                                          <  1:2155938‑M1/50‑9 (MQ=255)
        aTCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgta                                           >  1:5384723‑M1/1‑42 (MQ=255)
        aTCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgta                                           >  1:3966695‑M1/1‑42 (MQ=255)
         tCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtata                                         <  1:4321271‑M1/50‑10 (MQ=255)
         tCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtata                                         <  1:3801155‑M1/50‑10 (MQ=255)
         tCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtat                                          >  1:1783133‑M1/1‑41 (MQ=255)
          cGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtata                                         >  1:826212‑M1/1‑40 (MQ=255)
          cGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtata                                         >  1:2766055‑M1/1‑40 (MQ=255)
           gAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaa                                       <  1:6570824‑M1/50‑12 (MQ=255)
           gAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtata                                         >  1:6741536‑M1/1‑39 (MQ=255)
           gAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtata                                         >  1:7110945‑M1/1‑39 (MQ=255)
            aaCTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaag                                      <  1:1807153‑M1/50‑13 (MQ=255)
            aaCTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaag                                      <  1:4293637‑M1/50‑13 (MQ=255)
            aaCTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaag                                      <  1:5444827‑M1/50‑13 (MQ=255)
            aaCTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaag                                      >  1:2576127‑M1/1‑38 (MQ=255)
             aCTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagc                                     >  1:4327520‑M1/1‑37 (MQ=255)
             aCTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagc                                     >  1:4483931‑M1/1‑37 (MQ=255)
             aCTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagc                                     >  1:7379129‑M1/1‑37 (MQ=255)
              cTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagct                                    <  1:4735554‑M1/50‑15 (MQ=255)
              cTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagct                                    >  1:6581060‑M1/1‑36 (MQ=255)
              cTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagct                                    <  1:5316745‑M1/50‑15 (MQ=255)
               tttATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagctt                                   >  1:3183450‑M1/1‑35 (MQ=255)
               tttATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagctt                                   >  1:2176456‑M1/1‑35 (MQ=255)
                 tATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttc                                  <  1:6051554‑M1/49‑17 (MQ=255)
                 tATCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagctt                                   <  1:3612953‑M1/48‑16 (MQ=255)
                  aTCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttcc                                 <  1:2550569‑M1/49‑18 (MQ=255)
                   tCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccga                               >  1:5783579‑M1/1‑31 (MQ=255)
                   tCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccga                               >  1:1899849‑M1/1‑31 (MQ=255)
                   tCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccg                                <  1:6729016‑M1/49‑19 (MQ=255)
                   tCTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccg                                <  1:6613331‑M1/49‑19 (MQ=255)
                    cTGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgat                              >  1:3802843‑M1/1‑30 (MQ=255)
                     tGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgat                              >  1:3450617‑M1/1‑29 (MQ=255)
                     tGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgat                              >  1:4493896‑M1/1‑29 (MQ=255)
                     tGGTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgat                              >  1:6311966‑M1/1‑29 (MQ=255)
                      ggTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgat                              >  1:6108065‑M1/1‑28 (MQ=255)
                       gTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgag                           <  1:3785875‑M1/50‑24 (MQ=255)
                       gTTTTCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgag                           <  1:489194‑M1/50‑24 (MQ=255)
                          ttCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgt                        >  1:6900401‑M1/1‑24 (MQ=255)
                          ttCTCGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcg                         <  1:546444‑M1/49‑26 (MQ=255)
                           tctcGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtt                       >  1:5981721‑M1/1‑23 (MQ=255)
                           tctcGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgt                        <  1:223067‑M1/49‑27 (MQ=255)
                           tctcGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgt                        <  1:1228545‑M1/49‑27 (MQ=255)
                            ctcGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtt                       <  1:3503469‑M1/49‑28 (MQ=255)
                            ctcGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtt                       <  1:2302195‑M1/49‑28 (MQ=255)
                             tcGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgttta                     >  1:2203113‑M1/1‑21 (MQ=255)
                              cGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgttta                     >  1:5063545‑M1/1‑20 (MQ=255)
                              cGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgttta                     >  1:5438085‑M1/1‑20 (MQ=255)
                              cGTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgttta                     >  1:806291‑M1/1‑20 (MQ=255)
                               gTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttata                   <  1:3724052‑M1/50‑32 (MQ=255)
                               gTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgttta                     <  1:3893411‑M1/48‑30 (MQ=255)
                               gTTTCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgttt                      <  1:2840384‑M1/47‑29 (MQ=255)
                                tttCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatat                  >  1:2295121‑M1/1‑18 (MQ=255)
                                tttCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatat                  <  1:4690507‑M1/50‑33 (MQ=255)
                                tttCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatat                  <  1:1877967‑M1/50‑33 (MQ=255)
                                tttCACTAACCGAAGGAGatgcgtataaagcttccgatgagagtttatat                  >  1:3999532‑M1/1‑18 (MQ=255)
                                 ttCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatg                 >  1:1043578‑M1/1‑17 (MQ=255)
                                  tCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgt                >  1:1773448‑M1/1‑16 (MQ=255)
                                  tCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatg                 <  1:3890219‑M1/49‑34 (MQ=255)
                                  tCACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatg                 <  1:1059220‑M1/49‑34 (MQ=255)
                                   cACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgt                <  1:271704‑M1/49‑35 (MQ=255)
                                   cACTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgt                >  1:3343066‑M1/1‑15 (MQ=255)
                                    aCTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtca              <  1:6268149‑M1/50‑37 (MQ=255)
                                     cTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaa             <  1:2949900‑M1/50‑38 (MQ=255)
                                     cTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaa             <  1:3390793‑M1/50‑38 (MQ=255)
                                     cTAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaa             <  1:3498597‑M1/50‑38 (MQ=255)
                                      tAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaa             >  1:5577217‑M1/1‑12 (MQ=255)
                                      tAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaa             >  1:1606983‑M1/1‑12 (MQ=255)
                                      tAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaa             >  1:3446546‑M1/1‑12 (MQ=255)
                                      tAACCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaa             >  1:5162810‑M1/1‑12 (MQ=255)
                                       aaCCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaat            >  1:1414476‑M1/1‑11 (MQ=255)
                                       aaCCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaat            >  1:3405153‑M1/1‑11 (MQ=255)
                                        aCCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataa          <  1:1726567‑M1/50‑41 (MQ=255)
                                        aCCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataa          <  1:6637443‑M1/50‑41 (MQ=255)
                                        aCCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaata           <  1:7049109‑M1/49‑40 (MQ=255)
                                        aCCGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaat            >  1:752926‑M1/1‑10 (MQ=255)
                                         ccGAAGGAGatgcgtataaagcttccgttgagcgtttatatgtcaataag         >  1:6425499‑M1/1‑9 (MQ=255)
                                         ccGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataag         >  1:1158408‑M1/1‑9 (MQ=255)
                                         ccGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataag         <  1:628688‑M1/50‑42 (MQ=255)
                                          cGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataagc        <  1:2788461‑M1/50‑43 (MQ=255)
                                          cGAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataag         >  1:4540669‑M1/1‑8 (MQ=255)
                                           gAAGGAGatgtgtataaagcttccgatgagcgtttatatgtcaataag         >  1:5629426‑M1/1‑7 (MQ=255)
                                           gAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataagca       <  1:4551003‑M1/50‑44 (MQ=255)
                                           gAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataag         >  1:2060732‑M1/1‑7 (MQ=255)
                                           gAAGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataag         >  1:2906872‑M1/1‑7 (MQ=255)
                                             aGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataagcaga     <  1:2897216‑M1/50‑46 (MQ=255)
                                             aGGAGatgcgtataaagcttccgatgagcgtttatatgtcaataagcaga     <  1:3524014‑M1/50‑46 (MQ=255)
                                              ggAGatgcgtataaagcttccgatgagcgtttatatgtcaataagtagaa    >  1:2490784‑M1/1‑4 (MQ=255)
                                              ggAGatgcgtataaagcttccgatgagcgtttatatgtcaataagcagaa    >  1:6061762‑M1/1‑4 (MQ=255)
                                              ggAGatgcgtataaagcttccgatgagcgtttatatgtcaataagcagaa    <  1:1414527‑M1/50‑47 (MQ=255)
                                               gagatgcgtataaagcttccgatgagcgtttatatgtcaataagcagaac   <  1:582837‑M1/50‑48 (MQ=255)
                                               gagatgcgtataaagcttccgatgagcgtttatatgtcaataagcagaa    >  1:3215673‑M1/1‑3 (MQ=255)
                                                 gatgcgtataaagcttccgatgagcgtttatatgtcaataagcagaaca  >  1:3936866‑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 < 17 ≤ ATCG/ATCG < 32 ≤ 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.