Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,185,604 | C→G | 16.6% | L446L (CTC→CTG) | PP_3681 → | putative Helicase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,185,604 | 0 | C | G | 16.6% | 100.5 / 10.5 | 41 | L446L (CTC→CTG) | PP_3681 | putative Helicase |
Reads supporting (aligned to +/- strand): ref base C (16/18); new base G (4/3); total (20/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.97e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CTTATGAGGACTGGTGGAAGCTCACAAGCAGCTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCGGGGTCGAAACCCAAGCTTTCTGGG > NC_002947/4185459‑4185743 | cTTATGAGGACTGGTGGAAGCTCACAAGCAGCTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTcagca < 1:1351830/150‑1 (MQ=255) gtggAAGCTCACAAGCAGCTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGAc < 2:454055/147‑1 (MQ=255) gtggAAGCTCACAAGCAGCTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGAc > 1:454055/1‑147 (MQ=255) aaGCTCACAAGCAGCTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTccc > 1:1102835/1‑150 (MQ=255) aGCTCACAAGCAGCTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCg > 2:499171/1‑150 (MQ=255) gCTCACAAGCAGCTTCATCCGTGATGTAGGGTATCCAATCCTGACTGCGCTCTCAGGTGGCTACGAGTCGAAACGACGCTGGCGTGAGGAGGGTCGTAAGGGCAGCGATAGACGTCGAGCGCAGCTCACCGGATAAACGCCCATTccctg > 2:162528/1‑148 (MQ=255) cagcTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGCTCGAGATGGAGCGCTGCTCATCAGCGAAAAGACCATTCCCGGTGCCCTTGc < 1:743553/150‑1 (MQ=255) cATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCt < 1:179239/150‑1 (MQ=255) aTCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCtt > 2:363893/1‑150 (MQ=255) tCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCttt > 2:579279/1‑150 (MQ=255) cGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCa < 2:1102835/150‑1 (MQ=255) cGTGATGTAGATTCTCCCATCCTGACTGCGCTCTCAGGTGACTACGCGTCGACCAAACGCAGTTACCAATGCTTTCGTGCGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGATGCCCTTGCCCGGCTTTCa < 1:878621/150‑1 (MQ=255) tGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGGGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGa < 1:882018/150‑1 (MQ=255) ggTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTccc < 1:546761/120‑1 (MQ=255) ggTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTccc > 2:546761/1‑120 (MQ=255) aaTCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTcagcag < 2:875575/96‑1 (MQ=255) aaTCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTcagcag > 1:875575/1‑96 (MQ=255) aaTCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACg > 2:1017241/1‑150 (MQ=255) tCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCAtt > 1:800952/1‑107 (MQ=255) tCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCAtt < 2:800952/107‑1 (MQ=255) cTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTca < 1:579279/150‑1 (MQ=255) tctcAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTGGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATcc > 1:1357944/1‑150 (MQ=255) aGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCa < 2:1266178/123‑1 (MQ=255) aGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCa > 1:1266178/1‑123 (MQ=255) gTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTGAGCAGCgc < 2:1274158/79‑2 (MQ=255) gTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTGAGCAGCgc > 1:1274158/1‑78 (MQ=255) gTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACaaaa < 1:1017241/150‑1 (MQ=255) cGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAggg < 1:363893/150‑1 (MQ=255) gAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGa > 2:332675/1‑150 (MQ=255) gCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGt < 2:1357944/150‑1 (MQ=255) aCGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGtca > 1:575387/1‑150 (MQ=255) cGAATTCGTTAGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGtcat < 2:640108/150‑1 (MQ=255) gTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCAtt > 1:542815/1‑150 (MQ=255) cgtGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCggg < 1:499171/150‑1 (MQ=255) cAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCGGGGTCGAAAc < 2:542815/150‑1 (MQ=255) agagATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAg < 2:1125009/108‑1 (MQ=255) agagATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAg > 1:1125009/1‑108 (MQ=255) aGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCGGGGTCGAAACCCAAGCttt > 2:1414761/1‑150 (MQ=255) tcttttcgcTGCTGAGCAGCGCAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCGGGGTCGAAACCCAAGCTTTCt > 1:202435/7‑150 (MQ=255) cttttcgcTGCTGAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCGGGGTCGAAACCCAAGCTTTCTg > 2:980261/6‑150 (MQ=255) gggggctgctGAGCAGCGAAAAGACCATTTCCGGTGCGGTTCGCCTGCCTTCAGACGATCAGGCCGTACGGTTTCTCATACCCCACAAAAGCAAGGGGCGTGAGTTTCACACGGTCATCATTCTCGGGGTCGAAACCCAAGCTTTCTggg < 2:477552/146‑1 (MQ=255) gAGCTCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCGGGGTCGAAACCCAAGCTTTCTggg < 1:162528/150‑1 (MQ=255) ttcgcTGCTGAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAgg < 2:205882/60‑1 (MQ=255) ttcgcTGCTGAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAgg > 1:205882/3‑62 (MQ=255) | CTTATGAGGACTGGTGGAAGCTCACAAGCAGCTTCATCCGTGATGTAGGTTATCCAATCCTGACTGCGCTCTCAGGTGACTACGAGTCGAAACTACGCTTGCACGAAGTCGTTCGTGAGGTCAGAGATCGAGTTGGAGCGCTGCTCAGCAGCGAAAAGACCATTCCCGGTGCCCTTGCCCGGCTTTCAGACGATCAGGCCGTACGCTTTCTCACAATCCACAAAAGCAAGGGGCTTGAGTTTCACACGGTCATCATTCTCGGGGTCGAAACCCAAGCTTTCTGGG > NC_002947/4185459‑4185743 |
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: ‑ |