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

AGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACT  >  NZ_CP009273/4609797‑4609971
                                                                                     |                                                                                         
ctgggaaattgccggtgtcccggttctgtttgagtggggatatttctccgatcacgaaaacagcttcccgcaagagctgctggatCAACt                                                                                       <  2:162188‑M2/5‑1 (MQ=255)
ctgggaaattgccggtgtcccggttctgtttgagtggggatatttctccgatcacgaaaacagcttcccgcaagagctgctggatCAACt                                                                                       <  1:71450‑M2/5‑1 (MQ=255)
agagCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACt                                                                                       >  2:403766/1‑90 (MQ=255)
     aaattgccggtgtcccggttctgtttgagtggggatatttctccgatcacgaaaacagcttcccgcaagagctgctggatCAACTGGGCg                                                                                  >  2:154230‑M2/81‑90 (MQ=255)
        ttgccggtgtcccggttctgtttgagtggggatatttctccgatcacgaaaacagcttcccgcaagagctgctggatCAACTGGGCGAag                                                                               <  2:328704‑M2/13‑1 (MQ=255)
           ccggtgtcccggttctgtttgagtggggatatttctccgatcacgaaaacagcttcccgcaagagctgctggatCAACTGGGCGAAGAGc                                                                            <  2:336397‑M2/16‑1 (MQ=255)
                          tgtttgagtggggatatttctccgatcacgaaaacagcttcccgcaagagctgctggatCAACTGGGCGAAGAGCACATGAAAACCGGCa                                                             >  1:140462‑M2/60‑90 (MQ=255)
                          tgtttgagtggggatatttctccgatcacgaaaacagcttcccgcaagagctgctggatCAACTGGGCGAAGAGCACATGAAAACCGGCa                                                             >  2:326513‑M2/60‑90 (MQ=255)
                           gtttgagtggggatatttctccgatcacgaaaacagcttcccgcaagagctgctggatCAACTGGGCGAAGAGCACATGAAAACCGGCaa                                                            <  1:197165‑M2/32‑1 (MQ=255)
                              tgagtggggatatttctccgatcacgaaaacagcttcccgcaagagctgctggatCAACTGGGCGAAGAGCACATGAAAACCGGCAAGcc                                                         <  2:161506‑M2/35‑1 (MQ=255)
                              tgagtggggatatttctccgatcacgaaaacagcttcccgcaagagctgctggatCAACTGGGCGAAGAGCACATGAAAACCGGCAAGcc                                                         <  1:366359‑M2/35‑1 (MQ=255)
                              tgagtggggatatttctccgatcacgaaaacagcttcccgcaagagctgctggatCAACTGGGCGAAGAGCACATGAAAACCGGCAAGcc                                                         <  2:302264‑M2/35‑1 (MQ=255)
                                agtggggatatttctccgatcacgaaaacagcttcccgcaagagctgctggatCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGa                                                       >  2:370490‑M2/54‑90 (MQ=255)
                                 gtggggatatttctccgatcacgaaaacagcttcccgcaagagctgctggatCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGAt                                                      <  2:140462‑M2/38‑1 (MQ=255)
                                   ggggatatttctccgatcacgaaaacagcttcccgcgagagctgctggatCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGAttt                                                    <  2:304811‑M2/40‑1 (MQ=255)
                                     ggatatttctccgatcacgaaaacagcttcccgcaagagctgctggatCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTc                                                  >  2:66082‑M2/49‑90 (MQ=255)
                                                cgatcacgaaaacagcttcccgcaagagctgctggatCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGc                                       >  2:119825‑M2/38‑90 (MQ=255)
                                                      cgaaaacagcttcccgcaagagctgctggatCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTc                                 <  1:394235‑M2/59‑1 (MQ=255)
                                                        gCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCg                               <  1:403766/90‑1 (MQ=255)
                                                        aaaacagcttcccgcaagagctgctggatCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCg                               <  1:326513‑M2/61‑1 (MQ=255)
                                                              gcttcccgcaagagctgctggatCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTcc                         >  2:17476‑M2/24‑90 (MQ=255)
                                                              gcttcccgcaagagctgctggatCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTcc                         >  1:329113‑M2/24‑90 (MQ=255)
                                                                ttcccgcaagagctgctggatCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGAt                                                      >  2:274346‑M2/22‑59 (MQ=255)
                                                                ttcccgcaagagctgctggatCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGAt                                                      <  1:274346‑M2/38‑1 (MQ=255)
                                                                             tgctggatCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATg          >  2:245773‑M2/9‑90 (MQ=255)
                                                                                   atCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGaaa    >  1:240985‑M2/3‑90 (MQ=255)
                                                                                     cAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACt  >  2:148250/1‑90 (MQ=255)
                                                                                     |                                                                                         
AGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGAAGAGCACATGAAAACCGGCAAGCCGATTTTCTATACCTCCGCTGACTCCGTGTTCCAGATTGCCTGCCATGAAGAAACT  >  NZ_CP009273/4609797‑4609971

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