| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | CP009273 | 1,484,009 | T→C | 33.3% | pseudogene (213/252 nt) | gapC ← | pseudogene, GAP dehydrogenase; glyceraldehyde‑3‑phosphate dehydrogenase (second fragment) |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | CP009273 | 1,484,009 | 0 | T | C | 33.3% | 20.0 / 15.5 | 21 | pseudogene (213/252 nt) | gapC | pseudogene, GAP dehydrogenase; glyceraldehyde‑3‑phosphate dehydrogenase (second fragment) |
| Reads supporting (aligned to +/- strand): ref base T (8/6); new base C (3/4); total (11/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.59e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GCAGATGGTAAACATGGCTTGCATGAATATCTGCAGACCCAGGTGGTTTATTTACAGTCTTAATGAGTGAAAGAGGCGGAGGTTTTTTCCTCCGCCTGTGCGCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTT > CP009273/1483866‑1484148 | gcAGATGGTAAACATGGCTTGCATGAATATCTGCAGACCCAGGTGGTTTATTTACAGTCTTAATGAGTGAAAGAGGCGGAGGTTTTTTCCTCCGCCTGTGCGCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGTGTGACg < 2:124078/150‑1 (MQ=255) aCATGGCTTGCATGAATATCTGCAGACCCAGGTGGTTTATTTACAGTCTTAATGAGTGAAAGAGGCGGAGGTTTTTTCCTCCGCCTGTGCGCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGTGTGACGAAGCCATATTc > 2:147016/1‑150 (MQ=255) tGCATGAATATCTGCAGACCCAGGTGGTTTATTTACAGTCTTAATGAGTGAAAGAGGCGGAGGTTTTTTCCTCCGCCTGTGCGCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCg > 1:119270/1‑149 (MQ=255) atatCTGCAGACCCAGGTGGTTTATTTACAGTCTTAATGAGTGAAAGAGGCGGAGGTTTTTTCCTCCGCCTGTGCGCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGc < 2:10408/150‑1 (MQ=255) atCTGCAGACCCAGGTGGTTTATTTACAGTCTTAATGAGTGAAAGAGGCGGAGGTTTTTTCCTCCGCCTGTGCGCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGa > 2:157180/1‑150 (MQ=255) cttcattttatcaggtaaaaaaaagagcgaccgaagtcgctctttttagatcacagtgtcaCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGCGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAg < 1:124774‑M1/89‑1 (MQ=255) cattttatcaggtaaaaaaaagagcgaccgaagtcgctctttttagatcacagtgtcaCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGCGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGtt > 2:144073‑M1/59‑149 (MQ=255) cattttatcaggtaaaaaaaagagcgaccgaagtcgctctttttagatcacagtgtcaCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGCGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGt > 1:195267‑M1/59‑148 (MQ=255) gCGGAGGTTTTTTCCTCCGCCTGTGCGCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGc > 1:105828/1‑150 (MQ=255) ggTTTTTTCCTCCGCCTGTGCGCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTgg > 1:193960/1‑149 (MQ=255) gTTTTTTCCTCCGCCTGTGCGCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCa > 1:196072/1‑150 (MQ=255) gatcacagtgtcaCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGCGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAacac < 1:61672‑M1/137‑1 (MQ=255) gagTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCTGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGc > 2:143664/1‑149 (MQ=255) gCGAATTTTTCGAGGGTGCGAATAAGCTGCGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGata < 1:124975/150‑1 (MQ=255) gAATTTTTCGAGGGTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATc < 1:108452/150‑1 (MQ=255) gAATTTTTCGAGGGTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATc < 1:111091/150‑1 (MQ=255) gAGGGTGCGAATAAGCTGCGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAgaga < 2:61672/149‑1 (MQ=255) aGGGTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACt < 2:196072/150‑1 (MQ=255) gggTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTa < 1:147016/150‑1 (MQ=255) gTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTAtt > 2:7022/1‑150 (MQ=255) gAATAAGCTGCGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATttctt > 1:67507/1‑150 (MQ=255) | GCAGATGGTAAACATGGCTTGCATGAATATCTGCAGACCCAGGTGGTTTATTTACAGTCTTAATGAGTGAAAGAGGCGGAGGTTTTTTCCTCCGCCTGTGCGCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTT > CP009273/1483866‑1484148 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |