| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,566,599 | 0 | G | . | 40.4% | 132.2 / 79.8 | 57 | coding (351/1257 nt) | yfeO | predicted ion channel protein |
| Reads supporting (aligned to +/- strand): ref base G (20/14); new base . (0/23); total (20/37) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.67e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTG > minE/1566530‑1566662 | ggTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCATCCTgg > 1:3687457/1‑71 (MQ=255)ggTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTgg > 1:2487758/1‑71 (MQ=255)ggTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTgg > 1:1107723/1‑71 (MQ=255)ggTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTgg > 1:894242/1‑71 (MQ=255)ggTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTgg > 1:2337722/1‑71 (MQ=255)ggTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTgg > 1:1962994/1‑71 (MQ=255)ggTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTgg > 1:4211925/1‑71 (MQ=255)ggTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTgg > 1:4124144/1‑71 (MQ=255)ggTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTgg > 1:3420771/1‑71 (MQ=255) cgcTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATcc > 1:2733690/1‑70 (MQ=255) cgcTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATcc > 1:3827723/1‑70 (MQ=255) cgcTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATcc > 1:3508852/1‑70 (MQ=255) cTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGAt > 1:2435165/1‑71 (MQ=255) cTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGAt > 1:286869/1‑71 (MQ=255) cTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGAt > 1:3340117/1‑71 (MQ=255) cTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGAt > 1:4063707/1‑71 (MQ=255) cTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGAt > 1:137601/1‑71 (MQ=255) cTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGAt > 1:133075/1‑71 (MQ=255) cTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGAt > 1:1061106/1‑71 (MQ=255) cTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGTGCCGGAACATCCGATCATg < 1:4278097/71‑1 (MQ=255) cTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGCCGGCGTCAGCCTGGGGCCGGAACATCCGATCATg < 1:3051375/71‑1 (MQ=255) aaTTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCaa > 1:3806279/1‑60 (MQ=255) aaTTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCaa > 1:4296638/1‑60 (MQ=255) aaTTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCaa > 1:3113270/1‑60 (MQ=255) aaTTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCaa > 1:2647736/1‑60 (MQ=255) aaTTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCaa > 1:676367/1‑60 (MQ=255) aaTTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCaa > 1:2487382/1‑60 (MQ=255) aaTTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCaa > 1:89783/1‑60 (MQ=255) aaTTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCaa > 1:2000178/1‑60 (MQ=255) aaTTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCaa > 1:1955965/1‑60 (MQ=255) aaTTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCaa > 1:1805495/1‑60 (MQ=255) ttCTCGGTCTTGCTGGCGGCGTCAGCCT‑GGGCCGGAACATCCTATCATGACCGTCAATATCGCCCTTGc < 1:2825829/69‑1 (MQ=255) ttCTCGGTCTTGCTGGCGGCGTCAGCCT‑GGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGc < 1:3760095/69‑1 (MQ=255) ttCTCGGTCTTGCTGGCGGCGTCAGCCT‑GGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGc < 1:1001701/69‑1 (MQ=255) ttCTCGGTCTTGCTGGCGGCGTCAGCCT‑GGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGc < 1:599943/69‑1 (MQ=255) ttCTCGGTCTTGCTGGCGGCGTCAGCCT‑GGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGc < 1:4331613/69‑1 (MQ=255) ttCTCGGTCTTGCTGGCGGCGTCAGCCT‑GGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGc < 1:1378076/69‑1 (MQ=255) ttCTCGGTCTTGCTGGCGGCGTCAGCCT‑GGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGc < 1:1575041/69‑1 (MQ=255) ttCTCGGTCTTGCTGGCGGCGTCAGCCT‑GGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGc < 1:4228971/69‑1 (MQ=255) ttCTCGGTCTTGCTGGCGGCGTCAGCCT‑GGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGc < 1:1607748/69‑1 (MQ=255) ttCTCGGTCTTGCTGGCGGCGTCAGCCT‑GGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGc < 1:1612936/69‑1 (MQ=255) ttCTCGGTCTTGCTGGCGGCGTCAGCCT‑GGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGc < 1:4010270/69‑1 (MQ=255) ttCTCGGTCTTGCTGGCGGCGTCAGCCT‑GGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGc < 1:1777583/69‑1 (MQ=255) ttCTCGGTCTTGCTGGCGGCGTCAGCCT‑GGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGc < 1:1872522/69‑1 (MQ=255) ttCTCGGTCTTGCTGGCGGCGTCAGCCT‑GGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGc < 1:1912778/69‑1 (MQ=255) ttCTCGGTCTTGCTGGCGGCGTCAGCCT‑GGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGc < 1:341451/69‑1 (MQ=255) ttCTCGGTCTTGCTGGCGGCGTCAGCCT‑GGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGc < 1:3256692/69‑1 (MQ=255) ttCTCGGTCTTGCTGGCGGCGTCAGCCT‑GGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGc < 1:1981403/69‑1 (MQ=255) ttCTCGGTCTTGCTGGCGGCGTCAGCCT‑GGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGc < 1:2530867/69‑1 (MQ=255) ttCTCGGTCTTGCTGGCGGCGTCAGCCT‑GGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGc < 1:2873015/69‑1 (MQ=255) ttCTCGGTCTTGCTGGCGGCGTCAGCCT‑GGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGc < 1:2308522/69‑1 (MQ=255) ttCTCGGTCTTGCTGGCGGCGTCAGCCT‑GGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGc < 1:2665354/69‑1 (MQ=255) tctcGGTCTTGCTGGCGGCGTCAGCCT‑GGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGc < 1:1896961/68‑1 (MQ=255) tctcGGTCTTGCTGGCGGCGTCAGCCT‑GGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGc < 1:2402444/68‑1 (MQ=255) cgTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTc < 1:1247614/71‑1 (MQ=255) tCAGCCTGGGGCCGGAACATCGGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTg < 1:3556586/71‑1 (MQ=255) tCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTg < 1:4235362/71‑1 (MQ=255) tCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTg < 1:4258763/71‑1 (MQ=255) tCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTg < 1:2712854/71‑1 (MQ=255) tCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTg < 1:2274952/71‑1 (MQ=255) tCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTg < 1:2214318/71‑1 (MQ=255) tCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTg < 1:1626448/71‑1 (MQ=255) tCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTg < 1:620226/71‑1 (MQ=255) tCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTg < 1:3930521/71‑1 (MQ=255) tCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTg < 1:2423209/71‑1 (MQ=255) gcAGCCGGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTg < 1:1113175/70‑1 (MQ=255) | GGTTCCGCCCTCTGCGCTACCTGGACTTATCGTAGCATTAATTCTCGGTCTTGCTGGCGGCGTCAGCCTGGGGCCGGAACATCCGATCATGACCGTCAATATCGCCCTTGCGGTTGCCATTGGCGCTCGTCTG > minE/1566530‑1566662 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |