Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,848,210 | 0 | C | G | 37.0% | 37.9 / 14.7 | 27 | G344G (GGG→GGC) | hycC | hydrogenase 3, membrane subunit |
Reads supporting (aligned to +/- strand): ref base C (4/13); new base G (10/0); total (14/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.53e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCCAATACCACCGAGTTTTTCGATATCGCGATGACCGGTACGGAACCAGACGCTCCCCGCCCCGAGGAACAGTACGCTTTTGAACAGGCTATGGTTAAGCAGATGGTACAGACCACCGA > NC_000913/2848156‑2848274 | gCCAATACCACCGAGTTTTTCGATATCGCGATGACCGGTACGGAACCAGACGCTc > 2:809526/1‑55 (MQ=255) gCCAATACCACCGAGTTTTTCGATATCGCGATGACCGGTACGGAACCAGACGCTc < 1:809524/55‑1 (MQ=255) tACCACCGAGTTTTTCGATATCGCGATGACCGGTACGGAACCAGACGCTCCCCGCCCCGAGGAACAGt < 2:72673/68‑1 (MQ=255) tACCACCGAGTTTTTCGATATCGCGATGACCGGTACGGAACCAGACGCTCCCCGCCCCGAGGAACAGt < 2:122369/68‑1 (MQ=255) tACCACCGAGTTTTTCGATATCGCGATGACCGGTACGGAACCAGACGCTCCCCGCCCCGAGGAACAGt < 1:794917/68‑1 (MQ=255) ccaccGAGTTTTTCGATATCGCGATGACCGGTACGGAACCAGACGCTCCCCGCCCCGAGGa > 1:114317/1‑61 (MQ=255) ccaccGAGTTTTTCGATATCGCGATGACCGGTACGGAACCAGACGCTCCCCGCCCCGAGGa < 2:114317/61‑1 (MQ=255) ccGAGTTTTTCGATATCGCGATGACCGGCACGGAACCAGACGCTCCCCGCCCCGAGGAACAGTACGCt < 1:45805/68‑1 (MQ=255) cGAGTTTTTCGATATCGCGATGACCGGTACGGAACCAGACGCTGCCCGCCCCGAGGAACAGTACGCtt > 1:645714/1‑68 (MQ=255) gTTTTTCGATATCGCGATGACCGGTACGGAACCAGACGCTGCCCGCCCCGAGGAACAGTACGCTTTTg > 2:716031/1‑68 (MQ=255) tCGATATCGCGATGACCGGTACGGAACCAGACGCTCCCCGCCCCGAGGAACAGTACGCTTTTGAACAg < 2:639645/68‑1 (MQ=255) cGATATCGCGATGACCGGTACGGAACCAGACGCTCCCCGCCCCGAGGAACAGTACGCTTTTGAACAgg < 1:533239/68‑1 (MQ=255) cGATATCGCGATGACCGGTACGGAACCAGACGCTCCCCGCCCCGAGGAACAGTACGCTTTTGAACAgg < 1:771638/68‑1 (MQ=255) ccGGTACGGAACCAGACGCTGCCCGCCCCGAGGAACAGTACGCTTTTGAACAGGCTATGGTTAAGCAg > 1:53118/1‑68 (MQ=255) ccGGTACGGAACCAGACGCTGCCCGCCCCGAGGAACAGTACGCTTTTGAACAGGCTATGGTTAAGCAg > 2:337899/1‑68 (MQ=255) ccGGTACGGAACCAGACGCTCCCCGCCCCGAGGAACAGTACGCTTTTGAACAGGCTATGGTTAAGCAg < 1:137972/68‑1 (MQ=255) aCGGAACCAGACGCTCCCCGCCCCGAGGAACAGTACGCTTTTGAACAGGCTATGGTTAAGCAGATGGt < 1:40826/68‑1 (MQ=255) aCGGAACCAGACGCTCCCCGCCCCGAGGAACAGTACGCTTTTGAACAGGCTATGGTTAAGCAGATGGt < 1:24090/68‑1 (MQ=255) aCGGAACCAGACGCTCCCCGCCCCGAGGAACAGTACGCTTTTGAACAGGCTATGGTTAAGCAGATGGt < 2:762934/68‑1 (MQ=255) aCGGAACAAGACGCTGCCCGACCCGAGGAACAGTACGCATTTGAACAGGCAATGGTTAAGCAGAAGGt > 1:561961/1‑68 (MQ=255) aCCAGCCGCAGCCCACCCCGAGGAACAGAACGCATATGAACAGGCAATGGATTAGCAGATCGAACAGa > 1:798304/1‑68 (MQ=255) aCCAGACGCTGCCCGCCCCGAGGAACAGTACGCTTTTGGACAGGCTATGGTTAAGCAGATGGTACAGa > 1:165902/1‑68 (MQ=255) aCCAGACGCTCCCCGCCCCGAGGAACAGTACGCTTTTGAACAGGCTATg < 1:622334/49‑1 (MQ=255) aCCAGACGCTCCCCGCCCCGAGGAACAGTACGCTTTTGAACAGGCTATg > 2:622334/1‑49 (MQ=255) aCCAGACGCTCCCCGCCCCGAGGAACAGTACGCTTTTGAACAGGCTATGGTTAAGCAGATGGTACAgc > 2:233045/1‑67 (MQ=255) aCCAGACGCTCCCCGCCCCGAGGAACAGTACGCTTTTGAACAGGCTATGGTTAAGCAGATGGTACAgc > 2:713416/1‑67 (MQ=255) aCGCTGCCCGCCCCGAGGAACAGTACGCTTTTGAACAGGCTATGGTTAAGCAGATGGTACAGaccacc > 2:530041/1‑68 (MQ=255) aCGCTGCCCGCCCCGAGGAACAGTACGCTTTTGAACAGGCTATGGTTAAGCAGATGGTACAGaccacc > 2:399561/1‑68 (MQ=255) gCTGCCCGCCCCGAGGAACAGTACGCTTTTGAACAGGCTATGGTTAAGCAGATGGTACAGCCCACCgc > 1:105467/1‑67 (MQ=255) | GCCAATACCACCGAGTTTTTCGATATCGCGATGACCGGTACGGAACCAGACGCTCCCCGCCCCGAGGAACAGTACGCTTTTGAACAGGCTATGGTTAAGCAGATGGTACAGACCACCGA > NC_000913/2848156‑2848274 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |