| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 3,800,629 | G→A | R113C (CGT→TGT) | waaY ← | lipopolysaccharide core heptose (II) kinase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,800,629 | 0 | G | A | 93.3% | 33.9 / ‑3.6 | 15 | R113C (CGT→TGT) | waaY | lipopolysaccharide core heptose (II) kinase |
| Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (8/6); total (9/6) | |||||||||||
| 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.63e-01 | |||||||||||
CACCATTTTTTATAATGAAGTTACCACGATGGGGGTCGCCAGAAACCATGCCATGTTGATGTAAGGCATTAATTGATTGCTGAATTTTATTTTTTAGCGCATCATCAATATCGGGCATATCACACAACTCTATGCCATCAATATACTCGATGATCATGATATAAGTATGGACAAAACGTAAGGTTTTCCGTTCAGCCAATAAATAAAAGTCATTGAGTGTATTTAACCCTTCATTTCGTACTTTTTGGGTTTGCTCAAAAAGGCGTTCGTAATAATCACCTTTTAACAGAGACTTAAAGAAACGTTCGTTACGCTTAACTTTCGGAGAAAAAACCTTAAGAATCAATTTACCGTAATCGGTATCAATAAGCATAACTTTTGTATCATCAATAGAACGAAAAAC > NC_000913/3800452‑3800854 | cACCATTTTTTATAATGAAGTTACCACGATGGGGGTCGCCAGAAACCATGCCATGTTGCTGTAAGGCATTAATTGATTGCTGAATTTTATTTTTTAGCGCATCATCAATATCGGGCATATCACACAACTCTATGCCATCAATATACTAGATGATCATGATATAAGTATGGACAAAACATAAGGTTTTCCGTTCAGCCAATAAATAAAAGTCATTGAg > 1:53138/1‑217 (MQ=255)cACCATTTTTTATAATGAAGTTACCACGATGGGGGTCGCCAGAAACCATGCCATGTTGATGTAAGGCATTAATTGATTGCTGAATTTTATTTTTTAGCGCATCATCAATATCGGGCATATCACACAACTCTATGCCATCAATATACTCGATGATCATGATATAAGTATGGACAAAACATAAGGTTTTCCGTTCAGCCAATAAATAAAAGTCATTGAGTGTATTTAACCCTTCATTTCGTACTTTTTggg < 1:4946/249‑1 (MQ=255) ggTCGCCAGAAACCATGCCATGTTGATGTAAGGCATTAATTGCTTGCTGAATTTTATTTTTTAGCGCATCATCAATATCGGGCATATCACACAACTCTATGCCATCAATATACTCGATGATCATGATATAAGTATGGACAAAACATAAGGTTTTCCGTTCAGc > 1:230574/1‑163 (MQ=255) ggTCGCCAGAAACCATGCCATGTTGATGTAAGGCATTAATTGATTGCTGAATTTTATTTTTTAGCGCATCATCANTATNGGGCNTANCACACAACTCTATGCNATCAATATACTCGATGATCATGATATAAGTATGGACAAAACATANNGNNNTNCGTTCAGc < 2:230574/163‑1 (MQ=255) gTCGCCAGAAACCATGCCATGTTGCTGTAAGGCATTAATTGATTGCTGAATTTTATTTTTTCGCGCATCATCAATCTCGGGCATATCACACAACTCTATGCCATCAATATACTCGATGATCATGATATAAGTATGGACAAAACGTAAGGTTTTCCGTTCAGCCAATAAATAAAAGTCATTGAGTGTATTTAAACCTTCATTTCGTACttttt > 1:232406/1‑212 (MQ=255) ccatgTTGATGTAAGGCATTAATTGATTGCTGAATTTTATTTTTTAGCGCATCATCAATATCGGGCATATCACACAACTCTATGCCATCAATATACTCGATGATCATGATATAAGTATGGACAAAACATAAGGTTTTCCGTTCAGCCAATAAATAAAAGTCATTGAGTGTATTTAACCCTTCATTTCGTACTTTTTGGGTTTGCTCAAAAAGGCGTTCGTAATAATCACCTTTTAACAGAGACTTAAAGa > 1:72302/1‑250 (MQ=255) gattgCTGAATTTTATTTTTTAGCGCATCAGCAATATCGGGCATATCACACAACTCTATGCCATCAATATACTCGATGATCATGATATAAGTATGGACAAAACATAAGGTTTTCCGTTCAGCCAATAAATAAAAGTCATTGAGTGTATTTAAccc < 1:17008/155‑1 (MQ=255) ttgCTGAATTTTATTTTTTAGCGCATCATCAATATCGGGCATATCACACAACTCTATGCCATCAATATACTCGATGATCATGATATAAGTATGGACAAAACATAAGGTTTTCCGTTCAGCCAATAAATAAAAg > 1:134308/1‑133 (MQ=255) tgCTGAATTTTATTTTTTAGCGCATCATCAATATCGGGCATATCACACAACTCTATGCCATCAATATACTCGATGATCATGATATAAGTATGGACAAAACATAAGGTTTTCCGTTCAGCCAATAAATAAAAGTCATTGAGTGTATTTAACCCATCATTTCGTACTTTTTGGGTTTGCTCAAAAAGGCGGTCGTAATAATCACCTTTTAACAGAGACTTAAAGAAACGTTCGTTACGCTTAACTTTCGgag < 1:304636/250‑1 (MQ=255) ncgcATCATCAATATCGGGCATATCACACAACTCTATGCCATCAATATACTCGATGATCATGATATAAGTATGGACAAAACATAAGGTTTTCCGTTCAGCCAATAAATAAAAGTCATTGAGTGTATTTAACCCTTCATTTCGTACTTTTTgg > 1:84787/2‑152 (MQ=255) ncgcATCATCAATATCGGGCATATCACACAACTCTATGCCATCAATATACTCGATGATCATGATATAAGTATGGACAAAACATAAGGTTTTCCGTTCAGCCAATAAATAAAAGTCATTGAGTGTATTTAACCCTTCATTTCGTACTTTTTGGGTTTGCTCAAAAAGGCGTTCGTAATAATCAcc > 1:94207/2‑184 (MQ=255) gcgcATCATCAATATCGGGCATATCACACAACTCTATGCCATCAATATACTCGATGATCATGATATAAGTATGGACAAAACATAAGGTTTTCCGTTCAGCCAATAAATAAAAGTCATTGAGTGTATTTAACCCTTCATTTCGTACTTTTTGGGTTTGCTCAAAAAGGCGTTCGTAATAATCACCTTTTAACAGAGACTTAAAGAAACGTTCGTTACGCTTAACTTTCGGAGAAAAAACCTTAAGAATCAAt > 1:137360/1‑251 (MQ=255) cacaACTCTATGCCATCAATATACTCGATGATCAGTATAGAAGTATGGACAAAACATAAGGTTTTCCGTGCAGCCAATAAATAAAAGTCACTGTGTGTATTTATCCCTTCAGTGCGTGCTTTTTGGGTTTGCTCAAAAAGGCg < 1:269882/143‑1 (MQ=255) tanaCTNGNNGNNCATNATATAAGTATGGACAAAACATNNGGTGTTCCGNTCAGCCAATAAATAAAAGTCATTGAGTGTATTTAACCCTTCATTTCGTACTTTNTGGGTTGGCNNAANAAGGCNTNNNTNNNAANCACCTTTTAACAGAGACTTNAANNAACGTTNGTTACGCTTAACTTTCGGAGaaaaaa < 2:72302/192‑1 (MQ=255) atCATGATATAAGTATGGACAAAACATAAGGTTTTCCGTTCAGCCAATAAATAAAAGACATTGAGTGTCTTTAAACCTTCATTTCGTACTTTTTGGGTTTGCTCAACAAGGCGTTCGTAATAATCACCTTTTACCAGAGAATTAAAGAAACGTTCGTTACGCTTAACTTTCGGAGAAAAAACCTTAAGAATCAATTTACCGTAATCGGTATCACTAAGCATAACTTTTGTATCATCAATAGAACGAAAAAc > 1:151666/1‑251 (MQ=255) | CACCATTTTTTATAATGAAGTTACCACGATGGGGGTCGCCAGAAACCATGCCATGTTGATGTAAGGCATTAATTGATTGCTGAATTTTATTTTTTAGCGCATCATCAATATCGGGCATATCACACAACTCTATGCCATCAATATACTCGATGATCATGATATAAGTATGGACAAAACGTAAGGTTTTCCGTTCAGCCAATAAATAAAAGTCATTGAGTGTATTTAACCCTTCATTTCGTACTTTTTGGGTTTGCTCAAAAAGGCGTTCGTAATAATCACCTTTTAACAGAGACTTAAAGAAACGTTCGTTACGCTTAACTTTCGGAGAAAAAACCTTAAGAATCAATTTACCGTAATCGGTATCAATAAGCATAACTTTTGTATCATCAATAGAACGAAAAAC > NC_000913/3800452‑3800854 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |