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 52 [0] [0] 52 deoA thymidine phosphorylase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 46194770 (0.000)52 (1.040) 38/190 0.2 100% coding (1249/1323 nt) deoA thymidine phosphorylase
?NC_000913 4619487 = 0 (0.000)coding (1259/1323 nt) deoA thymidine phosphorylase

GACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTG  >  NC_000913/4619351‑4619588
                                                                                                                              |                                                                                                               
gACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttg                                                                                                     >  1:307208‑M1/1‑127 (MQ=255)
       tCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataa                                                                                              >  1:1076953‑M1/1‑120 (MQ=255)
         gATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaag                                                                                            <  2:307208‑M1/139‑22 (MQ=255)
                 cGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaa                                                                                    >  1:652259‑M1/1‑110 (MQ=255)
                    cGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagc                                                                                 >  2:1179516‑M1/1‑107 (MQ=255)
                      gCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcac                                                                               <  1:730123‑M1/139‑35 (MQ=255)
                         ttACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacacc                                                                            <  2:628937‑M1/139‑38 (MQ=255)
                             tGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccgtaaagcacaccaact                                                                        >  2:182646‑M1/1‑98 (MQ=255)
                                 atGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtct                                                                    <  1:741722‑M1/139‑46 (MQ=255)
                                 atGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgc                                                                                                    <  2:225535‑M1/107‑14 (MQ=255)
                                 atGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgc                                                                                                    >  1:225535‑M1/1‑94 (MQ=255)
                                   ggCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctat                                                                  <  2:361436‑M1/139‑48 (MQ=255)
                                    gcgcGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatc                                                                 >  1:1116315‑M1/1‑91 (MQ=255)
                                       cgTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgcc                                                              >  2:1077823‑M1/1‑88 (MQ=255)
                                       cgTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgcc                                                              >  1:58624‑M1/1‑88 (MQ=255)
                                        gTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccg                                                             <  1:1077823‑M1/139‑53 (MQ=255)
                                         tCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgt                                                            >  1:368623‑M1/1‑86 (MQ=255)
                                         tCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaa                                                                          >  2:188330‑M1/1‑86 (MQ=255)
                                         tCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaa                                                                          <  1:188330‑M1/125‑40 (MQ=255)
                                                gACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcg                                                     <  2:652259‑M1/139‑61 (MQ=255)
                                                  ccAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcggcgaa                                                   <  1:925102‑M1/139‑63 (MQ=255)
                                                  ccAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaa                                                   >  1:32033‑M1/1‑77 (MQ=255)
                                                  ccAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaa                                                   <  2:58624‑M1/139‑63 (MQ=255)
                                                     ggTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataa                                                <  2:32033‑M1/139‑66 (MQ=255)
                                                     ggTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataa                                                <  1:1311986‑M1/139‑66 (MQ=255)
                                                      gTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaa                                                   <  2:1255461‑M1/135‑63 (MQ=255)
                                                      gTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaa                                                   >  1:1255461‑M1/1‑73 (MQ=255)
                                                       tAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacg                                              <  1:1179516‑M1/139‑68 (MQ=255)
                                                       tAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacg                                              <  2:134290‑M1/139‑68 (MQ=255)
                                                         gACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggt                                            <  2:349523‑M1/139‑70 (MQ=255)
                                                            ggTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACTAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtata                                         >  1:719349‑M1/1‑67 (MQ=255)
                                                              tCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaata                                                 <  1:436820‑M1/129‑65 (MQ=255)
                                                              tCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaata                                                 >  2:436820‑M1/1‑65 (MQ=255)
                                                                aGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactga                                     >  1:1203584‑M1/1‑63 (MQ=255)
                                                                    tCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaa                                                   <  2:11905‑M1/121‑63 (MQ=255)
                                                                    tCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaa                                                   >  1:11905‑M1/1‑59 (MQ=255)
                                                                      cGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactgatctgat                               <  2:1241378‑M1/139‑83 (MQ=255)
                                                                        cTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactgatctgatca                             >  1:431573‑M1/1‑55 (MQ=255)
                                                                                      cgcgAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactgatctgatcatttaaatttgaagc               <  2:1116315‑M1/139‑99 (MQ=255)
                                                                                        cgAAAGACGAAAACCACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatctccgtatcagcgaataacggtatactgttctgatcatttaaattcgaagcac             <  1:823454‑M1/139‑101 (MQ=255)
                                                                                         gAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactgatctgatcatttaaatttgaagcact            >  2:1192280‑M1/1‑38 (MQ=255)
                                                                                         gAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaag                                                                                            <  1:710398‑M1/59‑22 (MQ=255)
                                                                                         gAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaag                                                                                            >  2:710398‑M1/1‑38 (MQ=255)
                                                                                            aGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactgatctgatcatttaaatttgaagcactgag         <  2:1076953‑M1/139‑105 (MQ=255)
                                                                                                 aaaaCAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacg                                              <  1:770636‑M1/97‑68 (MQ=255)
                                                                                                 aaaaCAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcac                                                                               <  2:539771‑M1/64‑35 (MQ=255)
                                                                                                 aaaaCAACTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcac                                                                               >  1:539771‑M1/1‑30 (MQ=255)
                                                                                                 aaaaCAACTGGCAGGAAACGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactatctatcgccgtatcagcgaataacg                                              >  2:770636‑M1/1‑30 (MQ=255)
                                                                                                   aacaacTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactgatctgatcatttaaatttgaagcactgagtacggag  <  2:368623‑M1/139‑112 (MQ=255)
                                                                                                   aacaacTGGCAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaataacggtatactgatctgatcatttaaatttgaagcactgagtacggag  <  2:27239‑M1/139‑112 (MQ=255)
                                                                                                            cAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaa                                                   >  1:718755‑M1/1‑19 (MQ=255)
                                                                                                            cAGGAAGCGGCGAAAGCGGcaattaaacttgccgataaagcaccggaaagcacaccaactgtctatcgccgtatcagcgaa                                                   <  2:718755‑M1/81‑63 (MQ=255)
                                                                                                                              |                                                                                                               
GACACCATCGATTACAGCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTG  >  NC_000913/4619351‑4619588

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