Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 3,400,225 | T→C | E288G (GAG→GGG) | mreB ← | dynamic cytoskeletal protein MreB |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,400,225 | 0 | T | C | 97.2% | 116.6 / ‑4.2 | 36 | E288G (GAG→GGG) | mreB | dynamic cytoskeletal protein MreB |
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (16/19); total (17/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.72e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01 |
GAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGG > NC_000913/3400148‑3400301 | gAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCCCgg < 2:644279/81‑1 (MQ=255) gAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCCCgg < 2:644277/81‑1 (MQ=255) gAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCCCg > 2:644333/1‑80 (MQ=255) aTGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGgag < 2:644281/81‑1 (MQ=255) aTGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGgag < 2:644280/81‑1 (MQ=255) ggTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTc < 1:644282/81‑1 (MQ=255) tttCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTCgg < 2:644283/81‑1 (MQ=255) cttcCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAg > 2:644334/1‑81 (MQ=255) ccATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGcc < 1:644284/80‑1 (MQ=255) cATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTTGGAAGCCAg < 1:644285/81‑1 (MQ=255) ttAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGtt < 2:644286/81‑1 (MQ=255) aaCAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGTTcc > 2:644332/1‑81 (MQ=255) aaCGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGTTCcggcg > 2:644348/1‑81 (MQ=255) cGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCggg < 2:644287/81‑1 (MQ=255) gTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCa > 1:644352/1‑81 (MQ=255) aaGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACt < 2:807660/80‑1 (MQ=255) ggTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGtt > 1:644335/1‑81 (MQ=255) gTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTc > 1:644337/1‑81 (MQ=255) gTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTc < 1:644288/81‑1 (MQ=255) ttACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTc > 2:644356/1‑80 (MQ=255) cagcagTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGc < 1:644289/81‑1 (MQ=255) gcagTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCaa < 1:644290/81‑1 (MQ=255) gcagTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCaa > 2:644359/1‑81 (MQ=255) gTGCGCCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCa > 1:644338/1‑81 (MQ=255) gcgcCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCAt > 1:644341/1‑80 (MQ=255) cgcCACCACCGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCAtt < 2:644291/80‑1 (MQ=255) accaccGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACcgc < 1:644292/81‑1 (MQ=255) caccGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCgcgc > 1:644355/1‑81 (MQ=255) ccGGTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTc < 2:644293/81‑1 (MQ=255) ggTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTca > 2:644343/1‑80 (MQ=255) gTGAGCACCATGCCGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTcac > 1:644364/1‑80 (MQ=255) cATGCCGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGt > 1:644336/1‑81 (MQ=255) aTGCCGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTc < 2:644294/81‑1 (MQ=255) ccGCGCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAg < 1:644295/80‑1 (MQ=255) cgcgCCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCg > 1:644340/1‑81 (MQ=255) gcgcCCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCgg > 1:644351/1‑81 (MQ=255) | GAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACCGGTGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGG > NC_000913/3400148‑3400301 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |