| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,879,536 | 0 | C | G | 39.1% | 32.3 / 10.1 | 25 | G79G (GGG→GGC) | ygbT | multifunctional endonuclease Cas1, CRISPR adaptation protein; DNA repair enzyme |
| Reads supporting (aligned to +/- strand): ref base C (3/11); new base G (9/0); total (14/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.37e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACAGCTGCATGCGAAACC > NC_000913/2879473‑2879596 | gagCAGCTAATCTGAACGCGCACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTCCccca < 1:779138/68‑1 (MQ=255) gcagcTTATCTGAACGCGCACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTCCcccacc < 2:409294/68‑1 (MQ=255) agcTTATCTGAACGCGCACCTCCAGGCTGACCAGAAGCATAAACACGCACGCCCGCTTCGccccccca > 2:757920/1‑68 (MQ=255) agcTTATCTGAACGCGCACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTCCcccaccca < 2:679984/68‑1 (MQ=255) agcTTATCTGAACGCGCACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTCCcccaccca < 1:680202/68‑1 (MQ=255) cTTATCTGAACGCGCACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTCCCCCACCCATa < 1:691091/68‑1 (MQ=255) cgcgCACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGt < 1:224021/68‑1 (MQ=255) gCTGACCAGAAGCATAAACACGCACGCCCGCTTCGCCCACCCGTACCAACAAAGTTCCAACATgcgcg > 1:664846/1‑67 (MQ=255) gCTGACCAGAAGCATAAACACGACCGCCCGCTTCGCCCACCCATACCACCAATGTTCCAACTTGCGCa > 1:921370/1‑68 (MQ=255) gCTGACCAGAAGCATAAACACGAACGCCCGCTTCGCCCACCCATACCAACAATGTTCCAACTTGCGCa > 1:401305/1‑68 (MQ=255) gCTGACCAGAAGCATAAACACGAACGCCCGCTTCGCCCACCCATACCAACAATGTTCCAACTTGCGCa > 1:564625/1‑68 (MQ=255) gCTGACCAGAAGCATAAACACGAACGCCCGCTTCGCCCACCCATACCAACAATGTTCCAACTTGCGCa > 2:376505/1‑68 (MQ=255) gCTGACCAGAAGCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCa > 1:223307/1‑68 (MQ=255) gCTGACCAGAAGCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCa > 2:20391/1‑68 (MQ=255) cTGGCCGGAAGCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAAAGTTCCAACTTGCGCAg > 1:145312/1‑68 (MQ=255) cTGCCCAGAAGCATAAACACGAACGCCCGCTTCACCCACCCATACCAACAATGTTCCAACATGCGCAg > 1:321939/1‑68 (MQ=255) ccAGAAGCATAAACACGAACGCCCGCTTCGCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAg > 2:291135/1‑68 (MQ=255) aGAAGCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGc < 2:308166/68‑1 (MQ=255) cacGAGCGCCCGCTTGGCCCACCCATACCAACAATGTTCCGACTTGCGCAGCCAGGCGTACAGCTGCa > 1:686068/1‑68 (MQ=255) cacGAACGCCCGCTTCGCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACAGCTGCa > 2:267526/1‑68 (MQ=255) cacGAACGCCCGCTTCACCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACAGCTGCa > 1:792345/1‑68 (MQ=255) cGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACAGCTGCATg < 2:346602/68‑1 (MQ=255) ccGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACAGCTGCATGCGAAAcc < 2:223307/68‑1 (MQ=255) ccGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACAGCTGCATGCGAAAcc < 2:251172/68‑1 (MQ=255) ccGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACAGCTGCATGCGAAAcc < 2:824424/68‑1 (MQ=255) | GAGCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACAGCTGCATGCGAAACC > NC_000913/2879473‑2879596 |
| 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: ‑ |