Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 645,910 | 0 | T | G | 53.8% | 9.9 / 13.9 | 13 | T13P (ACC→CCC) | citE | citrate (pro‑3S)‑lyase subunit beta |
Reads supporting (aligned to +/- strand): ref base T (0/6); new base G (7/0); total (7/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.83e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGC > NZ_CP009273/645848‑645990 | gATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGGGGTTTTTGTTGTTTCAGCGGAGc > 2:199345/1‑90 (MQ=255) gTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGGGGTTTTTGTTGTTTGAGCGAAGCGGGAAACAt > 2:270307/1‑90 (MQ=255) gTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGGTTTTCGTTGTTGCAGCGAAGCGGGAAACAt > 1:139504/1‑90 (MQ=255) tGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCa < 1:444395/90‑1 (MQ=255) tCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGGGGGGTTTTTGGTTTGGCAGCGAAACGGAAAACAATGGCAAacctc > 1:469753/1‑90 (MQ=255) gcgGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTccc < 1:81248/90‑1 (MQ=255) gcgGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGGGCGGGGTTTTCGTTGTTGCAGCGAAGCGGGAATCATTTGCAATCCTccc > 1:133530/1‑90 (MQ=255) aGGCACAAACAACATGCTGCGGCGGGGGGGAGGTTTTCGTTTTTGGAGCGAAGCGGGAAACATTGGCAAACCTCCCCTGGCAGAGCCggg > 1:330890/1‑90 (MQ=255) gCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGAt < 2:283325/90‑1 (MQ=255) acaacaTGCTGCGGCGGGTGCGTGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCt < 2:326545/90‑1 (MQ=255) acaacaTGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCt < 2:362559/90‑1 (MQ=255) aTGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCACATGGCAGAGCCGGGATACCGCTGgcgc < 2:133530/90‑1 (MQ=255) tgcGGCGGGGGCGGGGTTTTCGTTGTTGTAGCCGAACGGAAAAAATTGGCAATCCTCCCCTGGCAGAGGCGGGGTAGCGCCGGCGCGGGc > 1:362559/1‑90 (MQ=255) | GATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGGCACAAACAACATGCTGCGGCGGGTGCGAGTTTTACGTTGTTGCAGCGAAGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGC > NZ_CP009273/645848‑645990 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |