Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 4,620,828 43 bp→AATGGGTGCGGCTTGATGCCCTCACCCCGGCCCTCTCCCACAGGGAGAGGGTGAAAACATGTTGAT intergenic (+2/‑13) deoB → / → deoD phosphopentomutase/purine nucleoside phosphorylase

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4620828 4620870 43 38 [0] [0] 37 deoB/deoD phosphopentomutase/purine nucleoside phosphorylase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 46208270 (0.000)27 (0.810)
+AATGGGTGCGGCTTGATGCCCTCACCCCGGCCCTCTCCCACAGGGAGAGGGTGAAAACATGTTGAT
24/346 0.4 100% intergenic (+1/‑56) deoB/deoD phosphopentomutase/purine nucleoside phosphorylase
?NC_000913 4620871 = 0 (0.000)intergenic (+45/‑12) deoB/deoD phosphopentomutase/purine nucleoside phosphorylase

GGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACG  >  NC_000913/4620592‑4621072
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                     
ggAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttg                                                                                                                                                                                                                                        <  2:258980‑M1/251‑16 (MQ=255)
       aTGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgcccc                                                                                                                                                                                                                                 >  2:469188‑M1/1‑229 (MQ=255)
         gTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCGCGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgctgcttgatgccctc                                                                                                                                                                                                                                <  1:36075‑M1/250‑24 (MQ=255)
                         gacgacATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGACAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctccc                                                                                                                                                                                                                <  1:155883‑M1/250‑40 (MQ=255)
                             acagcctgatCCGCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagg                                                                                                                                                                                                           <  2:172214‑M1/247‑45 (MQ=255)
                                      tCCTCACTGCTGACCCCGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATGCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGTATATTGGTCAGACTCTGGCAAAATATTTTGGTACTGCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacatggagagggt                                                                                                                                                                                                   <  2:454689‑M1/250‑53 (MQ=255)
                                              gcgtgccacggtTGGGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGTGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTTGTCAGACTCGGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagatggtgaaaacatg                                                                                                                                                                                          <  2:83605‑M1/246‑62 (MQ=255)
                                                 ggcctcggtggtgatccgacctggaCCGGTACGGACCACACTCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGTTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGGTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacaggtagagggtgaaaacatgtt                                                                                                                                                                                        <  2:414099‑M1/238‑64 (MQ=255)
                                                                         aCCGGTACTGACCACACGCGGGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccc                                                                                                                                                               <  2:328637‑M1/251‑89 (MQ=255)
                                                                          ccGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttggtgccctcaccccggccctctcccaccgggagagggtgaaaacatgtttataaggataaaaaaatgggtacccc                                                                                                                                                              >  2:392486‑M1/1‑162 (MQ=255)
                                                                                  gACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacat                                                                                                                                                        <  1:79883‑M1/249‑96 (MQ=255)
                                                                                             tGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggagaaaacatgttggttagggtaaaacaaaggctacccccaacattaaagcagaaatt                                                                                                                                           >  1:31792‑M1/1‑143 (MQ=255)
                                                                                                        cGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCGTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgc                                                                                                                                <  2:436197‑M1/251‑120 (MQ=255)
                                                                                                            aCTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaac                                                                                                                                                                                             <  2:22993‑M1/186‑59 (MQ=255)
                                                                                                               ggtttacggcccgaaaGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgta                                                                                                                         <  1:222147‑M1/244‑127 (MQ=255)
                                                                                                                              aGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggc                                                                                                          <  1:418095‑M1/251‑142 (MQ=255)
                                                                                                                              aGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggc                                                                                                          <  1:178070‑M1/251‑142 (MQ=255)
                                                                                                                                           tCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatggcctcaccccggccctctcccacagggagagggtgaaaacatgttggtaagggtaaaaaaatggctaccccacacattaatgcagaaatgggcgatttcgctgaagtagatttgatgccaggcgacccgctgcgtg                                                                                             >  1:277862‑M1/1‑97 (MQ=255)
                                                                                                                                           tCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctccccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacaccttaatgcaaaaatgggcgatttcgctgacgg                                                                                                                          >  2:108491‑M1/1‑97 (MQ=255)
                                                                                                                                           tCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgt                                                                                                                          <  1:108491‑M1/222‑126 (MQ=255)
                                                                                                                                                 ggTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggggaaaacatgttgataaggataaaacaatggctaccccacacattaaagcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaggg                                                                                       >  2:11688‑M1/1‑91 (MQ=255)
                                                                                                                                                    caccgtgAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtat                                                                                     <  1:469188‑M1/247‑163 (MQ=255)
                                                                                                                                                     accgtGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacc                                                                                                      <  2:128286‑M1/230‑146 (MQ=255)
                                                                                                                                                     accgtGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcacccccgccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagtttttatgcccggcgacc                                                                                                      >  1:128286‑M1/3‑87 (MQ=255)
                                                                                                                                                          gAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTTCTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgattgcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtattttgct                                                                               <  1:93739‑M1/250‑169 (MQ=255)
                                                                                                                                                                              tcaGACTCTGGCAAAATATTTTGGTACTTCTGATATTGAATATGGCAAAGCCATGTTCTGATaatgggtgcggc                                                                                                                                                                                                                                           <  2:281492‑M1/73‑13 (MQ=255)
                                                                                                                                                                              tcaGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggc                                                                                                                                                                                                                                           >  1:281492‑M1/2‑62 (MQ=255)
                                                                                                                                                                                  aCTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcacccccgccctctcccccaaggagagggggaaaacatgtttatcaggataaaacaatggccaccccacccattaatgccaaaattggcggattcgctgacggagttttgattccaggcgacccgctgcgtgcggagtatattgctgaaactttcctttgagatgcccgtg                                                      >  2:402718‑M1/1‑58 (MQ=255)
                                                                                                                                                                                   ctctGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggggagggggaaaacctgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgattt                                                                                                                                   >  2:45157‑M1/1‑57 (MQ=255)
                                                                                                                                                                                   ctctGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgattt                                                                                                                                   <  1:45157‑M1/173‑117 (MQ=255)
                                                                                                                                                                                       ggCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaa                                                    <  1:459572‑M1/248‑196 (MQ=255)
                                                                                                                                                                                         cAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCAGGTTCTGAGaatgggtgtggcttgatgcccacaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtttgaagtatattgctgaatctttccttgaagatgcccgtgaagtgaa                                               <  1:961‑M1/251‑201 (MQ=255)
                                                                                                                                                                                                          tCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggtat                                <  1:392486‑M1/249‑216 (MQ=255)
                                                                                                                                                                                                                atGGAATATGGCAAAGCCATGGTCGGATaatgggtgcggcttgatgccctcaccccggccctcgcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgcggacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggtatgctgggc                         <  2:186092‑M1/250‑223 (MQ=255)
                                                                                                                                                                                                                              aGCCATGTTTGGATaatgggtgctgctcgttgccctcgccctggccctcgcccacagggagagggtgaaaacatgttgataatgataaaacaatggctaccccacgcattaatgcagaaatgggcgatttcgcttacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaatatgtccgtgaagtgaacaacgttcgcggtatgctgggcttcaccggtacttac          <  2:323095‑M1/251‑238 (MQ=255)
                                                                                                                                                                                                                                  aTGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagagaacaacgttcgcggtatgctgggcttcaccggtacttacaaag      <  2:80542‑M1/251‑242 (MQ=255)
                                                                                                                                                                                                                                    gTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggtatgctgggcttcaccggtacttacaaagg     <  1:46224‑M1/250‑243 (MQ=255)
                                                                                                                                                                                                                                      tCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggggaaaacatgttggtaaggataaaacaatggctaccccacacattaattcagaaattggcgattttgctgacgtagtttttatgccaagcgacccgctgcgtgcgaagtatatttctgaaactttccttgaagatgcccgtgaagtgacaaacgttcgcggtgtgctggggttcaccggtacttccaaaggccc  >  2:328793‑M1/1‑6 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                     
GGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACG  >  NC_000913/4620592‑4621072

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: