Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 3,459,986 | G→A | S57N (AGC→AAC) | gspF → | Type II secretion system protein GspF |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,459,986 | 0 | G | A | 96.0% | 74.3 / ‑3.9 | 25 | S57N (AGC→AAC) | gspF | Type II secretion system protein GspF |
Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (11/13); total (12/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.80e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.19e-01 |
CAAGGGATCATTGATGCTAACGATGAACGTCAGGCACGACTGCGGCTGCGTGAAGAAGGGCTTTTCCTGCTGGATATTCGCCCCCAAAAAAGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAGCCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCA > NC_000913/3459862‑3460115 | cAAGGGATCATTGATGCTAACGATGAACGTCAGGCACGACTGCGGCTGCGTGAAGAAGGGCTTTTCCTGCTGGATATTCGCCCCCAAAAAAGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAACCATAGTg > 1:20492/1‑133 (MQ=255) gggATCATTTGTGGTTACCGTGGACGTCAGGCACGGCTGGGGCTGCGGGAAGAAGGGCTTTTCCTGCTGGGTTGTCGCCCCCAAAAAAGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAACCATAGTGAAc < 2:92996/133‑1 (MQ=255) aTCATTGATGCTAACGATGAACGTCAGGCACGACTGCGGCTGCGTGAAGAAGGGCTTTTCCTGCTGGATATTCGCCCCCAAAAAAGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAACCATAGTGAACTGa < 1:49402/133‑1 (MQ=255) gCTAACGATGAACGTCAGGCACGACTGCGGCTGCGTGAAGAAGGGCTTTTCCTGCTGGATATTCGCCCCCAAAAAAGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAACCATAGTGAACTGACGCTTTTCa > 1:90314/1‑133 (MQ=255) aaCGATGAACGTCAGGCACGACTGCGGCTGCGTGAAGAAGGGCTTTTCCTGCTGGATATTCGCCCCCAAAAAAGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAACCATAGTGAACTGACGCTTTTCAccc > 2:57309/1‑133 (MQ=255) aaCGATGAACGTCAGGCACGACTGCGGCTGCGTGAAGAAGGGCTTTTCCTGCTGGATATTCGCCCCCAAAAAAGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAACCATAGTGAACTGACGCTTTTCAccc > 2:362080/1‑133 (MQ=255) aaCGTCAGGCACGACTGCGGCTGCGTGAAGAAGGGCTTTTCCTGCTGGATATTCGCCCCCAAAAAAGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAACCATAGTGAACTGACGCTTTTCACCCGGCAGtt > 2:76722/1‑133 (MQ=255) cGTCAGGCTCGACTGCGGGTGCCTGAAGAACGGCTTTTCCTGCTGGATATTCGCCCCCAAAAAAGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAACCATAGTGAACTGACGCTTTTCACCCGGCAGTTgg < 1:249509/133‑1 (MQ=255) tCAGGCACGACTGCGGCTGCGTGAAGAAGGGCTTTTCCTGCTGGATATTCGCCCCCAAAAAAGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAACCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCa < 1:106790/133‑1 (MQ=255) gACTGCGGCTGCGTGAAGAAGGGCTTTTCCTGCTGGATATTCGCCCCCAAAAAAGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAACCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAg < 1:407237/133‑1 (MQ=255) tGCGGCTGCGTGAAGAAGGGCTTTTCCTGCTGGATATTCGCCCCCAAAAAAGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAACCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAgcgc < 1:76722/133‑1 (MQ=255) gCGTGAAGAAGGGCTTTTCCTGCTGGATATTCGCCCCCAAAAAAGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAACCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCa > 2:376339/1‑133 (MQ=255) gTGAAGAAGGGCTTTTCCTGCTGGATATTCGCCCCCAAAAAAGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAACCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAATCTTAAGCGCAGCGGCAtt > 1:216233/1‑133 (MQ=255) gTGAAGAAGGGCTTTTCCTGCTGGATATTCGCCCCCAAAAAAGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAACCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCAtt < 2:216233/133‑1 (MQ=255) ccTGCTGGATATTCGCCCCCAAAAAAGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAACCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGc < 1:376339/133‑1 (MQ=255) ctgGATATTCGCCCCCAAAAAAGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAACCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTg > 1:266077/1‑133 (MQ=255) atatTCGCCCCCAAAAAAGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAACCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGt < 2:86788/133‑1 (MQ=255) cccccAAAAAAGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAACCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGt < 2:90314/133‑1 (MQ=255) ccAAAAAAGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAACCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTcaa < 1:362080/133‑1 (MQ=255) cAAAAAAGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAACCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTcaac > 1:322803/1‑133 (MQ=255) gTTCGGGAGTAAAAACACGTCGCCCGAGGATCAACCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAg < 1:57309/133‑1 (MQ=255) cGTCGCCCGAGGATCAACCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTg > 2:110696/1‑133 (MQ=255) tCGCCCGAGGATCAACCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGAc > 1:481217/1‑133 (MQ=255) aGGATCAACCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTaa < 2:266077/133‑1 (MQ=255) tCAGCCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCa > 2:389715/1‑133 (MQ=255) | CAAGGGATCATTGATGCTAACGATGAACGTCAGGCACGACTGCGGCTGCGTGAAGAAGGGCTTTTCCTGCTGGATATTCGCCCCCAAAAAAGTTCGGGAGTAAAAACACGTCGCCCGAGGATCAGCCATAGTGAACTGACGCTTTTCACCCGGCAGTTGGCAACCTTAAGCGCAGCGGCATTACCCCTGGAAGAGAGCCTTGCCGTAATCGGTCAACAAAGCAGTAATAAACGACTGGGTGACGTGTTAAATCA > NC_000913/3459862‑3460115 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |