Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,878,956 | 0 | C | T | 39.3% | 24.0 / 23.1 | 28 | A273T (GCC→ACC) | ygbT | multifunctional endonuclease Cas1, CRISPR adaptation protein; DNA repair enzyme |
Reads supporting (aligned to +/- strand): ref base C (10/7); new base T (4/7); total (14/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.40e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATAT > NC_000913/2878889‑2879015 | aCAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCACCAGCGGc > 1:804122/1‑68 (MQ=255) gCGGAATGGCAACAGGCTGTGCATCTTCAGGCGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACg < 2:169666/68‑1 (MQ=255) ggAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTc < 1:706174/68‑1 (MQ=255) gAAAGGAAAAAGGCTGTGAATCTACAGGTGGGGCCGGAGGTTGTATTTCTACAGCGGTAAGAACgaca > 1:471230/1‑65 (MQ=255) tGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGTAAGCACGTCctct < 1:17914/68‑1 (MQ=255) ggCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCgg > 1:1077841/1‑68 (MQ=255) ggCTGAGCATATTCAGGTGGGGCCGGCGGAAGTATTTCTCCAGCGGTAAGCACGTCCTCTATAAGCgg > 2:175367/1‑68 (MQ=255) gCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGTAAGCACGTCCTCTATAAGCGGa > 2:370935/1‑68 (MQ=255) gCTGTGCAGCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGa < 1:414856/68‑1 (MQ=255) gCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCaa < 1:791510/68‑1 (MQ=255) cATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAAt > 1:836064/1‑68 (MQ=255) ttCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGTAAGCGCGTCCTCTATAAGCGGAATCAATTTgg > 2:845099/1‑68 (MQ=255) ttCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTgg > 1:864511/1‑68 (MQ=255) gTGGGGCCGGCGGTTGTATTTCTCCAGCGGTAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAAt < 1:370935/68‑1 (MQ=255) gTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAAt > 2:730529/1‑68 (MQ=255) gTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAAt > 1:904434/1‑68 (MQ=255) gTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAAt < 1:942950/68‑1 (MQ=255) gTGGGGCCGGCGGGTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAAt > 2:317389/1‑68 (MQ=255) gggCCGGCGGTTGTATTTCTCCAGCGTTAATCACGTCCACTATAAGCTTAATCAATTTGTCTAAtttt < 1:175367/68‑1 (MQ=255) gggCCGGCGGTTGTATTTCTCCAGCGGTAAGCACGTCCTCTATAAGCTGATTCAATTTGGCTAATGtt < 1:153233/68‑1 (MQ=255) gggCCGGCGGTTGTATTTCTCCAGCGGTAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGtt < 1:584774/68‑1 (MQ=255) gggCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGtt < 2:395572/68‑1 (MQ=255) gggCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGtt < 2:434056/68‑1 (MQ=255) ccGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAAGGTTTta > 2:522161/1‑68 (MQ=255) tttttCTCCAGCGGTTTTCACGTCCTCTATAAGCGTAATCAATTTGGCTTATGTTTTACTACAGCGaa < 2:471230/66‑1 (MQ=255) tttCTCCAGCGGTAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGaaaa < 1:467439/68‑1 (MQ=255) ttCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGaaaaa > 2:966997/1‑68 (MQ=255) ctcCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAata > 2:298625/1‑68 (MQ=255) tcCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAatgt > 2:413706/1‑66 (MQ=255) | ACAGGAAGCGGAATGGCAACAGGCTGTGCATCTTCAGGTGGGGCCGGCGGTTGTATTTCTCCAGCGGCAAGCACGTCCTCTATAAGCGGAATCAATTTGGCTAATGTTTTACTACTGCGAAAAATAT > NC_000913/2878889‑2879015 |
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: ‑ |