Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 2,808,251 | T→G | K516N (AAA→AAC) | gshA ← | glutamate‑‑cysteine ligase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,808,251 | 0 | T | G | 90.2% | 130.5 / 6.8 | 41 | K516N (AAA→AAC) | gshA | glutamate‑‑cysteine ligase |
Reads supporting (aligned to +/- strand): ref base T (1/3); new base G (19/18); total (20/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.06e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01 |
TCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGCTCTTCACGCAGCAGA > NZ_CP009273/2808136‑2808380 | tcttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGGTTTTCCAGCCACACCGCaaa < 2:2694/136‑1 (MQ=255) tcttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGGTTTTCCAGCCACACCGCaaa < 1:709233/136‑1 (MQ=255) tcttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGGTTTTCCAGCCACACCGCAAAc < 1:37392/137‑1 (MQ=255) cccAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGGTTTTCCAGCCACACCGCAAACGGtt < 2:826583/137‑1 (MQ=255) cAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGGTTTTCCAGCCACACCGCAAACGGTTCgg > 2:730112/1‑138 (MQ=255) gTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGGTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCg > 1:495513/1‑139 (MQ=255) tCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGGTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCgg > 1:410512/1‑139 (MQ=255) tCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGGTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCgg > 1:653921/1‑139 (MQ=255) aGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTc > 1:856787/1‑139 (MQ=255) aGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGGTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGc < 1:18758/136‑1 (MQ=255) aGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGGTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCtt < 2:170991/138‑1 (MQ=255) tGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGGTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTc > 1:780194/1‑139 (MQ=255) aCGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGGTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCctgct > 2:466693/1‑138 (MQ=255) aCGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGGTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCctgct > 1:689519/1‑138 (MQ=255) aCGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGGTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCctgct > 1:143729/1‑138 (MQ=255) cgcgCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGGTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGa < 2:484451/139‑1 (MQ=255) cATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGc < 2:352965/139‑1 (MQ=255) atcatcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGGTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGcc > 1:712748/1‑139 (MQ=255) cccTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGGTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTcgcg > 2:119975/1‑139 (MQ=255) ccTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGGTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTcgcg < 2:679824/138‑1 (MQ=255) aaTTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGGTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCgg > 1:598669/1‑139 (MQ=255) tGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGGTTTTCCAGCCACACCGCaaa < 1:757572/70‑1 (MQ=255) tGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGGTTTTCCAGCCACACCGCaaa > 2:757572/1‑70 (MQ=255) gCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCgcgc < 2:850783/139‑1 (MQ=255) aCGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGGTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAAttt > 1:421748/1‑139 (MQ=255) gCCTTTTTCTTTTCTGTCAGGCGTGGTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCgg > 1:806606/1‑139 (MQ=255) ccTTTTTCTTTTCTGTCAGGCGTGGTTTTCCAGCCACACCGCaaa > 1:183786/1‑45 (MQ=255) ccTTTTTCTTTTCTGCCAGGCGTGGTTTTCCAGCCACACCGCaaa < 2:183786/45‑1 (MQ=38) cTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGct < 2:736527/139‑1 (MQ=255) tttttCTTTTCTGTCAGGCGTGGTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGctc > 2:685022/1‑139 (MQ=255) tttcttttctGTCAGGCGTGGTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGCTCtt > 2:384209/1‑139 (MQ=255) tttcttttctGTCAGGCGTGGTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGCTCtt > 1:261595/1‑139 (MQ=255) cttttctGTCAGGCGTGGTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGCTCTTCa < 1:670465/138‑1 (MQ=255) cttttctGTCAGGCGTGGTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGCTCTTCa < 1:466693/138‑1 (MQ=255) tttctGTCAGGCGTGGTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGCTCTTCACg > 2:281921/1‑138 (MQ=255) ttctGTCAGGCGTGGTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGCTCTTCACg < 2:806606/137‑1 (MQ=255) tctGTCAGGCGTGGTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGCTCTTCACgca < 2:717563/138‑1 (MQ=255) ctGTCAGGCGTGGTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGCTCTTCACgcag < 2:143729/138‑1 (MQ=255) ctGTCAGGCGTGGTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGCTCTTCACgcag < 1:385549/138‑1 (MQ=255) tCAGGCGTGGTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGCTCTTCACgcagca < 1:176903/137‑1 (MQ=255) tCAGGCGTGGTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGCTCTTCACGCAGCAGa < 2:821931/139‑1 (MQ=255) | TCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACCGCAAACGGTTCGGTATCAGCGGCTTCCATTTCCTGCTGACGGCGTTCAGACGCCTCGCGCTCGGCTACAAAATCCTCTTCGCGCAGAATTTCCAGCGGCTCTTCACGCAGCAGA > NZ_CP009273/2808136‑2808380 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |