Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 1,094,108 | A→G | E256G (GAA→GGA) | ghrA → | glyoxylate/hydroxypyruvate reductase GhrA |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,094,108 | 0 | A | G | 92.3% | 33.9 / ‑3.2 | 13 | E256G (GAA→GGA) | ghrA | glyoxylate/hydroxypyruvate reductase GhrA |
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (6/6); total (7/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.77e-01 |
GTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTTAAAGGCGCAATGTTGGATGTTTTTAATCGTGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTACCCGTCCCGCTGAA > NZ_CP009273/1094026‑1094193 | gtTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTTAAAGGCGCAATGTTGGATGTTTTTAATCGTGGACCCTTa < 1:110744/90‑1 (MQ=255) tgtgGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTTAAAGGCGCAATGTTGGATGTTTTTAATCGTGGACCCTTACCGCCTGa > 2:240033/1‑90 (MQ=255) tgtgGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTTAAAGGCGCAATGTTGGATGTTTTTAATCGTGGACCCTTACCGCCTGa > 2:51601/1‑90 (MQ=255) tgGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTTAAAGGCGCAATGTTGGATGTTTTTAATCGTGGACCCTTACCGCCTGaaa > 2:244874/1‑90 (MQ=255) aaGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTTAAAGGCGCAATGTTGGATGTTTTTAATCGTGGACCCTTACCGCCTGAAAGTc < 1:332102/90‑1 (MQ=255) gACCTGCTCGCGGCGCTGGATAGCGGCAAAGTTAAAGGCGCAATGTTGGATGTTTTTAATCGTGGACCCTTACCGCCTGAAAGTCCGctc < 1:55855/90‑1 (MQ=255) tGGATAGCGGCAAAGTTAAAGGCGCAATGTTGGATGTTTTTAATCGTGGACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCAcgcg < 1:248368/90‑1 (MQ=255) aGCGGCAAAGTTAAAGGCGCAATGTTGGATGTTTTTAATCGTGGACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACg < 2:292720/90‑1 (MQ=255) aaaGTTAAAGGCGCAATGTTGGATGTTTTTAATCGTGGACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAaca > 2:387096/1‑90 (MQ=255) aaaGTTAAAGGCGCAATGTTGGATGTTTTTAATCGTGGACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAaca > 2:84346/1‑90 (MQ=255) tGTTGGATGTTTTTAATCGTGGACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGa < 2:77788/90‑1 (MQ=255) tGGATGTTTTTAATCGTGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTa > 1:144439/1‑90 (MQ=255) cGTGGACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTACCCGTCCCGCTGaa > 2:342050/1‑90 (MQ=255) | GTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTTAAAGGCGCAATGTTGGATGTTTTTAATCGTGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTACCCGTCCCGCTGAA > NZ_CP009273/1094026‑1094193 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |