| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 254,727 | A→C | G330G (GGT→GGG) | pepD ← | aminoacyl‑histidine dipeptidase (peptidase D) |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 254,727 | 0 | A | C | 96.4% | 191.4 / ‑0.7 | 57 | G330G (GGT→GGG) | pepD | aminoacyl‑histidine dipeptidase (peptidase D) |
| Reads supporting (aligned to +/- strand): ref base A (2/0); new base C (31/23); total (33/24) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.07e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01 | |||||||||||
CAGTGAATTTCTACGTTATTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCGAGATTTCGCAATCAGGGCAGCTTTATCGTTCGCTACAGAGTCCAGCAACAAG > NC_000913/254632‑254820 | cAGTGAATTTCTACGTTATTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACCCCGt > 2:112127/1‑100 (MQ=255) tGAATTTCTACGTTATTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACCCCGTTCg > 1:650677/1‑100 (MQ=255) aTTTCTACGTTATTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACCCCGTTCgggg > 1:384042/1‑100 (MQ=255) tttCTACGTTATTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACGCCTTTGGCTACATCGGAGTTACGAATCACTCCGTTCgggg < 1:444860/99‑1 (MQ=255) ttCTACGTTATTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACCCCGTTCGGGGTg > 2:63707/1‑100 (MQ=255) ttCTACGTTATTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACCCCGTTCGGGGTg > 2:521549/1‑100 (MQ=255) tCTACGTTATTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACCCCGTTCGGGGTgg < 2:818824/100‑1 (MQ=255) tCTACGTTATTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACCCCGTTCGGGGTgg > 2:250728/1‑100 (MQ=255) tACGTTATTGGCAGTCATGGTCACCACACCGGCGTTCATGGACGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACCCCGTTCGGGGTGGCg < 2:169723/100‑1 (MQ=255) cGTTATTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACCCCGTTCGGGGTgattcg < 1:957448/100‑6 (MQ=255) cGTTATTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACCCCGTTCGGGGTGGCGtt > 1:350140/1‑100 (MQ=255) gTTATTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACCCCGTTCGGGGTGGCGTTc < 1:933258/100‑1 (MQ=255) gTTATTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTc > 2:146725/1‑100 (MQ=255) gTTATTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTc > 2:146760/1‑100 (MQ=255) aTTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACCCCGTTCGGGGTGGCGTTcagc < 1:1073912/100‑1 (MQ=255) cagtcaTGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACCCCGTTCGGGGTGGCGTTCAGCAGACg > 2:1184361/1‑100 (MQ=255) cagtcaTGGTCACCACACCGACGTTCAGGGAGGTTTCAACAACACCTTTGGCTACATCGGAGTTACGAATCACCCCGTTCGGGGTGGCGTTCAGCAGACg > 2:242571/1‑100 (MQ=255) gTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAgg < 2:434549/100‑1 (MQ=255) gTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAgg < 1:508592/100‑1 (MQ=255) tCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGt > 1:914547/1‑100 (MQ=255) caccacACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTa > 2:1273782/1‑100 (MQ=255) accacACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTAt > 1:532531/1‑100 (MQ=255) ccacACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATc > 2:490027/1‑100 (MQ=255) cacCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATcgcg < 1:37141/100‑1 (MQ=255) ccGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCgag < 2:134856/100‑1 (MQ=255) gTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCGAGATTTc < 2:1137139/100‑1 (MQ=255) tCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCGAGATTTCGc < 2:415224/100‑1 (MQ=255) tCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCGAGATTTCGc < 2:650673/100‑1 (MQ=255) tCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCGAGATTTCGc > 1:274845/1‑100 (MQ=255) aGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCGAGATTTCGCaa > 1:718011/1‑100 (MQ=255) gggAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCGAGATTTCGCAAt > 2:664940/1‑100 (MQ=255) gggAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCGAGATTTCGCAAt > 1:842997/1‑100 (MQ=255) gggAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCGAGATTTCGCAAt > 1:1029453/1‑100 (MQ=255) gggAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCGAGATTTCGCAAt > 1:51478/1‑100 (MQ=255) ggTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCGAGATTTCGCAATCAgg < 1:1188211/100‑1 (MQ=255) gTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCGAGATTTCGCAATCAggg > 2:237608/1‑100 (MQ=255) cacaCCTTTGGCTACATCGGAGTTACGAATCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCGAGATTTCGCAATCAGGGCAGAttt > 1:55026/1‑99 (MQ=255) acCTTTGGCTACATCGGAGTTACGAATCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCGAGATTTCGCAATCAGGGCAGCTTTATCg < 2:51477/100‑1 (MQ=255) acCTTTGGCTACATCGGAGTTACGAATCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCGAGATTTCGCAATCAGGGCAGCTTTATCg < 1:356238/100‑1 (MQ=255) ccTTTGGCTACATCGGAGTTACGAATCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCGAGATTTCGCAATCAGGGCAGCTTTATCGt > 1:548888/1‑100 (MQ=255) ggCTACATCGGAGTTACGAATCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCGAGATTTCGCAATCAGGGCAGCTTTATCGTTCGCt > 1:232660/1‑100 (MQ=255) tACATCGGAGTTACGAATCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCGAGATTTCGCAATCAGGGCAGCTTTATCGTTCGCTACa < 2:51395/100‑1 (MQ=255) aCATCGGAGTTACGAATCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCGATATTTCGCAATCAGGGCAGCTTTATCGTTCGCTACag < 1:920150/100‑1 (MQ=255) aCATCGGAGTTACGAATCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCGAGATTTCGCAATCAGGGCAGCTTTATCGTTCGCTACag > 1:1213612/1‑100 (MQ=255) cATCGGAGTTACGAATCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCGAGATTTCGCAATCAGGGCAGCTTTATCGTTCGCTACaga < 1:431365/100‑1 (MQ=255) gAGTTACGAATCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCGAGATTTCGCAATCAGGGCAGCTTTATCGTTCGCTACAGAGTCCa > 1:723462/1‑100 (MQ=255) gTTACGAATCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCGAGATTTCGCAATCAGGGCAGCTTTATCGTTCGCTACAGAGTCCAGc < 1:433847/100‑1 (MQ=255) gTTACGAATCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCGAGATTTCGCAATCAGGGCAGCTTTATCGTTCGCTACAGAGTCCAGc < 1:521553/100‑1 (MQ=255) gTTACGAATCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCGAGATTTCGCAATCAGGGCAGCTTTATCGTTCGCTACAGAGTCCAGc < 1:193491/100‑1 (MQ=255) ttACGAATCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCGAGATTTCGCAATCAGGGCAGCTTTATCGTTCGCTACAGAGTCCAGca < 1:1041032/100‑1 (MQ=255) tACGAATCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCGAGATTTCGCAATCAGGGCAGCTTTATCGTTCGCTACAGAGTCCAGcaa > 1:966147/1‑100 (MQ=255) ctaaTCACCCCGTTCGGGGTGGCGTTCAGAAGACGAATAAAGGTATCGCGAGATTTCGCAATCAGGGCAGCTTTATCGTTCGCTACAGAGTCCAGcaaca < 2:512653/98‑1 (MQ=255) ggaTCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCGAGATTTCGCAATCAGGGCAGCTTTATCGTTCGCTACAGAGTCCAGcaacaa > 1:880008/3‑100 (MQ=255) gAATCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCGAGATTTCGCAATCAGGGCAGCTTTATCGTTCGCTACAGAGTCCAGcaacaa > 1:1015956/1‑100 (MQ=255) aaTCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCGAGATTTCGCAATCAGGGCAGCTTTATCGTTCGCTACAGAGTCCAGCAACAAg > 1:26066/1‑100 (MQ=255) aaTCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCGAGATTTCGCAATCAGGGCAGCTTTATCGTTCGCTACAGAGTCCAGCAACAAg > 1:562138/1‑100 (MQ=255) aaTCACCCCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCGAGATTTCGCAATCAGGGCAGCTTTATCGTTCGCTACAGAGTCCAGCAACAAg > 1:156319/1‑100 (MQ=255) | CAGTGAATTTCTACGTTATTGTCAGTCATGGTCACCACACCGACGTTCAGGGAGGTTTCAACCACACCTTTGGCTACATCGGAGTTACGAATCACACCGTTCGGGGTGGCGTTCAGCAGACGAATAAAGGTATCGCGAGATTTCGCAATCAGGGCAGCTTTATCGTTCGCTACAGAGTCCAGCAACAAG > NC_000913/254632‑254820 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |