Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,533,863 | A→G | 11.7% | G220G (GGT→GGC) | PP_4867 ← | BraC‑like branched‑chain amino acid ABC transporter substrate‑binding protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,533,863 | 0 | A | G | 11.7% | 81.9 / 2.7 | 34 | G220G (GGT→GGC) | PP_4867 | BraC‑like branched‑chain amino acid ABC transporter substrate‑binding protein |
Reads supporting (aligned to +/- strand): ref base A (15/15); new base G (2/2); total (17/17) | |||||||||||
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.32e-01 |
ATTGCGCGCCGCCGGCGGTAGACACCAGTTCGTCGGTGACGATGCCATCATCGGACATGAACTTGACGTCCTTCAGGCCTTGCTCGCGCAGCTGGCGTACCAGCGGGCCGGCTTCCGGGTGCAGGCCGCCGAAGTAGACCACGTCGGCACCGGTAGAGCGGATCTTGGTGACCACCGCGCTGAAGTCTTTCTCGCCACGGGTCAGGCCTTCGTACAGCACAGGCTTCACGCCGCGTTTTTCCAGCTGCGCCTTGGTCGCGTCAGCCAGGCCCTGGCCGTAGGTGTCCTTGTCGTGCAA > NC_002947/5533715‑5534012 | aTTGCGCGCCGCCGGCGGTAGACACCAGTTCGTCGGTGACGATGCCATCATCGGACATGAACTTGACGTCCTTCAGGCCTTGCTCGCGCAGCTGGCGTACCAGCGGGCCGGCTTCCGGGTGCAGGCCGCCGAAGTAGACCACGTCGGCAc < 2:465366/150‑1 (MQ=255) gcgcgcCGCCGGCGGTAGACACCAGTTCGTCGGTGACGATGCCATCATCGGACATGAACTTGACGTCCTTCAGGCCTTGCTCGCGCAGCTGGCGTACCAGCGGGCCGGCTTCCGGGTGCAGGCCGCCGAAGTAGACCACGTCGGCACCgg > 2:757464/1‑150 (MQ=255) gccgGCGGTAGACACCAGTTCGTCGGTGACGATGCCATCATCGGACATGAACTTGACGTCCTTCAGGCCTTGCTCGCGCAGCTGGCGTACCAGCGGGCCGGCTTCCGGGTGCAGGCCGCCGAAGTAGACCACGTCGGCACCGGTAGAGCg < 2:721393/150‑1 (MQ=255) gTTCGTCGGTGACGATGCCATCATCGGACATGAACTTGACGTCCTTCAGGCCTTGCTCGCGCAGCTGGCGTACCAGCGGGCCGGCTTCCGGGTGCAGGCCGCCGAAGTAGACCACGTCGGCACCGGTAGAGCGGATCTTGGTGACCACcg > 2:32129/1‑150 (MQ=255) atcatcGGACATGAACTTGACGTCCTTCAGGCCTTACTCGCGCAGCTGACGTACCAGCGGTCCGGCTTCCGGGTGCAGGCCGCCGAAGTAGACCACGTCGGCACCGGTAGAGCGGCTCTTGGTGACCACCGCACTGGAGTCTTTCTCGcc > 2:128923/1‑150 (MQ=255) tcatcGGACATGAACTTGACGTCCTTCAGGCCTTGCTCGCGCAGCTGGCGTACCAGCGGGCCGGCTTCCGGGTGCAGGCCGCCGAAGTAGACCACGTCGGCACCGGTAGAGCGGATCTTGGTGACCACCGCGCTGAAGTCTTTCTCGCCa < 2:24475/150‑1 (MQ=255) cATGAACTTGACGTCCTTCAGGCCTTGCTCGCGCAGCTGGCGTACCAGCGGGCCGGCTTCCGGGTGCAGGCCGCCGAAGTAGACCACGTCGGCACCGGTAGAGCGGATCTTGGTGACCACCGCGCTGAAGTCTTTCTCGCCACGGGTCAg < 2:369445/150‑1 (MQ=255) cTTCAGGCCTTGCTCGCGCAGCTGGCGTACCAGCGGGCCGGCTTCCGGGTGCAGGCCGCCGAAGTAGACCACGTCGGCACCGGTAGAGCGGATCTTGGTGACCACCGCGCTGAAGTCTTTCTCGCCACGGGTCAGGCCTTCGTACAGcac > 1:218466/1‑150 (MQ=255) ttCAGGCCTTGCTCGCGCAGCTGGCGTACCAGCGGGCCGGCTTCCGGGTGCAGGCCGCCGAAGTAGACCACGTCGGCGCCGGTAGAGCGGATCTTGGTGACCACCGCGCTGAAGTCTTtctc > 2:487853/1‑122 (MQ=255) ttCAGGCCTTGCTCGCGCAGCTGGCGTACCAGCGGGCCGGCTTCCGGGTGCAGGCCGCCGAAGTAGACCACGTCGGCGCCGGTAGAGCGGATCTTGGTGACCACCGCGCTGAAGTCTTtctc < 1:487588/122‑1 (MQ=255) ttCAGGCCTTGCTCGCGCAGCTGGCGTACCAGCGGGCCGGCTTCCGGGTGCAGGCCGCCGAAGTAGACCACGTCGGCGCCGGTAGAGCGGATCTTGGTGACCACCGCGCTGAAGTCTTtctc < 1:487853/122‑1 (MQ=255) ttCAGGCCTTGCTCGCGCAGCTGGCGTACCAGCGGGCCGGCTTCCGGGTGCAGGCCGCCGAAGTAGACCACGTCGGCGCCGGTAGAGCGGATCTTGGTGACCACCGCGCTGAAGTCTTtctc > 2:487588/1‑122 (MQ=255) ttCAGGCCTTGCTCGCGCAGCTGGCGTACCAGCGGGCCGGCTTCCGGGTGCAGGCCGCCGAAGTAGACCACGTCGGCACCGGTAGAGCGGATCTTGGTGACCACCGCGCTGAAGTCTTTCTCGCCACGGGTCAGGCCTTCGTACAGcaca > 2:362176/1‑150 (MQ=255) gCCTTGCTCGCGCAGCTGGCGTACCAGCGGGCCGGCTTCCGGGTGCAGGCCGCCGAAGTAGACCACGTCGGCACCGGTAGAGCGGATCTTGGTGACCACCGCGCTGAAGTCTTTCTCGCCACGGGTCAGGCCTTCGTACAGCACAGGCtt > 2:715138/1‑150 (MQ=255) ttGCTCGCGCAGCTGGCGTACCAGCGGACCGGCTTCCGGGTGCAGGCCGCCGAAGTAGACCACGTCGGCACCGGTAGAGCGGATCTTGGTGACCACCGCGCTGAAGTCTTTCTCGCCACTGGTCAGGCCTTCGTACAGCACAGGCTTCAc < 1:748457/150‑1 (MQ=255) aGCTGGCGTACCAGCGGGCCGGCTTCCGGGTGCAGGCCGCCGAAGTAGACCACGTCGGCACCGGTAGAGCGGATCTTGGTGa < 2:101014/82‑1 (MQ=255) aGCTGGCGTACCAGCGGGCCGGCTTCCGGGTGCAGGCCGCCGAAGTAGACCACGTCGGCACCGGTAGAGCGGATCTTGGTGa > 1:101014/1‑82 (MQ=255) aGCTGGCGTACCAGCGGGCCGGCTTCCGGGTGCAGGCCGCCGAAGTAGACCACGTCGGCACCGGTAGAGCGGATCTTGGTGACCACCGCGCTGAAGTCTTTCTCGCCACGGGTCAGGCCTTCGTACAGCACAGGCTTCACGCCGCGtttt > 1:774412/1‑150 (MQ=255) aCCAGCGGGCCGGCTTCCGGGTGCAGGCCGCCGAAGTAGACCACGTCGGCACCGGTAGAGCGGATCTTGGTGACCACCGCGCTGAAGTCTTTCTCGCCACGGGTCAGGCCTTCGTACAGCACAGGCTTCACGCCGCGTTTTTCCAGCTgc < 1:757464/150‑1 (MQ=255) aGCGGGCCGGCTTCCGGGTGCAGGCCGCCGAAGTAGACCACGTCGGCACCGGTAGAGCGGATCTTGGTGACCACCGCGCTGAAGTCTTTCTCGCCACGGGTCAGGCCTTCGTACAGCACAGGCTTCACGCCGCGTTTTTCCAGCTGCGcc > 1:849761/1‑150 (MQ=255) gCCGGCTTCCGGGTGCAGGCCGCCGAAGTAGACCACGTCGGCACCGGTAGAGCGGATCTTGGTGACCACCGCGCTGAAGTCTTTCTCGCCACGGGTCAGGCCTTCGTACAGCACAGGCTTCACGCCGCGTTTTTCCAGCTGCGCCTTGGt > 1:789532/1‑150 (MQ=255) atCCGGGTGCAGGCCGCCAAAGGAGACAACGTAGGCACCGGTAGAGCGGATCTTGGTGACAACCGAGCAGAAGTATTTATCAACAAGGGTTAGGCCTTCATACAGCACAGGCTTCACTCCGCGTCTTTCCAGCTGCGCCATTGTCGCTTc > 2:364064/2‑150 (MQ=255) tCCGGGTGCAGGCCGCCGAAGTAGACCACGTCGGCACCGGTAGAGCGGATCTTGGTGACCACCGCGCTGAAGTCTTTCTCGCCACGGGTCAGGCCTTCGTACAGCACAGGCTTCACGCCGCGTTTTTCCAGCTGCGCCTTGGTcttttca < 2:849761/150‑7 (MQ=255) cGGGTGCAGGCCGCCGAAGTAGACCACGTCGGCACCGGTAGAGCGGATCTTGGTGACCACCGCGCTGAAGTCTTTCTCGCCACGGGTCAgg < 2:215696/91‑1 (MQ=255) cGGGTGCAGGCCGCCGAAGTAGACCACGTCGGCACCGGTAGAGCGGATCTTGGTGACCACCGCGCTGAAGTCTTTCTCGCCACGGGTCAgg > 1:215696/1‑91 (MQ=255) ggTGCAGGCCGCCGAAGTAGACCACGTCGGCACCGGTAGAGCGGATCTTGGTGACCACCGCGCTGAAGTCTTTCTCGCCACGGGTCAGGCCTTCGTACAGCACAGGCTTCACGCCGCGTTTTTCCAGCTGCGCCTTGGTCGCGTCAGCCa < 1:32129/150‑1 (MQ=255) gTGCAGGCCGCCGAAGTAGACCACGTCGGCACCGGTAGAGCGGATCTTGGTGACCACCGCGCTGAAGTCTTTCTCGCCACGGGTCAGGCCTTCGTACAGCACAGGCTTCACGCCGCGTTTTTCCAGCTGCGCCTTGGTCGCGTCAGCCAg < 2:789532/150‑1 (MQ=255) gccgccGAAGTAGACCACGTCGGCACCGGTAGAGCGGATCTTGGTGACCACCGCGCTGAAGTCTTTCTCGCCACGGGTCAGGCCTTCGTACAGCACAGGCTTCACGCCGCGTTTTTCCAGCTGCGCCTTGGTCGCGTCAGCCAGGCCCTg > 1:831326/1‑150 (MQ=255) aaGTAGACCACGTCGGCACCGGTAGAGCGGATCTTGGTGACCACCGCGCTGAAGTCtt < 2:215931/58‑1 (MQ=255) aaGTAGACCACGTCGGCACCGGTAGAGCGGATCTTGGTGACCACCGCGCTGAAGTCtt > 1:215931/1‑58 (MQ=255) aaGTAGACCACGTCGGCACCGGTAGAGCGGATCTTGGTGACCACCGCGCTGAAGTCTTTCTCGCCACGGGTCAGGCCTTCGTACAGCACAGGCTTCACGCCGCGTTTTTCCAGCTGCGCCTTGGTCGCGTCAGCCAGGCCCTGGCCGTAg < 1:362176/150‑1 (MQ=255) gACCACGTCGGCACCGGTAGAGCGGATCTTGGTGACCACCGCGCTGAAGTCTTTCTCGCCACGGGTCAGGCCTTCGTACAGCACAGGCTTCACGCCGCGTTTTTCCAGCTGCGCCTTGGTCGCGTCAGCCAGGCCCTGGCCGTAGGTGTc < 1:715138/150‑1 (MQ=255) aCGTCGGCACCGGTAGAGCGGATCTTGGTGACCACCGCGCTGAAGTCTTTCTCGCCACGGGTCAGGCCTTCGTACAGCACAGGCTTCACGCCGCGTTTTTACAGCTGCGCCTTGGTCGCGTCAGCCAGGCCCTGGCCGTAGGTGTCCTTg > 1:746146/1‑150 (MQ=255) gTCGGCACCGGTAGAGCGGATCTTGGTGACCACCGCGCTGAAGTCTTTCTCGCCACGGGTCAGGCCTTCGTACAGCACAGGCTTCACGCCGCGTTTTTCCAGCTGCGCCTTGGTCGCGTCAGCCAGGCCCTGGCCGTAGGTGTCCTTGTc < 2:218466/150‑1 (MQ=255) aCCGGTAGAGCGGATCTTGGTGACCACCGCGCTGAAGTCTTTCTCGCCACGGGTCAGGCCTTCGTACAGCACAGGCTTCACGCCGCGTTTTTCCAGCTGCGCCTTGGTCGCGTCAGCCAGGCCCTGGCCGTAGGTGTCCTTGTCGTGCaa > 1:445899/1‑150 (MQ=255) | ATTGCGCGCCGCCGGCGGTAGACACCAGTTCGTCGGTGACGATGCCATCATCGGACATGAACTTGACGTCCTTCAGGCCTTGCTCGCGCAGCTGGCGTACCAGCGGGCCGGCTTCCGGGTGCAGGCCGCCGAAGTAGACCACGTCGGCACCGGTAGAGCGGATCTTGGTGACCACCGCGCTGAAGTCTTTCTCGCCACGGGTCAGGCCTTCGTACAGCACAGGCTTCACGCCGCGTTTTTCCAGCTGCGCCTTGGTCGCGTCAGCCAGGCCCTGGCCGTAGGTGTCCTTGTCGTGCAA > NC_002947/5533715‑5534012 |
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: ‑ |