Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,235,452 | 0 | T | G | 25.0% | 64.5 / 10.5 | 32 | F85V (TTT→GTT) | yhbQ | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (4/20); new base G (8/0); total (12/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.71e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGTTACGGGCGGAATATCGCGTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCA‑GGGTTTGCAGAGCTGTTAAGCAGTCTG‑CAAACCCCGGAGATTAAAAGCGATTGAAATGCTCGTGATACTCAACCAGGC > REL606/3235378‑3235526 | cGTTACGGGCGGAATATCGCGTTAAGCAACTGACCAAACGACAAAAAGAGCGTCTGGTAGCGGAAGGCGCA‑GGGtt < 2:961889/76‑1 (MQ=255) gTTGCGGGCGGAATAGCGCGTTAAGCAACTGACCAAACGGCAAGAAGAGCGTCTGGTAGCGGAAGGCGCA‑GGGttt < 2:326265/76‑1 (MQ=255) gTTACGGGCGGAATATCGCGTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCA‑GGGttt < 1:1332035/76‑1 (MQ=255) gTTACGGGCGGAATATCGCGTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCA‑GGGttt < 2:61085/76‑1 (MQ=255) ttACGGGCGGAATATCGCGTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCA‑GGGGTTg > 2:1829386/1‑76 (MQ=255) tACGGGCGGAATATCGCGTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCA‑GGGTTTGc < 2:1078046/76‑1 (MQ=255) aCGGGCGGAATATCGCGTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCA‑GGGTTTGCa < 2:742616/76‑1 (MQ=255) cGGGCGGAATATCGCGTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCA‑GGGTTTGCag < 2:171743/76‑1 (MQ=255) ggAATATCGCGTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCA‑GGGTTTGCAGAGCTg < 1:1271957/76‑1 (MQ=255) ggAATATCGCGTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCA‑GGGTTTGCAGAGCTg < 2:1569517/76‑1 (MQ=255) ggAATATCGCGTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCA‑GGGGGTGCGGAGcgg > 1:1865812/1‑74 (MQ=255) aaTATCGCGTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCA‑GGGGTTGCAGAGCTGtt > 1:1671910/1‑76 (MQ=255) atatCGCGTTAAGCAACTGACCAAACGGCAAAATGAGCGTCTGGTAGCGGAAGGCGCA‑GGGTTTGCAGAGCTGTTa < 2:139992/76‑1 (MQ=255) atatCGCGTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCA‑GGGTTTGCAGAGCTGTTa < 2:1671779/76‑1 (MQ=255) atCGCGTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCA‑GGGTTTGCAGAGCTGTTaan < 1:500960/76‑2 (MQ=255) tCGCGTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCA‑GGGTTTGCAGAGCTGTTAAGc < 1:894468/76‑1 (MQ=255) cgcgTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCA‑GGGTTTGCAGAGCTGTTAAGCa < 1:1909343/76‑1 (MQ=255) cgTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCA‑GGGGGTGCGGAGCCGTTATGCAga > 2:1146882/1‑75 (MQ=255) ttAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCA‑GGGTTTGCAGAGCTGTTAAGCAGTCt < 1:1602132/76‑1 (MQ=255) ttAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCA‑GGGTTTGCAGAGCTGTTAAGCAGTCt < 1:910614/76‑1 (MQ=255) aaGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCAGGGGTTTGCAGGGC‑GTTGAGCGGTGTG‑c > 2:1471489/1‑76 (MQ=255) gCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCA‑GGGTTTGCAGAGCTGTTAAGCAGTCTG‑Caa < 2:228950/76‑1 (MQ=255) gCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCA‑GGGTTTGCAGAGCTGTTAAGCAGTCTG‑Caa < 2:27779/76‑1 (MQ=255) cAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCA‑GGGGGTGCAGAGCTGTTAAGCAGTCTG‑Caaa > 2:1981792/1‑76 (MQ=255) aaCTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCA‑GGGTTTGCAGAGCTGTTAAGCAGTCTG‑CAAAc < 2:1671910/76‑1 (MQ=255) acggACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCA‑GGGTTTGCAGAGCTGTTAAGCAGTCTG‑CAAAcc < 2:34360/73‑1 (MQ=255) tGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCA‑GGGGTTGCAGAGGTGTTTAAGAGGCTG‑CAAAaccc > 1:1719228/1‑76 (MQ=255) gACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCA‑GGGTTTGCAGAGCTGTTAAGCAGTCTG‑CAAACCCCg < 2:362160/76‑1 (MQ=255) gACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCA‑GGGTTTGCAGAGCTGTTAAGCAGTCTG‑CAAACCCCg < 2:49879/76‑1 (MQ=255) gACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCA‑GGGTTTGCAGAGCTGTTAAGCAGTCTG‑CAAACCCCg < 2:1966421/76‑1 (MQ=255) ccAAACGGAAAAAGGGGCTTCTGGTGGCGGAGGGCGCG‑GGTTTTGCAGAGTTTTAAAGCAGTCTG‑CAACCCCCGga < 2:853858/76‑1 (MQ=255) aaaCGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCA‑GGGGGTGCAGAGATGGTGAGGAGAATG‑GACACCCCggggg > 2:1579551/1‑73 (MQ=255) aaaCGGAAAAAAGGGCTTCGGGGGGCGGAAGGCGCA‑GGTTTTGCAGAGCTGTTAAGCAGTCTG‑CAACCCCCGgaga < 2:1069915/76‑1 (MQ=255) gaaaaaaGAGCGTCGGGTGGGGGAGGGCGCG‑GGGTTTGCAGAGTTTTAAGGCAGTCTG‑CAACCCCCGGAGATTaaa < 2:532339/74‑1 (MQ=255) cAAAAAGAGCGTCTGGTAGCGGAAGGCGCA‑GGGGTTGCAGAGCCGTTAAGGAGGCTG‑CAAAACCCGGGGAATtaaa > 2:1345571/1‑76 (MQ=255) aaaaaGAGCGTCTGGTAGCGGAAGGCGCA‑GGGGGTGTCGAGCTGTTTAGCAGTGTGCCAAACCCCGGGGGTTaaaa > 1:156453/1‑76 (MQ=255) aaaGAGCGTCTGGTAGCGGAAGGCGCA‑GGGGTTGCAGAGCTGTTTAGCAGTCTG‑CAAAACCCGGGGATTAAAAACg > 2:1868792/1‑76 (MQ=255) aaGAGCGTCTGGTAGCGGAAGGCGCA‑GGGGTTGCAGAGCTGTTAAGCAGTCTG‑CAAAACCCGGGGGTTAAAAgggg > 2:38354/1‑73 (MQ=255) acngCGTCTGGTAGCGGAAGGCGCA‑GGGGTTGCAGAGCTGTTAAGCAGTCTG‑CAAACCCCGGAGGTTAAAAGCGAt > 2:1406126/4‑76 (MQ=255) agCGTCTGGTAGCGGAAGGCGCA‑GGGGTTGCAGAGCTGTTAAGCAGTCTG‑CCAACCCCGGGGATTAAAAGGGAATg > 2:205758/1‑76 (MQ=255) cGTGTGGTAGCGGAGGGCGCG‑GGGTTTGCAGAGCTGTTAAGCAGTCTG‑CAAACCCCGGAGATTAAAAGCGATTGaa < 2:1903816/76‑1 (MQ=255) cTGGTAGCGGAAGGCGCA‑GGGGTTGCAGAGGTGTTTAGCAGGCTG‑CAAACCCCCGGGGTTAAAAACGATTGaaaat > 2:619189/1‑74 (MQ=255) cTGGTAGCGGAAGGCGCA‑GGGGTTGCAGAGCTGTTAAGCAGTCTG‑CAAACCCCGGAGATTAAAAACGATTGAAATg > 2:1446616/1‑76 (MQ=255) cTGGTAGCGGAAGGCGCA‑GGGGGTGCAGAGCTGTTAAGCAGTCTG‑CAAACCCCGGAGATTAAAAACGATTGaaaag > 2:1873855/1‑74 (MQ=255) tAGCGGAAGGCGCA‑GGGTTTGCAGAGCTGTTAAGCAGTCTG‑CAGACCCCGGAGATTAAAAGCGATTGAAATGCTCg > 1:393910/1‑76 (MQ=255) tAGCGGAAGGCGCA‑GGGGTTGCAGAGCTGTTAAACAGTCTG‑CAAACCCCCGAGATTAAAAACGAATGAAAAGATCg > 2:1921236/1‑76 (MQ=255) gagggCGCG‑GGGTTGGCAGGTCTGTAAAGCAGTCTG‑CAAACCCCGGAGATTAAAAGCGATTGAAATGCTCGTGATa < 2:1841002/73‑1 (MQ=255) aaGGCGCA‑GGGGCTGCAGAGGCGTTTAGCAACCTG‑CAAACCCCGGGGGTTAAAAACGATTGAAAAGCTTGTGGTAc > 2:387844/1‑76 (MQ=255) cA‑GGGTTTGCAGAGCTGTTAAGCAGTCTG‑CAAACCCCGGAGATTAAAAGCGATTGAAATGCTCGTGATACTCAAcc > 2:1969226/1‑76 (MQ=255) a‑GGGTTTGCAGAGCTGTTAAGCAGTCTG‑CAAACCCCGGAGATTAAAAGCGATTGAAATGCTCGTGATACTCAACCa > 2:1822868/1‑76 (MQ=255) a‑GGGTTTGCAGAGCTGTTAAGCAGTCTG‑CAAACCCCGGAGATTAAAAGCGATTGAAATGCTCGTGATACTCAACCa > 2:1337899/1‑76 (MQ=255) a‑GGGGTCGCAGAGCTGTTAAGCAGTCTG‑CAAACCCCGGGGATTAAAAGCGGTTGAAATGCCCGTGATCCTCAATCa > 2:1502386/1‑76 (MQ=255) ggTTGGAAGAGTTGTAAAGCAGTCTA‑CAACCCCCGGAGATTAAAAGCGATTGAAATGCTCGTGATACTCAACCAgg < 1:1576820/76‑1 (MQ=255) gTTTGCAGAGCTGTTAAGCAGTCTG‑CAGACCCCGGAGATTAAAAGCGATTGAAATGCTCGTGATACTCAACCAGGc < 2:393910/76‑1 (MQ=255) | CGTTACGGGCGGAATATCGCGTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCA‑GGGTTTGCAGAGCTGTTAAGCAGTCTG‑CAAACCCCGGAGATTAAAAGCGATTGAAATGCTCGTGATACTCAACCAGGC > REL606/3235378‑3235526 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |