Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 2,231,892 | 0 | A | C | 38.2% | 52.1 / 10.2 | 34 | T25P (ACC→CCC) | yejH | predicted ATP‑dependet helicase |
Reads supporting (aligned to +/- strand): ref base A (16/5); new base C (0/13); total (16/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.84e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCATGATTTTTACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCACGGCTGGCTC > REL606/2231818‑2231964 | ccATGATTTTTACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAAc > 2:9770/1‑76 (MQ=255) aTTTTTACACTTCGCTCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTg < 1:953787/76‑1 (MQ=255) ttttACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGcc > 2:696641/1‑76 (MQ=255) tttACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCg > 1:1029525/1‑76 (MQ=255) ttCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGt < 1:49663/76‑1 (MQ=255) gCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGtgct > 1:16585/1‑76 (MQ=255) cccATAACAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGtgctg < 1:64156/76‑1 (MQ=255) tatCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCa < 1:565364/76‑1 (MQ=255) naGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCg > 1:238181/2‑76 (MQ=255) aGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCgg < 1:798937/76‑1 (MQ=255) gCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCTTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGgc < 1:237771/76‑1 (MQ=255) cAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGgcg < 1:341265/76‑1 (MQ=255) aataagCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGgcgc < 2:548206/73‑1 (MQ=255) aagaagCGGGGGAGGCAACGCTCAGCCGTTGGCGTCGTCCTAAAACCCCTGCCTTTATCGTGCTGCCCACCGgcgc < 2:17406/76‑1 (MQ=255) aagaagCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGgcgc > 2:870216/1‑76 (MQ=255) aagaagCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCAACAGgcgc > 2:1227230/1‑76 (MQ=255) cggTGGAAGCCGCGCTCAACCATTTCCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTaa < 1:1181993/74‑1 (MQ=255) ccGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTaa > 2:1051180/1‑76 (MQ=255) ccGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTaa > 2:396384/1‑76 (MQ=255) gTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTaaaa > 1:960209/1‑76 (MQ=255) tGGATGCCACGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAg < 2:1063054/76‑1 (MQ=255) ccACGCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGt < 2:16585/76‑1 (MQ=255) cACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTg > 2:864413/1‑76 (MQ=255) gCTCAACCATTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATa < 2:1298322/76‑1 (MQ=255) cTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAg > 1:873664/1‑76 (MQ=255) cTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAg > 2:55140/1‑76 (MQ=255) aTTTTCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACt < 1:674908/76‑1 (MQ=255) ttCGTCTTCATAAACCCCCTGCCGTTATCGTGCTGACCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGc < 1:550033/76‑1 (MQ=255) ttCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGc < 2:1502071/76‑1 (MQ=255) ttCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGc < 2:1502086/76‑1 (MQ=255) ttCGTCGTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGc < 1:355442/76‑1 (MQ=255) tcgtcgTCATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCa < 2:699570/76‑1 (MQ=255) tcATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCACggctg > 1:1325934/1‑76 (MQ=255) cATAAACCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCACggctgg < 1:279581/76‑1 (MQ=255) cATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCACggctgg > 1:107553/1‑76 (MQ=255) aTAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCACggctggc > 1:999998/1‑76 (MQ=255) aTAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCACggctggc > 1:1221240/1‑76 (MQ=255) aaacCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCACGGCTGGCTc < 2:25058/76‑1 (MQ=255) aaaaCCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCACGGCCGGCTc > 1:862175/1‑76 (MQ=255) | CCATGATTTTTACACTTCGCCCATATCAGCAAGAAGCCGTGGATGCCACGCTCAACCATTTTCGTCGTCATAAAACCCCTGCCGTTATCGTGCTGCCCACCGGCGCAGGTAAAAGCCTGGTGATAGCGGAACTGGCACGGCTGGCTC > REL606/2231818‑2231964 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |