Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 3,012,599 | T→C | 7.5% | intergenic (+206/+15) | yqeA → / ← yqeB | putative amino acid kinase/XdhC‑CoxI family protein with NAD(P)‑binding Rossman fold |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,012,599 | 0 | T | C | 7.5% | 172.2 / 3.6 | 53 | intergenic (+206/+15) | yqeA/yqeB | putative amino acid kinase/XdhC‑CoxI family protein with NAD(P)‑binding Rossman fold |
Reads supporting (aligned to +/- strand): ref base T (25/24); new base C (2/2); total (27/26) | |||||||||||
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.88e-01 |
ATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGCCCGGGCTTTATCAGAAACGCTGGTGAAA > NC_000913/3012508‑3012728 | aTAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATc < 4:2214/93‑1 (MQ=255) aTAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATc > 3:2214/1‑93 (MQ=255) tGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCActt > 1:198625/1‑120 (MQ=255) tGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCActt < 2:198625/120‑1 (MQ=255) aTGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCActtct > 4:266877/1‑118 (MQ=255) aTGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCActtct < 3:266877/118‑1 (MQ=255) tGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAAca > 3:986194/1‑148 (MQ=255) aCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAAcatca > 2:871039/1‑148 (MQ=255) tGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTa < 1:871039/148‑1 (MQ=255) ccTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCAt > 4:282293/1‑83 (MQ=255) ccTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCAt < 3:282293/83‑1 (MQ=255) ccTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCAt > 4:171527/1‑83 (MQ=255) ccTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCAt < 3:171527/83‑1 (MQ=255) tGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCAt > 4:1039998/1‑67 (MQ=255) tGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCAt < 3:1039998/67‑1 (MQ=255) aTAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCg > 2:67985/1‑108 (MQ=255) aTAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCg < 1:67985/108‑1 (MQ=255) aTAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACg > 3:1078550/1‑116 (MQ=255) aTAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACg < 4:1078550/116‑1 (MQ=255) aaGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTg > 2:1127931/1‑103 (MQ=255) aaGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTg < 1:1127931/103‑1 (MQ=255) cgcATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTa > 3:1114656/1‑72 (MQ=255) cgcATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTa < 4:1114656/72‑1 (MQ=255) aTCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCAt > 3:1122983/1‑114 (MQ=255) aTCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCAt < 4:1122983/114‑1 (MQ=255) tCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTAcgccg < 2:633797/141‑1 (MQ=255) tCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTAcgccg > 1:633797/1‑141 (MQ=255) gCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGccc > 2:873087/1‑101 (MQ=255) gCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGccc < 1:873087/101‑1 (MQ=255) cAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTg < 4:418249/102‑1 (MQ=255) cAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTg > 3:418249/1‑102 (MQ=255) cAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGcc > 3:1121917/1‑148 (MQ=255) cAATTTGGCGTTTGCCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGcc > 2:227486/1‑148 (MQ=255) ttGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGCCCGGGCTTTATc < 4:1121917/148‑1 (MQ=255) tcaGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTAcg < 1:123281/116‑1 (MQ=255) tcaGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTAcg > 2:123281/1‑116 (MQ=255) gTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGCCCGGGCTTTATCAGAAACGCt < 1:227486/148‑1 (MQ=255) cAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGCCCGGGCTTTATCAGAAACGCt < 4:898652/144‑1 (MQ=255) cAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGCCCGGGCTTTATCAGAAACGCt > 3:898652/1‑144 (MQ=255) cAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGCCCGGGCTTTATCAGAAACGCt > 3:356022/1‑144 (MQ=255) cAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGCCCGGGCTTTATCAGAAACGCt < 4:356022/144‑1 (MQ=255) ggCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGCCCGGGCTTTATCAGAAACGCTGGTGaa > 1:502820/1‑148 (MQ=255) gCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGCCCGGGCTTTATCAGAAACGCTGGTGaaa < 2:502820/148‑1 (MQ=255) gggtgggtAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGCCCGGGCTTTATCAGAAACg < 3:1053941/136‑1 (MQ=255) gggtgggtAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGCCCGGGCTTTATCAGAAACg > 4:1053941/1‑136 (MQ=255) ttacccACCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGCCCggg < 2:290594/110‑1 (MQ=255) ttacccACCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGCCCggg > 1:290594/5‑114 (MQ=255) ttacccACCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGCCCggg > 1:310125/5‑114 (MQ=255) ttacccACCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGCCCggg < 2:310125/110‑1 (MQ=255) aaTACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTa > 4:682893/1‑93 (MQ=255) aaTACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTa < 3:682893/93‑1 (MQ=255) tACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGCCCGGGCTTTATCAGAAACGCTGGTGaa < 3:603215/134‑1 (MQ=255) tACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGCCCGGGCTTTATCAGAAACGCTGGTGaa > 4:603215/1‑134 (MQ=255) ccATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTg > 3:1032041/1‑60 (MQ=255) ccATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTg < 4:1032041/60‑1 (MQ=255) | ATAGATTGCAATGCACTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTGGCGTTTGTCATCAGTCTCAGGCCGGGTGGGTAATACCATCCGGCCATTTTCATTTAAGCCACTTCCAGCACTTCTTTTGTCGCTTTCACGCCCTGATGCATCAACATCATTAACGCCTCAAGTACGCCGCCGCCAATCGCCCGGGCTTTATCAGAAACGCTGGTGAAA > NC_000913/3012508‑3012728 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |