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

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 46208271 (0.020)22 (0.500)
+AATGGGTGCGGCTTGATGCCCTCACCCCGGCCCTCTCCCACAGGGAGAGGGTGAAAACATGTTGAT
17/350 1.2 96.8% intergenic (+1/‑56) deoB/deoD phosphopentomutase/purine nucleoside phosphorylase
?NC_000913 4620871 = 1 (0.020)intergenic (+45/‑12) deoB/deoD phosphopentomutase/purine nucleoside phosphorylase

GCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGG  >  NC_000913/4620589‑4621067
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                
gCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcg                                                                                                                                                                                                                                        <  1:60062‑M1/249‑11 (MQ=255)
         gATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccc                                                                                                                                                                                                                             >  1:144726‑M1/1‑230 (MQ=255)
                    tgcGCGATGACGACATCCTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctccccccggcc                                                                                                                                                                                                                  >  1:233616‑M1/1‑219 (MQ=255)
                                    ccttatcctCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTCACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggaga                                                                                                                                                                                                  <  2:159945‑M1/247‑49 (MQ=255)
                                    ccTGATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggag                                                                                                                                                                                                   <  1:404521‑M1/250‑48 (MQ=255)
                                       gATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagg                                                                                                                                                                                                <  1:86529‑M1/250‑51 (MQ=255)
                                       gATCCTCACTGCTGACCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATGCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACGGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagag                                                                                                                                                                                                 <  2:562635‑M1/249‑50 (MQ=255)
                                        aTCCTCACTGCTGACCACTGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggt                                                                                                                                                                                              <  2:652716‑M1/251‑53 (MQ=255)
                                             cactgctgaCCACGGTTGTGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaa                                                                                                                                                                                         <  1:180261‑M1/247‑58 (MQ=255)
                                                                    cGACATGGACCGGTAATGACCACACGAGTGAAAACATTCCGGTAATTGTTTACGGCCCGAAAGTAAAACCGGGCAAACTGGGTAAACGTGAAACCTTCGCGGATATTGGTCAGACTCGGGCAAAATATTTAGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcacaccggcactctcccacagggagagggtgaaaacatgttgataaggataaaacaat                                                                                                                                                                  <  2:513447‑M1/251‑81 (MQ=255)
                                                                         tGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagag                                                                                                                                                                                                 <  2:534130‑M1/215‑50 (MQ=255)
                                                                         tGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagag                                                                                                                                                                                                 >  1:534130‑M1/1‑166 (MQ=255)
                                                                            aCCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccccagggagagggtgaaaacatgttgataaggataaaacaatggctaccc                                                                                                                                                          >  1:378312‑M1/1‑163 (MQ=255)
                                                                                                       aTTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGTCTCACTGGGTCTCCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATTGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataagacaatggctaccccacacattaatgcagaaatgggcgatt                                                                                                                               <  2:573123‑M1/251‑116 (MQ=255)
                                                                                                       aTTCCGGTACTGGTTTACGGCCCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatt                                                                                                                               <  1:666678‑M1/251‑116 (MQ=255)
                                                                                                                       acggccCGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACCTTCGCGGATATTGGACAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggggcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttt                                                                                                                <  2:408635‑M1/248‑131 (MQ=255)
                                                                                                                           cccGAAAGTAAAACCGGGCTCACTGGGTCACCGTGAAACTGTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatg                                                                                                           <  1:330344‑M1/251‑136 (MQ=255)
                                                                                                                                     aaaCCGGGCTCACTGGGGCACCGTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatacaaaaatgagcgatttcgctgaagtagtttttaagaccgacgaac                                                                                                 >  2:441669‑M1/1‑106 (MQ=255)
                                                                                                                                                       caccgtgAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctccccccagccctctcccacagggggggggtgaaaacatgttgataaggataaaacaattgataccccacacattaatgcagaaatgggcgatttccctgccgtcgtttttatgccaggcggcccgctgcggggggaatat                                                                                >  1:39396‑M1/4‑88 (MQ=255)
                                                                                                                                                       caccgtgAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtat                                                                                <  1:505988‑M1/247‑163 (MQ=255)
                                                                                                                                                       caccgtgAAACCTTCGCGGATATTGGTCAGACTCTGGCAAAATAATTTGGTACTTCTGATATTGAATATGGCAAAGCCATGTTCTGATaataggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtat                                                                                <  2:546951‑M1/247‑163 (MQ=255)
                                                                                                                                                         ccgTGAAACCTTCGCGGATATTGGTCAGACTCTGGCAGAATATTTTGGTACTTCTGATATGTAATATGGCAAAGCCATGTTCTGATaatgggttcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggcttccccacacattaatgcagaaatgggagatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatatt                                                                             <  2:511092‑M1/250‑166 (MQ=255)
                                                                                                                                                                  cTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagggggtgaaaacatgttgataaggataaaacaatggctaccccacacatta                                                                                                                                                 >  1:174499‑M1/1‑77 (MQ=255)
                                                                                                                                                                  cTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacatta                                                                                                                                                 <  2:174499‑M1/174‑98 (MQ=255)
                                                                                                                                                                   ttCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggggagggggaaaacatgttgg                                                                                                                                                                                 >  1:220296‑M1/1‑76 (MQ=255)
                                                                                                                                                                   ttCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttga                                                                                                                                                                                 <  2:220296‑M1/141‑66 (MQ=255)
                                                                                                                                                                     cgtggataTTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaacttt                                                                  <  2:291151‑M1/247‑177 (MQ=255)
                                                                                                                                                                     cgcgGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaacttt                                                                  <  1:276264‑M1/250‑177 (MQ=255)
                                                                                                                                                                              tggtcagacTCTGGTAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaaga                                                         <  1:644309‑M1/246‑186 (MQ=255)
                                                                                                                                                                                ggcagACTCTGTCAAAATATTTTGGTGCTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaacccggccctctcccacagggagagggtgaaaacatgttgataaggataaaaccatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgc                                                      <  2:378312‑M1/249‑189 (MQ=255)
                                                                                                                                                                                                      tGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAAAAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGTGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGGTGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTAc                                >  1:659433/1‑251 (MQ=255)
                                                                                                                                                                                                        gTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcgg                              <  1:601797‑M1/251‑213 (MQ=255)
                                                                                                                                                                                                                     ggAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggtatgctgggctt                  <  1:619381‑M1/250‑225 (MQ=255)
                                                                                                                                                                                                                          atGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggtatgctgggcttcaccg             <  1:535314‑M1/250‑230 (MQ=255)
                                                                                                                                                                                                                                    cATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggtatgctgggcttcaccggtacttacaaa  <  1:575972‑M1/251‑241 (MQ=255)
                                                                                                                                                                                                                                    cATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggtatgctgggcttcaccggtacttacaaa  <  1:356465‑M1/251‑241 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                
GCCGGAGCTGATGTCTCTGCTGCGCGATGACGACATCCTGATCCTCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGG  >  NC_000913/4620589‑4621067

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