| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 2,879,257 | C→T | K172K (AAG→AAA) | cas1 ← | multifunctional nuclease Cas1 |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,879,257 | 0 | C | T | 76.5% | 57.7 / 16.9 | 34 | K172K (AAG→AAA) | cas1 | multifunctional nuclease Cas1 |
| Reads supporting (aligned to +/- strand): ref base C (4/4); new base T (16/10); total (20/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.89e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.36e-01 | |||||||||||
AGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACGCCGTATTGCTTCGCCAGAAGTGCG > NC_000913/2879179‑2879326 | aGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCttt > 2:497176/1‑81 (MQ=255) tGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCTTTTTc > 1:497183/1‑81 (MQ=255) aaGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCTTTTTCCCa < 1:497149/81‑1 (MQ=255) gccgcTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCTTTTTCCCAGTCTTTc > 2:497181/1‑81 (MQ=255) cTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCTTTTTCCCAGTCTTTTGGAt > 2:497177/1‑81 (MQ=255) cTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGAt < 1:497150/81‑1 (MQ=255) ttCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCTTTTTCCCAGTCTTTCGGATc > 2:497193/1‑81 (MQ=255) cAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCg > 2:497178/1‑80 (MQ=255) aGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCTTTTTCCCAGTCTTTCGGATCg > 1:497189/1‑79 (MQ=255) aGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTa < 2:497151/81‑1 (MQ=255) gTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCTTTTTCCCAGTCTTTCGGATCGTAg < 2:497152/81‑1 (MQ=255) ttACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCTTTTTCCCAGTCTTTCGGATCGTAGc < 2:497153/81‑1 (MQ=255) aCGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCTTTTTCCCAGTCTTTCGGATCGTAGCg < 1:497155/80‑1 (MQ=255) aCGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCTTTTTCCCAGTCTTTCGGATCGTAGCg < 1:497154/80‑1 (MQ=255) cGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGAc > 1:497186/1‑81 (MQ=255) gCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCTTTTTCCCAGTCTTTCGGATCGTAGCGACg > 1:497199/1‑81 (MQ=255) cGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCTTTTTCCCAGTCTTTCGGATCGTAGCGACGtc > 2:497180/1‑81 (MQ=255) tataAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCTTTTTCCCAGTCTTTCGGATCGTAGCGACGtcca > 1:497182/1‑81 (MQ=255) aaaCAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCTTTTTCCCAGTCTTTCGGATCGTAGCGACGtccattc > 2:497179/1‑81 (MQ=255) aaaCAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCTTTTTCCCAGTCTTTCGGATCGTAGCGACGtccattc > 2:497202/1‑81 (MQ=255) aCAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGACGtccattcca > 1:497184/1‑81 (MQ=255) gcgcTAATGCATTGGTTGATCGTATCGCCTTTTTCCCAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACGCCGTAt < 1:497156/81‑1 (MQ=255) gcTAATGCATTGGTTGATCGTATCGCCTTTTTCCCAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACGCCGTATTg < 1:497157/81‑1 (MQ=255) cTAATGCATTGGTTGATCGTATCGCCTTTTTCCCAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACGCCGTATTGc < 1:497158/81‑1 (MQ=255) aaTGCATTGGTTGATCGTATCGCCTTTTTCCCAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACGCCGTATTGCtt > 2:497196/1‑81 (MQ=255) aTGCATTGGTTGATCGTATCGCCTTTTTCCCAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACGCCGTATTGCTTc > 1:497187/1‑81 (MQ=255) cATTGGTTGATCGTATCGCCTTTTTCCCAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACGCCGTATTGCTTCGc > 2:497190/1‑80 (MQ=255) ggTTGATCGTATCGCCTTTTTCCCAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACGCCGTATTGCTTCGCCAGa < 1:497160/80‑1 (MQ=255) ggTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACGCCGTATTGCTTCGCCAGa < 1:497159/80‑1 (MQ=255) ttGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACGCCGTATTGCTTCGCCAGAAGt < 2:497161/81‑1 (MQ=255) gATCGTATCGCCTTTTTCCCAGTCTTTCGGCTCGTAGCGACGTCCATTCCATGTCACGCCGTATTGCTTCGCCAGAAGTGc > 1:497192/1‑81 (MQ=255) gATCGTATCGCCTTTTTCCCAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACGCCGTATTGCTTCGCCAGAAGTg < 1:497162/80‑1 (MQ=255) gATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACGCCGTATTGCTTCGCCAGAAGTGc > 2:497195/1‑81 (MQ=255) cGTATCGCCTTTTTCCCAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACGCCGTATTGCTTCGCCAGAAGTGCg > 1:497188/1‑79 (MQ=255) cGTATCGCCTTTTTCCCAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACGCCGTATTGCTTCGCCAGAAGTGCg > 2:497198/1‑79 (MQ=255) | AGCTGCAAGTATCGCCGCTTCAGTTACGCCGTATAAACAGGAAGTTGCAGCGCTAATGCATTGGTTGATCGTATCGCCCTTTTCCCAGTCTTTCGGATCGTAGCGACGTCCATTCCATGTCACGCCGTATTGCTTCGCCAGAAGTGCG > NC_000913/2879179‑2879326 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |