Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 3,908,697 | T→C | N248S (AAC→AGC) | pstA ← | phosphate ABC transporter membrane subunit PstA |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,908,697 | 0 | T | C | 97.3% | 121.2 / ‑4.2 | 37 | N248S (AAC→AGC) | pstA | phosphate ABC transporter membrane subunit PstA |
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (24/12); total (24/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.51e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01 |
TCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAACGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCGGTCATAATCCCGGACACCGAC > NC_000913/3908566‑3908825 | tcGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAACGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGCTGGCGATc > 2:7205/1‑140 (MQ=255) cGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCACACGGGCTCATCGGAAACTTAAAGATGGTCATCGGCAGGCTGGCGATCg > 2:281850/1‑140 (MQ=255) cGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAACGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGCTGGCGATCg > 2:412352/1‑140 (MQ=255) cGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAACGGGCTCATAGCAAACTCAAAGATGGTCACCAGCAGGCTGGCGATCg > 1:348165/1‑140 (MQ=255) cAAAAACAACGCGCGTCAGAATCTTCAGCCGCAGTACACACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAACGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGCTGGCGATCGGc < 1:344042/140‑1 (MQ=255) cAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAACGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGc < 2:421733/140‑1 (MQ=255) cAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAACGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGCTGGCGATCGGc < 1:773534/140‑1 (MQ=255) cgcCAGAATGTTCAGCACCAGTCCGCACAGGTTAAAGATCAATACCCCTGCCCAGGCCAATTGCTGCCATTCCGCAAACGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGCTGGCGATCGGCTGCATCATGTCCg < 2:217987/140‑1 (MQ=255) ccAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAACGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGCTGGCGATCGGCTGCATCATGTCCGTg < 1:703906/140‑1 (MQ=255) gTTCAGCAGCAGTACGCACAGGGTAATGATCAATAGCCCTGCCCAGGCGAATTGCTGCCATTCCGCAAACGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGCTGGCGATCGGCTGCATCATGTCCGTGCTCCAGa < 2:348165/140‑1 (MQ=255) gTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAACGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGCTGGCGATCGGCTGCATCATGTCCGTGCTCCAGa < 1:7205/140‑1 (MQ=255) tCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAACGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGCTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAAc > 1:461032/1‑140 (MQ=255) cagcagTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAACGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGCTGGCGATCGGCTGCATCATGTCCGTGCTCCAGa > 2:622875/1‑134 (MQ=255) cagcagTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAACGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGCTGGCGATCGGCTGCATCATGTCCGTGCTCCAGa < 1:622875/134‑1 (MQ=255) aGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAACGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGCTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAgcg > 1:455378/1‑136 (MQ=255) aGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAACGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGCTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAgcg < 2:455378/136‑1 (MQ=255) aaTACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAACGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGCTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAgcgg > 2:601674/1‑140 (MQ=255) aaTACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAACGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGCTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAgcgg > 1:86561/1‑140 (MQ=255) cGGCCCAGGCCAATTGCTGCCATTCCGCAAACGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGCTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGtt > 1:628032/1‑140 (MQ=255) cGGCCCAGGCCAATTGCTGCCATTCCGCAAACGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGCTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGtt > 2:447146/1‑140 (MQ=255) gCCCAGGCCAATTGCTGCCATTCCGCAAACGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGCTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTc > 1:1807/1‑140 (MQ=255) cccAGGCCAATTGCTGCCATTCCGCAAACGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGCTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCa < 2:86561/140‑1 (MQ=255) aTTGCTGCCATTCCGCAAACGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGCTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATa > 1:788762/1‑140 (MQ=255) ttGCTGCCATTCCGCAAACGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGCTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATAc > 1:158572/1‑140 (MQ=255) gcCATTCCGCAAACGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGCTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCa > 2:746593/1‑140 (MQ=255) gCTCATCGCAAACTTAAAGATCGTCACCGGCAGGCTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATAc < 1:412352/140‑1 (MQ=255) aaGATCGTCACCGGCAGGCTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCGGTCATAATCCCGGa > 1:180208/1‑140 (MQ=255) aaGATCGTCACCGGCAGGCTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCGGTCATAATCCCGGa < 2:788762/140‑1 (MQ=255) aTCGTCACCGGCAGGCTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCTATCACCAGCAGGCTACCGGTCATAATCCCGGacac > 2:341033/1‑140 (MQ=255) tCGTCACCGGCAGGCTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCGGTCATAATCCCGGACAcc < 1:114374/140‑1 (MQ=255) cGTCACCGGCAGGCTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCTATCGCCATCAGTATACCGGTCCTAATCCCGGACACCg > 1:262222/1‑140 (MQ=255) gTCACCGGCAGGCTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCGGTCATAATCCCGGACACCGa > 1:815462/1‑140 (MQ=255) tCACCGGCAGGCTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCGGTCATAATCCCGGACACCGAc > 1:74685/1‑140 (MQ=255) tCACCGGCAGGCTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCGGTCATAATCCCGGACACCGAc > 1:547250/1‑140 (MQ=255) tCACCGGCAGGCTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCGGTCATAATCCCGGACACCGAc > 1:522982/1‑140 (MQ=255) tCACCGGCAGGCTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCGGTCATAATCCCGGACACCGAc > 1:174740/1‑140 (MQ=255) tCACCGGCAGGCTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCGGTCATAATCCCGGACACCGAc > 1:115300/1‑140 (MQ=255) | TCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAACGGGCTCATCGCAAACTTAAAGATCGTCACCGGCAGGTTGGCGATCGGCTGCATCATGTCCGTGCTCCAGAACTGGTTGGAGAGCGCGGTAAACAGCAGCGGCGCGGTTTCACCAGCAATACGGGCAATCGCCAGCAGGATACCGGTCATAATCCCGGACACCGAC > NC_000913/3908566‑3908825 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |