Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP009273 | 1,484,072 | G→A | 39.1% | pseudogene (150/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,072 | 0 | G | A | 39.1% | 9.3 / 23.2 | 23 | pseudogene (150/252 nt) | gapC | pseudogene, GAP dehydrogenase; glyceraldehyde‑3‑phosphate dehydrogenase (second fragment) |
Reads supporting (aligned to +/- strand): ref base G (9/5); new base A (3/6); total (12/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.14e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GCGGAGGTTTTTTCCTCCGCCTGTGCGCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGG > CP009273/1483941‑1484214 | 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) cgtgACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGtata < 1:164270/149‑1 (MQ=255) gACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAAc > 2:91425/1‑149 (MQ=255) cGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCaa < 2:124774/150‑1 (MQ=255) cATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACt > 1:130723/1‑150 (MQ=255) gTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTAtt < 2:67507/150‑1 (MQ=255) gTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTggtggt > 1:87715/1‑150 (MQ=255) tttCACCAGTTGTAAATCGCCCACAGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGc > 1:76112/1‑150 (MQ=255) cACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGtt > 2:144453/1‑150 (MQ=255) tCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCg < 1:7022/150‑1 (MQ=255) cccACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCgg > 2:178992/1‑148 (MQ=255) | GCGGAGGTTTTTTCCTCCGCCTGTGCGCGTCAGAGTTTAGCGAATTTTTCGAGGGTGCGAATAAGCTGTGTGACGAAGCCATATTCGTTATCGTACCAGGCGACCGTTTTCACCAGTTGTAAATCGCCCACGGCGGTAATTTCCGTTTGCGTGGCATCAAACACCGAACCGAAATGGCTGCCAATGATATCGGAAGAGACTATTTCTTCATCGGTATAACCAAATGACTCGTTATTGGTGGTTGCTTGTTTAAGTGCGTTATTCACCTCTTCGG > CP009273/1483941‑1484214 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |