Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 3,910,589 | C→T | A313T (GCG→ACG) | pstS ← | phosphate ABC transporter periplasmic binding protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,910,589 | 0 | C | T | 96.8% | 83.6 / ‑6.1 | 31 | A313T (GCG→ACG) | pstS | phosphate ABC transporter periplasmic binding protein |
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (15/15); total (15/16) | |||||||||||
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.81e-01 |
CCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTAATATTGGTCTTCCACGCAGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCG‑CCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACC > NC_000913/3910355‑3910811 | cccTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTAATATTGGTCTTCCACGCAGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCG‑TCTGTTTAGCCCCGGtt > 2:286165/1‑251 (MQ=255) cccGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTAATATTGGTCTTCCACGCAGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCG‑TCTGTTTAGCCCCGGTTTTGTAc > 1:129037/1‑248 (MQ=255) cccGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTAATATTGGTCTTCCACGCAGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCG‑TCTGTTTAGCCCCGGTTTTGTAc < 2:129037/248‑1 (MQ=255) cccGGTGGGTTAAAAGCAGGCTTGGTGGCAGCCATAAGTTACTCTTCAGTTAAACGCGGTGACAAAGTCATGTTGTAAGGCCGGATGCGGCGTAAAACACCGTGCCCGGCCTGGAGTTTTATTAGTACAGCGGCGTACCGCTACTGTCTTTAATATCGGCCTTCCACGCAGCGCTAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTGCC‑TCTGTTTAGCCCCGTTTCTGTACGcc < 2:192543/251‑1 (MQ=255) ccGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTAATATTGGTCTTCCACGCAGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCG‑TCTGTTTAGCCCCGGTTTTGTACGccc < 1:276924/251‑1 (MQ=255) gggTTAACAGCAGGCGTTGTTTCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACAGCGTACCCGGCTTGGAGGTTTATTAGTACAGCGGCTTACCGTTACTGTCTTGAATACTGGTCTTCCACGCAGCGCGAGCCTGTTCATCGACTCTAGCCGGGAGGCTGTCGTAATCCAGGTCGTGCT‑TCTGTTTAGCACCGGTTTTGTACGCCCAGTc < 2:316230/250‑1 (MQ=255) aGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTAATATTGGGCTTCCACGCAGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCG‑TCTGTTTAGCCCCGGGTTTGTACGCCCAGTCGAAGAATTTCAg < 2:30798/251‑1 (MQ=255) ggCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAACCGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTAATATTGGTCTTCCACGCAGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCG‑TCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGc > 1:192490/1‑251 (MQ=255) cTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTATACATAGTCATTTTGTAATGCCGGATGCGGCAACAAACACCGTACCCGGCCTGGAGTTTTATTCGTAAAGCGGCTTACCGCTACTGTCTTTAATATTGGTCTTCCACGCAACGCGAACCTGTTCACCTACACTATCCGGCAGGCTGGCGAAATCCCGGTCGTTCG‑TCTGTTTAGCCCCGGTTTTGTAAGCCCCGTCGAAGAATTTCAgacc > 2:24362/1‑248 (MQ=255) tgggtGCAGCCATCAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTAATATTGGTCTTCCACGCAGCGCGCACCTGTTCAACTACACTATCCGGCAGGCTGGCGAAATCCAGGTCGTTCG‑TCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGACTTTCAGCACtt > 2:164281/5‑251 (MQ=255) aTAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTAATATTGGTCTTCCACGCAGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCG‑TCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGt < 1:286165/251‑1 (MQ=255) tctTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTAATATTGGTCTTCCACGCAGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCG‑TCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTttc > 1:195089/1‑251 (MQ=255) tctTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTAATATTGGTCTTCCACGCAGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCG‑TCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTttc > 1:280330/1‑251 (MQ=255) ctTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCTGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTCCAGCGGCTTACCGCTACTGTCTTTAATATTGGTCTTCCCCGCAGCGCGAACCTGTTCAACTTCACTATCCGTCAGGCTGGCGTAATCCAGGTCGTTCG‑TCTGTTTAGCCCCGGTTTTGGACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGtt < 2:201363/248‑1 (MQ=255) aaaGTCTTTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTAATATTGGTCTTCCACGCAGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCG‑TCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATc > 1:85543/1‑248 (MQ=255) aaGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTAATATTGGTCTTCCACGCAGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCG‑TCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGaa < 2:80623/251‑1 (MQ=255) cagtttGTAATGCCGGATGCGGCGTAAAACACCGGACCCGGCCTGGAGTTTTATTAGTACAGCAGCTTACCGCTACTGTCTTTAATATTGGTCTTCCACTCAGCGCGAACCCGTTCAACTACACTAGCCGGCAGGCTGGCGTAATCCAGGTCGTTCTTCCT‑TTTAGCCCCGGGTTTGTACGCCCAGTCGAAGAATTTCAGCACGTCTGTGCCGTGTGCTGGTTTCTTCTGATCTTTGTGGATCAGAATGa < 1:41480/247‑1 (MQ=255) tttGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTAATATTGGTCTTCCACGCAGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCG‑TCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCt < 1:22689/204‑1 (MQ=255) aTGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTAATATTGGTCTTCCACGCAGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCG‑TCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGt > 1:56226/1‑248 (MQ=255) gAGTTTTATTAGTACAGCGGCTTCCCGCTACTGTCTTTAATATTGGTCTTCCACGCAGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCG‑TCTGTTTAGCCCCGGTTTTGTACGCCCATTTGAAGAATTTCAGCGCTTCGGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCtggttggt < 1:54057/251‑1 (MQ=255) atagtGGTCCTCCACGCAGCTCGAACCTGTTCAACTACACTATCCGGCATGCTGGCGTAATCCATGTCGGTCG‑TCGGTTTAGCCCTGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGGTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATATGCCATGCATCTTCGCCTTTCTTGTTGGTCAGATCCTGAGCCAAGGTTTTTCTCCAGTCTGCCCCttt < 2:133061/247‑1 (MQ=255) atTGGTCTGCCAAGCAGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTAG‑TCTGTTTAGCCCAGGTTTTGAGCGACCAGTCGACGAATTTCAGCACTTATGTGCGTTGTTCTAGTTTCTGCTGGTCTTTGTGGATAAGATTGCCCGTGGACGAGGTAATAGGTGATGCATCTTCGCCTTTCTGGTTGGTCTGATCCTGAGCGATGGTTTTGCTCCAGTCTGCACCTTTtg > 2:110116/1‑251 (MQ=255) gAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAGTCCAGGACGTTCG‑TCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCg > 1:40888/1‑203 (MQ=255) gAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCG‑TCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTttct > 2:265636/1‑251 (MQ=255) gAACCTGTTCAACTACACTACCCGGCAGGCTGGCGTAATCCAGGTCGTTCG‑TCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCGGCGCCTTTCTGGTTGTTCAGTTCCTGAGCg < 2:40888/203‑1 (MQ=255) ccTGTTCAAATACACTATCCGGCAGGCTGGCGTAATCCAGGTAGTTCG‑TCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGCATTTCAGCACTGCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAATGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTttcttcc < 2:233172/251‑2 (MQ=255) tCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCG‑TCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTtcggtcgg > 2:237096/1‑251 (MQ=255) acacTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCG‑TCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAAc > 2:45040/1‑251 (MQ=255) agcccgCAGGCCGGCGTAATCCAGGGCGTTCG‑TCTGGTTAGCCCCGGTTTTGTACGCCTAGTCGAAGAATTGCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTTGTAGAGGTAATAGGCTATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTGGCTGCATTAGCGAAGTTTTCTGCGGTCGGACTAACCGGtt < 2:258713/246‑1 (MQ=255) ggCAGGCTGGCGTAATCCAGGTCGTTCG‑TCTGTTTAGCCCAGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGCTCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTAcc > 2:190395/1‑251 (MQ=255) aaTCCAGGTCGTTCG‑TCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTtcgg < 2:86762/219‑1 (MQ=255) | CCCTTTTTACCCGGTGGGTTAAAAGCAGGCTTGGTTGCAGCCATAAGTTACTCTTCAGTTAAACGCGTTTACAAAGTCATTTTGTAATGCCGGATGCGGCGTAAAACACCGTACCCGGCCTGGAGTTTTATTAGTACAGCGGCTTACCGCTACTGTCTTTAATATTGGTCTTCCACGCAGCGCGAACCTGTTCAACTACACTATCCGGCAGGCTGGCGTAATCCAGGTCGTTCG‑CCTGTTTAGCCCCGGTTTTGTACGCCCAGTCGAAGAATTTCAGCACTTCTGTGCCTTGTTCTGGTTTCTTCTGATCTTTGTGGATCAGAATGAACGTGGTAGAGGTAATAGGCCATGCATCTTCGCCTTTCTGGTTGGTCAGATCCTGAGCGAAGGTTTTGCTCCAGTCTGCACCTTTTGCTGCATTAGCGAAGTTTTCTTCGGTCGGACTAACCGGTTTACC > NC_000913/3910355‑3910811 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |