Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 4,626,183 (GCG)3→2 coding (272‑274/1383 nt) radA → DNA recombination protein

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4626176 4626178 3 15 [0] [0] 16 radA DNA recombination protein

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 4626179 =0 (0.000)12 (0.380) 12/82 0.7 100% coding (268/1383 nt) radA DNA recombination protein
?NC_000913 = 4626175 0 (0.000)coding (264/1383 nt) radA DNA recombination protein

CCCGCACCAGGGTTACCGCCAATCAGAATGGCACTTCCTGGCACCACGCCGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4626231‑4626179
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgccgccTAGTACGCGGTCGAACTCTTTAAATCCGGTGGAAAAACGCGGCAGCTC  <  NC_000913/4626175‑4626128
                                                                                                     
CCCGCACCAGGGTTACCGCCAATCAGAATGGCACTTCCTGGCACCACGC                                                      <  1:1004187/49‑1
 CCGCACCAGGGTTACCGCCAATCAGAATGGCACTTCCTGGCACCACGCC                                                     <  1:2677807/49‑1
     ACCAGGGTTACCGCCAATCAGAATGGCACTTCCTGGCACCACGCCGCC                                                  >  1:1766317/1‑48
           GTTACCGCCAATCAGAATGGCACTTCCTGGCACCACGCCGCCTAGTACG                                           <  1:910303/49‑1
            TTACCGCCAATCAGAATGGCACTTCCTGGCACCACGCCGCCTAGTACGC                                          <  1:1559052/49‑1
            TTACCGCCAATCAGAATGGCACTTCCTGGCACCACGCCGCCTAGTACGC                                          >  1:2458392/1‑49
             TACCGCCAATCAGAATGGCACTTCCTGGCACCACGCCGCCTAGTACGCG                                         >  1:1035980/1‑49
                     ATCAGAATGGCACTTCCTGGCACCACGCCGCCTAGTACGCGGTCGAACT                                 >  1:338042/1‑49
                         GAATGGCACTTCCTGGCACCACGCCGCCTAGTACGCGGTCGAACTCTTT                             <  1:203770/49‑1
                              GCACTTCCTGGCACCACGCCGCCTAGTACGCGGTCGAACTCTTTAAATC                        >  1:2231991/1‑49
                                ACTTCCTGGCACCACGCCGCCTAGTACGCGGTCGAACTCTTTAAATCCG                      >  1:2329626/1‑49
                                ACTTCCTGGCACCACGCCGCCTAGTACGCGGTCGAACTCTTTAAATCCG                      <  1:1449795/49‑1
                                    CCTGGCACCACGCCGCCTAGTACGCGGTCGAACTCTTTAAATCCGGT                    <  1:1780533/47‑1
                                      TGGCACCACGCCGCCTAGTACGCGGTCGAACTCTTTAAATCCGGTGGA                 >  1:2724526/1‑48
                                          ACCACGCCGCCTAGTACGCGGTCGAACTCTTTAAATCCGGTGGAAAAAC            <  1:1838247/49‑1
                                                CCGCCTAGTACGCGGTCGAACTCTTTAAATCCGGTGGAAAAACGCGGCA      <  1:2715088/49‑1
                                                  GCCTAGTACGCGGTCGAACTCTTTAAATCCGGTGGAAAAACGCGGCAGC    >  1:2353126/1‑49
                                                    CTAGTACGCGGTCGAACTCTTTAAATCCGGTGGAAAAACGCGGCAGCTC  <  1:1279227/49‑1
                                                                                                     
CCCGCACCAGGGTTACCGCCAATCAGAATGGCACTTCCTGGCACCACGCCGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4626231‑4626179
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgccgccTAGTACGCGGTCGAACTCTTTAAATCCGGTGGAAAAACGCGGCAGCTC  <  NC_000913/4626175‑4626128

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