Predicted mutation
evidence seq id position mutation annotation gene description
JC NZ_CP009273 4,609,810 Δ72 bp coding (391‑462/1224 nt) deoB → phosphopentomutase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 46098093 (0.190)23 (1.540) 17/156 0.1 90.6% coding (390/1224 nt) deoB phosphopentomutase
?NZ_CP009273 4609882 = 2 (0.130)coding (463/1224 nt) deoB phosphopentomutase

GGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGC  >  NZ_CP009273/4609720‑4609897
                                                                                         |                                                                                        
ggTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGAT                                                                                          >  2:313030‑M1/1‑90 (MQ=255)
ggTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGAT                                                                                          >  2:131358‑M1/1‑90 (MQ=255)
     cTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaact                                                                                     <  2:162188‑M1/90‑6 (MQ=255)
     cTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaact                                                                                     <  1:71450‑M1/90‑6 (MQ=255)
      tGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTg                                                                                    >  1:327636/1‑90 (MQ=255)
          aaaTTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcg                                                                                >  2:154230‑M1/1‑80 (MQ=255)
             ttGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaag                                                                             <  2:328704‑M1/90‑14 (MQ=255)
                ccGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagc                                                                          <  2:336397‑M1/90‑17 (MQ=255)
                               tGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggca                                                           >  1:140462‑M1/1‑59 (MQ=255)
                               tGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggca                                                           >  2:326513‑M1/1‑59 (MQ=255)
                                gTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaa                                                          <  1:197165‑M1/90‑33 (MQ=255)
                                   tGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagcc                                                       <  1:366359‑M1/90‑36 (MQ=255)
                                   tGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagcc                                                       <  2:161506‑M1/90‑36 (MQ=255)
                                   tGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagcc                                                       <  2:302264‑M1/90‑36 (MQ=255)
                                     aGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccga                                                     >  2:370490‑M1/1‑53 (MQ=255)
                                      gTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccgat                                                    <  2:140462‑M1/90‑39 (MQ=255)
                                        ggggATATTTCTCCGATCACGAAAACAGCTTCCCGCGAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccgattt                                                  <  2:304811‑M1/90‑41 (MQ=255)
                                          ggATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccgattttc                                                >  2:66082‑M1/1‑48 (MQ=255)
                                                     cGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccgattttctatacctccgc                                     >  2:119825‑M1/1‑37 (MQ=255)
                                                           cGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccgattttctatacctccgctgactc                               <  1:394235‑M1/90‑60 (MQ=255)
                                                             aaaaCAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccgattttctatacctccgctgactccg                             <  1:326513‑M1/90‑62 (MQ=255)
                                                                   gCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccgattttctatacctccgctgactccgtgttcc                       >  2:17476‑M1/1‑23 (MQ=255)
                                                                   gCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccgattttctatacctccgctgactccgtgttcc                       >  1:329113‑M1/1‑23 (MQ=255)
                                                                    cTTCCCGCAAGAGCTGCTGGATAACCTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCt                      <  1:251958/90‑1 (MQ=255)
                                                                     ttCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccgat                                                    >  2:274346‑M1/1‑21 (MQ=255)
                                                                     ttCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccgat                                                    <  1:274346‑M1/59‑39 (MQ=255)
                                                                             agagCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACt             >  2:403766/1‑90 (MQ=255)
                                                                                  tgctgGATcaactgggcgaagagcacatgaaaaccggcaagccgattttctatacctccgctgactccgtgttccagattgcctgccatg        >  2:245773‑M1/1‑8 (MQ=255)
                                                                                        aTcaactgggcgaagagcacatgaaaaccggcaagccgattttctatacctccgctgactccgtgttccagattgcctgccatgaagaaa  >  1:240985‑M1/1‑2 (MQ=255)
                                                                                         |                                                                                        
GGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGC  >  NZ_CP009273/4609720‑4609897

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