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)45 (1.450) 32/190 0.1 100% coding (1249/1323 nt) deoA thymidine phosphorylase
?NC_000913 4619487 = 0 (0.000)coding (1259/1323 nt) deoA thymidine phosphorylase

GCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTA  >  NC_000913/4619367‑4619621
                                                                                                                        |                                                                                                                                      
tcgattacagcgtcggctttactgatatggcgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATaa                                                                                                                      >  2:615959‑M2/121‑139 (MQ=255)
          cgtcggctttactgatatggcgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGaa                                                                                                            >  2:265964‑M2/111‑139 (MQ=255)
          cgtcggctttactgatatggcgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGaa                                                                                                            >  1:356991‑M2/111‑139 (MQ=255)
             cggctttactgatatggcgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGc                                                                                                         >  2:28475‑M2/108‑139 (MQ=255)
              ggctttactgatatggcgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGca                                                                                                        >  2:406409‑M2/107‑139 (MQ=255)
                    actgatatggcgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGca                                                                                                        >  1:762338‑M2/101‑133 (MQ=255)
                    actgatatggcgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGca                                                                                                        <  2:762338‑M2/33‑1 (MQ=255)
                          atggcgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCt                                                                                            >  1:108932‑M2/95‑139 (MQ=255)
                             gcgcgtctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATc                                                                                         >  1:72815‑M2/92‑139 (MQ=255)
                                  tctgggcgaccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGt                                                                                    >  2:787665‑M2/87‑139 (MQ=255)
                                           ccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGaa                                                                           <  2:108932‑M2/62‑1 (MQ=255)
                                           ccaggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGaa                                                                           <  1:265964‑M2/62‑1 (MQ=255)
                                             aggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa                                                                         <  2:209520‑M2/64‑1 (MQ=255)
                                              ggtagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATaa                                                                        <  1:663379‑M2/65‑1 (MQ=255)
                                                tagacggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACg                                                                      <  2:137691‑M2/67‑1 (MQ=255)
                                                   acggtcagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGta                                                                   <  2:293924‑M2/70‑1 (MQ=255)
                                                        cagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATaa                                                                        >  2:291281‑M2/65‑129 (MQ=255)
                                                        cagcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATaa                                                                        <  1:291281‑M2/65‑1 (MQ=255)
                                                         agcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCaa                                                                                                  <  2:98948‑M2/39‑1 (MQ=255)
                                                         agcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACg                                                                      <  2:405440‑M2/67‑1 (MQ=255)
                                                         agcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACg                                                                      >  1:405440‑M2/64‑130 (MQ=255)
                                                         agcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATActga                                                             >  1:314086‑M2/64‑139 (MQ=255)
                                                         agcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATActga                                                             >  1:38569‑M2/64‑139 (MQ=255)
                                                         agcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGAACACCaa                                                                                                  >  1:98948‑M2/64‑102 (MQ=255)
                                                          gcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATActgat                                                            <  1:787665‑M2/77‑1 (MQ=255)
                                                          gcgtccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATActgat                                                            <  1:615959‑M2/77‑1 (MQ=255)
                                                             tccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATActgatctg                                                         >  2:255837‑M2/60‑139 (MQ=255)
                                                              ccgctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggagaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATActgatctga                                                        <  1:255837‑M2/81‑1 (MQ=255)
                                                                 ctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGaa                                                                           <  2:474424‑M2/62‑1 (MQ=255)
                                                                 ctggcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGaa                                                                           >  1:474424‑M2/56‑117 (MQ=255)
                                                                    gcggttatccacgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCAttt                                                  <  2:310583‑M2/87‑1 (MQ=255)
                                                                               cgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGaa                                                                           <  1:330474‑M2/62‑1 (MQ=255)
                                                                               cgcgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGaa                                                                           >  2:330474‑M2/42‑103 (MQ=255)
                                                                                 cgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGt                                                                                              >  2:694927‑M2/40‑82 (MQ=255)
                                                                                 cgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGt                                                                                              <  1:694927‑M2/43‑1 (MQ=255)
                                                                                 cgaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCAc                                     <  2:356991‑M2/100‑1 (MQ=255)
                                                                                  gaaagacgaaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACt                                    <  2:314086‑M2/101‑1 (MQ=255)
                                                                                          aaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACg                                                                      <  1:452160‑M2/67‑1 (MQ=255)
                                                                                          aaaacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACg                                                                      >  2:452160‑M2/31‑97 (MQ=255)
                                                                                            aacaactggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGgag                          <  1:28475‑M2/111‑1 (MQ=255)
                                                                                                 ctggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGgag                          >  2:662337‑M2/24‑134 (MQ=255)
                                                                                                 ctggcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGgag                          <  1:662337‑M2/111‑1 (MQ=255)
                                                                                                    gcaggaagcggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATg                  >  2:740830‑M2/21‑139 (MQ=255)
                                                                                                            cggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACg                                                                      <  1:478923‑M2/67‑1 (MQ=255)
                                                                                                            cggcgaaagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACg                                                                      >  2:478923‑M2/13‑79 (MQ=255)
                                                                                                                  aagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAAt                                                                          <  1:555411‑M2/63‑1 (MQ=255)
                                                                                                                  aagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACg                                                                      <  1:769342‑M2/67‑1 (MQ=255)
                                                                                                                  aagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACg                                                                      >  2:769342‑M2/7‑73 (MQ=255)
                                                                                                                  aagcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCAGTATCAGCGAAt                                                                          >  2:555411‑M2/7‑69 (MQ=255)
                                                                                                                    gcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATaa                                                                        <  2:602547‑M2/65‑1 (MQ=255)
                                                                                                                    gcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATaa                                                                        >  1:602547‑M2/5‑69 (MQ=255)
                                                                                                                    gcggCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATttatta  <  2:72815‑M2/135‑1 (MQ=255)
                                                                                                                        |                                                                                                                                      
GCGTCGGCTTTACTGATATGGCGCGTCTGGGCGACCAGGTAGACGGTCAGCGTCCGCTGGCGGTTATCCACGCGAAAGACGAAAACAACTGGCAGGAAGCGGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAACGGTATACTGATCTGATCATTTAAATTTGAAGCACTGAGTACGGAGAACATATGAAACGTGCATTTATTA  >  NC_000913/4619367‑4619621

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: