Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,115,751 | T→G | 31.1% | H28P (CAC→CCC) | PP_4500 ← | TrpH family protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,115,751 | 0 | T | G | 31.1% | 17.4 / 11.5 | 19 | H28P (CAC→CCC) | PP_4500 | TrpH family protein |
Reads supporting (aligned to +/- strand): ref base T (6/7); new base G (3/3); total (9/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 = 1.00e+00 |
ATGGTCGCGCCGCCCCACGTGCAGGACAGCTCCACCCCGCTGACCCAGCGCATACCCAGGTCAACGCAGGCCTGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCATAGAGGAGCGCTTTCGCTGGATTCGATGTTTGTTATTATGCCGTCACATCCCGATTCTGG > NC_002947/5115616‑5115892 | atggTCGCGCCGCCCCACGTGCAGGACAGCTCCACCCCGCTGACCCAGCGCATACCCAGGTCAACGCAGGCCTGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGGGCTCATGAGCACgg < 2:202725/150‑1 (MQ=255) cACGTGCAGGACAGCTCCACCCCGCTGACCCAGCGCATACCCAGGTCAACGCAGGCCTGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGaa < 1:439271/150‑1 (MQ=255) aGGACAGCTCCACCCCGCTGACCCAGCGCATACCCAGGTCAACGCAGGCCTGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGGGCTCATGAGCACggggt > 2:28071/1‑129 (MQ=255) aGGACAGCTCCACCCCGCTGACCCAGCGCATACCCAGGTCAACGCAGGCCTGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGGGCTCATGAGCACggggt < 1:28071/131‑3 (MQ=255) tCCACCCCGCTGACCCAGCGCATACCCAGGTCAACGCAGGCCTGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGGGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCg > 2:444318/1‑150 (MQ=255) cccGCTGACCCAGCGCATACCCAGGTCAACGCAGGCCTGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGa > 1:584765/1‑150 (MQ=255) aCCCAGCGCATACCCAGGTCAACGCAGGCCTGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTg > 1:81511/1‑150 (MQ=255) ccAGCGCATACCCAGGTCAACGCAGGCCTGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAAc > 1:23692/1‑111 (MQ=255) ccAGCGCATACCCAGGTCAACGCAGGCCTGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAAc < 2:23692/111‑1 (MQ=255) aGCGCATACCCAGGTCAACGCAGGCCTGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCtgt < 1:609658/150‑1 (MQ=255) aGCGCATACCCAGGTCAACGCAGGCCTGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCtgt > 1:482103/1‑150 (MQ=255) gTCAACGCAGGCCTGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATc < 2:584765/150‑1 (MQ=255) aGGCCTGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGGGCTCATGAGc > 1:419238/1‑79 (MQ=255) aGGCCTGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGGGCTCATGAGc < 2:419238/79‑1 (MQ=255) cTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCa > 2:247837/1‑126 (MQ=255) cTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCa < 1:247837/126‑1 (MQ=255) ggggtCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCATAGAGGAGCGCTTTCGCTGGACTCGAtgttt < 2:182615/147‑1 (MQ=255) ggTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCATAGAGGAGCGCTTTCGCTGGATTCGAtgttt < 2:489513/150‑1 (MQ=255) gCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCATAGAGGAGCGCTTTCGCTGGATTCGATGTTTGTTATTATGCCGTCACATCCCGATTCTgg > 1:129389/1‑150 (MQ=255) | ATGGTCGCGCCGCCCCACGTGCAGGACAGCTCCACCCCGCTGACCCAGCGCATACCCAGGTCAACGCAGGCCTGGCGCGCTTCGGGCAGGCCTTCGAGGGTGTCATGGTCGGTCAGTGCCAGCGCTTGCACCCCGTGCTCATGAGCCCGGGCAACCAGTACCGAAGGCGACAGGGCGCCGTCGGAGGCCGTGCTGTGACAGTGCAGATCAACATTCATAGAGGAGCGCTTTCGCTGGATTCGATGTTTGTTATTATGCCGTCACATCCCGATTCTGG > NC_002947/5115616‑5115892 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |