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 83 [0] [0] 81 [dgcJ] [dgcJ]

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

GGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAA  >  NC_000913/1873418‑1873515
                                                 |                                                
ggccgcaatcgaactttatctggttttctcgtttcactaaccgaaggaga                                                  <  1:2402140‑M2/1‑1 (MQ=255)
 gccgcaatcgaactttatctggttttctcgtttcactaaccgaaggagAt                                                 >  1:40177‑M2/49‑50 (MQ=255)
     caatcgaactttatctggttttctcgtttcactaaccgaaggagATGCg                                              >  1:3683136‑M2/45‑49 (MQ=255)
     caatcgaactttatctggttttctcgtttcactaaccgaaggagATGCg                                              >  1:1858251‑M2/45‑49 (MQ=255)
      aatcgaactttatctggttttctcgtttcactaaccgaaggagATGCGta                                            >  1:1027273‑M2/44‑50 (MQ=255)
      aatcgaactttatctggttttctcgtttcactaaccgaaggagATGCGta                                            <  1:1969446‑M2/7‑1 (MQ=255)
       atcgaattttatctggttttctcgtttcactaaccgaaggagATGCGta                                            >  1:5347675‑M2/43‑49 (MQ=255)
       atcgaactttatctggttttctcgtttcactaaccgaaggagATGCGta                                            >  1:4948734‑M2/43‑49 (MQ=255)
        tcgaactttatctggttttctcgtttcactaaccgaaggagATGCGtata                                          <  1:4431267‑M2/9‑1 (MQ=255)
        tcgaactttatctggttttctcgtttcactaaccgaaggagATGCGtat                                           >  1:3739478‑M2/42‑49 (MQ=255)
         cgaactttatctggttttctcgtttcactaaccgaaggagATGCGtata                                          >  1:95401‑M2/41‑49 (MQ=255)
         cgaactttatctggttttctcgtttcactaaccgaaggagATGCGtata                                          >  1:1463273‑M2/41‑49 (MQ=255)
          gaactttatctggttttctcgtttcactaaccgaaggagATGCGtata                                          >  1:6630942‑M2/40‑48 (MQ=255)
          gaactttatctggttttctcgtttcactaaccgaaggagATGCGtata                                          >  1:5415523‑M2/40‑48 (MQ=255)
          gaactttatctggttttctcgtttcactaaccgaaggagATGCGtata                                          >  1:4116855‑M2/40‑48 (MQ=255)
          gaactttatctggttttctcgtttcactaaccgaaggagATGCGTATaaa                                        <  1:2621997‑M2/11‑1 (MQ=255)
          gaactttatctggttttctcgtttcactaaccgaaggagATGCGTATaaa                                        <  1:2801921‑M2/11‑1 (MQ=255)
          gaactttatctggttttctcgtttcactaaccgaaggagATGCGTATaa                                         <  1:5557649‑M2/10‑1 (MQ=255)
           aactttatctggttttctcgtttcactaaccgaaggagATGCGTATAAAg                                       >  1:5935168‑M2/39‑50 (MQ=255)
            actttatctggttttctcgtttcactaaccgaaggagATGCGTATAAAGc                                      >  1:1706460‑M2/38‑50 (MQ=255)
             ctttatctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCt                                     >  1:3614279‑M2/37‑50 (MQ=255)
              tttatctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCtt                                    >  1:3985099‑M2/36‑50 (MQ=255)
               ttatctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTc                                   >  1:4863842‑M2/35‑50 (MQ=255)
                tatctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTcc                                  >  1:1014901‑M2/34‑50 (MQ=255)
                tatctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTcc                                  >  1:2581988‑M2/34‑50 (MQ=255)
                 atctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTcc                                  <  1:4291999‑M2/17‑1 (MQ=255)
                  tctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCg                                 <  1:2471218‑M2/18‑1 (MQ=255)
                  tctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGa                                >  1:2178965‑M2/32‑50 (MQ=255)
                   ctggttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGAt                               >  1:938651‑M2/31‑50 (MQ=255)
                      gttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATg                              <  1:4570793‑M2/21‑1 (MQ=255)
                      gttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGa                             <  1:2887339‑M2/22‑1 (MQ=255)
                      gttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGa                             <  1:718476‑M2/22‑1 (MQ=255)
                      gttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGa                             <  1:4656272‑M2/22‑1 (MQ=255)
                      gttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAg                            <  1:73072‑M2/23‑1 (MQ=255)
                      gttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAg                            <  1:271135‑M2/23‑1 (MQ=255)
                       ttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGc                           <  1:6194448‑M2/24‑1 (MQ=255)
                       ttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGc                           <  1:1608077‑M2/24‑1 (MQ=255)
                       ttttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGc                           >  1:4811111‑M2/27‑50 (MQ=255)
                        tttctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCg                          >  1:2043729‑M2/26‑50 (MQ=255)
                          tctcgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGt                         <  1:136643‑M2/26‑1 (MQ=255)
                             cgtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTta                      >  1:3795901‑M2/21‑49 (MQ=255)
                              gtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGttt                       <  1:2612803‑M2/28‑1 (MQ=255)
                              gtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGttt                       <  1:3710113‑M2/28‑1 (MQ=255)
                              gtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTtata                    >  1:4075315‑M2/20‑50 (MQ=255)
                              gtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTtata                    >  1:2125748‑M2/20‑50 (MQ=255)
                              gtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTtat                     <  1:3449298‑M2/30‑1 (MQ=255)
                              gtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTtat                     <  1:1369489‑M2/30‑1 (MQ=255)
                              gtttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTta                      <  1:4400672‑M2/29‑1 (MQ=255)
                               tttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTatat                   >  1:3035447‑M2/19‑50 (MQ=255)
                               tttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTatat                   >  1:4064018‑M2/19‑50 (MQ=255)
                                ttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATg                  >  1:6740898‑M2/18‑50 (MQ=255)
                                ttcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATg                  >  1:6556834‑M2/18‑50 (MQ=255)
                                 tcactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATg                  <  1:6390497‑M2/33‑1 (MQ=255)
                                  cactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGt                 <  1:286889‑M2/34‑1 (MQ=255)
                                  cactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGt                 <  1:6521258‑M2/34‑1 (MQ=255)
                                  cactaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGt                 <  1:4909995‑M2/34‑1 (MQ=255)
                                   actaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTc                >  1:3453353‑M2/15‑49 (MQ=255)
                                   actaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCa               <  1:4098317‑M2/36‑1 (MQ=255)
                                   actaaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCa               <  1:3947494‑M2/36‑1 (MQ=255)
                                     taaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCaa              >  1:1218123‑M2/13‑49 (MQ=255)
                                      aaccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAAt             >  1:3091066‑M2/12‑49 (MQ=255)
                                       accgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAAt             <  1:4752950‑M2/38‑1 (MQ=255)
                                       accgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATa            <  1:6060904‑M2/39‑1 (MQ=255)
                                        ccgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAg          >  1:6739573‑M2/10‑50 (MQ=255)
                                         cgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAg          >  1:3391561‑M2/9‑49 (MQ=255)
                                         cgaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGc         <  1:1164471‑M2/42‑1 (MQ=255)
                                          gaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAg          >  1:6644627‑M2/8‑48 (MQ=255)
                                          gaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAg          >  1:3036669‑M2/8‑48 (MQ=255)
                                          gaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAg          >  1:4912758‑M2/8‑48 (MQ=255)
                                          gaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAg          >  1:3462864‑M2/8‑48 (MQ=255)
                                          gaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCa        <  1:4258008‑M2/43‑1 (MQ=255)
                                          gaaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCa        <  1:415582‑M2/43‑1 (MQ=255)
                                           aaggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAg          >  1:6241465‑M2/7‑47 (MQ=255)
                                            aggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAg       >  1:2984046‑M2/6‑49 (MQ=255)
                                            aggagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGa      <  1:5652472‑M2/45‑1 (MQ=255)
                                              gagATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAAc    <  1:309868‑M2/47‑1 (MQ=255)
                                               agATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAAc    >  1:2313239‑M2/3‑49 (MQ=255)
                                               agATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAAc    >  1:1193271‑M2/3‑49 (MQ=255)
                                               agATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACa   <  1:4843301‑M2/48‑1 (MQ=255)
                                                gATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACaa  <  1:2958527‑M2/49‑1 (MQ=255)
                                                gATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACa   >  1:4772979‑M2/2‑49 (MQ=255)
                                                 |                                                
GGCTTCTCTTCCGGTATTTACGCGATGAAAGAAAACGATACGTTACATGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAACAA  >  NC_000913/1873418‑1873515

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