| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 1,825,134 | G→C | 8.8% | A54A (GCC→GCG) | mutS ← | DNA mismatch repair protein MutS |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 1,825,134 | 0 | G | C | 8.8% | 130.3 / 2.7 | 45 | A54A (GCC→GCG) | mutS | DNA mismatch repair protein MutS |
| Reads supporting (aligned to +/- strand): ref base G (21/20); new base C (2/2); total (23/22) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01 | |||||||||||
CACCGGGCCCTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAG‑‑AAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCGGGGTGTGT > NC_002947/1824987‑1825278 | cACCGGGCCCTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCgcgggc > 2:1072511/1‑150 (MQ=255) aCCGGGCCCTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTTACCgcgggcg > 1:725697/1‑150 (MQ=255) ccTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCCCGCGGTCAgtc < 1:319892/150‑3 (MQ=255) ttGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCC‑‑CTGACCGCCCGCGGTCAgtcggcc < 1:162715/150‑7 (MQ=255) tGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGa < 2:1017116/150‑1 (MQ=255) tGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGAATGACCGCGGGCGGTCAGGGTGa < 2:72518/150‑1 (MQ=255) tggtggCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGAGCATGTGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGatat > 1:618326/1‑150 (MQ=255) ggCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCcag < 1:510483/150‑1 (MQ=255) cTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCttttt > 1:415121/1‑150 (MQ=255) tCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTc < 2:585405/150‑1 (MQ=255) aTCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCcag > 2:866032/1‑125 (MQ=255) aTCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGAACGCGGGCGGTCAGGGTGATATCcag < 1:866032/125‑1 (MQ=255) caccGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTa < 2:725697/150‑1 (MQ=255) accGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTag > 2:418688/1‑150 (MQ=255) aCTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTag‑‑aaga > 1:590513/1‑150 (MQ=255) gAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAG‑‑AAGATCTCGTa > 2:324717/1‑150 (MQ=255) cTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCAGCTTCGTAG‑‑AAGATCTCGTAGaa > 2:559137/1‑150 (MQ=255) ccAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAG‑‑AAGATCTCGTAGAAGTCGc > 1:678885/1‑150 (MQ=255) ccAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAACCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAG‑‑AAGATCTCGTAGAAGTCGc < 2:618326/150‑1 (MQ=255) cAGCTTGGCCAGGTAGCCTTCCAGCGAATGGCACGGAAGCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAG‑‑AAGATCTCGTAGAAGTCGcc > 1:213523/1‑150 (MQ=255) ggCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAG‑‑AAGATCTCGTAGAAGTCGCCCATGCg < 2:415121/150‑1 (MQ=255) gCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAG‑‑AAGATCTCGTAGAAGTCGCCCATGCgg > 1:728566/1‑150 (MQ=255) gTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAG‑‑AAGATCTCGTAGAAGTCGCCCATGCGGTAGaa > 1:760574/1‑150 (MQ=255) gTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAG‑‑AAGATCTCGTAGAAGTCGCCCATGCGGTAGaa < 2:590513/150‑1 (MQ=255) gTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAG‑‑AAGATCTCGTAGAAGTCGCCCATGCGGTAGaa < 1:1065980/150‑1 (MQ=255) cTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGCCTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAG‑‑AAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCa < 2:678885/150‑1 (MQ=255) tCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAG‑‑AAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGc < 2:728566/150‑1 (MQ=255) aCGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAG‑‑AAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGgcgcc > 1:381376/1‑149 (MQ=255) cGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAG‑‑AAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTg < 2:760574/150‑1 (MQ=255) ggAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAG‑‑AAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTgg < 1:324717/150‑1 (MQ=255) aaTCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAG‑‑AAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGtt > 2:251838/1‑150 (MQ=255) tCCCGCACATGGGGATAGCCTGGCCAGCCTACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATATTCGTAGGCGAGATCTTATTTTACTCG‑CCGT‑CGGGAGTACATCAGCTGGTCCGGGTGCTGGcttt > 2:242850/1‑150 (MQ=255) ccGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAG‑‑AAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTc < 1:955931/150‑1 (MQ=255) ccGCACATGGGGATGGACTGGAAGGCCGACTGACCGCGGGCGGTCAGGGTAATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTAGTAG‑‑AAGATCTCGTAAAAGTCGACCATGACGTAGCACATAAGCTGGTCCGGGTGCTGGttttta > 2:785291/1‑149 (MQ=255) cGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAG‑‑AAGATCTCg > 1:734615/1‑98 (MQ=255) cGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAG‑‑AAGATCTCg < 2:734615/98‑1 (MQ=255) aTGGGGATGGACTGGCCGGCCGACTGACCGCGGGAGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAG‑‑AAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTc > 1:214052/1‑150 (MQ=255) ccTGACCGCCCGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAG‑‑AAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCgggg > 2:256871/2‑150 (MQ=255) ccTGACCGCCCGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAG‑‑AAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCgggg > 2:178448/2‑150 (MQ=255) aCTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAG‑‑AAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCgggg < 1:200523/150‑1 (MQ=255) cTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAG‑‑AAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCGGGgt > 1:231547/1‑150 (MQ=255) gACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAG‑‑AAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCGGGgtgt < 1:418688/150‑1 (MQ=255) gACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATATTCGTAG‑‑AAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCGGGgtgt > 1:923565/1‑150 (MQ=255) aCCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAG‑‑AAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCGGGgtgtg > 1:480535/1‑150 (MQ=255) ccGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAG‑‑AAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCGGGgtgtgt < 1:695561/150‑1 (MQ=255) cgcgGGCGGTCAGGGTGATATCCGGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAG‑‑AAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTg < 1:975988/109‑1 (MQ=255) cgcgGGCGGTCAGGGTGATATCCGGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAG‑‑AAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTg > 2:975988/1‑109 (MQ=255) | CACCGGGCCCTTGCTGGTGGCGGGGTCGCCGATCTGCTCACAGATCACCACCGACTCGCCGAGCTTGACCAGCTTGGCCAGGTAGCCTTCCAGCGAATGGAACGGAATCCCGCACATGGGGATGGACTGGCCGGCCGACTGACCGCGGGCGGTCAGGGTGATATCCAGCAGTTTCGCGGCCTTTTTCGCATCTTCGTAG‑‑AAGATCTCGTAGAAGTCGCCCATGCGGTAGAACATCAGCTGGTCCGGGTGCTGGTTTTTCAGCTTCCAGTACTGCTGCATCATCGGGGTGTGT > NC_002947/1824987‑1825278 |
| 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: ‑ |