| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 709,862 | C→G | 18.4% | R202G (CGC→GGC) | ribF → | bifunctional riboflavin kinase/FMN adenylyltransferase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 709,862 | 0 | C | G | 18.4% | 58.9 / 10.2 | 28 | R202G (CGC→GGC) | ribF | bifunctional riboflavin kinase/FMN adenylyltransferase |
| Reads supporting (aligned to +/- strand): ref base C (13/9); new base G (3/2); total (16/12) | |||||||||||
| 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.64e-01 | |||||||||||
GCCAACACGGTCATCCAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTGCCGGTGATGGG > NC_002947/709715‑710011 | gccAACACGGTCATCCAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGc > 2:171133/1‑150 (MQ=255) ccTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCaa > 2:116805/1‑150 (MQ=255) cagcaCCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCaa < 1:116805/150‑1 (MQ=255) cagcaCCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCaa < 2:750316/150‑1 (MQ=255) gcaCCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAgc < 1:670979/150‑1 (MQ=255) cAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGcc < 2:299342/150‑1 (MQ=255) aaGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGGAGCTGGCGGGCCAGCttc < 1:485465/106‑3 (MQ=255) aaGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGGAGCTGGCGGGCCAGCttc > 2:485465/1‑104 (MQ=255) gACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTg > 2:751325/1‑150 (MQ=255) aCTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTgg < 1:259479/150‑1 (MQ=255) cTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATc < 2:707191/150‑1 (MQ=255) ctgctgGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCCACCGCCAACATCCAGCTCAAGAGCCGCCGCGTGCCGCTCAACGGTGTATACCTGGCCAGCATCGTGATCGACGGc > 2:384001/1‑150 (MQ=255) ctgctgGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGc > 1:30570/1‑150 (MQ=255) tgGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAg > 1:69748/1‑150 (MQ=255) gccCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATc > 1:156607/1‑134 (MQ=255) gccCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATc < 2:156607/134‑1 (MQ=255) gCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCATCGAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCg < 1:384001/150‑1 (MQ=255) tgcACGGCCAGAAGCTGGCCCGCCAGCTTGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCggg < 1:55804/83‑1 (MQ=255) tgcACGGCCAGAAGCTGGCCCGCCAGCTTGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCggg > 2:55804/1‑83 (MQ=255) gcACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGccc > 2:560832/1‑150 (MQ=255) aCGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACAGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGAGCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCa > 1:609989/1‑150 (MQ=255) cgaGCTGGCGGGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTGc > 2:53643/3‑150 (MQ=255) gaGCTGGCGGGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCg > 2:509075/2‑70 (MQ=255) gaGCTGGCGGGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCg < 1:509075/69‑1 (MQ=255) aaGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTg > 1:577550/1‑93 (MQ=255) aaGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTg < 2:577550/93‑1 (MQ=255) cTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCTCACCGTTGCCGGt > 2:128569/1‑150 (MQ=255) cTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTGCCGGt > 2:183213/1‑150 (MQ=255) cGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTGCCGGTGATggg > 1:481271/1‑150 (MQ=255) | GCCAACACGGTCATCCAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTGCCGGTGATGGG > NC_002947/709715‑710011 |
| 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: ‑ |