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)30 (1.020) 27/258 0.2 100% noncoding (768/768 nt) IS1 repeat region
?NC_000913 = 4620788 0 (0.000)coding (1186/1224 nt) deoB phosphopentomutase

TCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGA  >  NC_000913/4620633‑4620926
                                                                                                                                                   |                                                                                                                                                  
tCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCaaaa                                                                                                                                                   >  5:27527/1‑149 (MQ=255)
      cTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATAtttt                                                                                                                                             <  5:71798/149‑1 (MQ=255)
       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                                                                                                                                          >  5:137913‑M2/1‑136 (MQ=255)
                    gATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGgg                                                                                                                                          <  6:137913‑M2/138‑3 (MQ=255)
                       ccGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGg                                                                                                                                           <  2:97923‑M2/134‑2 (MQ=255)
                       ccGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGg                                                                                                                                           >  1:97923‑M2/1‑133 (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)
                                                         acaTTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATAtttt                                                                                                                                             >  1:42429/1‑98 (MQ=255)
                                                         acaTTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATAtttt                                                                                                                                             <  2:42429/98‑1 (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                                                                                   <  1:128843‑M2/66‑1 (MQ=255)
                                                                                                    ctctaaaacaccatcatacactaaatcagtaagttggcagcatcaccAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGtt                                                                                   >  2:128843‑M2/48‑113 (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                                                                                 <  7:41997‑M2/107‑60 (MQ=255)
                                                                                                            aTCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgt                                                                                 >  8:41997‑M2/1‑48 (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                                     <  3:73747‑M2/144‑104 (MQ=255)
                                                                                                                   aaCCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaacgacagcgacttccgtccc                                     >  4:73747‑M2/1‑41 (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)
                                                                                                                                                   |                                                                                                                                                  
TCACCGCTGACCACGGTTGCGATCCGACCTGGACCGGTACTGACCACACGCGTGAACACATTCCGGTACTGGTATATGGCCCGAAAGTAAAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGA  >  NC_000913/4620633‑4620926

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.