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

ATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCAAATGGATTTTATTGTTTTACATCAACTTATGCGGGTGTGAAATTTTACCAATTTACATTTTTTTGCACTCGTTTAAGTCTAAAAAATGAGCATGATT‑‑TTGTTCTGTAGAAAGAAGCAGTTAAGCTA‑‑GGCGGATTGAAGATTCGCAGGAGAGCGAGATGGCTAACGAATTAACCTGGCATGACGTGCTGGCTGAAGAGAAGCAGCAACCCTATTTTCTTAATACCCTTCAGACCGTCGCCAGCGAGCGGCAGTCCGGCGTCACTATCTACCCACCACAAAAAGATGTCTTTAACGCGTTCCGCTTTACAGAGTTGGGTGACGTTAAAGTGGTGATTCTCGGCCAG  >  NC_010473/2806340‑2806729
                                                                                                                                                                                                                 |                                                                                                                                                                                            
aTCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCAAATGGATTTTATTGTTTTACATCAACTTATGCGGGTGTGAAATTTTACCAATTTACATTTTTTTGCACTCGTTTAAGTCTAAAAAATGAGCATGATT‑‑TTGTTCTGTAGAAAGAAGCAGTTAGG‑‑A‑‑GGTGGATTGAAGATTCGCAGGAGAGCgaga                                                                                                                                                                                              <  1:1812S2/200‑1 (MQ=255)
           aaaTCAATAAAAACCACACATCAAGTATGGTCGCAAATGGATTTTATTGTTTTCATCAACTTATGCGGGTGTGAAATTTTACCAATTTACATTTTTTTGCACTCGTTTAAGTCTAAAAAATGAGCATGATT‑‑TTGTTCTGTAGAAAGAAGCAGTTAAGCTA‑‑GGCGGATTGAAGATTCGCAGGAGAGCGAGAgagcggccgc                                                                                                                                                                                    <  1:6781S7‑M1/200‑11 (MQ=255)
                     aaaCCACACATCAAGTATGGTCGCAAATGGATTTTATTGTTTTACATCAACTTATGCGGGTGTGAAATTTTACCGAATTTACATTTTTTTGCACTCGTTTAAGTCTAAAAAATGAGCATGATT‑‑TTGTTCTGTAGAAAGAAGCAGTTAAGCTA‑‑GGCGGATTGAAGATTCGCAGGAGAGCGAGAgagcggccgccaccgcgg                                                                                                                                                                            <  1:27835S207‑M1/200‑19 (MQ=255)
                                      tGGTCGCAAATGGATTTTATTGTTTTACATCAACTTATGCGGGTGTGAAATTTTACCAATTTACATTTTTTTGCACTCGTTTAAGTCTAAAAAATGAGCATGATT‑‑TTGTTCTGTAGAAAGAAGCAGTTAAGCTA‑‑GGCGGATTGAAGATTCGCAGGAGAGCGAGAgagcggccgccaccgcgggaagttcctatactttct                                                                                                                                                          >  1:24699S6‑M1/1‑164 (MQ=255)
                                       ggTCGCAAATGGATTTTATTGTTTTACATCAACTTATGCGGGTGTGAAATTTTACCAATTTACATTTTTTTGCACTCGTTTAAGTCTAAAAAATGACATGATT‑‑TTGTTCTGTAGAAAGAAGCAGTTAAGCTA‑‑GGCGGATTGAAGATTCGCAGGAGAGCGAGAgagcggccgccaccgcgggaagttcctatactttctag                                                                                                                                                        >  1:29904S18‑M1/1‑162 (MQ=255)
                                               aTGGATTTTATTGTTTTACATCAACTTATGCGGGTGTGAAATTTTACCAATTTACATTTTTTTGCACTCGTTTAAGTCTAAAAAATGAGCATGATT‑‑TTGTTCTGTAGAAAGAAGCAGTTAAGCTACGGAGCGGATTGAAGATTCGCAGGAGAGCGAGAgagcggccgccaccgcgggaagttcctatactttctagaata                                                                                                                                                    >  1:12855S14‑M1/1‑158 (MQ=255)
                                                                              gggTGTGAAATTTTACCAATTTACATTTTTTTGCACTCGTTTAAGTCTAAAAAATGAGCACTATT‑‑TTGTTCTGTAGAAAGAAGCAGTTAAGCTA‑‑GGCGGATTGAAGATTCGCAGGAGAGCGAGAgagcggccgccaccgcgggaagttcctatactttctagagaataggaacttcgcatgcacgcagcatatgcatttg                                                                                                                  <  1:12743S67‑M1/200‑77 (MQ=255)
                                                                                         tttACCAATTTACATTTTTTTGCACTCGTTTAAGTCTAAAAAATGAGCATGATTTGTTTTTCTTGTAGAAAGAAGCAGTTAAGCTA‑‑GGCGGATTGAAGATTCGCAGGAGAGCGAGAgagcggccgccaccgcgggagtttctctgatactttctagagaataggaacttcgcatgcacgcagcatatgcatttgcggtaa                                                                                                          <  1:10546S61‑M1/200‑85 (MQ=255)
                                                                                         tttACCAATTTACATTTTTTTGCACTCGTTTAAGTCTAAAAAATGAGCATGATT‑‑TTGTTCTGTAGAAAGAAGCAGTTAAGCTA‑‑GGCGGATTGAAGATTCGCAGGAGAGCGAGAgagcggccgccaccgcgggaagttcctatactttctagagaataggaacttcgcatgcacgcagcatatgcatttgcggggaaatgc                                                                                                       >  1:1907S78‑M1/1‑113 (MQ=255)
                                                                                          ttACCAATTTACATTTTTTTGCACTCGTTTAAGTCTAAAAAATGAGCATGATT‑‑TTGTTCTGTAGAAAGAAGCAGTTAAGCTA‑‑GGCGGATTGAAGATTCGCAGGAGAGCGAGAgagcggccgccaccgcgggaagttcctatactttctagagaataggaacttcgcatgcacgcagcatatgcatttgcgggaaaatgcc                                                                                                      <  1:32087S78‑M1/200‑89 (MQ=255)
                                                                                                     aCATTTTTTTGCACTCGTTTAAGTCTAAAAAATGAGCATGATT‑‑TTGTTCTGTAGAAAGAAGCAGTTAAGCTA‑‑GGCGGATTGAAGATTCGCAGGAGAGCGAGAgagcggccgccaccgcgggaagttcctatactttctagagaataggaacttcgcatgcacgcagcatatgcatttgcggcatgcccggatggcagagt                                                                                            <  1:14235S5‑M1/200‑99 (MQ=255)
                                                                                                                cACTCGTTTAAGTCTAAAAAATGAGCATGATT‑‑TTGTTCTGTAGAAAGAAGCAGTTAAGCTA‑‑GGCGGATTGAAGATTCGCAGGAGAGCGAGAgagcggccgccaccgcgggaagttctatactttctagagaataggaacttcgcatgcacgcagcatatgcatttgcggggaaatgccggatggcagagttgccacccgg                                                                                 >  1:5567S2‑M1/1‑91 (MQ=255)
                                                                                                                     gTTTAAGTCTAAAAAATGAGCATGATT‑‑TTGTTCTGTAGAAAGAAGCAGTTAAGCTA‑‑GGCGGATTGAAGATTCGCAGGAGAGCGAGAgagcggccgccaccgcgggaagttcctatactttctagagaataggaacttcgcatgcacgcagcatatgcatttgcgggggaaatgccggatggcagagttgccacccggctg                                                                            <  1:2149S1‑M1/200‑115 (MQ=255)
                                                                                                                                                     ttCTGTAGAAAGAAGCAGTTAAGCTA‑‑GGCGGATTGAAGATTCGCAGGAGAGCGAGAgagcggccgccaccgcgggaagttcctatactttctagagaataggaacttcgcatgcacgcagcatatgcatttgcggggaaatgccggatggcagagttgccacccggctgatttatcaggctttattctgacgccaccatt                                              <  1:3103S50‑M1/200‑145 (MQ=255)
                                                                                                                                                         tGTAGAAAGAAGCAGTTAAGCTA‑‑GGCGGATTGAGAGATTCGCAGGAGAGCGAGAgagcggccgccaccgcgggaagttcctatactttctagagaataggaacttcgcatgcacgcagcatatgcatttgcggggaaatgccggatggcaaactgccacccggctgatttatcaggctttattctgacgccaccattcac                                            <  1:9708S93‑M1/200‑147 (MQ=255)
                                                                                                                                                               aagaagCAGTTAAGCTA‑‑GGCGGATTGAAGATTCGCAGGAGAGCGAGAgagctggccgccaccgcgggaagttcctatactttctagagaataggaacttcgcatgcacgcagcatatgcatttgcggggaaatgccggatggcagagttgccacccggctgatttatcaggctttattctgacgccaccattcaccaagc                                     <  1:1894S17‑M1/200‑154 (MQ=255)
                                                                                                                                                               aagaagCAGTTAAGCTA‑‑GGCGGATTGAAGATTCGCAGGAGAGCGAGAgagcggccgccaccgcgggaagttcctatactttctaaatagaactcgcatgcacgcagcatatgcatttgcggggaaatgccggatggcagagttgccacccggctgatttatcaggctttattctgacgccaccattcaccaagcaaaacg                                     >  1:23878S24‑M1/1‑47 (MQ=255)
                                                                                                                                                                      aGTTAAGCTA‑‑GGCGGATTGAAGATTCGCAGGAGAGCGAGAgagcggccgccaccgcgggaagttcctatactttctagagaataggaacttcgcatgcacgcagcatatgcatttgcggggaaatgccggatggcagagttgccacccggctgatttatcaggctttattctgacgccaccattcaccaagcaaaacgcc                              <  1:18994S60‑M1/200‑161 (MQ=255)
                                                                                                                                                                               a‑‑GGCGGATTGAAGATTCGCAGGAGAGCGAGAgagcggccgccaccgcgggaagttcctatactttctagagaataggaacttcgcataggcatattgcgtttgcggggaaatgccggatggcagagttgccacccggctgatttatcaggctttattctgacgccaccattcaccaagcaaaacgccggttgcgacagag                     >  1:14233S10‑M1/1‑31 (MQ=255)
                                                                                                                                                                                  ggCGGATTGAAGATTCGCAGGAGAGCGAGAgagcggccgccaccgcgggaagttcctatactttctagagaataggaacttcgcgtgcacgcagcatatgcatttgcggggaaatgccggatggcagagttgccacccggctgatttatcaggctttattctgacgccaccattcaccaagcaaaacgccggttgcgaca                    >  1:723S1‑M1/1‑30 (MQ=255)
                                                                                                                                                                                                      aGGAGAGCGAGAgagcggccgccaccgcgggagttcaatactttctagagaataggaacttcgcatgcacgcagcatatgcatttgcggggaaatgccggatggcagagttgccacccggctgatttatcaggctttattctgacgccaccattcaccaagcaaaacgccggttgcgacagagatattcagcccggcaac  <  1:29180S16‑M1/200‑189 (MQ=255)
                                                                                                                                                                                                                 |                                                                                                                                                                                            
ATCTTTGATTTAAATCAATAAAAACCACACATCAAGTATGGTCGCAAATGGATTTTATTGTTTTACATCAACTTATGCGGGTGTGAAATTTTACCAATTTACATTTTTTTGCACTCGTTTAAGTCTAAAAAATGAGCATGATT‑‑TTGTTCTGTAGAAAGAAGCAGTTAAGCTA‑‑GGCGGATTGAAGATTCGCAGGAGAGCGAGATGGCTAACGAATTAACCTGGCATGACGTGCTGGCTGAAGAGAAGCAGCAACCCTATTTTCTTAATACCCTTCAGACCGTCGCCAGCGAGCGGCAGTCCGGCGTCACTATCTACCCACCACAAAAAGATGTCTTTAACGCGTTCCGCTTTACAGAGTTGGGTGACGTTAAAGTGGTGATTCTCGGCCAG  >  NC_010473/2806340‑2806729

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