| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_007779 | 3,300,073 | 0 | T | G | 46.2% | 40.8 / 18.5 | 40 | F85V (TTT→GTT) | yhbQ | endonuclease |
| Reads supporting (aligned to +/- strand): ref base T (5/16); new base G (18/0); total (24/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 9.98e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CGCCAGTGGGCGATCGTTCGCTGGCGTTACGGGCGGAATATCGCGTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCAGGGTTTGCAGAGCTGTTAAGCAGTC‑TGCAAACCCCGGAGATTAAAAGCGATTGAAATGCTCGTGATACTCAACCAGGCCGGTAACGCCATTCAGCGCGTCATCC > NC_007779/3299975‑3300173 | cgcCAGTGGGCGATCGTTCGCTGGCGTTACGGGCGGAATATCGCGTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCCGGGttt > 1:1318321/1‑101 (MQ=255) gcCAGTGGGCGATCGTTCGCTGGCGTTACGGGCGGAATATCGCGTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCAGGGTTTg < 1:857876/101‑1 (MQ=255) gcCAGTGGGCGATCGTTCGCTGGCGTTACGGGCGGAATATCGCGTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCAGGGGTTg > 3:956229/1‑101 (MQ=255) aGTGGGCGATCGTTCGCTGGCGTTACGGGCGGAATATCGCGTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCAGGGGTTGCag > 3:408201/1‑101 (MQ=255) tGGGCGATCGTTCGCTGGCGTTACGGGCGGAATATCGCGTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCAGGGTTTGCagag < 1:3345/101‑1 (MQ=255) gATCGTTCGCTGGCGTTACGGGCGGAATATCGCGTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCAGGGGTTGCAGAGCTGtt > 1:37939/1‑101 (MQ=255) gATCGTTCGCTGGCGTTACGGGCGGAATATCGCGTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCAGGGGTTGCAGAGCTGtt > 1:154625/1‑101 (MQ=255) gATCGTTCGCTGGCGTTACGGGCGGAATATCGCGTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCAGGGGGTGCAGAGCTGtt > 1:1422501/1‑101 (MQ=255) aTCGTTCGCTGGCGTTACGGGCGGAATATCGCGTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCAGGGGTTGCAGAGCTGTTa > 2:138415/1‑101 (MQ=255) tCGCTGGCGTTACGGGCGGAATATCGCGTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCAGGGTTTGCAGAGCTGTTAAGCAg < 3:1216115/101‑1 (MQ=255) tCGCTGGCGTTACGGGCGGAATATCGCGTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCAGGGTTTGCAGAGCTGTTAAGCAg < 3:1265777/101‑1 (MQ=255) cTGGCGTTACGGGCGGAATATCGCGTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCAGGGTTTGCAGAGCTGTTAAGCAGTC‑t < 1:1289593/101‑1 (MQ=255) tGGCGTTACGGGCGGAATATCGCGTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCAGGGTTTGCAGAGCTGTTAAGCAGTC‑Tg < 1:114593/101‑1 (MQ=255) gTTACGGGCGGAATATCGCGTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCAGGGGTTGCAGAGCTGTTTAGCAGTC‑TGCaaa > 1:1352101/1‑101 (MQ=255) ttACGGGCGGAATATCGCGTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCAGGGGTTGCAGAGCTGTTTAGCAGTC‑TGCAAAc > 1:1341020/1‑101 (MQ=255) tACGGGCGGAATATCGCGTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCAGGGGTTGCAGAGCTGTTAAGCAGTC‑TGCcaaac > 3:1414121/1‑99 (MQ=255) gCGGAATATCGCGTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCAGGGGTTGCAGAGCTGTTAAGCAGTC‑TGCAAACCCCggg > 1:705555/1‑100 (MQ=255) cGGAATATCGCGTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCAGGGTTTGCAGAGCTGTTAAGCAGTC‑TGCAAACCCCGgag < 1:660680/101‑1 (MQ=255) ggAATATCGCGTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCAGGGTTTGCAGAGCTGTTAAGCAGTC‑TGCAAACCCCGgaga < 1:141529/101‑1 (MQ=255) ggAATATCGCGTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCAGGGGTTGCCGAGCTGTTTAGCAGTCTTCCAAACCCCggga > 3:18800/1‑99 (MQ=255) aaTATCGCGTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCAGGGGTTGCAGAGCTGTTAAGCAGTC‑TGCAAACCCCGGAGAtt > 1:1151824/1‑101 (MQ=255) cgcgTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCAGGGTTTGCAGAGCTGTTAAGCAGTC‑TGCAAACCCCGGAGATTAAAAg > 1:1285474/1‑101 (MQ=255) cgTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCAGGGTCTGCCGAACTGTTAAGCAGGC‑TGCAAACCCCCGAGACTAAAAACg > 1:902432/1‑101 (MQ=255) cgTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCAGGGGTTGCAGAGCTGTTAAGCAGTC‑TGCAAACCCCCGAGAGTAAAAGCg > 1:1039296/1‑101 (MQ=255) tAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCAGGGTTTGCAGAGCTGTTAAGCAGTC‑TGCAAACCCCGGAGATTAAAAGCGAtt < 1:1552242/101‑1 (MQ=255) aaaGCAACTGACCAAACGGAAAAAAGAGCGTCTGGTAGCGGAAGGCGCGGGGTTTGCAGAGCTGTTAAGCAGTC‑TGCAACCCCCGGAGATTAAAAGCGAtt < 1:563503/100‑1 (MQ=255) aaGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCAGGGGTTGCAGAGCTGTTAAGCAGTC‑TGCAAACCCCGGAGATTAAAAGCGATTg > 1:240955/1‑101 (MQ=255) cAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCAGGGTTTGCAGAGCTGTTAAGCAGTC‑TGCAAACCCCGGAGATTAAAAGCGATTGaaa < 1:1289562/101‑1 (MQ=255) aaCTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCAGGGGTTGCAGAGCTGTTAAGCAGTC‑TGCAAACCCCGGAGATTAAAAGCGATTGAAAt > 1:154014/1‑101 (MQ=255) cAAAAAGAGCGTCTGGTAGCGGAAGGCGCAGGAACTTCAGAGCTGTTAAGCACTC‑CGTAAACCCCGGAGCTTAAAAGCCATTCAAATGCTCGTGAtctcca > 3:1307224/1‑96 (MQ=255) aaaaaGAGCGTCTGGTAGCGGAAGGCGCAGGGGTTGCAGAGCTGTTAAGCAGTC‑TGCAAACCCCGGAGATTAAAAGCGATTGAAATGCTCGTGATACTCaa > 1:371904/1‑101 (MQ=255) aaGAGCGTCTGGTAGCGGAAGGCGCAGGGTTTGCAGAGCTGTTAAGCAGTC‑TGCAAACCCCGGAGATTAAAAGCGATTGAAATGCTCGTGATACTCAACCa < 1:946806/101‑1 (MQ=255) agCGTGGGGTAGGGGAAGGCGGGGGGTTTGCAGAGCTGTTAAGCAGTC‑TGAAAACCCCGGAGATTAAAAGCGATTGAAATGCTCGTGATACTCAACCAGGc < 1:325368/101‑1 (MQ=255) cGTCTGGTAGCGGAAGGCGCAGGGGTTGCAGAGCTGTTAAGCAGTC‑TGCAAACCCCGGAGATTAAAAGCGATTGAAATGCTCGTGATACTCAACCAGGCCg > 1:684233/1‑101 (MQ=255) gTCTGGTAGCGGAAGGCGCAGGGGTTGCAGAGCTGTTTAGCAGTC‑TGCAAACCCCGGAGATTAAAAGCGATTGAAATGCTCGTGATACTCAACCAGGCCgg > 3:1392404/1‑101 (MQ=255) gTCGGGGAGCGGAGGGCGCAGGGTTTGCAGAGCTGTTAAGCAGTC‑TGCAAACCCCGGAGATTAAAAGCGATTGAAATGCTCGTGATACTCAACCAGGCCgg < 3:1268015/101‑1 (MQ=255) tGGTAGCGGAAGGCGCAGGGTTTGCAGAGCTGTTAAGCAGTC‑TGCAAACCCCGGAGATTAAAAGCGATTGAAATGCTCGTGATACTCAACCAGGCCGGTaa > 1:452467/1‑101 (MQ=255) tAGCGGAAGGCGCAGGGTTTGCAGAGCTGTTAAGCAGTC‑TGCAAACCCCGGAGATTAAAAGCGATTGAAATGCTCGTGATACTCAACCAGGCCGGGAACGc > 1:365449/1‑101 (MQ=255) ggAAGGCGCAGGGTTTGCAGAGCTGTTAAGCAGTC‑TGCAAACCCCGGAGATTAAAAGCGATTGAAATGCTCGTGATACTCAACCAGGCCGGTAACGCCAtt > 1:1187576/1‑101 (MQ=255) gAAGGCGCAGGGTTTGCAGAGCTGTTAAGCAGTC‑TGCAAACCCCGGAGATTAAAAGCGATTGAAATGCTCGTGATACTCAACCAGGCCGGTAACGCCATTc < 1:220496/101‑1 (MQ=255) gcAGGGTTTGCAGAGCTGTTAAGCAGTC‑TGCAATCCCCGGAGATTAAAAGCGATTGAAATGCTCGTGATACTCAACCAGGCCGGTAACGCCATTCAGcgcg < 1:508197/101‑1 (MQ=255) tttGCAGAGCTGTTAAGCAGTC‑TGCAAACCCCGGAGATTAAAAGCGATTGAAATGCTCGTGATACTCAAACAAGCCGGGAACGCCAATCAGcgcgcccccc > 1:1311550/1‑95 (MQ=255) | CGCCAGTGGGCGATCGTTCGCTGGCGTTACGGGCGGAATATCGCGTTAAGCAACTGACCAAACGGCAAAAAGAGCGTCTGGTAGCGGAAGGCGCAGGGTTTGCAGAGCTGTTAAGCAGTC‑TGCAAACCCCGGAGATTAAAAGCGATTGAAATGCTCGTGATACTCAACCAGGCCGGTAACGCCATTCAGCGCGTCATCC > NC_007779/3299975‑3300173 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |