| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 3,639,879 | (T)8→7 | intergenic (‑159/‑232) | uspB ← / → uspA | universal stress protein B/universal stress protein A |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,639,879 | 0 | T | . | 85.7% | 86.4 / 13.4 | 28 | intergenic (‑159/‑232) | uspB/uspA | universal stress protein B/universal stress protein A |
| Reads supporting (aligned to +/- strand): ref base T (2/2); new base . (10/14); total (12/16) | |||||||||||
| 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.53e-01 | |||||||||||
AACGCGCGTAGTGATGAGAAATAGCGCGCCATGTTAACAATGCAAACGACACATAAAGCCCAAAATAATGCGACGGTGCTTATCATACCTCCTCCCCGGCGACCTGCCCGCGGAGTTCCACCCCGGGGCTACCGCTCCCGATACGCTGCCAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCGGGTCGCTATAAGGTAAGGATGGTCTTAACACTGAATCTTTACGGCTGGGTTAGCCCCGCGCACGTAGTTCGCAGGACGCGGGTGACGTAACGGCACAAGAAACGCTAGCTGGCCAGTCATCGACAACTTTATGGAAGGAGTAACACTATGGCTTA > NC_000913/3639635‑3640118 | aacGCGCGTAGTGATGAGAAATAGCGCGCCATGTTAACAATGCAAACGACACATAAAGCCCAAAATAATGCGACGGTGCTTATCATACCTCCTCCCCGGCGACCTGCCCGCGGAGTTCCACCCCGGGGCTACCGCTCCCGATACGCTGCCAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTT‑CCAGAAg < 2:80733/251‑1 (MQ=255)aacGCGCGTAGTGATGAGAAATAGCGCGCCATGTTAACAATGCAAACGACACATAAAGCCCAAAATAATGCGACGGTGCTTATCATACCTCCTCCCCGGCGACCTGCCCGCGGAGTTCCACCCCGGGGCTACCGCTCCCGATACGCTGCCAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTAGAGATAGTGTAGGAGATCAGGTTGTTTTTTT‑CCAGAAg < 1:166943/251‑1 (MQ=255) aaTAGCGCGCCATGTTAACAATGCAAACGACACATAAATCCCAAAATAATGCGACGGTGCTTAGCATACCTCCTCCCCGGCGACCTGCCCGCGGAGTTCCACCCCGGTGCTACCGCTCCCGATACGCTGCCAATCAGTTAACACCAGGTCTTGGAGAAACCGCGGTTGTGGTGACCAACATACGAGCGGCTCTAGAGATAGTGTAGGAGATCAGGT‑GGTTTTTTTCCAGAATTTTAACCACTATCAatata < 1:171797/251‑1 (MQ=255) aGCGCGCCATGTTAACAATGCAAACGACACATAAAGCCCAAAATAATGCGACGGTGCTTATCATACCTCCTCCCCGGCGACCTGCCCGCGGAGTTCCACCCCGGGGCTACCGCTCCCGATACGCTGCCAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTT‑CCAGAAGGTTAACCACTATCAATATATTc > 2:60759/1‑251 (MQ=255) cgcgcCATGTTAACAATGCAAACGACACATAAAGCCCAAAATAATGCGACGGTGCTTATCATACCTCCTCCCCGGCGACCTGCCCGCGGAGTTCCACCCCGGGGCTACCGCTCCCGATACGCTGCCAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTT‑CCAGAAGGTTAACCACTATCAATATATTCAt > 1:150528/1‑251 (MQ=255) aTGTTAACAATGCAAACGACACATAAAGCCCAAAATAATGCGACGGTGCTTATCATACCTCCTCCCCGGCGACCTGCCCGCGGAGTTCCACCCCGGGGCTACCGCTCCCGATACGCTGCCAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCg > 2:64756/1‑250 (MQ=255) aaTGCAAACGACACATAAAGCCCAAAATAATGCGACGGTGCTTATCATACCTCCTCCCCGGCGACCTGCCCGCGGAGTTCCACCCCGGGGCTACCGCTCCCGATACGCTGCCAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTT‑CCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGtt > 2:21469/1‑251 (MQ=255) gCAAACGACACATAAAGCCCAAAATAATGCGACGGTGCTTATCATACCTCCTCCCCGGCGACCTGCCCGCGGAGTTCCACCCCGGGGCTACCGCTCCCGATACGCTGCCAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGttt > 1:198587/1‑250 (MQ=255) gCAAACGACACATAAAGCCCAAAATAATGCGACGGTGCTTATCATACCTCCTCCCCGGCGACCTGCCCGCGGAGTTCCACCCCGGGGCTACCGCTCCCGATACGCTGCCAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTT‑CCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTAt > 2:263833/1‑251 (MQ=255) ccAAAATAATGCGACGGTGCTTATCATACCTCCTCCCCGGCGACCTGCCCGCGGAGTTCCACCCCGGGGCTACCGCTCCCGATACGCTGCCAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTT‑CCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGc > 2:115955/1‑251 (MQ=255) gCGACGGTGCTTATCATACCTCCTCCCCGGCGACCTGCCCGCGGAGTTCCACCCCGGGGCTACCGCTCCCGATACGCTGCCAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTT‑CCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACgg < 1:151762/251‑1 (MQ=255) cccGGCGACCTGCCCGCGGAGTTCCACCCCGGGGCTACCGCTCCCGATACGCTGCCAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTT‑CCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCGGGTCGCTAtaaggtaag < 1:34669/251‑1 (MQ=255) gACCTGCCCGCGGAGTTCCACCCCGGGGCTACCGCTCCCGATACGCTGCCAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTGTTT‑CCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGATTAAGCAAGGCGGATTGACTGATCATCCGGGTCGCTAGAAGGTAAGGATGGt < 1:175467/251‑1 (MQ=255) gCCCGCGGAGTTCCACCCCGGGGCTACCGCTCCCGATACGCTGCCAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTT‑CCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCTGGTCGCTATAAGGTAATGATGGTCtt < 1:51486/249‑1 (MQ=255) ccACCCCGGGGCTACCGCTCCCGATACGCTGCCAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTT‑CCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCGGGTCGCTATAAGGTAAGGATGGTCTTAACACTGAATCttt < 2:95624/251‑1 (MQ=255) aCCGCTCCCGATACGCTGCCAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTT‑CCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCGGGTCGCTATAAGGTAAGGATGGTCTTAACACTGAATCTTTACGGCTGGGtt < 1:226782/249‑1 (MQ=255) cAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGGTTTTTT‑CCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCGGGTCGCTATAAGGTAAGGATGGTCTTAACACTGAATCGTTACGGCGGGGTTAGCCCCGCGCACGTAGTTCGCAgg < 2:256559/250‑1 (MQ=255) acCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCGGGTCGCTATAAGGTAAGGATGGTCTTAACACTGAATCTTTACGGCTGGGTTAGCCCCGCGCACGTAGTTCGCAGGACGCGGGt < 1:64756/251‑1 (MQ=255) ggTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTT‑CCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCGGGTCGCTATAAGGTAAGGATGGTCTTAACACTGAATCTTTACGGCTGGGTTAGCCCCGCGCACGTAGTTCGCAGGACGCGGGTGACGt > 1:229890/1‑251 (MQ=255) cGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTT‑CCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCGGGTCGCTAtaaggtaag > 1:55437/1‑166 (MQ=255) gtggtgACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTT‑CCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCGGGTCTCTATAAGGTAAGGATGGTCTTAACGCTGAATCTTTACGGCTGGGTTAGCCCCGCGCTCGTAGTTCGCAGGACGCGGGTGACGTAACGGCACAAGAAACGCTAg < 1:42638/250‑1 (MQ=255) gACCAACATACGAGCGGCTCTATAGATAGTTTAGGAGATCAGGTTGTTTTTTT‑CCATAAGGTTAACCACTGTCAATATATTCATGTCGAAAATTTGGTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCGGGTCGCTAGAAGGTAAGGATGGTCTTAACACGGAATCTTTGCGGCTGGGTTAGCCCCGCGGACGGAGTGCGCAGGACGCGTGTGACGTAACGGCACAAGATACGCTAGCTGGc < 1:168436/250‑1 (MQ=255) gACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTT‑CCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCGGGTCGCTATAAGGTAAGGATGGTCTTAACACTGAATCTTTACGGCTGGGTTAGCCCCGCGCACGTAGTTCGCAGGACGCGGGTGACGTAACGGCACAAGAAACTCTAGCTGGcc < 2:102338/251‑1 (MQ=255) aTACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTT‑CCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCGGGTCGCTATAAGGTAAGGATGGTCTTAACACTGAATCTTTACGGCTGGGTTAGCCCCGCGCACGTAGTTCGCAGGACTCGGGTGACGTAACGGCACAAGAAACGCTAGCTGGCCAGTCATc < 1:136507/251‑1 (MQ=255) tACGAGCCTCTCTCTAGATAGTGTAGGAGATCCGGTTGTTTTTTT‑CCAGCAGGTTAACCACTATCAATACATTAATGTCGCAACTTTGTTTATATAACGAGTAAGCAAGGCGGATTGACGGATCATCCGGGTCGCTATAAGGTAAGGATGGTCTTAACACTGAATCTTTACGGCTGGGTTAGCCCCGCGCACGTAGTTCGCCGGACGCGGGTGACGTAAAGGCACCAGAAACGCTAGCTGGCCAGCCATCg > 2:277197/1‑251 (MQ=255) tagTGTAGGAGATCAGGTTGTTTTTTT‑CCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCGGGTCGCTATAAGGTAAGGATGGTCTTAACACTGAATCTTTACGGCTGGGTTAGCCCCGCGCACGTAGTTCGCAGGACGCGGGTGACGTAACGGCACCAGAAACGCTAGCTGGCCAGTCATCGACAACTTTATGGCAgg > 2:233821/1‑249 (MQ=255) ggAGATCAGGTTGTTTTTTT‑CCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCGGGTCGCTATAAGGTAAGGATGGTCTTAACACTGAATCTTTACGGCTGGGTTAGCCCCGCGCACGTAGTTCGCAGGACGCGGGTGACGTAACGGCACAAGAAACGCTAGCTGGCCAGTCATCGACAACTTTATGGAAGGAGTAACACt > 1:250513/1‑251 (MQ=255) ggTTGTTTTTTT‑CCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCGGGTCGCTATAAGGTAAGGATGGTCTTAACACTGAATCTTTACGGCTGGGTTAGCCCCGCGCACGTAGTTCGCAGGACGCGGGTGACGTAACGGCACAAGAAACGCTAGCTGGCCAGTCATCGACAACTTTATGGAAGGAGTAACACTATGGCTta < 1:21469/251‑1 (MQ=255) | AACGCGCGTAGTGATGAGAAATAGCGCGCCATGTTAACAATGCAAACGACACATAAAGCCCAAAATAATGCGACGGTGCTTATCATACCTCCTCCCCGGCGACCTGCCCGCGGAGTTCCACCCCGGGGCTACCGCTCCCGATACGCTGCCAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCGGGTCGCTATAAGGTAAGGATGGTCTTAACACTGAATCTTTACGGCTGGGTTAGCCCCGCGCACGTAGTTCGCAGGACGCGGGTGACGTAACGGCACAAGAAACGCTAGCTGGCCAGTCATCGACAACTTTATGGAAGGAGTAACACTATGGCTTA > NC_000913/3639635‑3640118 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |