Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,310,453 | G→C | 9.1% | H362Q (CAC→CAG) | PP_1144 ← | membrane protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,310,453 | 0 | G | C | 9.1% | 165.6 / 7.4 | 55 | H362Q (CAC→CAG) | PP_1144 | membrane protein |
Reads supporting (aligned to +/- strand): ref base G (24/26); new base C (3/2); total (27/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.69e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01 |
GGTGGTCATGACCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTGGCCTGATGTTCGA > NC_002947/1310306‑1310601 | ggtggtCATGACCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTg < 1:887751/150‑1 (MQ=255) ggtggtCATGACCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTg < 1:847450/150‑1 (MQ=255) tCATGACCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGt > 1:946825/1‑144 (MQ=255) tCATGACCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGt < 2:946825/144‑1 (MQ=255) aTGACCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCgg > 2:231508/1‑150 (MQ=255) cATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTTCGGGTAGc > 2:364/1‑150 (MQ=255) cATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCgg > 2:259143/1‑145 (MQ=255) cATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCgg < 1:259143/145‑1 (MQ=255) cATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGc > 1:1078600/1‑150 (MQ=255) aTGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGATGGGCCAGGTGGCGGGCGTCGTTGCGGCTGTTGAGGCTGGTGTCCACAAAGAAGTCGTCGGGACCAAAGCGCGCAACCAGGGCGGGGTTGCGGGTAGCg > 2:222002/1‑150 (MQ=255) gcgcCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATTCGCTTGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGctgt < 1:1342467/150‑4 (MQ=255) gcgcCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGc < 1:272140/150‑1 (MQ=255) gcgcCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGc < 2:678917/150‑1 (MQ=255) gccgccCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCACGACCTGGTTGCgcgcttt < 2:123305/145‑7 (MQ=255) gccgccCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCACGACCTGGTTGCgcgcttt > 1:123305/1‑139 (MQ=255) cgcccgccCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTg > 1:633424/1‑150 (MQ=255) ccAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGt > 2:1003424/1‑150 (MQ=255) aGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCAGGGCCCGGAGGACCGCGTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTg > 1:602938/1‑150 (MQ=255) ggggTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCAc < 1:389329/150‑1 (MQ=255) tttCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGag < 1:1394626/150‑1 (MQ=255) aTCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGagaagc > 2:176681/1‑150 (MQ=255) tCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGGCCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTg > 2:1266243/1‑103 (MQ=255) tCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCGTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTg < 1:1266243/103‑1 (MQ=255) ccGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTc < 1:1003424/150‑1 (MQ=255) tGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTgg < 2:633424/150‑1 (MQ=255) aGGTGGCGGCCTTCGTCGCGGCTGTTCACGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAgg < 2:678373/150‑1 (MQ=255) ttCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGAt > 2:1164945/1‑150 (MQ=255) gTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGtt > 2:983203/1‑135 (MQ=255) gTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGtt < 1:983203/135‑1 (MQ=255) ttCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATgg > 1:270345/1‑150 (MQ=255) ggCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGa < 2:1008527/150‑1 (MQ=255) aCCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAgg > 1:1418209/1‑150 (MQ=255) aCCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAgg < 1:222002/150‑1 (MQ=255) cacaGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGa < 2:595399/120‑1 (MQ=255) cacaGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGa > 1:595399/1‑120 (MQ=255) aGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGa < 2:868469/150‑1 (MQ=255) cTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGc > 1:111487/1‑150 (MQ=255) cTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGc < 1:231508/150‑1 (MQ=255) gtcgCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAg > 2:1072433/1‑150 (MQ=255) gtcgCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAg > 2:687240/1‑150 (MQ=255) cgCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGc < 2:270345/150‑1 (MQ=255) accaAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGtt > 1:355497/1‑150 (MQ=255) accaAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGtt < 1:687240/150‑1 (MQ=255) accaAAGCGCGCAACCAGGTCGTGGATGTTGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGtt > 1:750031/1‑150 (MQ=255) aaaGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTc > 2:210484/1‑150 (MQ=255) aaaGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGATGGATATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTc < 2:750031/150‑1 (MQ=255) tacccgcAACCACGACCTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGc > 2:158182/5‑150 (MQ=255) aGCGCGCAACCAGGTCGTGGTTGAGGAAAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGc < 1:364/150‑1 (MQ=255) ccgcAACCACGACCTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGaaa < 1:249042/106‑1 (MQ=255) ccgcAACCACGACCTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGaaa > 2:249042/2‑107 (MQ=255) aaCCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTGGCCTg > 1:821028/1‑150 (MQ=255) gTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTg < 2:122536/94‑1 (MQ=255) gTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTg > 1:122536/1‑94 (MQ=255) gTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTGGCCTGATGTTc < 2:1078600/150‑1 (MQ=255) cGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTGGCCTGATGTTCGa < 2:602938/150‑1 (MQ=255) | GGTGGTCATGACCATGCGCCGCCCGCCCAGGTCGATGGGGTCTTTCATCCGTTGCATGATGCGCTGGGCCAGGTGGCGGGCTTCGTCGCGGCTGTTCAGGCTGGTGACCACACAGAACTCGTCGCCACCAAAGCGCGCAACCAGGTCGTGGTTGCGGGTAGCGGCCTTGATGTGGTTGGCGATCACCTTGAGCAGCTCGTCACCGGCGTCGTGGCCAAGGCTGTCGTTGATACGCTTGAAATGGTCGATGTCGAGGAACATGACAGCCAGGCGGTTTTCGCTGGCCTGATGTTCGA > NC_002947/1310306‑1310601 |
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: ‑ |