Predicted mutation
evidence seq id position mutation annotation gene description
JC JC NC_000913 4,620,780 IS1 (–) +9 bp coding (1178‑1186/1224 nt) deoB → phosphopentomutase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3583428 =NA (NA)30 (1.020) 26/260 0.3 100% noncoding (1/768 nt) IS1 repeat region
?NC_000913 4620780 = 0 (0.000)coding (1178/1224 nt) deoB phosphopentomutase
* ? NC_000913 = 3584195NA (NA)31 (1.060) 27/258 0.2 100% noncoding (768/768 nt) IS1 repeat region
?NC_000913 = 4620788 0 (0.000)coding (1186/1224 nt) deoB phosphopentomutase

TGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTT  >  NC_000913/4620640‑4620933
                                                                                                                                                    |                                                                                                                                                 
gctgatgtccggcggtgcttttgccgttacgcaccaccccgtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccAATATTTTg                                                                                                                                                   <  6:39459‑M2/9‑1 (MQ=255)
    cACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGggt                                                                                                                                                <  7:64296‑M2/148‑4 (MQ=255)
     tgtccggcggtgcttttgccgttacgcaccaccccgtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccAATATTTTGGTACt                                                                                                                                              >  8:85705‑M2/136‑149 (MQ=255)
       tccggcggtgcttttgccgttaagcaccaccccgtcagtagctgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccAATATTTTGGTACTTc                                                                                                                                            <  2:135271‑M2/16‑1 (MQ=255)
             gATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGggtaatgactcc                                                                                                                                       <  8:3380‑M2/148‑13 (MQ=255)
             gATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGgg                                                                                                                                                 <  6:137913‑M2/138‑3 (MQ=255)
             gATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGgg                                                                                                                                                 >  5:137913‑M2/1‑136 (MQ=255)
                ccGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGg                                                                                                                                                  >  1:97923‑M2/1‑133 (MQ=255)
                ccGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGg                                                                                                                                                  <  2:97923‑M2/134‑2 (MQ=255)
                                 tGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGggtaatgactccaacttattgatagtgttttat                                                                                                                  <  7:58542‑M2/149‑34 (MQ=255)
                                     cacaCGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGggtaatgactccaacttattgatagtgttttatgttc                                                                                                              >  6:27527‑M2/1‑112 (MQ=255)
                                      acacGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGggtaatgactccaacttattgatagtgttttatgttca                                                                                                             <  2:146722‑M2/149‑39 (MQ=255)
                                       cacGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGggtaatgactccaacttattgatagtgttttatgttcag                                                                                                            >  8:58542‑M2/1‑110 (MQ=255)
                                          gcgTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGggtaatgactccaacttattgatagtgttttatgttcagat                                                                                                          <  7:17346‑M2/148‑42 (MQ=255)
                                                 tgaacaggagggacagctgatagaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGc                                                                                                  >  4:98822‑M2/92‑149 (MQ=255)
                                                   caTTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgat                                                                                                >  1:146722‑M2/1‑98 (MQ=255)
                                                    aTTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgat                                                                                                <  5:67767‑M2/148‑52 (MQ=255)
                                                         ggTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgacttt                                                                                          >  7:3380‑M2/1‑92 (MQ=255)
                                                                      agaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccAATATTTTGGTACTTCTGat                                                                                                                                        >  2:17264‑M2/71‑90 (MQ=255)
                                                                      agaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccAATATTTTGGTACTTCTGat                                                                                                                                        <  1:17264‑M2/20‑1 (MQ=255)
                                                                       gaaacagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTg                                                                             <  3:98822‑M2/79‑1 (MQ=255)
                                                                        cccGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGggtaatgactccaacttattgatagtgttttatgttcaga                                                                                                           >  2:81331‑M2/1‑77 (MQ=255)
                                                                        cccGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGggtaatgactccaacttattgatagtgttttatgttcaga                                                                                                           <  1:81331‑M2/117‑41 (MQ=255)
                                                                          acagaagccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTAt                                                                         >  2:111663‑M2/67‑149 (MQ=255)
                                                                               gTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATAACGGCCAGACTCTGGCAAAATATTTTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgatttt                                                                    <  2:29882‑M2/149‑80 (MQ=255)
                                                                                gccactggagcacctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGc                                                                                                  <  8:44795‑M2/58‑1 (MQ=255)
                                                                                    aCCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaa                                                               >  1:29882‑M2/1‑65 (MQ=255)
                                                                                      cGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGAAAAATATTTTGggtaatgactccaacttattgatactgtttaatgttcagataatgaccgatgacttagtcattaagctccaccgatattgagaacg                                                             <  7:24265‑M2/149‑87 (MQ=255)
                                                                                            cctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATaaaa                                                       <  7:85705‑M2/101‑1 (MQ=255)
                                                                                             ctctaaaacaccatcatacactaaatcagtaagttggcagcatcaccAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGtt                                                                                          >  2:128843‑M2/48‑113 (MQ=255)
                                                                                             ctctaaaacaccatcatacactaaatcagtaagttggcagcatcaccAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGtt                                                                                          <  1:128843‑M2/66‑1 (MQ=255)
                                                                                             ctcaaaaacaccatcatacactaaatcagtaagttggcagcatcaccAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATaaaa                                                       >  1:65575‑M2/48‑148 (MQ=255)
                                                                                                 ggTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgactt                                                   >  6:67767‑M2/1‑52 (MQ=255)
                                                                                                  aaacaccatcatacactaaatcagtaagttggcagcatcaccAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAAc                                                      <  8:1539‑M2/102‑1 (MQ=255)
                                                                                                  aaacaccatcatacactaaatcagtaagttggcagcatcaccAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAAc                                                      >  7:1539‑M2/43‑144 (MQ=255)
                                                                                                     aTCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgt                                                                                        >  8:41997‑M2/1‑48 (MQ=255)
                                                                                                     aTCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgt                                                                                        <  7:41997‑M2/107‑60 (MQ=255)
                                                                                                      accatcatacactaaatcagtaagttggcagcatcaccAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGa                                                            >  5:139226‑M2/39‑134 (MQ=255)
                                                                                                      accatcatacactaaatcagtaagttggcagcatcaccAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGa                                                            <  6:139226‑M2/96‑1 (MQ=255)
                                                                                                          tcatacactaaatcagtaagttggcagcatcaccAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCca                                         >  7:84468‑M2/35‑149 (MQ=255)
                                                                                                            aaCCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtccc                                            >  4:73747‑M2/1‑41 (MQ=255)
                                                                                                            aaCCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtccc                                            <  3:73747‑M2/144‑104 (MQ=255)
                                                                                                            aaCCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGggaaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccg                                       <  4:30930‑M2/149‑109 (MQ=255)
                                                                                                             tacactaaatcagtaagttggcagcatcaccAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGa                                                                                                             >  2:114000‑M2/32‑78 (MQ=255)
                                                                                                             tacactaaatcagtaagttggcagcatcaccAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGa                                                                                                             <  1:114000‑M2/47‑1 (MQ=255)
                                                                                                                   gcgGATATCGGCCAGACTCTGGCAAAATATTTTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccag                                 >  8:17346‑M2/1‑34 (MQ=255)
                                                                                                                    aatcagtaagttggcagcatcaccAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGa                                                            >  2:94882‑M2/25‑120 (MQ=255)
                                                                                                                    aatcagtaagttggcagcatcaccAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGa                                                            <  1:94882‑M2/96‑1 (MQ=255)
                                                                                                                          taagttggcagcatcaccAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACaa                                                    <  8:72314‑M2/104‑1 (MQ=255)
                                                                                                                          taagttggcagcatcaccAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACaa                                                    >  7:72314‑M2/19‑122 (MQ=255)
                                                                                                                                  aCTCTGGCAAAATATTTTGggtaatgactccatcttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcaga                    >  2:121279‑M2/1‑19 (MQ=255)
                                                                                                                                  aCTCTGGCAAAATATTTTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgacccagccgtgccaggtgctgcctcaga                    <  1:121279‑M2/146‑128 (MQ=255)
                                                                                                                                    gcatcaccAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAg                                                              <  7:2543‑M2/94‑1 (MQ=255)
                                                                                                                                    gcatcaccAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAg                                                              >  8:2543‑M2/9‑102 (MQ=255)
                                                                                                                                    gcatcaccAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATg                                                  >  8:43721‑M2/9‑114 (MQ=255)
                                                                                                                                       tcaccAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGc            >  5:32860‑M2/6‑149 (MQ=255)
                                                                                                                                       tcaccAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGc            >  3:25543‑M2/6‑149 (MQ=255)
                                                                                                                                        caccAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCt           >  7:57538‑M2/5‑149 (MQ=255)
                                                                                                                                          aaaaTATTTTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatg         <  6:145645‑M2/149‑139 (MQ=255)
                                                                                                                                          aaaaTATTTTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgt                                      <  6:105342‑M2/120‑110 (MQ=255)
                                                                                                                                           aaaTATTTTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggt             >  6:14046‑M2/1‑10 (MQ=255)
                                                                                                                                           aaaTATTTTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggt             <  5:14046‑M2/144‑135 (MQ=255)
                                                                                                                                                  ccgggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtcccagccgtgccaggtgctgcctcagattcaggttatgccgctca  >  3:11562‑M2/3‑3 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTT  >  NC_000913/4620640‑4620933

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