Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,879,536 | 0 | C | G | 40.7% | 35.4 / 12.9 | 28 | G79G (GGG→GGC) | ygbT | multifunctional endonuclease Cas1, CRISPR adaptation protein; DNA repair enzyme |
Reads supporting (aligned to +/- strand): ref base C (2/14); new base G (11/0); total (14/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.98e-06 | |||||||||||
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 | gagCAGCTTATCTGAACGCGCACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTCCccca < 1:488005/68‑1 (MQ=255) agcagcTTCTCTGAACGCGCACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTCCcccac < 2:112832/68‑1 (MQ=255) agcagcTTATCTGAACGCGCACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTCCcccac < 1:77262/68‑1 (MQ=255) agcTTATCTGAACGCGCACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTCCcccaccca < 1:715439/68‑1 (MQ=255) agcTTATCTGAACGCGCACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTCCcccaccca < 1:813266/68‑1 (MQ=255) aTCTGAACGCGCACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTCCCCCACCCATACca < 2:345344/68‑1 (MQ=255) tGAACGCGCACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTGCCCCACCCATACcaaca > 1:95597/1‑68 (MQ=255) tGAAAGCGCACCTCCAGGCTGACCAGAAGCATAAACACGAACGACCGCTTCACCCACCCATACcaaca > 2:214275/1‑68 (MQ=255) aCGCGCACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTCGCCCACCCATACCAACAATg > 1:705605/1‑68 (MQ=255) cACCTCCAGGCTGACCAGAAGCATAAACACGAACGCCCGCTTCGCCCACCCAAACCAACAAAGTTcaa > 2:605670/1‑66 (MQ=255) gCTGACCAGAAGCATAAACACGCACGCCCGCTTCGCCCGCCCATACCAACAATGTTCCAACTTGCGCa > 1:534874/1‑68 (MQ=255) gCTGACCAGAAGCATAAACACGAACGCCCGCTTGGCCCACCCATACCAACAATGTTCCAACTTGCGCa > 1:425763/1‑68 (MQ=255) gCTGACCAGAAGCATAAACACGAACGCCCGCTTCGCCCGCCCATACCAACAATGTTCCAACTTGCGCa > 1:70800/1‑68 (MQ=255) gCTGACCAGAAGCATAAACACGAACGCCCGCTTAGCCCAACCATACCAACAATGTTAAAACATGCGCa > 2:763646/1‑68 (MQ=255) aCCAGAAGCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCa < 1:223514/68‑1 (MQ=255) ccAGAAGCATAAACACGAACGCCCGCTTCGCCCCCACATACCAACAATGTTCCAACTTGCGCAGCCAg > 2:106176/1‑68 (MQ=255) cAGAAGCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAgg < 2:129957/68‑1 (MQ=255) aGAAGCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGc < 1:218027/68‑1 (MQ=255) aGAAGCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGc < 2:755238/68‑1 (MQ=255) aaGCATAAACACGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGt < 2:154575/68‑1 (MQ=255) aGCATAAACACGAACGCCCGCTTAGCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTa > 1:175816/1‑68 (MQ=255) cATAAACACGAACGCCCGCTTGGCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACa > 2:584488/1‑68 (MQ=255) cacGAACGCCCGCTTCGCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACAGCTGCa > 2:560488/1‑68 (MQ=255) cacGAACGCCCGCTTCCCCCACCCATACAAACAATGTTCCAACTTGCGCAGCCAGGCGTACAGCTGCa > 2:227571/1‑68 (MQ=255) acGAACGCCCGCTTCGCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACAGCTGCAt > 2:458390/1‑68 (MQ=255) cGAACGCCCGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACAGCTGCATg < 2:234086/68‑1 (MQ=255) agAACGCCCTCTTCCCCCACCCATACCTACTATGTTCCAACTTGCGCAGCCAGGCGTACAGCTGCATg < 1:789792/67‑1 (MQ=255) ccGCTTCCCCCACCCATACCAACAATGTTCCAACTTGCGCAGCCAGGCGTACAGCTGCATGCGAAAcc < 2:705607/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: ‑ |