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 16 [0] [0] 17 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)13 (0.410) 13/84 0.8 100% coding (268/1383 nt) radA DNA recombination protein
?NC_000913 = 4626175 0 (0.000)coding (264/1383 nt) radA DNA recombination protein

TCCCCGCACCAGGGTTACCGCCAATCAGAATGGCACTTCCTGGCACCACGCCGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4626233‑4626179
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgccgccTAGTACGCGGTCGAACTCTTTAAATCCGGTGGAAAAACGCGGCAGCTCT  <  NC_000913/4626175‑4626127
                                                                                                        
TCCCCGCACCAGGGTTACCGCCAATCAGAATGGCACTTCCTGGCACCAC                                                         <  1:687606/49‑1
  CCCGCACCAGGGTTACCGCCAATCAGAATGGCACTTCCTGGCACCACGCC                                                      <  1:507554/50‑1
      CACCAGGGTTACCGCCAATCAGAATGGCACTTCCTGGCACCACGCCGCCT                                                  <  1:951181/50‑1
      CACCAGGGTTACCGCCAATCAGAATGGCACTTCCTGGCACCACGCCGCC                                                   >  1:2738279/1‑49
       ACCAGGGTTACCGCCAATCAGAATGGCACTTCCTGGCACCACGCCGCC                                                   >  1:1070336/1‑48
        CCAGGGTTACCGCCAATCAGAATGGCACTTCCTGGCACCACGCCGCCTAG                                                <  1:2529787/50‑1
              TTACCGCCAATCAGAATGGCACTTCCTGGCACCACGCCGCCTAGTACGCG                                          <  1:694767/50‑1
                  CGCCAATCAGAATGGCACTTCCTGGCACCACGCCGCCTAGTACGCGGTCG                                      <  1:1721795/50‑1
                   GCCAATCAGAATGGCACTTCCTGGCACCACGCCGCCTAGTACGCGGTCGA                                     <  1:1541803/50‑1
                        TCAGAATGGCACTTCCTGGCACCACGCCGCCTAGTACGCGGTCGAACTCT                                >  1:803973/1‑50
                          AGAATGGCACTTCCTGGCACCACGCCGCCTAGTACGCGGTCGAACTCTTT                              <  1:1763232/50‑1
                             ATGGCACTTCCTGGCACCACGCCGCCTAGTACGCGGTCGAACTCTTTAAA                           >  1:124438/1‑50
                             ATGGCACTTCCTGGCACCACGCCGCCTAGTACGCGGTCGAACTCTTTAAA                           <  1:470454/50‑1
                              TGGCACTTCCTGGCACCACGCCGCCTAGTACGCGGTCGAACTCTTTAAAT                          <  1:808404/50‑1
                                   CTTCCTGGCACCACGCCGCCTAGTACGCGGTCGAACTCTTTAAATCCGGT                     >  1:207479/1‑50
                                      CCTGGCACCACGCCGCCTAGTACGCGGTCGAACTCTTTAAATCCGGTG                    <  1:2434432/48‑1
                                            ACCACGCCGCCTAGTACGCGGTCGAACTCTTTAAATCCGGTGGAAAAACG            >  1:313511/1‑50
                                                  CCGCCTAGTACGCGGTCGAACTCTTTAAATCCGGTGGAAAAACGCGGCAG      >  1:2531849/1‑50
                                                      CTAGTACGCGGTCGAACTCTTTAAATCCGGTGGAAAAACGCGGCAGCTCT  >  1:997706/1‑50
                                                                                                        
TCCCCGCACCAGGGTTACCGCCAATCAGAATGGCACTTCCTGGCACCACGCCGCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/4626233‑4626179
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑cgccgccTAGTACGCGGTCGAACTCTTTAAATCCGGTGGAAAAACGCGGCAGCTCT  <  NC_000913/4626175‑4626127

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