| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 825,102 | G→A | A336V (GCG→GTG) | ybhR ← | putative ABC transporter permease |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 825,102 | 0 | G | A | 95.1% | 113.4 / ‑2.0 | 41 | A336V (GCG→GTG) | ybhR | putative ABC transporter permease |
| Reads supporting (aligned to +/- strand): ref base G (1/1); new base A (20/19); total (21/20) | |||||||||||
| 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.14e-01 | |||||||||||
CGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAA > NC_000913/824970‑825237 | cggTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCACATCCt < 1:494672/139‑1 (MQ=255) gtgtCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCACATCCTTCaaa < 1:575442/139‑1 (MQ=255) tGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCACATCCTTCAAATAAATCTGCt < 1:13087/139‑1 (MQ=255) ggCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCACATCCTTCaaataaat > 2:34445/1‑133 (MQ=255) ggCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCACATCCTTCaaataaat < 1:34445/133‑1 (MQ=255) aCCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCACATCCTTCAAATAAATCTGCTTGGTAATGTCCGTaaa > 2:418913/1‑131 (MQ=255) aCCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCACATCCTTCAAATAAATCTGCTTGGTAATGTCCGTaaa < 1:418913/131‑1 (MQ=255) ccTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCATTAGCGGCCACAAACTATTCCACACAATATCCAGACTCACATCCTTCAAATAAATCTGCGTGGTAATGTCCGTAAAGTGGCGAAt > 1:591868/1‑139 (MQ=255) ccTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCACATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAAt < 2:1025622/139‑1 (MQ=255) ccTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCACATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAAt < 1:353173/139‑1 (MQ=255) ttACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCACATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAg > 2:970715/1‑139 (MQ=255) tACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCACATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGaa > 1:978086/1‑135 (MQ=255) tACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCACATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGaa < 2:978086/135‑1 (MQ=255) cGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCACATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGt > 1:241334/1‑139 (MQ=255) tCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCACATCCTTCaaataa > 2:99355/1‑92 (MQ=255) tCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCACATCCTTCaaataa < 1:99355/92‑1 (MQ=255) gcgTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCACATCCt > 1:282383/1‑82 (MQ=255) gcgTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCACATCCt < 2:282383/82‑1 (MQ=255) cgTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCACATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGt < 1:67049/126‑1 (MQ=255) cgTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCACATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGt > 2:67049/1‑126 (MQ=255) gTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCACATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAgg > 2:793800/1‑139 (MQ=255) tACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCACATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCg > 1:214596/1‑116 (MQ=255) tACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCACATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCg < 2:214596/116‑1 (MQ=255) gTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCACATCCTTCaaataaat > 1:969499/1‑74 (MQ=255) gTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCACATCCTTCaaataaat < 2:969499/74‑1 (MQ=255) ggCCACAAACTATTCCACACAATATCCAGACTCACATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGt < 1:970715/139‑1 (MQ=255) gCCACAAACTATTCCACACAATATCGAGACTCACATCCCTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCAACGTCAGGTTTTGCAGCAATACCGGCATGTTTTCCACCGGAGAAACGTa > 2:701462/1‑139 (MQ=255) gCCACAAACTATTCCACACAATATCCGGACTCACATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTa > 2:473123/1‑139 (MQ=255) cacaAACTATTCCACACAATATCCGGACTCACATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAAc < 1:473123/139‑1 (MQ=255) cacaAACTATTCCACACAATATCCAGACTCACATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAg > 2:701511/1‑95 (MQ=255) cacaAACTATTCCACACAATATCCAGACTCACATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAAc > 2:290210/1‑139 (MQ=255) cacaAACTATTCCACACAATATCCAGACTCACATCCTTCAAATAAATCTGCTAGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAg < 1:701511/95‑1 (MQ=255) tCCACACAATATCCAGACTCACATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATg < 2:349212/107‑1 (MQ=255) tCCACACAATATCCAGACTCACATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATg > 1:349212/1‑107 (MQ=255) tatCCAGACTCACATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCg < 2:650186/93‑1 (MQ=255) tatCCAGACTCACATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCg > 1:650186/1‑93 (MQ=255) atCCAGACTCACATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGc < 1:328793/139‑1 (MQ=255) ccAGACTCACATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTAAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGcat > 1:888448/1‑139 (MQ=255) cTCACATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATaa > 2:243767/1‑139 (MQ=255) cGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTc < 2:124461/93‑1 (MQ=255) cGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTc > 1:124461/1‑93 (MQ=255) | CGGTAGTGTCTTTGCTGGCGAAAGATAAGAAGTTACATCACCTTACGTCTAAACATCGCGTACGCCGCTGACCCTGTCGTGGCCGTTATCACCAGTAGCGGCCACAAACTATTCCACACAATATCCAGACTCGCATCCTTCAAATAAATCTGCTTGGTAATGTCCGTAAAGTGGCGAATAGGGTTAATCCACGTCAGGTTTTGCAGCCATACCGGCATGTTTTCCACCGGAGAAACGTAACCGGAAAGGAGAATGGCGGGCATCATAA > NC_000913/824970‑825237 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |