Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,070,462 | 0 | T | C | 23.8% | 127.7 / 42.0 | 63 | intergenic (‑61/+53) | yihR/yihS | putative aldose 1‑epimerase YihR/sulfoquinovose isomerase |
Reads supporting (aligned to +/- strand): ref base T (30/18); new base C (9/6); total (39/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.87e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTT > NC_000913/4070413‑4070506 | tGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCt < 1:3485448/50‑1 (MQ=255) tGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCt > 1:3486421/1‑50 (MQ=255) tGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCt > 1:4335674/1‑50 (MQ=255) aGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCttt > 1:2450375/1‑50 (MQ=255) aGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCttt < 1:497065/50‑1 (MQ=255) aGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCtt < 1:1119647/49‑1 (MQ=255) aGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCtt < 1:4985142/49‑1 (MQ=255) ggctgcaagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTg > 1:1319998‑M2/29‑50 (MQ=255) ggCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg > 1:5581563/1‑50 (MQ=255) ggCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg < 1:5278314/50‑1 (MQ=255) ggCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTg < 1:994994/50‑1 (MQ=255) gctgcaagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTgg > 1:3514537‑M2/28‑50 (MQ=255) gCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCtttag > 1:83249/1‑48 (MQ=255) gCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgg < 1:2035663/50‑1 (MQ=255) gCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTgg < 1:3687370/50‑1 (MQ=255) tGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg > 1:1557791/1‑49 (MQ=255) tGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTggg > 1:3919760/1‑49 (MQ=255) caagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTggg > 1:3445342‑M2/24‑47 (MQ=255) caagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTggg > 1:258996‑M2/24‑47 (MQ=255) tGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGgaga < 1:2416291/50‑1 (MQ=255) aagttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGga > 1:3711977‑M2/23‑48 (MQ=255) agttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGagag > 1:1957106‑M2/22‑50 (MQ=255) gttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAgg > 1:6473029‑M2/21‑50 (MQ=255) ttgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAggg < 1:1193157‑M2/31‑1 (MQ=255) tgcagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGt > 1:2736283‑M2/19‑50 (MQ=255) aaaCCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAggg < 1:2201200/49‑1 (MQ=255) aaaCCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAggg < 1:2767262/49‑1 (MQ=255) aaaCCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGt > 1:3775917/1‑50 (MQ=255) cagaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGtt < 1:6581203‑M2/33‑1 (MQ=255) aCCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGt > 1:2970277/1‑48 (MQ=255) aCCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTa < 1:5222571/50‑1 (MQ=255) cccAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAg < 1:4288368/50‑1 (MQ=255) agaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAg < 1:265911‑M2/35‑1 (MQ=255) gaccgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAgg < 1:6096446‑M2/36‑1 (MQ=255) ccAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg > 1:4325891/1‑50 (MQ=255) cAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAggg < 1:5098670/50‑1 (MQ=255) cAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAggg < 1:2116434/50‑1 (MQ=255) cAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAggg < 1:2507810/50‑1 (MQ=255) accgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAggg < 1:3650776‑M2/37‑1 (MQ=255) accgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAggg < 1:4497640‑M2/37‑1 (MQ=255) accgttacgagatGGACGGTCCCCTCGCCCCCTTGGGGAGAGGGTTAggg > 1:4328052‑M2/14‑50 (MQ=255) aaaaCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGgag < 1:3100375/50‑1 (MQ=255) aaaaCTGCACAGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGgag < 1:4978853/50‑1 (MQ=255) aaaCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGgaga < 1:1733066/50‑1 (MQ=255) aaaCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGgaga < 1:4647474/50‑1 (MQ=255) aCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAgg < 1:2470685/50‑1 (MQ=255) aCTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAgg > 1:3218284/1‑50 (MQ=255) cTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAggg > 1:3975081/1‑50 (MQ=255) cTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAggg > 1:3625391/1‑50 (MQ=255) cACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGa > 1:3248125/1‑49 (MQ=255) ccGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGaaaa > 1:4148579/1‑50 (MQ=255) ggACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAcc > 1:3503128/1‑50 (MQ=255) ggACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAcc > 1:1867221/1‑50 (MQ=255) gACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAAccc > 1:4096055/1‑50 (MQ=255) ggTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCa > 1:1933812/1‑50 (MQ=255) ggTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCa > 1:4180839/1‑50 (MQ=255) ggTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCa > 1:5987693/1‑50 (MQ=255) tCCCCTCGCCCGTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCAcc < 1:6122447/50‑1 (MQ=255) tCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCAcc > 1:6686380/1‑50 (MQ=255) ccccTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCAccc > 1:2167312/1‑50 (MQ=255) cTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAt > 1:4595174/1‑50 (MQ=255) cTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAt > 1:2510982/1‑50 (MQ=255) cTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAt > 1:5640823/1‑50 (MQ=255) gCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtttt > 1:281061/1‑50 (MQ=255) gCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtttt > 1:1454265/1‑50 (MQ=255) gCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtttt > 1:2608200/1‑50 (MQ=255) ccccTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtttt > 1:4384351/1‑49 (MQ=255) ccTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtttt > 1:6528431/1‑47 (MQ=255) ccTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtttt > 1:2604943/1‑47 (MQ=255) ccTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtttt > 1:3655193/1‑47 (MQ=255) ccTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTAtttt > 1:2830837/1‑47 (MQ=255) | TGAGGCTGATGTCAAACCCAAAAACTGCACCGGACGGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAAAACCCGCACCCTATTTT > NC_000913/4070413‑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: ‑ |