Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 3,090,347 Δ708 bp coding (1‑708/708 nt) endA → DNA‑specific endonuclease I

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 3090347 3091054 708 26 [0] [0] 25 endA DNA‑specific endonuclease I

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 30903460 (0.000)24 (0.790) 22/98 0.1 100% intergenic (+94/‑1) yggI/endA conserved protein YggI/DNA‑specific endonuclease I
?NC_000913 3091055 = 0 (0.000)intergenic (+1/‑79) endA/rsmE DNA‑specific endonuclease I/16S rRNA m(3)U1498 methyltransferase

AATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACC  >  NC_000913/3091007‑3091104
                                                |                                                 
caaaacagctttcgctacgttgctggctcgttttaacacggagaagtgcc                                                  >  1:1385854‑M2/49‑50 (MQ=255)
  aaacagctttcgctacgttgctggctcgttttaacacggagaagtgCCTa                                                <  1:2650877‑M2/4‑1 (MQ=255)
   aacagctttcgctacgttgctggctcgttttaacacggagaagtgCCTa                                                >  1:1519790‑M2/46‑49 (MQ=255)
     cagctttcgctacgttgctggctcgttttaacacggagaagtgCCTacac                                             >  1:1937881‑M2/44‑50 (MQ=255)
      agctttcgctacgttgctggctcgttttaacacggagaagtgCCTACACt                                            <  1:220887‑M2/8‑1 (MQ=255)
        ctttcgctacgttgctggctcgttttaacacggagaagtgCCTACACTAg                                          <  1:58697‑M2/10‑1 (MQ=255)
            cgctacgttgctggctcgttttaacacggagaagtgCCTACACTAGCgg                                       >  1:877734‑M2/37‑49 (MQ=255)
              ctacgttgctggctcgttttaacacggagaagtgCCTACACTAGCGGGAt                                    <  1:198590‑M2/16‑1 (MQ=255)
                 cgttgctggctcgttttaacacggagaagtgCCTACACTAGCGGGATTCt                                 >  1:1416272‑M2/32‑50 (MQ=255)
                  gttgctggctcgttttaacacggagaagtgCCTACACTAGCGGGATTCtt                                >  1:1426068‑M2/31‑50 (MQ=255)
                   ttgctggctcgttttaacacggagaagtgCCTACACTAGCGGGATTCttt                               <  1:2061011‑M2/21‑1 (MQ=255)
                   ttgctggctcgttttaacacggagaagtgCCTACACTAGCGGGATTCttt                               >  1:2521843‑M2/30‑50 (MQ=255)
                     gctggctcgttttaacacggagaagtgCCTACACTAGCGGGATTCtttt                              <  1:907123‑M2/22‑1 (MQ=255)
                         gctcgttttaacacggagaagtgCCTACACTAGCGGGATTCTTTTTGtt                          >  1:2099428‑M2/24‑49 (MQ=255)
                         gctcgttttaacacggagaagtgCCTACACTAGCGGGATTCTTTTTGtt                          <  1:1691872‑M2/26‑1 (MQ=255)
                         gctcgttttaacacggagaagtgCCTACACTAGCGGGATTCTTTTTGt                           <  1:1859657‑M2/25‑1 (MQ=255)
                             gttttaacacggagaagtgCCTACACTAGCGGGATTCTTTTTGTTAAccc                     >  1:517117‑M2/20‑50 (MQ=255)
                              ttttaacacggagaagtgCCTACACTAGCGGGATTCTTTTTGTTAaccct                    >  1:760155‑M2/19‑49 (MQ=255)
                              ttttaacacggagaagtgCCTACACTAGCGGGATTCTTTTTGTTAAcccc                    <  1:239187‑M2/32‑1 (MQ=255)
                               tttaacacggagaagtgCCTACACTAGCGGGATTCTTTTTGTTAAcccc                    >  1:2149456‑M2/18‑49 (MQ=255)
                                   acacggagaagtgCCTACACTAGCGGGATTCTTTTTGTTAACCCCTAccc               <  1:2187623‑M2/37‑1 (MQ=255)
                                         agaagtgCCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCAcgcg         >  1:871117‑M2/8‑50 (MQ=255)
                                           aagtgCCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTa       <  1:321280‑M2/45‑1 (MQ=255)
                                               gCCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAAc   >  1:2468744‑M2/2‑50 (MQ=255)
                                                ccTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAAcc  >  1:1375104/1‑50 (MQ=255)
                                                |                                                 
AATCATAACCCGTATGTGCAACGCGCTTGCCAGGCGCGAAAGAGCTAACCTACACTAGCGGGATTCTTTTTGTTAACCCCTACCCCACGCGTACAACC  >  NC_000913/3091007‑3091104

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