Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 3,328,104 | A→G | V292A (GTG→GCG) | mreB ← | cell wall structural complex MreBCD, actin‑like component MreB |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,328,104 | 0 | A | G | 97.6% | 146.8 / ‑3.7 | 41 | V292A (GTG→GCG) | mreB | cell wall structural complex MreBCD, actin‑like component MreB |
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (30/10); total (31/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.03e-01 |
GACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACC‑GGTGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTC > REL606/3328034‑3328175 | gACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACC‑GGTGAGCGCCATg < 2:833725/76‑1 (MQ=255) gACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACC‑GGTGAGCGCCATg > 2:222058/1‑76 (MQ=255) gACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACC‑GGTGAGCGCCATg < 2:1881700/76‑1 (MQ=255) aCTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACC‑GGTGAGCGCCATGc > 2:700989/1‑76 (MQ=255) aCTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACC‑GGGGAGCGCCACGc > 1:948639/1‑76 (MQ=255) cTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACC‑GGTGAGCGCCATGcc > 2:105969/1‑76 (MQ=255) aTGCCGGTTTCTCCCATACACAAACGGCCAAGTTTCCGCGGCAGTGCCCCCCCGCC‑GGTGAGCGCCATGCCGCGCt < 2:603796/76‑1 (MQ=255) gCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACC‑GGTGAGCGCCATGCCGCGCTCg > 1:1703315/1‑76 (MQ=255) gTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACC‑GGTGAGCGCCATGCCGCGCTCGgaga > 2:1102326/1‑76 (MQ=255) tttCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACC‑GGTGAGCGCCATGCCGCGCTCGGAGAt > 2:1843569/1‑76 (MQ=255) tttCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACC‑GGTGAGCGCCATGCCGCGCTCGGAGAt > 1:1392224/1‑76 (MQ=255) tcttcCATAAACAACCGGTCAAGGTTACGCAGCAGTGCGCCCCCCCC‑GGTGAGCGCCATGCCGCGCTCGGAGATGt < 2:1222138/76‑1 (MQ=255) cttcCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACC‑GGTGAGCGCCATGCCGCGCTCGGAGATGTc > 1:1444018/1‑76 (MQ=255) cccATAAACAACCGCTGAGGGTTACGCACGTGTGCGCCCCCACG‑GGTGAGCGCCATGCCGCGCTCGGAGATGTCgg < 1:1197248/75‑1 (MQ=255) aTTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACC‑GGTGAGCGCCATGCCGCGCTCGGAGATGTCGGAAg > 1:1332102/1‑76 (MQ=255) ttAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACC‑GGTGAGCGCCATGCCGCGCTCGGAGATGTCGGAAGc > 2:1337841/1‑76 (MQ=255) tAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACC‑GGTGAGCGCCATGCCGCGCTCGGAGATGTCGGAAGcc > 2:625103/1‑76 (MQ=255) aCAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACC‑GGTGAGCGCCATGCCGCGCTCGGAGATGTCGGAAGCCAg > 2:1019745/1‑76 (MQ=255) cAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACC‑GGTGAGCGCCATGCCGCGCTCGGAGATGTCGGAAGCCAGt > 2:1371009/1‑76 (MQ=255) aaccGGTCAAGGTTACGCAGCAGTGCGCCACCACC‑GGTGAGCGCCATGCCGCGCTCGGAGATGTCGGAAGCCAGtt < 1:1787046/73‑1 (MQ=255) aaccGGTCAAGGTTACGCAGCAGTGCGCCACCACC‑GGTGAGCGCCATGCCGCGCTCGGAGATGTCGGAAGCCAGtt < 2:1469888/73‑1 (MQ=255) aaaCGGTCAAGGTTACGCAGCAGTGCGCCACCACC‑GGTGAGCGCCATGCCGCGCTCGGAGCTGTCGGAAGCCAGtt > 2:877051/1‑76 (MQ=255) aaaCGGTCAAGGTTACGCAGCAGTGCGCCACCACC‑GGTGAGCGCCATGCCGCGCTCGGAGATGTGGGTAGCCGGtt > 2:1354505/1‑76 (MQ=255) agcGGTCAGGGTTACGCGGCGGTGCGCCCCCCCC‑GGTGAGCGCCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTc < 2:1296728/74‑1 (MQ=255) gTCAAGGTTACGCAGCAGTGCGCCACCACC‑GGTGAGCGCCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCcggc > 2:230993/1‑76 (MQ=255) tCAAGGTTACGCACCAGTGCCCCCCCACC‑GGTGAGCGCCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCcggcg < 1:928888/76‑1 (MQ=255) cAAGTTTACGCAGCAGTGCGCCCCCACC‑GGTGAGCGCCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCcggcgg < 1:1909390/76‑1 (MQ=255) cAAGGTTACGCAGCAGTGCGCCACCACC‑GGTGAGCGCCATGCCGCGCTCGGAGATGGCGGAAGCCAGTTCcggcgg > 2:1237874/1‑76 (MQ=255) aaGGTTACGCAGCAGTGCGCCACCACC‑GGTGAGCGCCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCggg > 1:136305/1‑76 (MQ=255) aaGGTTACGCAGCAGTGCGCCACCACC‑GGTGAGCGCCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCggg > 2:1511892/1‑76 (MQ=255) aaGGTTACGCAGCAGCGCGCCACCACC‑GGTGAGCGCCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCggg < 2:696178/76‑1 (MQ=255) aGGTTACGCAGCAGTGCGCCACCACC‑GGTGAGCGCCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGc > 2:1603991/1‑76 (MQ=255) ggTTACGCAGCAGTGCGCCACCACC‑GGTGAGCGCCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGcc > 1:1002037/1‑75 (MQ=255) ttACGCAGCGGTGCGCCACCGCC‑GGTGGGCGCAATGCCGCGCGCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACt < 2:1464575/76‑1 (MQ=255) ttACGCAGCAGTGCGCCACCACC‑GGTGAGCGCCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGccat > 2:1615159/1‑73 (MQ=255) caCGCAGCGGGGCGCCCCCGCG‑GGTGGGCGCCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTg < 1:1546808/75‑1 (MQ=255) cagcagTGCGCCACCACC‑GGTGAGCGCCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTcc > 2:1043607/1‑76 (MQ=255) agcagGGCGCCACCCCC‑GGTGAGCGCCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCa < 2:1209450/76‑1 (MQ=255) cagTGCGCCCCCCCC‑GGTGAGCGCCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGt < 2:1281002/76‑1 (MQ=255) caggcgcCCCCCCCGGGTGCGCGCCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGt < 2:460607/73‑1 (MQ=255) ggcgcCACCACC‑GGTGACCGCCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCa < 2:1669443/75‑1 (MQ=255) gcgcCACCACC‑GGTGAGCGCCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCaa > 2:55460/1‑76 (MQ=255) gcCACCACC‑GGTGAGCGCCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACCGTTCCCGTGCaaac > 2:1587692/1‑74 (MQ=255) caccacC‑GGTGAGCGCCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCAt > 2:1449430/1‑76 (MQ=255) ccccc‑GGTGAGCGCCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTa < 2:1551708/73‑1 (MQ=255) c‑GGTGAGCGCCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTAcccc > 2:328811/1‑74 (MQ=255) gTGAGCGCCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCCGTGCAACCATTACCgcgc > 2:45322/1‑76 (MQ=255) gTGAGCGCCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCgcgc < 2:1408391/76‑1 (MQ=255) tGAGCGCCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCt > 2:37610/1‑76 (MQ=255) gAGCGCCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTc > 2:616009/1‑76 (MQ=255) gAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGGTc > 1:59271/1‑76 (MQ=255) | GACTGGAATGCCGGTTTCTTCCATTAACAAACGGTCAAGGTTACGCAGCAGTGCGCCACCACC‑GGTGAGCACCATGCCGCGCTCGGAGATGTCGGAAGCCAGTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTC > REL606/3328034‑3328175 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |