| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 1,761,970 | A→T | 10.4% | N359I (AAC→ATC) ‡ | PP_1567 → | capsid protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 1,761,970 | 0 | A | T | 10.4% | 133.2 / 8.2 | 48 | N359I (AAC→ATC) ‡ | PP_1567 | capsid protein |
| Reads supporting (aligned to +/- strand): ref base A (25/18); new base T (3/2); total (28/20) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01 | |||||||||||
CTGGGCGGACATCGAGCTGACCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCAGAGTCGTTCGTTCACGGC > NC_002947/1761824‑1762118 | cTGGGCGGACATCGAGCTGACCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCGGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGAGAGTTCCTGGTAGGCGCGATGAACAt > 1:1044624/1‑150 (MQ=255) ggCGGACATCGAGCTGACCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATcgc > 1:1286534/1‑150 (MQ=255) gCGGACATCGAGCTGACCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAGGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGTTCATcgcg < 2:218219/150‑1 (MQ=255) gCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGaa > 2:917745/1‑150 (MQ=255) gCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGATGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGaa > 1:72148/1‑150 (MQ=255) ccAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCgg > 1:121077/1‑126 (MQ=255) ccAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCgg < 2:121077/126‑1 (MQ=255) cAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGc > 2:66370/1‑150 (MQ=255) ccTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCg > 2:1005787/1‑150 (MQ=255) cTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGa < 2:1286534/150‑1 (MQ=255) agcTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCGGGGGGTGGACAGCAACGCGGTACCAGAAGGCGAGGTCCTGGTAAGCGCGATGAACATCGCGGCCCAGGTGTTCGGCCGCGAGGAAGCGGCTGTCGAGGtc > 2:606494/3‑150 (MQ=255) tCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGtct > 2:500715/1‑150 (MQ=255) tGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAt > 1:1221881/1‑104 (MQ=255) tGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACAt < 2:1221881/104‑1 (MQ=255) aaTGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATcgcg < 2:965668/102‑1 (MQ=255) aaTGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATcgcg > 1:965668/1‑102 (MQ=255) aaTGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGATGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGaa < 2:72148/150‑1 (MQ=255) tCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGGCGCGGTCTCCACCGAAGACg > 2:698712/1‑150 (MQ=255) cAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAg < 1:527274/132‑1 (MQ=255) cAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAg > 2:527274/1‑132 (MQ=255) aGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTg > 2:180644/1‑150 (MQ=255) ggAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGa < 2:1102243/150‑1 (MQ=255) gcgTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTcc > 2:294515/1‑150 (MQ=255) aCAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCAc < 1:180644/150‑1 (MQ=255) gcgcCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCaa > 1:1224127/1‑150 (MQ=255) cTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATg < 2:1224127/150‑1 (MQ=255) tgGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTg > 2:70054/1‑150 (MQ=255) ggcaGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTACGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGa < 1:316092/147‑1 (MQ=255) aGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCa > 2:938392/1‑150 (MQ=255) tGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTc < 1:202808/77‑1 (MQ=255) tGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTc > 2:202808/1‑77 (MQ=255) tGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATcc > 2:1029742/1‑150 (MQ=255) aGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAgcgg > 2:965569/1‑149 (MQ=255) aGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTa < 1:282024/150‑1 (MQ=255) gAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTc < 2:634846/42‑1 (MQ=255) gAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTc > 1:634846/1‑42 (MQ=255) gAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGc > 2:1117449/1‑150 (MQ=255) gAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGCGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGc > 1:1045392/1‑150 (MQ=255) ttCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAACGGCTGTCGTGTTCTCTTCCGAAACCGCTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTAGACCGccag < 1:414370/150‑3 (MQ=255) ccTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTAAAGTGTTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCaga < 1:978817/150‑1 (MQ=255) gTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGa > 1:572816/1‑100 (MQ=255) gTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGa < 2:572816/100‑1 (MQ=255) ggCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCAGAGtcgttc < 1:500715/150‑1 (MQ=255) gccgcgATGTTCATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCAGAGtcgttc > 1:1209419/3‑150 (MQ=255) ccgcgATGTTCATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTc < 2:1025641/49‑1 (MQ=255) ccgcgATGTTCATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTc > 1:1025641/2‑50 (MQ=255) cgcgATGTTCATCGCGGCCCAGGTGTTCGACCGCATGAAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCAGAGtcgttcgt > 1:1066002/1‑150 (MQ=255) tGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCAGAGTCGTTCGTTCACg < 1:1005787/150‑1 (MQ=255) aaCATCGCGGCCCAGGTGTTCGAACGAGAGGAGGGAGTATACCCGTCCAGCTCCTAAGACGGTGACAACTTACGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCAGAGTCGTTCGTTCACGGc < 1:1249024/150‑1 (MQ=255) | CTGGGCGGACATCGAGCTGACCAAGGACGCCAACGGCTCCTACATCTGGGTGAATGTCCAGGAAGGCGGCGTACAGCGCCTGTGGAAGCTGCCGGTGGTGGACAGCAACGCGGTGCCAGAAGGCGAGTTCCTGGTAGGCGCGATGAACATCGCGGCCCAGGTGTTCGACCGCGAGGAAGCGGCTGTCGAGGTCTCCACCGAAGACGGTGACAACTTCCGCACCAACATGGTGACCATCCGCGCCGAGGAGCGCCTGGCGCTGGCGGTCTACCGCCCAGAGTCGTTCGTTCACGGC > NC_002947/1761824‑1762118 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |