| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 2,406,948 | 1 | . | C | 46.4% | 169.8 / 236.3 | 97 | coding (513/2547 nt) | cyaA | adenylate cyclase |
| Reads supporting (aligned to +/- strand): ref base . (34/18); new base C (27/18); total (61/36) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.75e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CGGCGGTACGATAAAATTCGTCAAGCAGCAGTATATGCTGGGTGGAGCCACAATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCTGACTTC > minE/2406890‑2407013 | cggcggTACGATAAAATTCGTCAAGCAGCAGTATATGCTGGGTGGAGCCACAATCTTCGCCCCCCCAGGCt > 1:4357716/1‑71 (MQ=255)cggcggTACGATAAAATTCGTCAAGCAGCAGTATATGCTGGGTGGAGCCACAATCTTCGCCCCCCCAGGCt > 1:3977547/1‑71 (MQ=255) gTACGATAAAATTCGTCAAGCAGCAGTATATGCTGGGTGGAGCCACAATCTTCGCCCCCCCAGGCTGCCGc < 1:500364/71‑1 (MQ=255) gTACGATAAAATTCGTCAAGCAGCAGTATATGCTGGGTGGAGCCACAATCTTCGCCCCCCCAGGCTGCCGc < 1:1863999/71‑1 (MQ=255) gTACGATAAAATTCGTCAAGCAGCAGTATATGCTGGGTGGAGCCACAATCTTCGCCCCCCCAGGCTGCCGc < 1:2411/71‑1 (MQ=255) gTACGATAAAATTCGTCAAGCAGCAGTATATGCTGGGTGGAGCCACAATCTTCGCCCCCCCAGGCTGCCGc < 1:2746246/71‑1 (MQ=255) gTACGATAAAATTCGTCAAGCAGCAGTATATGCTGGGTGGAGCCACAATCTTCGCCCCCCCAGGCTGCCGc < 1:2859271/71‑1 (MQ=255) gTACGATAAAATTCGTCAAGCAGCAGTATATGCTGGGTGGAGCCACAATCTTCGCCCCCCCAGGCTGCCGc < 1:475180/71‑1 (MQ=255) gTACGATAAAATTCGTCAAGCAGCAGTATATGCTGGGTGGAGCCACAATCTTCG‑CCCCCCAGGCTGCCGCt > 1:118875/1‑71 (MQ=255) tACGATAAAATTCGTCAAGCAGCAGTATATGCTGGGTGGAGCCACAATCTTCGCCCCCCCAGGCTGCCGc < 1:2925267/70‑1 (MQ=255) aCGATAAAATTCGTCAAGCAGCAGTATATGCTGGGTGGAGCCACAATCTTCGCCCCCCCAGGCTGCCGCtt > 1:1285841/1‑71 (MQ=255) aCGATAAAATTCGTCAAGCAGCAGTATATGCTGGGTGGAGCCACAATCTTCGCCCCCCCAGGCTGCCGCtt > 1:393441/1‑71 (MQ=255) aCGATAAAATTCGTCAAGCAGCAGTATATGCTGGGTGGAGCCACAATCTTCGCCCCCCCAGGCTGCCGCtt > 1:4009204/1‑71 (MQ=255) aCGATAAAATTCGTCAAGCAGCAGTATATGCTGGGTGGAGCCACAATCTTCGCCCCCCCAGGCTGCCGCtt > 1:2274875/1‑71 (MQ=255) aCGATAAAATTCGTCAAGCAGCAGTATATGCTGGGTGGAGCCACAATCTTCGCCCCCCCAGGCTGCCGCtt > 1:2812527/1‑71 (MQ=255) aCGATAAAATTCGTCAAGCAGCAGTATATGCTGGGTGGAGCCACAATCTTCGCCCCCCCAGGCTGCCGCtt > 1:2977320/1‑71 (MQ=255) aCGATAAAATTCGTCAAGCAGCAGTATATGCTGGGTGGAGCCACAATCTTCGCCCCCCCAGGCTGCCGCtt > 1:2312141/1‑71 (MQ=255) aaTTCGTCAAGCAGCAGTATATGCTGGGTGGAGCCACAATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTAt < 1:118729/70‑1 (MQ=255) ttCGTCAAGCAGCAGTATATGCTGGGTGGAGCCACAATCTTCGCCCCCCCAGGCTGCCGCtt > 1:2704115/1‑62 (MQ=255) ttCGTCAAGCAGCAGTATATGCTGGGTGGAGCCACAATCTTCGCCCCCCCAGGCTGCCGCTTTCATTATGa > 1:1684441/1‑71 (MQ=255) ttCGTCAAGCAGCAGTATATGCTGGGTGGAGCCACAATCTTCGCCCCCCCAGGCTGCCGCTTTCATTATGa > 1:3625413/1‑71 (MQ=255) ttCGTCAAGCAGCAGTATATGCTGGGTGGAGCCACAATCTTCGCCCCCCCAGGCTGCCGCTTTCATTATGa > 1:1566993/1‑71 (MQ=255) ttCGTCAAGCAGCAGTATATGCTGGGTGGAGCCACAATCTTCGCCCCCCCAGGCTGCCGCTTTCATTATGa > 1:3520895/1‑71 (MQ=255) cagcagTATATGCTGGGTGGAGCCACAATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGt > 1:4102856/1‑71 (MQ=255) cagcagTATATGCTGGGTGGAGCCACAATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGt > 1:1942518/1‑71 (MQ=255) cagcagTATATGCTGGGTGGAGCCACAATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGt > 1:2090367/1‑71 (MQ=255) cagcagTATATGCTGGGTGGAGCCACAATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGt > 1:2998445/1‑71 (MQ=255) cagTATATGCTGGGTGGAGCCACAATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGttt > 1:3423287/1‑70 (MQ=255) cagTATATGCTGGGTGGAGCCACAATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGttt > 1:3993093/1‑70 (MQ=255) tatGCTGGGTGGAGCCACAATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACg > 1:2621733/1‑55 (MQ=255) tatGCTGGGTGGAGCCACAATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTtcatc > 1:3255918/1‑70 (MQ=255) tatGCTGGGTGGAGCCACAATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTtcatc > 1:2632989/1‑70 (MQ=255) tatGCTGGGTGGAGCCACAATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTtcatc > 1:322044/1‑70 (MQ=255) tatGCTGGGTGGAGCCACAATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTtcatc > 1:1498829/1‑70 (MQ=255) tatGCTGGGTGGAGCCACAATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTtcatc > 1:3740381/1‑70 (MQ=255) tatGCTGGGTGGAGCCACAATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTtcatc > 1:167264/1‑70 (MQ=255) gggTGGAGCCACAATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAg > 1:3515407/1‑70 (MQ=255) gggTGGAGCCACAATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAg > 1:1597702/1‑70 (MQ=255) gggTGGAGCCACAATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAg > 1:184248/1‑70 (MQ=255) gggTGGAGCCACAATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAg > 1:1813548/1‑70 (MQ=255) gggTGGAGCCACAATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAg > 1:1766577/1‑70 (MQ=255) gAGCCACAATCTTCGCCCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGgaaga < 1:939227/71‑1 (MQ=255) aGCCACAATCTTCGCCCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGgaagaa < 1:1073745/71‑1 (MQ=255) aGCCACAATCTTCGCCCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGgaagaa < 1:2196176/71‑1 (MQ=255) gCCACAATCTTCGCCCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGaagaag > 1:3243813/1‑71 (MQ=255) gCCACAATCTTCGCCCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGaagaag > 1:482273/1‑71 (MQ=255) ccACAATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCt > 1:3968035/1‑71 (MQ=255) ccACAATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCt < 1:3016279/71‑1 (MQ=255) ccACAATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCt < 1:2986446/71‑1 (MQ=255) ccACAATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCt < 1:114576/71‑1 (MQ=255) ccACAATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCt < 1:1576266/71‑1 (MQ=255) ccACAATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCt < 1:2955268/71‑1 (MQ=255) ccACAATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCt < 1:1489536/71‑1 (MQ=255) ccACAATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCt < 1:4033958/71‑1 (MQ=255) ccACAATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCt < 1:2532838/71‑1 (MQ=255) ccACAATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCt > 1:1858877/1‑71 (MQ=255) ctcaATCTTCGCCCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCt < 1:2583399/69‑1 (MQ=255) cacaATCTTCGCCCCCCCAGGCTGCCGCTTTCAtt > 1:3895876/1‑35 (MQ=38) cacaATCTTCGCCCCCCCAGGCTGCCGCTTTCATTATGACGGAAg > 1:4306468/1‑45 (MQ=255) cacaATCTTCGCCCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCt > 1:3115098/1‑71 (MQ=255) cacaATCTTCGCCCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCt > 1:3088546/1‑71 (MQ=255) cacaATCTTCGCCCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCt < 1:410250/71‑1 (MQ=255) cacaATCTTCGCCCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCt > 1:4263319/1‑71 (MQ=255) cacaATCTTCGCCCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCt < 1:2495885/71‑1 (MQ=255) cacaATCTTCGCCCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCt > 1:2292888/1‑71 (MQ=255) cacaATCTTCGCCCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCt > 1:2229171/1‑71 (MQ=255) cacaATCTTCGCCCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCt > 1:1624673/1‑71 (MQ=255) cacaATCTTCGCCCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCt > 1:1588028/1‑71 (MQ=255) cacaATCTTCGCCCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCt > 1:557525/1‑71 (MQ=255) cacaATCTTCGCCCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCt > 1:107000/1‑71 (MQ=255) caATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGgaaga > 1:3668149/1‑64 (MQ=255) caATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGATGAAGCTGAc > 1:944413/1‑71 (MQ=255) caATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCt < 1:1385020/68‑1 (MQ=255) caATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCTGAc > 1:3969128/1‑71 (MQ=255) caATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCTGAc > 1:325196/1‑71 (MQ=255) caATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCTGAc > 1:938725/1‑71 (MQ=255) caATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCTGAc > 1:808942/1‑71 (MQ=255) caATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCTGAc > 1:690247/1‑71 (MQ=255) aaTCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCTGACt < 1:2991419/71‑1 (MQ=255) aaTCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCTGACt < 1:3957592/71‑1 (MQ=255) aaTCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCTGACt < 1:758111/71‑1 (MQ=255) aaTCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCTGACt < 1:622433/71‑1 (MQ=255) aTCTTCGCCCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCTGACt < 1:1004718/71‑1 (MQ=255) aTCTTCGCCCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCTGACt < 1:2130239/71‑1 (MQ=255) aTCTTCGCCCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCTGACt < 1:2719907/71‑1 (MQ=255) aTCTTCGCCCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCTGACt < 1:4010776/71‑1 (MQ=255) aTCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCTGACtt > 1:2330446/1‑71 (MQ=255) aTCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCTGACtt > 1:2286502/1‑71 (MQ=255) aTCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCTGACtt > 1:137426/1‑71 (MQ=255) aTCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCTGACt > 1:4238271/1‑70 (MQ=255) aTCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCTGACt > 1:401863/1‑70 (MQ=255) tCTTCGCCCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCTGACt < 1:2837394/70‑1 (MQ=255) cTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCTGACTTc > 1:2508488/1‑70 (MQ=255) cG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCTGACTTc < 1:2512932/67‑1 (MQ=255) cG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCTGACTTc < 1:2584361/67‑1 (MQ=255) cG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCTGACTTc < 1:1152315/67‑1 (MQ=255) cG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCTGACTTc < 1:2856723/67‑1 (MQ=255) | CGGCGGTACGATAAAATTCGTCAAGCAGCAGTATATGCTGGGTGGAGCCACAATCTTCG‑CCCCCCAGGCTGCCGCTTTCATTATGACGGAAGCGGTTTTCATCAATCAGGAAGAAGCTGACTTC > minE/2406890‑2407013 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |