Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 4,619,478 Δ9 bp coding (1250‑1258/1323 nt) deoA → thymidine phosphorylase

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4619478 4619486 9 66 [0] [0] 65 deoA thymidine phosphorylase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 46194770 (0.000)65 (1.320) 40/202 0.1 100% coding (1249/1323 nt) deoA thymidine phosphorylase
?NC_000913 4619487 = 0 (0.000)coding (1259/1323 nt) deoA thymidine phosphorylase

CGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTG  >  NC_000913/4619339‑4619588
                                                                                                                                          |                                                                                                               
cGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGG                                                                                                                 >  2:227522‑M1/1‑139 (MQ=255)
   cAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaa                                                                                                              <  2:928383‑M1/139‑4 (MQ=255)
          cTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaact                                                                                                       >  2:548375‑M1/1‑129 (MQ=255)
              cacCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgcc                                                                                                   >  1:1007800‑M1/1‑125 (MQ=255)
                  aTCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgata                                                                                               <  2:986443‑M1/139‑19 (MQ=255)
                    cGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaa                                                                                             >  2:1168960‑M1/1‑119 (MQ=255)
                             cGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaa                                                                                    >  2:843083‑M1/1‑110 (MQ=255)
                                 ggCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagca                                                                                <  2:103253‑M1/139‑34 (MQ=255)
                                 ggCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagca                                                                                >  1:1175430‑M1/1‑106 (MQ=255)
                                             atGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtct                                                                    <  1:1168960‑M1/139‑46 (MQ=255)
                                             atGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaa                                                                          >  2:708269‑M1/1‑94 (MQ=255)
                                             atGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaa                                                                          <  1:708269‑M1/133‑40 (MQ=255)
                                                   cgTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgcc                                                              >  1:1106080‑M1/1‑88 (MQ=255)
                                                     tCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgt                                                            >  1:411899‑M1/1‑86 (MQ=255)
                                                     tCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgt                                                            >  2:79345‑M1/1‑86 (MQ=255)
                                                      cTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgta                                                           <  2:418355‑M1/139‑55 (MQ=255)
                                                       tGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtat                                                          <  1:682384‑M1/139‑56 (MQ=255)
                                                         ggCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatca                                                        >  1:1162893‑M1/1‑82 (MQ=255)
                                                         ggCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagca                                                                                <  2:591115‑M1/115‑34 (MQ=255)
                                                         ggCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagca                                                                                >  1:591115‑M1/1‑82 (MQ=255)
                                                              ccAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaa                                                   <  2:1101605‑M1/139‑63 (MQ=255)
                                                              ccAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaa                                                   <  2:277337‑M1/139‑63 (MQ=255)
                                                              ccAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaa                                                   <  2:411899‑M1/139‑63 (MQ=255)
                                                              ccAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaa                                                   <  2:1162893‑M1/139‑63 (MQ=255)
                                                              ccAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaa                                                   <  2:688759‑M1/139‑63 (MQ=255)
                                                              ccAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaa                                                   <  1:536507‑M1/139‑63 (MQ=255)
                                                              ccAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaa                                                   <  1:265579‑M1/139‑63 (MQ=255)
                                                              ccAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaa                                                   <  1:1058207‑M1/139‑63 (MQ=255)
                                                              ccAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcacc                                                                                        <  1:1196176‑M1/102‑26 (MQ=255)
                                                              ccAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcacc                                                                                        <  1:967113‑M1/102‑26 (MQ=255)
                                                              ccAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcacc                                                                                        >  2:1196176‑M1/1‑77 (MQ=255)
                                                              ccAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcacc                                                                                        >  2:967113‑M1/1‑77 (MQ=255)
                                                               cAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaat                                                  <  2:1007800‑M1/139‑64 (MQ=255)
                                                                aGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaata                                                 <  2:675612‑M1/139‑65 (MQ=255)
                                                                 ggTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataa                                                <  1:102326‑M1/139‑66 (MQ=255)
                                                                   tAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacg                                              <  1:548375‑M1/139‑68 (MQ=255)
                                                                   tAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacg                                              <  2:823284‑M1/139‑68 (MQ=255)
                                                                   tAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacg                                              <  1:79345‑M1/139‑68 (MQ=255)
                                                                   tAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcaccgtatcagcgaataacg                                              <  1:324547‑M1/139‑68 (MQ=255)
                                                                      aCGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggta                                           <  1:843083‑M1/139‑71 (MQ=255)
                                                                        ggTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacg                                              <  1:970356‑M1/134‑68 (MQ=255)
                                                                        ggTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacg                                              >  2:970356‑M1/1‑67 (MQ=255)
                                                                          tCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatact                                       >  1:76510‑M1/1‑65 (MQ=255)
                                                                            aGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaa                                                   <  2:250857‑M1/125‑63 (MQ=255)
                                                                            aGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAGCTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaa                                                   >  1:250857‑M1/1‑63 (MQ=255)
                                                                            aGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacacaaactgtctatcgccgtatcagcgaataacggtatactga                                     >  2:1203812‑M1/1‑63 (MQ=255)
                                                                            aGCGTCCGCTGGCGGTTAACCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactga                                     >  1:676083‑M1/1‑63 (MQ=255)
                                                                                tCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactgatctg                                 <  2:971887‑M1/139‑81 (MQ=255)
                                                                                tCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataa                                                >  1:1158702‑M1/1‑59 (MQ=255)
                                                                                tCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataa                                                <  2:1158702‑M1/124‑66 (MQ=255)
                                                                                      ggCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactgatctgatcatt                           >  2:930483‑M1/1‑53 (MQ=255)
                                                                                        cGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataa                                                >  2:63492‑M1/1‑51 (MQ=255)
                                                                                        cGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataa                                                <  1:63492‑M1/116‑66 (MQ=255)
                                                                                        cGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaa                                                   >  2:1149516‑M1/1‑51 (MQ=255)
                                                                                        cGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaa                                                   <  1:1149516‑M1/113‑63 (MQ=255)
                                                                                         ggTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactgatctgatcatttaa                        <  2:598603‑M1/139‑90 (MQ=255)
                                                                                                 aCGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactgatctgatcatttaaatttgaag                >  1:655788‑M1/1‑42 (MQ=255)
                                                                                                  cgcgAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactgatctgatcatttaaatttgaagc               <  1:1203812‑M1/139‑99 (MQ=255)
                                                                                                    cgAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactgatctgatcatttaaatttgaagcac             <  2:76510‑M1/139‑101 (MQ=255)
                                                                                                       aaGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactgatctgatcatttaaatttgaagcactga          <  2:939883‑M1/139‑104 (MQ=255)
                                                                                                            gAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactgatctgatcatttaaatttgaagcactgagtacg     <  1:227522‑M1/139‑109 (MQ=255)
                                                                                                            gAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactgatctgatcatttaaatttgaagcactgagtacg     <  1:930483‑M1/139‑109 (MQ=255)
                                                                                                               aacaacTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactgatctgatcatttaaatttgaagcactgagtacggag  <  2:676083‑M1/139‑112 (MQ=255)
                                                                                                               aacaacTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccatatcagcgaataacggtatactgatctgatcatttaaatttgaagcactgtgtacggag  >  1:739622‑M1/1‑28 (MQ=255)
                                                                                                                    cTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactgatctgatcatttaaatttgaagc               <  1:1124913‑M1/121‑99 (MQ=255)
                                                                                                                    cTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactgatctgatcatttaaatttgaagc               >  2:1124913‑M1/1‑23 (MQ=255)
                                                                                                                                          |                                                                                                               
CGTCAGGCATCTGACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTG  >  NC_000913/4619339‑4619588

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