| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 1,842,635 | C→A | 16.7% | R734L (CGC→CTC) | gacS ← | sensor protein GacS |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 1,842,635 | 0 | C | A | 16.7% | 51.8 / 4.1 | 24 | R734L (CGC→CTC) | gacS | sensor protein GacS |
| Reads supporting (aligned to +/- strand): ref base C (8/12); new base A (2/2); total (10/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.59e-01 | |||||||||||
CACCTGGGCCAGTTGCCGTTCGCTGATCGGTTTGGTCAGGTAGTCGTCCATGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTCTGGGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAACAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGG > NC_002947/1842487‑1842749 | caccTGGGCCAGTTGCCGTTCGCTGATCGGTTTGGTCAGGTAGTCGTCCATGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTCTGGGTGTTCTCCCACAGGCg < 1:653648/150‑1 (MQ=255) ccAGTTGCCGTTCGCTGATCGGTTTGGTCAGGTAGTCGTCCATGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTCTGGGTGTTCTCCCACAGGCGGATCTGtt < 2:329622/150‑1 (MQ=255) tGCCGTTCGCTGATCGGTTTGGTCAGGTAGTCGTCCATGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTCTGGGTGTTCTCCCACAGGCGGATCTGTTCGGTa < 2:653301/150‑1 (MQ=255) ccGTTCGCTGATCGGTTTGGTCAGGTAGTCGTCCATGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTCTGGGTGTTCTCCCACAGGCGGATCTGTTCGGTACa > 1:487819/1‑150 (MQ=255) gTTTGGTCAGGTAGTCGTCCATGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTCTGGGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCa > 1:115863/1‑150 (MQ=255) aGTCGTCCATGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTCTGGGTGTTCTCCCACAGGCGGATCTGTTCGGTAAAGGCGCGGCCGTCCATGCCGGGCATCc < 1:312036‑M1/150‑2 (MQ=255) tgCAACAACGCGCGTTTTTCGTTGGCCATGGTATGGATCCTGAGGGCGACGATCGGCAGCGGGTTGCCGCTCTGGGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCcgcataaccattgttgac < 1:209701‑M1/150‑19 (MQ=255) gcgcgTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTCTGGGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCcgcataaccattgttgacggccagtac > 1:230252‑M1/1‑123 (MQ=255) gcgcgTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTCTGGGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCcgcataaccattgttgacggccagtac > 1:720362‑M1/1‑123 (MQ=255) cGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTCTGGGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCcgcataaccattgttgacggccagtacctcggcaccc < 2:188066‑M1/150‑38 (MQ=255) ttGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTCTGGGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGc < 2:81106/94‑1 (MQ=255) ttGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTCTGGGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGc > 1:81106/1‑94 (MQ=255) catggcatggGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTCTGGGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCcgcataaccattgttgacggccagtacctcggcacccaggtctt > 2:295109‑M1/1‑106 (MQ=255) atggGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTCTGGGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCcgcataaccattgttgacggccagtacctcggcacccaggtcttcgagca < 2:115863‑M1/150‑51 (MQ=255) tggGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTCTGGGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCcgcataaccattgttgacggccagtacctcggcacccaggtcttcgagcag < 1:102024‑M1/150‑52 (MQ=255) ggCGACGATCGGCAGCGGGTTGCCGCTCTGGGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCcgcataaccattgttgacggccagtacctcggcacccaggtcttcgagcagggtttttacca < 2:515248‑M1/150‑63 (MQ=255) gCGACGATCGGCAGCGGGTTGCCGCTCTGGGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCcgcataaccattgttgacggccagtacctcggcacccaggtcttcgagcagggtttttaccag < 2:230252‑M1/150‑64 (MQ=255) gCAGCGGGTTGCCGCTCTGGGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCcgcataaccattgttgacggccagtacctcggcacccaggtcttcgagcagggtttttaccagcagcaggttg > 1:394239‑M1/1‑77 (MQ=255) gtgtTCTCCCACAGGAGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCcgcataaccattgttgacggccagtacctcggcacccaggtcttcgagcagg > 2:156816‑M1/1‑57 (MQ=255) gtgtTCTCCCACAGGAGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCcgcataaccattgttgacggccagtacctcggcacccaggtcttcgagcagg > 2:157033‑M1/1‑57 (MQ=255) gtgtTCTCCCACAGGAGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCcgcataaccattgttgacggccagtacctcggcacccaggtcttcgagcagg < 1:157033‑M1/109‑53 (MQ=255) gtgtTCTCCCACAGGAGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCcgcataaccattgttgacggccagtacctcggcacccaggtcttcgagcagg < 1:156816‑M1/109‑53 (MQ=255) ccACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCcgcataaccattgttgacggccagtacctcggcacccaggtcttcga < 1:667944‑M1/96‑48 (MQ=255) ccACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCcgcataaccattgttgacggccagtacctcggcacccaggtcttcga > 2:667944‑M1/1‑49 (MQ=255) | CACCTGGGCCAGTTGCCGTTCGCTGATCGGTTTGGTCAGGTAGTCGTCCATGCCGCTGTGCAACAACGCGCGTTTTTCGTTGGCCATGGCATGGGCAGTGAGGGCGACGATCGGCAGCGGGTTGCCGCTCTGGGTGTTCTCCCACAGGCGGATCTGTTCGGTACAGGCGCGGCCGTCCATGCCGGGCATCTGCACATCCATCAGCACCAGGTCGAACAATTCTTCCTGCACGGCCTGTACCGCCGCATAACCATTGTTGACGG > NC_002947/1842487‑1842749 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |