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 25 [0] [0] 26 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)24 (0.940) 22/84 0.0 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:1935129/49‑1
 CCCCGCACCAGGGTTACCGCCAATCAGAATGGCACTTCCTGGCACCACGC                                                       >  1:795161/1‑50
   CCGCACCAGGGTTACCGCCAATCAGAATGGCACTTCCTGGCACCACGCCG                                                     >  1:2083582/1‑50
    CGCACCAGGGTTACCGCCAATCAGAATGGCACTTCCTGGCACCACGCCGC                                                    >  1:423686/1‑50
     GCACCAGGGTTACCGCCAATCAGAATGGCACTTCCTGGCACCACGCCGCC                                                   <  1:1198664/50‑1
       ACCAGGGTTACCGCCAATCAGAATGGCACTTCCTGGCACCACGCCGCCTA                                                 <  1:2009688/50‑1
          AGGGTTACCGCCAATCAGAATGGCACTTCCTGGCACCACGCCGCCTAGTA                                              <  1:1243847/50‑1
           GGGTTACCGCCAATCAGAATGGCACTTCCTGGCACCACGCCGCCTAGTAC                                             >  1:638159/1‑50
            GGTTACCGCCAATCAGAATGGCACTTCCTGGCACCACGCCGCCTAGTACG                                            >  1:527867/1‑50
              TTACCGCCAATCAGAATGGCACTTCCTGGCACCACGCCGCCTAGTACGCG                                          <  1:472817/50‑1
                 CCGCCAATCAGAATGGCACTTCCTGGCACCACGCCGCCTAGTACGCGGT                                        <  1:2145087/49‑1
                   GCCAATCAGAATGGCACTTCCTGGCACCACGCCGCCTAGTACGCGGTCG                                      >  1:1872862/1‑49
                    CCAATCAGAATGGCACTTCCTGGCACCACGCCGCCTAGTACGCGGTCGAA                                    >  1:917382/1‑50
                      AATCAGAATGGCACTTCCTGGCACCACGCCGCCTAGTACGCGGTCGAACT                                  <  1:261078/50‑1
                      AATCAGAATGGCACTTCCTGGCACCACGCCGCCTAGTACGCGGTCGAACT                                  >  1:783523/1‑50
                       ATCAGAATGGCACTTCCTGGCACCACGCCGCCTAGTACGCGGTCGAACTC                                 <  1:365925/50‑1
                        TCAGAATGGCACTTCCTGGCACCACGCCGCCTAGTACGCGGTCGAACTCT                                >  1:1435858/1‑50
                          AGAATGGCACTTCCTGGCACCACGCCGCCTAGTACGCGGTCGAACTCTTT                              <  1:1366021/50‑1
                          AGAATGGCACTTCCTGGCACCACGCCGCCTAGTACGCGGTCGAACTCTT                               >  1:752072/1‑49
                           GAATGGCACTTCCTGGCACCACGCCGCCTAGTACGCGGTCGAACTCTTTA                             <  1:1219134/50‑1
                           GAATGGCACTTCCTGGCACCACGCCGCCTAGTACGCGGTCGAACTCTTTA                             <  1:1647223/50‑1
                            AATGGCACTTCCTGGCACCACGCCGCCTAGTACGCGGTCGAACTCTTTAA                            <  1:2000361/50‑1
                              TGGCACTTCCTGGCACCACGCCGCCTAGTACGCGGTCGAACTCTTTAAAT                          <  1:1033690/50‑1
                              TGGCACTTCCTGGCACCACGCCGCCTAGTACGCGGTCGAACTCTTTAAAT                          <  1:2152048/50‑1
                               GGCACTTCCTGGCACCACGCCGCCTAGTACGCGGTCGAACTCTTTAAATC                         >  1:317790/1‑50
                                GCACTTCCTGGCACCACGCCGCCTAGTACGCGGTCGAACTCTTTAAATCC                        <  1:1854323/50‑1
                                        TGGCACCACGCCGCCTAGTACGCGGTCGAACTCTTTAAATCCGGTGGAA                 >  1:2029829/1‑49
                                         GGCACCACGCCGCCTAGTACGCGGTCGAACTCTTTAAATCCGGTGGAAA                >  1:1826794/1‑49
                                          GCACCACGCCGCCTAGTACGCGGTCGAACTCTTTAAATCCGGTGGAAA                <  1:69291/48‑1
                                                      CTAGTACGCGGTCGAACTCTTTAAATCCGGTGGAAAAACGCGGCAGCTCT  >  1:420611/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 < 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.