| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 4,504,966 | G→C | 31.6% | V10L (GTG→CTG) | tusE → | sulfur transfer protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 4,504,966 | 0 | G | C | 31.6% | 21.3 / 14.1 | 19 | V10L (GTG→CTG) | tusE | sulfur transfer protein |
| Reads supporting (aligned to +/- strand): ref base G (5/8); new base C (4/2); total (9/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.50e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01 | |||||||||||
CGCCCTGGACGAGGATGTACAGGCCCGCGCGGTCGGCAACGACCTGGCAAAGGCCGTGGACTACGCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGAATAC > NC_002947/4504825‑4505115 | cGCCCTGGACGAGGATGTACAGGCCCGCGCGGTCGGCAACGACCTGGCAAAGGCCGTGGACTACGCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGcc < 1:121965/150‑1 (MQ=255) taGGCCCGCGCGGTCGGCAACGACCTGGCAAAGGCCGTGGACTACGCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTcc < 1:263460/149‑1 (MQ=255) ccGCGCGGTCGGCAACGACCTGGCAAAGGCCGTGGACTACGCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTc > 2:579584/1‑150 (MQ=255) aCGACCTGGCAAAGGCCGTGGACTACGCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTg < 2:430865/150‑1 (MQ=255) aCCTGGCAAAGGCCGTGGACTACGCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGtc > 2:337018/1‑150 (MQ=255) aaaGGCCGTGGACTACGCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCacgac > 2:526316/1‑150 (MQ=255) cGCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGATGCCGCCGAGGCTTTg > 2:7102/1‑150 (MQ=255) gTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAg < 1:579584/150‑1 (MQ=255) tCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATcc < 1:397629/150‑1 (MQ=255) tGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGtt < 1:385892/150‑1 (MQ=255) tGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCACCTGATCGccaa < 2:345749/66‑3 (MQ=38) tGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCACCTGATCGccaa > 1:345749/1‑64 (MQ=38) cAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTgg > 1:100653/1‑150 (MQ=255) aaCAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTggg < 2:515151/150‑1 (MQ=255) cagggCGATCACCTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTc > 1:466603/4‑150 (MQ=255) agggCGATCACCTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGAcggc > 1:464228/3‑112 (MQ=255) agggCGATCACCTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGAcggc < 2:464228/110‑1 (MQ=255) gggCGATCACCTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTa > 1:513293/2‑150 (MQ=255) gATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAg < 1:235725/150‑1 (MQ=255) gTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCac < 1:44776/54‑1 (MQ=255) gTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCac > 2:44776/1‑54 (MQ=255) gTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGAATAc > 1:65604/1‑150 (MQ=255) | CGCCCTGGACGAGGATGTACAGGCCCGCGCGGTCGGCAACGACCTGGCAAAGGCCGTGGACTACGCCGGGTTCGTCGAACTGTCGCTGCAATACGACAAGGTCAACAGCTGGCTATGAATACGCTCACCGTTGGCGATCAGGTGATCGCCCTGGACAAGGATGGCTTCCTGGTCGATCTGCAAGATTGGTCTCACGACGCCGCCGAGGCTTTGGCCAAGCGCGAGGGTATCCCGTTGACGGCGGACCACTGGGAAATCCTCGAACTGCTCCGTCAGTTCTACCAGGAATAC > NC_002947/4504825‑4505115 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |