Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_010473 2,806,542 689 bp→GAGCGGCCGCCACCGCGGGAAGTTCCTATACTTTCTAGAGAATAGGAACTTCGCATGCACGCAGCATATGC coding (2‑690/690 nt) ung → uracil‑DNA glycosylase

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_010473 2806542 2807230 689 21 [0] [0] 19 ung uracil‑DNA glycosylase

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_010473 = 28065410 (0.000)15 (0.610)
+GAGCGGCCGCCACCGCGGGAAGTTCCTATACTTTCTAGAGAATAGGAACTTCGCATGCACGCAGCATATGC
13/256 1.2 100% coding (1/690 nt) ung uracil‑DNA glycosylase
?NC_010473 2807231 = 0 (0.000)intergenic (+1/+47) ung/yfiF uracil‑DNA glycosylase/tRNA/rRNA methyltransferase

TCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGAAAGCACCGCATCCGTCGCCGCTTTCGGCGCATCGTGGATTCTTTGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATGCCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGATGGCAGAGTTGCCACCCGGCTGATTTATCAGGCTTTATTCTGACGCCACCATTCACCAAGCAAAACGCCGGTTGCGACAGAGATATTCAGCCCGGCAACGTTGCCCGTACCGTCAATCTTCACGCGCAGATCGTTCGGATCGCGTGCGGCATC  >  NC_010473/2807036‑2807403
                                                                                                                                                                                                   |                                                                                                                                                                             
gggtgtgaaattttaccaatttacatttttttgcactcgtttaagtctaaaaaatgagcactattttgttctgtagaaagaagcagttaagctaggcggattgaagattcgcaggagagcgagagagcggccgccaccgcgggaagttcctatactttctagagaataggaacttcgcatgcacgcagcatatgcATTTg                                                                                                                                                                           <  1:12743S67‑M2/5‑1 (MQ=255)
      tttaccaatttacatttttttgcactcgtttaagtctaaaaaatgagcatgatttgtttttcttgtagaaagaagcagttaagctaggcggattgaagattcgcaggagagcgagagagcggccgccaccgcgggagtttctctgatactttctagagaataggaacttcgcatgcacgcagcatatgcATTtgcggtaa                                                                                                                                                                    <  1:10546S61‑M2/11‑4 (MQ=255)
           tttaccaatttacatttttttgcactcgtttaagtctaaaaaatgagcatgattttgttctgtagaaagaagcagttaagctaggcggattgaagattcgcaggagagcgagagagcggccgccaccgcgggaagttcctatactttctagagaataggaacttcgcatgcacgcagcatatgcATTTGCGGGGAAATGc                                                                                                                                                               >  1:1907S78‑M2/185‑200 (MQ=255)
            ttaccaatttacatttttttgcactcgtttaagtctaaaaaatgagcatgattttgttctgtagaaagaagcagttaagctaggcggattgaagattcgcaggagagcgagagagcggccgccaccgcgggaagttcctatactttctagagaataggaacttcgcatgcacgcagcatatgcATTTGCGGGAAAATGcc                                                                                                                                                              <  1:32087S78‑M2/17‑1 (MQ=255)
                      acatttttttgcactcgtttaagtctaaaaaatgagcatgattttgttctgtagaaagaagcagttaagctaggcggattgaagattcgcaggagagcgagagagcggccgccaccgcgggaagttcctatactttctagagaataggaacttcgcatgcacgcagcatatgcatttgcggcatgcccggatggcagagt                                                                                                                                                    <  1:14235S5‑M2/27‑20 (MQ=255)
                                  cactcgtttaagtctaaaaaatgagcatgattttgttctgtagaaagaagcagttaagctaggcggattgaagattcgcaggagagcgagagagcggccgccaccgcgggaagttctatactttctagagaataggaacttcgcatgcacgcagcatatgcATTTGCGGGGAAATGCCGGATGGCAGAGTTGCCACCCgg                                                                                                                                        >  1:5567S2‑M2/162‑200 (MQ=255)
                                      gtttaagtctaaaaaatgagcatgattttgttctgtagaaagaagcagttaagctaggcggattgaagattcgcaggagagcgagagagcggccgccaccgcgggaagttcctatactttctagagaataggaacttcgcatgcacgcagcatatgcATTTGCGGGGGAAATGCCGGATGGCAGAGTTGCCACCCGGCTg                                                                                                                                     <  1:2149S1‑M2/43‑1 (MQ=255)
                                                                    ttctgtagaaagaagcagttaagctaggcggattgaagattcgcaggagagcgagagagcggccgccaccgcgggaagttcctatactttctagagaataggaacttcgcatgcacgcagcatatgcATTTGCGGGGAAATGCCGGATGGCAGAGTTGCCACCCGGCTGATTTATCAGGCTTTATTCTGACGCCACCAtt                                                                                                      <  1:3103S50‑M2/73‑1 (MQ=255)
                                                                      tgtagaaagaagcagttaagctaggcggattgagagattcgcaggagagcgagagagcggccgccaccgcgggaagttcctatactttctagagaataggaacttcgcatgcacgcagcatatgcATTTGCGGGGAAATGCCGGATGGCAAACTGCCACCCGGCTGATTTATCAGGCTTTATTCTGACGCCACCATTCAc                                                                                                   <  1:9708S93‑M2/75‑1 (MQ=255)
                                                                            aagaagcagttaagctaggcggattgaagattcgcaggagagcgagagagctggccgccaccgcgggaagttcctatactttctagagaataggaacttcgcatgcacgcagcatatgcATTTGCGGGGAAATGCCGGATGGCAGAGTTGCCACCCGGCTGATTTATCAGGCTTTATTCTGACGCCACCATTCACCAAGc                                                                                              <  1:1894S17‑M2/81‑1 (MQ=255)
                                                                                  aagaagcagttaagctaggcggattgaagattcgcaggagagcgagagagcggccgccaccgcgggaagttcctatactttctaaatagaactcgcatgcacgcagcatatgcATTTGCGGGGAAATGCCGGATGGCAGAGTTGCCACCCGGCTGATTTATCAGGCTTTATTCTGACGCCACCATTCACCAAGCAAAACg                                                                                        >  1:23878S24‑M2/114‑200 (MQ=255)
                                                                                    agttaagctaggcggattgaagattcgcaggagagcgagagagcggccgccaccgcgggaagttcctatactttctagagaataggaacttcgcatgcacgcagcatatgcATTTGCGGGGAAATGCCGGATGGCAGAGTTGCCACCCGGCTGATTTATCAGGCTTTATTCTGACGCCACCATTCACCAAGCAAAACGcc                                                                                      <  1:18994S60‑M2/89‑1 (MQ=255)
                                                                                              ggcggattgaagattcgcaggagagcgagagagcggccgccaccgcgggaagttcctatactttctagagaataggaacttcgcgtgcacgcagcatatgcATTTGCGGGGAAATGCCGGATGGCAGAGTTGCCACCCGGCTGATTTATCAGGCTTTATTCTGACGCCACCATTCACCAAGCAAAACGCCGGTTGCGACa                                                                            >  1:723S1‑M2/102‑200 (MQ=255)
                                                                                                 aggcggattgaagattcgcaggagagcgagagagcggccgccaccgcgggaagttcctatactttctagagaataggaacttcgcataggcatattgcGTTTGCGGGGAAATGCCGGATGGCAGAGTTGCCACCCGGCTGATTTATCAGGCTTTATTCTGACGCCACCATTCACCAAGCAAAACGCCGGTTGCGACagag                                                                         >  1:14233S10‑M2/99‑200 (MQ=255)
                                                                                                                  aggagagcgagagagcggccgccaccgcgggagttcaatactttctagagaataggaacttcgcatgcacgcagcatatgcATTTGCGGGGAAATGCCGGATGGCAGAGTTGCCACCCGGCTGATTTATCAGGCTTTATTCTGACGCCACCATTCACCAAGCAAAACGCCGGTTGCGACAGAGATATTCAGCCCGGCAAc                                                        <  1:29180S16‑M2/119‑1 (MQ=255)
                                                                                                                                       caccgcggaagttcctatactttctagagaataggaacttcgcatgcacgcagcatatgcATTTGCGGGGAAATGCCGGATGGCAGAGTTGCCACCCGGCTGATTTATCAGGCTTTATTCTGACGCCACCATTCACCAAGCAAAACGCCGGTTGCGACAGAGATATTCAGCCCGGCAACGTTGCCCGTACCGTCAATCtt                                   <  1:6781S6‑M2/140‑1 (MQ=255)
                                                                                                                                               aagttcctatactttctagagaataggaacttcgcatgcacgcagcatatgcATTTGCGGGGAAATGCCGGATGGCAGAGTTGCCACCCGGCTGATTTATCAGGCTTTATTCTGACGCCACCATTCACCAAGCAAAACGCCGGTTGCGACAGAGATATTCAGCCCGGCAACGTTGCCCGTACCGTCAATCTTCACGCGCa                           <  1:27835S206‑M2/148‑1 (MQ=255)
                                                                                                                                                                agagaataggaacttcgcatgcacgcagcatatgcATTTGCGGGGAAATGCCGGATGGCAGAGTTGCCACCCGGCTGATTTATCAGGCTTTATTCTGACGCCACCATTCACCAAGCAAAACGCCGGTTGCGACAGAGATATTCAGCCCGGCAACGTTGCCCGTACCGTCAATCTTCACGCGCAGATCGTTCGGATCGCGt          >  1:24699S7‑M2/36‑200 (MQ=255)
                                                                                                                                                                        ggaacttcgcatgcacgcagcatatgcATTTGCGGGGAAATGCCGGATGGCAGAGTTGCCACCCGGCTGATTTATCAGGCTTTATTCTGACGCCACCATTCACCAAGCAAAACGCCGGTTGCGACAGAGATATTCAGCCCGGCAACGTTGCCCGTACCGTCAATCTTCACGCGCAGATCGTTCGGATCGCGTGCGGCATc  >  1:12855S15‑M2/28‑200 (MQ=255)
                                                                                                                                                                                                   |                                                                                                                                                                             
TCGCATGCGCAAAAGAAAGGGGCGATTATAGATAAGCAACGCCATCATGTACTGAAAGCACCGCATCCGTCGCCGCTTTCGGCGCATCGTGGATTCTTTGGCTGCAACCATTTTGTGCTGGCAAATCAGTGGCTGGAACAACGTGGCGAGACGCCGATTGACTGGATGCCAGTATTACCGGCAGAGAGTGAGTAAATTTGCGGGGAAATGCCGGATGGCAGAGTTGCCACCCGGCTGATTTATCAGGCTTTATTCTGACGCCACCATTCACCAAGCAAAACGCCGGTTGCGACAGAGATATTCAGCCCGGCAACGTTGCCCGTACCGTCAATCTTCACGCGCAGATCGTTCGGATCGCGTGCGGCATC  >  NC_010473/2807036‑2807403

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 48 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: