| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 5,378,931 | C→G | 11.8% | A156A (GCC→GCG) | recN → | DNA repair/recombination protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 5,378,931 | 0 | C | G | 11.8% | 100.5 / 5.3 | 34 | A156A (GCC→GCG) | recN | DNA repair/recombination protein |
| Reads supporting (aligned to +/- strand): ref base C (17/13); new base G (2/2); total (19/15) | |||||||||||
| 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.69e-01 | |||||||||||
ATCAACGGCACGCCCTGCCCGCTCGGCGACCTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAACCTCGGCCTGGGCGA > NC_002947/5378782‑5379074 | aTCAACGGCACGCCCTGCCCGCTCGGCGACCTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGcc > 1:785349/1‑150 (MQ=255)aTCAACGGCACGCCCTGCCCGCTCGGCGACCTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGcc > 2:260977/1‑150 (MQ=255) ggCACGCCCTGCCCGCTCGGCGACCTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAg > 2:301118/1‑150 (MQ=255) cACGCCCTGCCCGCTCGGCGACCTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCGGGCCAGGt < 2:790582/150‑1 (MQ=255) aCGCCCTGCCCGCTCGGCGACCTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTa < 1:301118/150‑1 (MQ=255) gCCCTGCCCGCTCGGCGACCTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCCGGCCCGCCAGGTACa > 1:621439/1‑150 (MQ=255) ccGCTCGGCGACCTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCa < 1:260977/150‑1 (MQ=255) gCTCGGCGACCTCAAGGCGCTGGGCGAACTGCTCAGCGTTATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGc > 1:265338/1‑150 (MQ=255) gtctgCGACCTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGCCACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGcc > 2:137085/5‑150 (MQ=255) cGGCGACCTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCCGGCCCGCCAGGTACAGCTTGCAGCCaa < 2:621439/150‑1 (MQ=255) gCGACCTCAAGGCGCTGGGCGAACTGCTGATTGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAgc < 1:481718/150‑1 (MQ=255) tcccatcTGGGCTAACAGCTGATCTATAGCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCACCGTCGCCTGCATGACGAATATGCCCGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAAt < 2:39166/144‑1 (MQ=255) aaGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAAt > 2:759245/1‑150 (MQ=255) tGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCACCGTCGCCTGCTTGACGAGTATGCCGTCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGAccc > 1:407646/1‑150 (MQ=255) ggCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGAAATCAGatcgga > 1:39166/1‑145 (MQ=255) aaCTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCctgg > 2:420761/1‑150 (MQ=255) tgATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCctggagctgg > 1:210430/1‑150 (MQ=255) aTCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCctggagctggag > 1:312118/1‑150 (MQ=255) cGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAgcg > 2:103007/1‑150 (MQ=255) acaGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCa > 2:249388/1‑150 (MQ=255) acaGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCa > 1:156060/1‑150 (MQ=255) cagccagcACGAGCACCAGTCCCTGCTCAAGACCGACACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCaa < 1:249388/150‑1 (MQ=255) cagcACGAGCACCAGTCCCTGCTCAAGACCGACACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCg < 1:759245/150‑1 (MQ=255) gcACGAGCACCAGTCCCTGCTCAAGACCGACACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCggg < 1:420761/150‑1 (MQ=255) tCCCTGCTCAAGACCGACACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAgcgc > 2:528453/1‑150 (MQ=255) ccTGCTCAAGACCGACACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCGGGCCAGGtcga < 2:341217/75‑4 (MQ=255) ccTGCTCAAGACCGACACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCGGGCCAGGtcga > 1:341217/1‑72 (MQ=255) cccaccGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCa < 2:210430/150‑1 (MQ=255) cccaccGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCa < 1:528453/150‑1 (MQ=255) cccaccGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCa < 1:103007/150‑1 (MQ=255) caccGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCa > 1:671920/1‑80 (MQ=255) caccGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCa < 2:671920/80‑1 (MQ=255) cGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACt < 2:156060/150‑1 (MQ=255) gAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTg > 1:724047/1‑150 (MQ=255) ctgtaCCTGGCGGGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAACCTCGGCCTg > 1:707171/5‑150 (MQ=255) ccTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTAGGGCCGCTAACTGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAACCTCGGCCTGGGCGa > 1:412989/1‑150 (MQ=255) | ATCAACGGCACGCCCTGCCCGCTCGGCGACCTCAAGGCGCTGGGCGAACTGCTGATCGATATCCACAGCCAGCACGAGCACCAGTCCCTGCTCAAGACCGACACCCACCGTCGCCTGCTTGACGAGTATGCCGGCGCCGTCGACCTGGCCCGCCAGGTACAGCTTGCAGCCAAGCGCTGGAATCAGACCCGCCTGGAGCTGGAGCGCCTCTCCAACTCGGGCGACGAGCAGCGCGCCCGCCACCAGCTGCTCAGTTACCAGCTGGAAGAACTGGACAACCTCGGCCTGGGCGA > NC_002947/5378782‑5379074 |
| 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: ‑ |