Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,879,536 | 0 | C | G | 44.0% | 10.5 / 22.5 | 25 | G79G (GGG→GGC) | ygbT | multifunctional endonuclease Cas1, CRISPR adaptation protein; DNA repair enzyme |
Reads supporting (aligned to +/- strand): ref base C (5/9); new base G (11/0); total (16/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.01e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACAGCTGCATG > NC_000913/2879469‑2879589 | gATAGAGCAGCTTATCTGAACGCGCACCTTCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTcc < 2:248775/68‑1 (MQ=255) gATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTcc < 2:561900/68‑1 (MQ=255) gATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTcc < 2:936889/68‑1 (MQ=255) gagCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTCCccca < 1:591324/68‑1 (MQ=255) agcTTATCTGAACGCGCACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTCCcccaccca < 2:720877/68‑1 (MQ=255) agcTTATCTGAACGCGCACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTCCcccaccca < 2:46865/68‑1 (MQ=255) cTTATCTGAACGCTCACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTCCCCCACCCATa < 1:515846/68‑1 (MQ=255) ttATCTGAACGCGCACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTCGCCCACCCATAc > 1:475117/1‑68 (MQ=255) ttATCTGAACGCGCACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTCGCCCACCCATAc > 2:970549/1‑68 (MQ=255) tGAACGCGCTCCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTCGCCCACACATACcaaca > 1:664661/1‑68 (MQ=255) tGAACGCGCACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTCCCCCACCCATACcaaca > 2:477690/1‑68 (MQ=255) aCGCGCACCTCCAGGCTGCCCAGAAGCATAAACACGAACGCACGCTTCGCCCACACATACCAACAATg > 1:539179/1‑68 (MQ=255) cgcgCACCTCCAGGCTGACCAGAAGCATAAACACGCACGCCCGCTTCCCCCACCCATACCAACAATGt > 2:817256/1‑68 (MQ=255) cACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTCGCCCACCCATACCAACAATGTTCCa > 1:335611/1‑68 (MQ=255) gCTGACCAGAATCATAAACACGAACGCCCGCTTACGCCAACAAGACCAACAATGTTCCAACTTGCGCa > 2:825933/1‑68 (MQ=255) gCTGACCAGAAGCATAAACACGAACGCCCGCTTGGCCCACCCATACCAACAATGTTCCAACTTgcgcg > 1:387979/1‑67 (MQ=255) gCTGACCAGAAGCATAAACACGAACGCCCGCTTCCCCCACACATACCAACAATGTTCCAACTTGCGCa > 1:224114/1‑68 (MQ=255) cTGACCAGAAGCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGATCCAACATGCGCAg > 2:8650/1‑68 (MQ=255) ccAGAAGCATAAACACGAACGCCCGCTTCGCCCACCCATACCAACAATGTTCCACCTTGCGCCGCCAg > 1:255716/1‑68 (MQ=255) aGAAGCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGc < 1:708888/68‑1 (MQ=255) aGAAGCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGc < 1:707882/68‑1 (MQ=255) aaGCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGt < 1:585362/68‑1 (MQ=255) aGCATAAACACGAACGCCCGCTTCGCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCCGGc > 1:829752/1‑65 (MQ=255) aGCATAAACACGAACGCCCGCTTCGCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTa > 2:572850/1‑68 (MQ=255) aGCATAAACACGAACGCCCGCTTCGCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTa > 1:392970/1‑68 (MQ=255) aGCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGc < 2:829749/65‑1 (MQ=255) cacGAACGCCCGCTTCGCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACAGCTGCa > 2:98484/1‑68 (MQ=255) cGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACAGCTGCATg < 2:167144/68‑1 (MQ=255) | GATAGAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACAGCTGCATG > NC_000913/2879469‑2879589 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |