| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,070,462 | 0 | T | C | 45.3% | 21.7 / 114.6 | 75 | intergenic (‑61/+53) | yihR/yihS | putative aldose 1‑epimerase YihR/sulfoquinovose isomerase |
| Reads supporting (aligned to +/- strand): ref base T (25/16); new base C (20/14); total (45/30) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.83e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTT > NC_000913/4070414‑4070506 | gAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCtt > 1:2915491/1‑50 (MQ=255)gAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCtt > 1:2173589/1‑50 (MQ=255)gAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCtt > 1:6396046/1‑50 (MQ=255) gggctgcaagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCctt < 1:1200794‑M2/21‑1 (MQ=255) gggctgcaagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCctt < 1:4584572‑M2/21‑1 (MQ=255) gggctgcaagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCctt > 1:48043‑M2/30‑50 (MQ=255) aGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCttt < 1:5398001/50‑1 (MQ=255) aGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCtt < 1:3265364/49‑1 (MQ=255) ggctgcaagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTg < 1:7359314‑M2/22‑1 (MQ=255) ggCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg > 1:689571/1‑50 (MQ=255) gctgcaagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTgg > 1:772411‑M2/28‑50 (MQ=255) gCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgg < 1:2471020/50‑1 (MQ=255) ctgcaagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTggg > 1:4984802‑M2/27‑50 (MQ=255) ctgcaagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTggg < 1:1545630‑M2/24‑1 (MQ=255) tgcaagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTggg > 1:5237263‑M2/26‑49 (MQ=255) tgcaagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTggg > 1:290395‑M2/26‑49 (MQ=255) gcaagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTggg > 1:3927504‑M2/25‑48 (MQ=255) caagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTggg > 1:5195936‑M2/24‑47 (MQ=255) aTGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgggg < 1:6708372/48‑1 (MQ=255) aTGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGgag < 1:329462/50‑1 (MQ=255) aTGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGgag < 1:2687180/50‑1 (MQ=255) aTGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGga < 1:5646763/49‑1 (MQ=255) tGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGga > 1:4500529/1‑48 (MQ=255) tGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGga > 1:4379483/1‑48 (MQ=255) tGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGga > 1:7099450/1‑48 (MQ=255) tGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGga > 1:4299336/1‑48 (MQ=255) gTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGagag > 1:5523629/1‑50 (MQ=255) agttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGagag < 1:1363769‑M2/29‑1 (MQ=255) tCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGagag < 1:7000243/49‑1 (MQ=255) tCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGagag < 1:1970908/49‑1 (MQ=255) tCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAgg < 1:6647518/50‑1 (MQ=255) tCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAgg < 1:646855/50‑1 (MQ=255) tCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAgg < 1:6138075/50‑1 (MQ=255) gttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAgg < 1:1636607‑M2/30‑1 (MQ=255) gttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAgg > 1:3958873‑M2/21‑50 (MQ=255) ttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAggg < 1:3753799‑M2/31‑1 (MQ=255) aaaCCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGt > 1:305440/1‑50 (MQ=255) aaaCCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGt < 1:6261206/50‑1 (MQ=255) aaaCCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGt > 1:4669719/1‑50 (MQ=255) gcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGt < 1:1192533‑M2/32‑1 (MQ=255) gcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGt < 1:3942551‑M2/32‑1 (MQ=255) gcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGt < 1:3441159‑M2/32‑1 (MQ=255) gcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGt > 1:265250‑M2/18‑49 (MQ=255) aaCCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGtt < 1:5015094/50‑1 (MQ=255) cagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGtt < 1:2483290‑M2/33‑1 (MQ=255) cagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGtt < 1:7148696‑M2/33‑1 (MQ=255) cagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGt > 1:3962845‑M2/17‑48 (MQ=255) aCCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGt > 1:4258510/1‑48 (MQ=255) aCCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTa < 1:1449132/50‑1 (MQ=255) agaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAg > 1:1663198‑M2/16‑50 (MQ=255) gaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAgg < 1:1306253‑M2/36‑1 (MQ=255) gaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAgg > 1:1646681‑M2/15‑50 (MQ=255) cAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAggg < 1:2195225/50‑1 (MQ=255) accgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAggg < 1:5685784‑M2/37‑1 (MQ=255) ccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTagggt > 1:6038259‑M2/13‑49 (MQ=255) ccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTagggt > 1:5883367‑M2/13‑49 (MQ=255) ccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTagggt > 1:1359252‑M2/13‑49 (MQ=255) ccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTagggt > 1:1027140‑M2/13‑49 (MQ=255) cgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTtagggtg < 1:3737321‑M2/39‑3 (MQ=255) cgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTtagggtg < 1:6075393‑M2/39‑3 (MQ=255) cgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTagggt > 1:650048‑M2/12‑48 (MQ=255) aaaaCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg > 1:3016654/1‑47 (MQ=255) gttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTagggt > 1:1807116‑M2/11‑47 (MQ=255) gttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTagggt > 1:3066625‑M2/11‑47 (MQ=255) gttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTagggt > 1:7280603‑M2/11‑47 (MQ=255) gttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTagggt > 1:4513753‑M2/11‑47 (MQ=255) aaCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGagag < 1:4534016/50‑1 (MQ=255) aCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAgg < 1:6026078/50‑1 (MQ=255) aCCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGaa > 1:2683026/1‑49 (MQ=255) ccGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGaaaa > 1:4756129/1‑50 (MQ=255) ggACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAcc > 1:6534813/1‑50 (MQ=255) ggACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAcc > 1:5742893/1‑50 (MQ=255) gACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAccc > 1:7369775/1‑50 (MQ=255) gACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAccc > 1:7002579/1‑50 (MQ=255) cccTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCAccc > 1:5989078/1‑49 (MQ=255) cTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAt > 1:2944668/1‑50 (MQ=255) tCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtt > 1:3384250/1‑50 (MQ=255) gCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtttt > 1:1950765/1‑50 (MQ=255) ccccTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtttt > 1:2304627/1‑49 (MQ=255) ccccTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtttt > 1:692534/1‑49 (MQ=255) cccTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtttt > 1:3260395/1‑48 (MQ=255) ccTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtttt > 1:2119159/1‑47 (MQ=255) cTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtttt > 1:5523829/1‑46 (MQ=255) | GAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTT > NC_000913/4070414‑4070506 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |