Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 3,579,168 | T→C | E339G (GAA→GGA) | ggt ← | gamma‑glutamyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,579,168 | 0 | T | C | 90.9% | 82.9 / 2.5 | 33 | E339G (GAA→GGA) | ggt | gamma‑glutamyltransferase |
Reads supporting (aligned to +/- strand): ref base T (1/2); new base C (15/15); total (16/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01 |
CAATTTGATCGGCAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATAT‑TCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGCCGTATTTCTTC > NZ_CP009273/3579082‑3579248 | cAATTTGATCGGCAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATAT‑CCCg < 1:119847/90‑1 (MQ=255) gATCGGCAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATAT‑CCCGAGCGGt > 1:644871/1‑90 (MQ=255) agcGGCAATAGATTTGGCATCGGCTTTACTGGGCCGCGCCTGCCACGGGACTATGACAAAATCCGGGTCGCCAAGATAT‑CCCGAGCGGTc < 2:783600/88‑1 (MQ=255) accGGCAATAGATTTGGCATAGGCTTTCTTGGTCAGCGCCTGCCACGGTACTTTGCCAAAATCCGGGTCGCCAAGATAT‑CCCGAGCGGTc < 2:544987/88‑1 (MQ=255) aTCGGCAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATAT‑CCCGAGCGGTc < 1:526972/90‑1 (MQ=255) aaTAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATAT‑CCCGAGCGGTCggcgta < 1:758588/90‑1 (MQ=255) tGGCATAGGCTTTATTGGTCAGCGCCTGCTACGGTACTTTGACAAAATCCGGGTCGCCAAGATAT‑CCCGAGCGGTCGGCGTAGGCGTAtt > 2:264963/1‑90 (MQ=255) tGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATAT‑CCCGAGCGGTCGGCGTAGGCGTAtt > 1:88969/1‑90 (MQ=255) tGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATAT‑CCCGAGCGGTCGGCGTAGGCGTAtt > 2:88592/1‑90 (MQ=255) tGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATAT‑CCCGAGCGGTCGGCGTAGGCGTAtt > 2:21527/1‑90 (MQ=255) tGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATAT‑CCCGAGCGGTCGGCGTAGGCGTAtt > 1:682749/1‑90 (MQ=255) tGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTATGACAAAATCCGGGTCGCCAAGATAT‑CCCGAGCGGTCGGCGTAGGCGTAtt > 1:512063/1‑90 (MQ=255) aTAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATAT‑CCCGAGCGGTCGGCGTAGGCGTATTtctc > 2:346723/1‑90 (MQ=255) tAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATAT‑CCCGAGCGGTCGGCGTAGGCGTATTTCTcc < 1:531231/90‑1 (MQ=255) tttATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATAT‑CCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTc < 1:399028/90‑1 (MQ=255) aTTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATAT‑CCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGc < 2:512063/90‑1 (MQ=255) tGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATAT‑CCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCa < 2:301210/90‑1 (MQ=255) ggTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATAT‑CCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCAt > 1:482852/1‑90 (MQ=255) tCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATAT‑CCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGa < 2:644871/90‑1 (MQ=255) cAGCGCCTGCCACGGTACTTTGGCAAAATCCGGGTCGCCAAGATAT‑CCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGAt > 1:238363/1‑90 (MQ=255) tGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATAT‑CCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTgcatc > 2:73207/1‑90 (MQ=255) aCTTTGACAAAATCCGGGTCGCCAAGATAT‑CCCGAGCGGTCGGGGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGgcg < 2:521968/90‑1 (MQ=255) aCTTTGACAAAATCCGGGTCGCCAAGATAT‑CCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGgcg < 2:174689/90‑1 (MQ=255) cTTTGACAAAATCCGGGTCGCCAAGATAT‑CCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGgcgc > 1:549200/1‑90 (MQ=255) gACAAAATCCGGGTCGCCAAGATAT‑CCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGcc > 1:707119/1‑90 (MQ=255) aCAAAATCCGGGTCGCCAAGATAT‑CCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCa < 1:531772/90‑1 (MQ=255) aaTCCGGGTCGCCAAGATAT‑TCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTc < 1:828354/55‑1 (MQ=255) aaTCCGGGTCGCCAAGATAT‑TCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTc > 2:828354/1‑55 (MQ=255) aaTCCGGGTCGCCAAGATAT‑TCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGc < 2:346186/90‑1 (MQ=255) aaTCCGGGTCGCCAAGATAT‑CCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGc > 2:494425/1‑90 (MQ=255) aaTCCGGGTCGCCAAGATAT‑CCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGc > 2:202777/1‑90 (MQ=255) ccAAGATATCCCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGCCGTATttctt < 1:629014/90‑1 (MQ=255) ccAAGATAT‑CCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGCCGTATttcttc < 1:734460/90‑1 (MQ=255) | CAATTTGATCGGCAATAGATTTGGCATAGGCTTTATTGGTCAGCGCCTGCCACGGTACTTTGACAAAATCCGGGTCGCCAAGATAT‑TCCGAGCGGTCGGCGTAGGCGTATTTCTCCGCTTCTGCCATGATTTGCATCGCATCGGCGCTGCCAAAGCCGTATTTCTTC > NZ_CP009273/3579082‑3579248 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |