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 22 [0] [0] 23 [dgcJ] [dgcJ]

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 18720290 (0.000)21 (0.710) 20/94 0.2 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

CGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/1871984‑1872029
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAAC  >  NC_000913/1873467‑1873513
                                                                                             
CGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGATG                                              >  1:1676297/1‑49
 GCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGATGC                                             <  1:2312553/49‑1
     TCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGATGCGTAT                                         <  1:528825/49‑1
     TCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGATGCGTA                                          >  1:837425/1‑48
       GAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGATGCGTATA                                        >  1:1289430/1‑48
        AACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGATGCGTATAAA                                      <  1:1347699/49‑1
         ACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGATGCGTATAAAG                                     >  1:1315392/1‑49
          CTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGATGCGTATAAAGC                                    >  1:1599602/1‑49
          CTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAGATGCGTATAAAGC                                    >  1:570285/1‑49
             TATCTGGTTTTCTCGTTTCACTAACCGAAGGAGATGCGTATAAAGCT                                   <  1:15409/47‑1
               TCTGGTTTTCTCGTTTCACTAACCGAAGGAGATGCGTATAAAGCTTCC                                <  1:63606/48‑1
                   GTTTTCTCGTTTCACTAACCGAAGGAGATGCGTATAAAGCTTCCGAT                             <  1:1843848/47‑1
                   GTTTTCTCGTTTCACTAACCGAAGGAGATGCGTATAAAGCTTCCG                               <  1:2232181/45‑1
                     TTTCTCGTTTCACTAACCGAAGGAGATGCGTATAAAGCTTCCGATGAGC                         >  1:1753137/1‑49
                       TCTCGTTTCACTAACCGAAGGAGATGCGTATAAAGCTTCCGATGAGCGT                       >  1:436661/1‑49
                           GTTTCACTAACCGAAGGAGATGCGTATAAAGCTTCCGATGAGCGTTTA                    >  1:975834/1‑48
                           GTTTCACTAACCGAAGGAGATGCGTATAAAGCTTCCGATGAGCGTTT                     <  1:1058804/47‑1
                                    ACCGAAGGAGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATA          <  1:1778781/49‑1
                                       GAAGGAGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAG        >  1:1375411/1‑48
                                          GGAGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAG     >  1:2342494/1‑48
                                           GAGATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAA   >  1:1167306/1‑49
                                             GATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAAC  >  1:942844/1‑48
                                                                                             
CGCAATCGAACTTTATCTGGTTTTCTCGTTTCACTAACCGAAGGAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/1871984‑1872029
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gATGCGTATAAAGCTTCCGATGAGCGTTTATATGTCAATAAGCAGAAC  >  NC_000913/1873467‑1873513

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