Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,623,312 | A→T | 19.2% | H159Q (CAT→CAA) | PP_4941 ← | hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,623,312 | 0 | A | T | 19.2% | 89.7 / 18.7 | 43 | H159Q (CAT→CAA) | PP_4941 | hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (17/17); new base T (3/5); total (21/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.00e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.26e-01 |
CACCACCCAGCTTTCGTCGGCCCGGCTGTACAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCT‑‑GCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCG > NC_002947/5623165‑5623461 | caccacCCAGCTTTCGTCGGCCCGGCTGTACAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCT‑‑GCAGCCGGCAATg < 2:118395/150‑1 (MQ=255) accaccCAGCTTTCGTCGGCCCGGCTGTACAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCT‑‑GCAGCCGGCAATgt > 2:691357/1‑149 (MQ=255) caccCAGCTTTCGTCGGCCCGGCTGTACAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCT‑‑GCAGCCGGCATTgccg > 1:236840/1‑150 (MQ=255) caccCAGCTTTCGTCGGCCCGGCTGTACAAACCATAAGCGGCCAACTCGGCGGCATCGGCATGTGGCGCGAGCACCACCCATCGAACCCGCTTCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCT‑‑GCAGCCGGCAATgccg > 1:778956/1‑150 (MQ=255) tttCGTCGGCCCGGCTGTACAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCT‑‑GCAGCCGGCATTgccggctgcagg < 2:236840/150‑7 (MQ=255) tttCGTCGGCCCGGCTGTACAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCT‑‑GCAGCCGGCAATGCCGGCCGCGCa < 2:763334/150‑1 (MQ=255) cgtcgGCCCGGCTGTACAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCGCGGCCGGCATTgccggctgcaggg < 1:600271/150‑8 (MQ=255) ggCCCGGCTGTACAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGTCCCT‑‑GCAGCCGGCAATGCCGGCCGCGCAGTCGCaa < 1:188227/150‑1 (MQ=255) ggCCCGGCTGTACAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCT‑‑GCAGCCGGCAATGTCGGCCGCGCAGTCGCaa < 1:691357/150‑1 (MQ=255) cTGTACAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCT‑‑GCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCg < 2:754841/150‑1 (MQ=255) tACAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCT‑‑GCAGCCGGCATTgccggctgca > 1:710244/1‑127 (MQ=255) tACAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGACGGGGTTATCCACAGGCGCAACTGACCCT‑‑GCAGCCGGCATTGCCGGCtgca < 2:710244/131‑2 (MQ=255) aaaCCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCT‑‑GCAGCCGGCAATGCCTGCCGCGCAGTCGCAACGTACCGGCCTGc > 2:640383/1‑150 (MQ=255) aCCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCT‑‑GCGGCCGGCATTgccggctgcagg > 1:319022/1‑122 (MQ=255) aCCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCT‑‑GCGGCCGGCATTgccggctgcagg < 2:319022/128‑7 (MQ=255) aCCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCT‑‑GCAGCCGGCaa > 1:396147/1‑115 (MQ=255) aCCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCT‑‑GCAGCCGGCaa < 2:396147/115‑1 (MQ=255) aTACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCCCCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGCCCCT‑‑GCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTgg < 1:436398/150‑1 (MQ=255) tACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCT‑‑GCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGgc < 2:549228/150‑1 (MQ=255) gcgGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGTCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCT‑‑GCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCgca > 1:581849/1‑150 (MQ=255) aaCTCGGCGTCATCGGCATGTGGAGCGAACCCCATCACACGGAGCCGGTGCAGAGCGACCAAGGGGGTTATACACAGGCGCAACTGAGGCT‑‑GCGGCCGACAGGGGAGGCAGGGCAGTCGCAACGTAACGGCCTGCAGTGGCCTAGCAGGa > 2:124752/1‑150 (MQ=255) cGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCT‑‑GCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCa < 2:77033/150‑1 (MQ=255) cGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCT‑‑GCAGCCGGCAATGTCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGttta < 2:493459/150‑3 (MQ=255) ggCATGTGGCGCGATCACCCTCCCCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCT‑‑GCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAg < 2:581849/150‑1 (MQ=255) ggCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCT‑‑GCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAgcgc > 1:480648/1‑150 (MQ=255) accatcaccCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCT‑‑GCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCg > 2:790437/1‑150 (MQ=255) catcaccCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCT‑‑GCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGAc > 1:18517/1‑150 (MQ=255) ccGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCT‑‑GCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACg > 2:87930/1‑150 (MQ=255) gcTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCT‑‑GCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTAc > 1:90766/1‑81 (MQ=255) gcTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCT‑‑GCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTAc < 2:90766/81‑1 (MQ=255) tCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCT‑‑GCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGa > 1:151234/1‑150 (MQ=255) gACCGAGGGGGTTATCCACAGGCGCAACTGACCCT‑‑GCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATc > 2:333570/1‑150 (MQ=255) gACCGAGGGGGTTATCCACAGGCGCAACTGACCCT‑‑GCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATc > 1:485302/1‑150 (MQ=255) aGGGGGTTATCCACAGGCGCAACTGACCCT‑‑GCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCTGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCt < 2:421818/150‑1 (MQ=255) ggTTATCCACAGGCGCAACTGACCCT‑‑GCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTc < 2:169063/150‑1 (MQ=255) ttATCCACAGGCGCAACTGACCCT‑‑GCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCgc > 1:605484/1‑150 (MQ=255) acaGGCGCAACTGACCCT‑‑GCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCgccagcca < 2:605484/150‑1 (MQ=255) acaGGCGCAACTGACCCT‑‑GCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCgccagcca > 1:226521/1‑150 (MQ=255) acaGGCGCAACTGACCCT‑‑GCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCgcca > 1:199201/1‑146 (MQ=255) acaGGCGCAACTGACCCT‑‑GCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCgcca < 2:199201/146‑1 (MQ=255) ggCGCAACTGACCCT‑‑GCAGCCGGCAATGCCGGCCGCCCGGCCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCcagca < 1:640383/150‑1 (MQ=255) ggACTCT‑‑GCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTCCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGAcc > 1:619009/2‑150 (MQ=255) gACCCT‑‑GCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGAccc > 2:136707/1‑150 (MQ=255) gcggCCGGCATTGCCGGCTGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGt < 1:370981/147‑1 (MQ=255) gCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCga < 1:333570/150‑1 (MQ=255) ccgagaaTGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGGTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCTCGGGTCgag < 1:124752/145‑1 (MQ=255) aTGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCTCGGGTCgagctgctg < 1:87930/150‑7 (MQ=255) | CACCACCCAGCTTTCGTCGGCCCGGCTGTACAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCT‑‑GCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGCCACGCTCGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGAAAGCG > NC_002947/5623165‑5623461 |
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: ‑ |