Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,980,951 | 1 | . | G | 34.4% | 31.8 / 59.8 | 32 | coding (65/1386 nt) | yifK | putative transporter YifK |
Reads supporting (aligned to +/- strand): ref base . (14/7); new base G (4/7); total (18/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.42e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.23e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCT > NC_000913/3980902‑3980999 | ccAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCt < 1:2997252/50‑1 (MQ=255) ccAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCt < 1:4165163/50‑1 (MQ=255) agagCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑gg > 1:2813980/1‑50 (MQ=255) gagCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑ggg < 1:111567/50‑1 (MQ=255) gagCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑ggg > 1:2066642/1‑50 (MQ=255) agCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑gggg > 1:26040/1‑50 (MQ=255) cTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑gggggg > 1:3077550/1‑50 (MQ=255) aCAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑ggggggg > 1:710525/1‑48 (MQ=255) cAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑ggggggg < 1:2862790/48‑2 (MQ=255) aGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑GGGGGGCAcc < 1:4165347/50‑1 (MQ=255) cGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑GGGGGGCACCAt > 1:2996537/1‑50 (MQ=255) gTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑GGGGGGCACCAtt > 1:1715260/1‑50 (MQ=255) gTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑GGGGGGCACCAtt < 1:2124228/50‑1 (MQ=255) tGGGCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCAtt < 1:578035/50‑1 (MQ=255) ggCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTgg < 1:1233190/50‑1 (MQ=255) ggCTGGAAGCTCGACATATCGAACTCATCGCCCT‑GGGGGGCACCATTgg < 1:1601902/49‑1 (MQ=255) gCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGc < 1:2904420/50‑1 (MQ=255) gCTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGc < 1:3325132/50‑1 (MQ=255) cTGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCg > 1:677988/1‑50 (MQ=255) tGGAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGt > 1:3606621/1‑50 (MQ=255) tGGAAGCTCGACATATCGAACTCATCGCCCT‑GGGGGGCACCATTGGCGTc > 1:1931126/1‑50 (MQ=255) tGGAAGCTCGACATATCGAACTCATCGCCCT‑GGGGGGCACCATTGGCGTc > 1:2664445/1‑50 (MQ=255) gAAGCTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTc < 1:1230975/49‑1 (MQ=255) gAAGCTCGACATATCGAACTCATCGCCCT‑GGGGGGCACCATTGGCGTCgg < 1:559491/50‑1 (MQ=255) gCTCGACATATCGAACTCATCGCCCT‑GGGGGGCACCATTGGCGTCGGCCt > 1:2521297/1‑50 (MQ=255) cTCGACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCt < 1:3834438/50‑1 (MQ=255) cTCGACATATCGAACTCATCGCCCT‑GGGGGGCACCATTGGCGTCGGCCTg > 1:3855199/1‑50 (MQ=255) gACATATCGAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGtt > 1:2217167/1‑50 (MQ=255) gACATATCGAACTCATCGCCCT‑GGGGGGCACCATTGGCGTCGGCCTGtt > 1:3419080/1‑49 (MQ=255) gACATATCGAACTCATCGCCCT‑GGGGGGCACCATTGGCGTCGGCCTGtt > 1:647756/1‑49 (MQ=255) tatCGAACTCATCGCCCT‑GGGGGGCACCATTGGCGTCGGCCTGTTTATgg > 1:3713543/1‑50 (MQ=255) cGAACTCATCGCCCT‑GGGGGGCACCATTGGCGTCGGCCTGTTTATggggg > 1:3409059/1‑50 (MQ=255) gAACTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATgggg < 1:3571751/49‑1 (MQ=39) gAACTCATCGCCCT‑GGGGGGCACCATTGGCGTCGGCCTGTTTATggggg < 1:3378793/49‑1 (MQ=255) aCTCATCGCCCTGGGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgcc > 1:3849970/1‑50 (MQ=255) cATCGCCCT‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGgccgc > 1:1579990/1‑48 (MQ=255) tCGCCCT‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGt < 1:1730887/50‑1 (MQ=255) cccT‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTAcc < 1:1661033/50‑1 (MQ=255) cccT‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTAcc > 1:843673/1‑50 (MQ=255) ccT‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTAccc > 1:2948622/1‑50 (MQ=255) ccT‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTAccc > 1:3473267/1‑50 (MQ=255) t‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCt > 1:1298297/1‑49 (MQ=255) | CCAGAGCTACAGCGTGGGCTGGAAGCTCGACATATCGAACTCATCGCCCT‑GGGGGGCACCATTGGCGTCGGCCTGTTTATGGGGGCCGCCAGTACCCT > NC_000913/3980902‑3980999 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |