| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NZ_CP009273 | 258,435 | A→G | I408V (ATC→GTC) | proA → | glutamate‑5‑semialdehyde dehydrogenase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 258,435 | 0 | A | G | 83.9% | 73.1 / 6.7 | 31 | I408V (ATC→GTC) | proA | glutamate‑5‑semialdehyde dehydrogenase |
| Reads supporting (aligned to +/- strand): ref base A (1/4); new base G (14/12); total (15/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.33e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01 | |||||||||||
GGGTGCGGAAGTGGCGGTAAGCACACAAAAACTCCACGCGCGTGGCCCAATGGGGCTGGAAGCACTGACCACTTACAAGTGGATCGGCATTGGTGATTACACCATTCGTGCGTAAATAAAACCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACG > NZ_CP009273/258353‑258518 | gggTGCGGAAGTGGCGGTAAGCACACAAAAACTCCACGCGCGTGGCCCAATGGGGCTGGAAGCACTGACCACTTACAAGTGGGTCGGCAt < 2:354710/90‑1 (MQ=255)gggTGCGGAAGTGGCGGTAAGCACACAAAAACTCCACGCGCGTGGCCCAATGGGGCTGGAAGCACTGACCACTTACAAGTGGGTCGGCAt < 1:993683/90‑1 (MQ=255) ggAAGTGGCGGTAAGCACCCAAAACCTCCACGCGGGTGGCCCAATGGGGCTGGAAGCACTGACCCCTTACAAGTCGGTCGGAATGGGTGa < 2:631910/90‑1 (MQ=255) aaGTGGCGGTAAGCACACAAAAACTCCACGCGCGTGGCCCAATGGGGCTGGAAGCACTGACCACTTACAAGTGGGTCGGCATTGGTGAtt < 1:672454/90‑1 (MQ=255) ggCGGTAAGCACACAAAAACTCCACGCGCGTGGCCCAATGGGGGTGGAAGCACTGACCACTTACAAGTGGGTCGGCATTGGGGATTacac > 1:573723/1‑90 (MQ=255) ggCGGTAAGCACACAAAAACTCCACGCGCGTGGCCCAATGGGGCTGGAAGCACTGACCACTTACAAGTGGGTCGGCATTGGTGATTacac < 1:1060569/90‑1 (MQ=255) tAAGCACACAAAAACTCCACGCGCGTGGCCCAATGGGGCTGGAAGCACTGACCACTTACAAATGGGTCGGCATTGGTGATTACACCATTc > 1:11413/1‑90 (MQ=255) gCACACAAAAACTCCACGCGCGTGGCCCAATGGGGCTGGAAGCACTGACCACTTACAAGTGGGTCGGCATTGGTGATTACACCATTCGTg > 1:150586/1‑90 (MQ=255) acacaAAAACTCCACGCGCGTGTCCCAATGGGGCTGGAAGCACTGACCACTTACAAGTGGGTCGGCATTGGTGATTACACCATTCGTGCg > 2:234870/1‑90 (MQ=255) acacaAAAACTCCACGCGCGTGGCCCAATGGGGCTGGAAGCACTGACCACTTACAAGTGGGTCGGCATTGGTGATTACACCATTCGTGCg < 1:770406/90‑1 (MQ=255) acacaAAAACTCCACGCGCGTGGCCCAATGGGGCTGGAAGCACTGACCACTTACAAGTGGATCGGCATTGGTGATTACACCATTCGTGCg < 2:375331/90‑1 (MQ=255) acacaAAAACTCCACGCGCGTGGCCCAATGGGGCTGGAAGCACTGACCACTTACAAGTGGATCGGCATTGGTGATTACACCATTCGTGCg < 1:61497/90‑1 (MQ=255) cacaAAAACTCCACGCGCGTGGCCCAATGGGGCTGGAAGCACTGACCACTTACAAGTGGGTCGGCATTGGTGATTACACCATTCGTGCGt < 1:506394/90‑1 (MQ=255) cTCCACGCGCGTGGCCCAATGGGGCTGGAAGCACTGACCACTTACAAGTGGATCGGCATTGGTGATTACACCATTCGTGCGTAAATaaaa < 2:974374/90‑1 (MQ=255) ccACGCGCGTGGCCCAATGGGGCTGGAAGCACTGACCACTTACAAGTGGGTCGGCATTGGTGATTACACCATTCGTGCGTAAATAAAAcc > 1:1043687/1‑90 (MQ=255) gcgTGGCCCAATGGGGCTGGAAGCACTGACCACTTACAAGTGGATCGGCAt < 2:1032288/51‑1 (MQ=255) gcgTGGCCCAATGGGGCTGGAAGCACTGACCACTTACAAGTGGATCGGCAt > 1:1032288/1‑51 (MQ=255) gTGGCCCAATGGGGCTGGAAGCACTGACCACTTACAAGTGGGTCGGCATTGGTGATTACACCATTCGTGCGTAAATAAAACCGGGTGATg > 2:1043101/1‑90 (MQ=255) aaTGGGGCTGGAAGCACTGACCACTTACAAGTGGGTCGGCATTGGTGATTACACCATTCGTGCGTAAATAAAACCGGGTGATGCAAAAGt > 1:587356/1‑90 (MQ=255) ggggCTGGAAGCACTGACCACTTACAAGTGGGTCGGCATTGGTGATTACACCATTCGTGCGTAAATAAAACCGGGTGATGCAAAAGTAGc < 2:587356/90‑1 (MQ=255) gCTGGAAGCACTGACCACTTACAAGTGGGTCGGCATTGGTGATTACACCATTCGTGCGTAAATAAAACCGGGt > 1:910279/1‑73 (MQ=255) gCTGGAAGCACTGACCACTTACAAGTGGGTCGGCATTGGTGATTACACCATTCGTGCGTAAATAAAACCGGGt < 2:910279/73‑1 (MQ=255) gCTGGAAGCACTGACCACTTACAAGTGGGTCGGCATTGGTGATTACACCATTCGTGCGTAAATAAAACCGGGTGATGCAAAAGTAGCCAt > 1:181046/1‑90 (MQ=255) aaGCACTGACCACTTACAAGTGGGTCGGCATTGGTGATTACACCATTCGTGCGTAAATAAAACCGGGTGATGCAAAAGTAGCCATTTGAt > 1:899354/1‑90 (MQ=255) aaGCACTGACCACTTACAAGTGGGTCGGCATTGGTGATTACACCATTCGTGCGTAAATAAAACCGGGTGATGCAAAAGTAGCCATTTGAt < 1:584541/90‑1 (MQ=255) aCTGACCACTTACAAGTGGGTCGGCATTGGTGATTACACCATTCGTGCGTAAATAAAACCGGGGGATGCAAAAGTAGCCATTTGATTcac > 2:734100/1‑90 (MQ=255) tGACCACTTACAAGTGGGTCGGCATTGGTGATTACACCATTCGTGCGTAAATAAAACCGGGTGATGCAAAAGTAGCCATTTGATTCACaa < 2:1028591/90‑1 (MQ=255) ccACTTACAAGTGGGTCGGCATTGGTGATTACACCATTCGTGCGTAAATAAAACCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGc > 2:511948/1‑90 (MQ=255) tACAAGTGGGTCGGCATTGGTGATTACACCATTCGTGCGTAAATAAAACCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTg < 1:511948/90‑1 (MQ=255) aCAAGTGGGTCGGCATTGGTGATTACACCATTCGTGCGTAAATAAAACCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGa > 2:323062/1‑90 (MQ=255) aaGTGGGTCGGCATTGGTGATTACACCATTCGTGCGTAAATAAAACCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACg > 1:413896/1‑90 (MQ=255) | GGGTGCGGAAGTGGCGGTAAGCACACAAAAACTCCACGCGCGTGGCCCAATGGGGCTGGAAGCACTGACCACTTACAAGTGGATCGGCATTGGTGATTACACCATTCGTGCGTAAATAAAACCGGGTGATGCAAAAGTAGCCATTTGATTCACAAGGCCATTGACG > NZ_CP009273/258353‑258518 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |