Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,777,043 | C→G | 8.5% | A169P (GCA→CCA) | PP_2429 ← | membrane protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,777,043 | 0 | C | G | 8.5% | 144.1 / 4.7 | 47 | A169P (GCA→CCA) | PP_2429 | membrane protein |
Reads supporting (aligned to +/- strand): ref base C (20/23); new base G (2/2); total (22/25) | |||||||||||
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.39e-01 |
GAGGCGTTGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAAC > NC_002947/2776906‑2777185 | gAGGCGTTGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑g > 1:526824/1‑150 (MQ=255) aGGCGTTGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑gg < 2:146934/150‑1 (MQ=255) gACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCG‑CCCGGTGGCTGCAGGCACC‑ggggcggttg < 1:572030/150‑8 (MQ=255) acGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCa > 1:286081/1‑150 (MQ=255) acGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGATATTTTAGAGGCTCAGGCCAATACTGCAGCCCATGCCAATATAAGAAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCa > 1:462028/1‑150 (MQ=255) gTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCAcc > 2:509935/1‑150 (MQ=255) ggCCCCTTGTCTGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAggcgg < 1:797297/150‑1 (MQ=255) gCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑gg > 2:92382/1‑118 (MQ=255) gCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑gg < 1:92382/118‑1 (MQ=255) ggTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGc > 1:92450/1‑150 (MQ=255) aaGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCAGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCagag > 1:233892/1‑150 (MQ=255) tAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGgcgc < 2:162984/105‑1 (MQ=255) tAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGgcgc > 1:162984/1‑105 (MQ=255) agagGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGAt < 1:197910/150‑1 (MQ=255) agGCTCAGGCCAATACTGCAGCCCATGCCGATATCTACACATCCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCa < 2:462028/150‑1 (MQ=255) tCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCaa < 2:278341/150‑1 (MQ=255) tCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCAAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCaa < 2:286081/150‑1 (MQ=255) gCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGcc > 1:565249/1‑123 (MQ=255) gCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGcc < 2:565249/123‑1 (MQ=255) ccAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAaca < 1:502123/150‑1 (MQ=255) aaTACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACAcc > 1:455924/1‑150 (MQ=255) aCTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAgc < 1:509935/150‑1 (MQ=255) cccATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCt > 2:798060/1‑150 (MQ=255) agCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCAGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCca < 2:233892/149‑1 (MQ=255) cGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCacac > 2:531458/1‑150 (MQ=255) tCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAggcgg < 2:194436/96‑1 (MQ=255) tCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAggcgg > 1:194436/1‑96 (MQ=255) aacaaAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCAt > 2:54662/1‑98 (MQ=255) aacaaAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCAt < 1:54662/98‑1 (MQ=255) aacaaAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCTGCCGAGcatggc < 2:427714/115‑6 (MQ=255) aacaaAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCTGCCGAGcatggc > 1:427714/1‑110 (MQ=255) aaaaCCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAacgac > 2:273131/1‑150 (MQ=255) cccGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGcgcg > 1:205593/1‑150 (MQ=255) cAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCaaa < 2:205593/150‑1 (MQ=255) cTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAgg > 2:341200/1‑72 (MQ=255) cTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAgg < 1:341200/72‑1 (MQ=255) ttGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCaa < 1:724035/60‑1 (MQ=255) ttGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCaa > 2:724035/1‑60 (MQ=255) ttGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAacga < 1:88587/132‑1 (MQ=255) ttGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAacga > 2:88587/1‑132 (MQ=255) gTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGTCAATGACGGCGCGGCAAATACCCCGac < 1:419936/150‑1 (MQ=255) gTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGac < 2:526824/150‑1 (MQ=255) gcccGGTGGCTGCAGGCACCGGGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCa > 2:143709/2‑150 (MQ=255) ccccGGTGCCTACAGCCACC‑GGGCGCTCCCCACCCCCAACACCGTCAAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCACCAGCTGCCACACTGGCAACGACGACGCGGCAAATACCCCCACCACCGCCACGGCCATCaa < 2:77557/150‑1 (MQ=255) cccGGTGGCTGCAGGCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGc < 2:351288/119‑1 (MQ=255) cccGGTGGCTGCAGGCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGc > 1:351288/1‑119 (MQ=255) cccGGTGCCTGCAGCCACC‑GGGCGCGCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAAc < 1:491792/150‑1 (MQ=255) | GAGGCGTTGACACGTTGAGCCGAGGCCGCTCAGGCCCCTTGTCGGTCAAGTAGAGGCTCAGGCCAATACTGCAGCCCATGCCGATATCAACAAAACCCCCGCCAACACTTGGTTGAACACCCGCAACCGCCCCGGTGCCTGCAGCCACC‑GGGCGCTCCCCACCCCCAACACCGCCCAGGCGGCCATGCTCGGCAGAGCGATCAGCAAGAACACCAGCGCCAGCTGCCACACTGGCAACGACGGCGCGGCAAATACCCCGACCACCGCCACGGCCATCAAC > NC_002947/2776906‑2777185 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |