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 = 46098096 (0.140)57 (1.400) 34/156 0.1 93.8% coding (390/1224 nt) deoB phosphopentomutase
?NZ_CP009273 4609882 = 2 (0.050)coding (463/1224 nt) deoB phosphopentomutase

GGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGA  >  NZ_CP009273/4609720‑4609892
                                                                                         |                                                                                   
ggTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGAT                                                                                     >  1:1008234‑M1/1‑90 (MQ=255)
ggTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGAT                                                                                     >  1:4993‑M1/1‑90 (MQ=255)
ggTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGAT                                                                                     >  1:301304‑M1/1‑90 (MQ=255)
     cTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaact                                                                                >  1:1045918‑M1/1‑85 (MQ=255)
     cTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaact                                                                                >  1:1178679‑M1/1‑85 (MQ=255)
     cTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaact                                                                                <  2:1136253‑M1/90‑6 (MQ=255)
     cTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaact                                                                                <  1:96345‑M1/90‑6 (MQ=255)
     cTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaact                                                                                <  1:884731‑M1/90‑6 (MQ=255)
     cTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaact                                                                                <  1:442773‑M1/90‑6 (MQ=255)
       gggAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgg                                                                              >  2:557818‑M1/1‑83 (MQ=255)
           aaTTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcga                                                                          <  1:129487‑M1/90‑12 (MQ=255)
            aTTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGaa                                                                         <  2:592112/90‑1 (MQ=255)
             ttGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaag                                                                        <  1:444233‑M1/90‑14 (MQ=255)
               gCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagag                                                                      >  2:332902‑M1/1‑75 (MQ=255)
                ccGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagc                                                                     <  1:666219‑M1/90‑17 (MQ=255)
                ccGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagc                                                                     <  2:96774‑M1/90‑17 (MQ=255)
                ccGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagc                                                                     <  1:1034581‑M1/90‑17 (MQ=255)
                ccGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagc                                                                     <  2:991336‑M1/90‑17 (MQ=255)
                 cGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagca                                                                    >  1:1225430‑M1/1‑73 (MQ=255)
                             tCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccgg                                                        >  1:177‑M1/1‑61 (MQ=255)
                               tGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggca                                                      >  1:702956‑M1/1‑59 (MQ=255)
                               tGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggca                                                      >  1:841337‑M1/1‑59 (MQ=255)
                               tGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaact                                                                                >  2:1255595‑M1/1‑59 (MQ=255)
                               tGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaact                                                                                <  1:1255595‑M1/64‑6 (MQ=255)
                               tGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGtt                                                      <  1:264619/90‑1 (MQ=255)
                                   tGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagcc                                                  <  2:1095918‑M1/90‑36 (MQ=255)
                                    gAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccg                                                 >  2:111477‑M1/1‑54 (MQ=255)
                                    gAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccg                                                 >  1:752032‑M1/1‑54 (MQ=255)
                                     aGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccga                                                <  2:923647‑M1/90‑38 (MQ=255)
                                      gTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccgat                                               <  1:638897‑M1/90‑39 (MQ=255)
                                      gTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccgat                                               <  1:301934‑M1/90‑39 (MQ=255)
                                        ggggATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccgattt                                             >  1:141313‑M1/1‑50 (MQ=255)
                                             tatTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccgattttctat                                        <  2:207357‑M1/90‑46 (MQ=255)
                                               tttCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccgattttctatac                                      <  2:292284‑M1/90‑48 (MQ=255)
                                               tttCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccgattttctatac                                      <  2:113127‑M1/90‑48 (MQ=255)
                                                 tctcCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccgattttctatacct                                    <  2:551662‑M1/90‑50 (MQ=255)
                                                  ctcCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATc                                                                                    >  2:732970‑M1/1‑40 (MQ=255)
                                                  ctcCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATc                                                                                    <  1:732970‑M1/41‑2 (MQ=255)
                                                   tcCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccgattttctatacctcc                                  >  1:523328‑M1/1‑39 (MQ=255)
                                                      gATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccgattttctatacctccgct                               >  2:1019321‑M1/1‑36 (MQ=255)
                                                      gATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccgattttctatacctccgct                               >  2:1112778‑M1/1‑36 (MQ=255)
                                                      gATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccgattttctatacctccgct                               >  2:535283‑M1/1‑36 (MQ=255)
                                                      gATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTcc                               <  2:1131366/90‑1 (MQ=255)
                                                       aTCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCg                              >  1:393105/1‑90 (MQ=255)
                                                        tCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccgattttctatacctccgctga                             <  1:657895‑M1/90‑57 (MQ=255)
                                                        tCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccgattttctatacctccgctga                             <  2:177‑M1/90‑57 (MQ=255)
                                                         cACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccgattttctatacctccgctgac                            <  2:524384‑M1/90‑58 (MQ=255)
                                                         cACGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccgattttctatacctccgctgac                            >  1:399890‑M1/1‑33 (MQ=255)
                                                           cGAAAACAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccgattttctatacctccgctgactc                          <  1:758336‑M1/90‑60 (MQ=255)
                                                           cGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTAc                          <  2:393105/90‑1 (MQ=255)
                                                             aaacCAGCTTCCCGCAAGAGCTGCTGGACcaactgggcgaagagcacatgaaaaccggcaagccgattttctatacctccgctgactccg                        <  1:746141‑M1/90‑62 (MQ=255)
                                                             aaaaCAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccgattttctatacctccgctgactccg                        <  2:412319‑M1/90‑62 (MQ=255)
                                                             aaaaCAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccgattttctatacctccgctgactccg                        <  2:479597‑M1/90‑62 (MQ=255)
                                                               aaCAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccgattttctatacctccgctgactccgtg                      <  2:219487‑M1/90‑64 (MQ=255)
                                                               aaCAGCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccgattttctatacctccgctgactccgtg                      <  2:515967‑M1/90‑64 (MQ=255)
                                                                   gCTTCCCGCAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccgattttctatacctccgctgactccgtgttcc                  >  1:1024709‑M1/1‑23 (MQ=255)
                                                                           cAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccgat                                               <  2:254359‑M1/53‑39 (MQ=255)
                                                                           cAAGAGCTGCTGGATcaactgggcgaagagcacatgaaaaccggcaagccgat                                               >  1:254359‑M1/1‑15 (MQ=255)
                                                                             agagCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACt        <  2:1175472/90‑1 (MQ=255)
                                                                                gctgctGGATcaactgggcgaagagcacatgaaaaccggcaagccgattttctatacctccgctgactccgtgttccagattgcctgcca     <  1:920827‑M1/90‑81 (MQ=255)
                                                                                gctgctGGATcaactgggcgaagagcacatgaaaaccggcaagccgat                                               <  2:320016‑M1/48‑39 (MQ=255)
                                                                                gctgctGGATcaactgggcgaagagcacatgaaaaccggcaagccgat                                               >  2:795846‑M1/1‑10 (MQ=255)
                                                                                gctgctGGATcaactgggcgaagagcacatgaaaaccggcaagccgat                                               <  1:795846‑M1/48‑39 (MQ=255)
                                                                                gctgctGGATcaactgggcgaagagcacatgaaaaccggcaagccgat                                               >  1:320016‑M1/1‑10 (MQ=255)
                                                                                  tgctgGATcaactgggcgaagagcacatgaaaaccggcaagccgattttctatacctccgctgactccgtgttccagattgcctgccatg   >  2:534959‑M1/1‑8 (MQ=255)
                                                                                   gctgGATcaactgggcgaagagcacatgaaaaccggcaagccgattttctatacctccgctgactccgtgttccagattgcctgccatga  >  1:941948‑M1/1‑7 (MQ=255)
                                                                                         |                                                                                   
GGTCACTGGGAAATTGCCGGTGTCCCGGTTCTGTTTGAGTGGGGATATTTCTCCGATCACGAAAACAGCTTCCCGCAAGAGCTGCTGGATAAACTGGTCGAACGCGCTAATCTGCCGGGTTACCTCGGTAACTGCCACTCTTCCGGTACGGTCATTCTGGATCAACTGGGCGA  >  NZ_CP009273/4609720‑4609892

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