Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 966,701 | G→T | 9.1% | R7S (CGC→AGC) | PP_0828 ← | membrane protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 966,701 | 0 | G | T | 9.1% | 127.6 / 3.3 | 45 | R7S (CGC→AGC) | PP_0828 | membrane protein |
Reads supporting (aligned to +/- strand): ref base G (19/21); new base T (2/2); total (22/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.77e-01 |
CCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATGCAGGATGGGACGTCGAGGCTCGTCGAAGGTTCGAGGGTTGAATGTTGCAAGGGATTTTGTAGCTGAAACGAAGGGAAAAATCCGACCGGCCTATCAGGT > NC_002947/966553‑966846 | cccAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCgg < 2:486107/150‑1 (MQ=255) tCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCAcgc > 2:223091/1‑150 (MQ=255) tGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCgg > 1:571533/1‑127 (MQ=255) tGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCgg < 2:571533/127‑1 (MQ=255) gCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTACCCGGGATGACGTCTTTCACGCGAGGGt > 2:835493/1‑150 (MQ=255) ggTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCt > 1:580752/1‑150 (MQ=255) tttcatGCTGTTCTTGAACAGCAAGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCGTTCACGCGAGAGGACATCGa > 2:361208/6‑150 (MQ=255) ttcatGCTGTTCTTGAACAGCAAGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGaa < 1:361208/146‑1 (MQ=255) ttCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCggggg < 1:223091/150‑1 (MQ=255) tCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGa < 1:22718/150‑1 (MQ=255) tgAACAGCAAGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATGc < 1:533452/149‑1 (MQ=255) aCAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATGCAgg > 2:321892/1‑150 (MQ=255) cATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATGCAGGATgg < 1:436378/150‑1 (MQ=255) tGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATGCAGGATGGGa > 1:324753/1‑150 (MQ=255) aaTCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATGCAGGATGGGACGt < 1:313605/150‑1 (MQ=255) aaTCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATGCAGGATGGGACGt > 2:313605/1‑150 (MQ=255) aTCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATGCAGGATGGGACGTc > 2:376959/1‑150 (MQ=255) aTCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACAAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGAGGGATGCAGGGTGTGCCGTc > 2:607101/1‑150 (MQ=255) tCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAAcgcg < 2:721798/128‑1 (MQ=255) tCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAAcgcg < 2:780192/128‑1 (MQ=255) tCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAAcgcg > 1:721798/1‑128 (MQ=255) tCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAAcgcg > 1:780192/1‑128 (MQ=255) gATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATGCAGGATGGGACGTCGAg < 1:276268/150‑1 (MQ=255) cAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATGCAGGa < 2:592121/129‑1 (MQ=255) cAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATGCAGGa > 1:592121/1‑129 (MQ=255) gCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATGCAGGATGGGACGTCGAGGCTCGTCGAAg < 2:580752/150‑1 (MQ=255) cccAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGCAACGCGGGGGATGCAGGATGGGACGTCGAGGCTCGTCGAAGGTTCGa < 1:835493/150‑1 (MQ=255) taagcagcCAGAACAGCGACCAGCCGATATTGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTGTTTCACGCGAGGGTACCTTGAAAAGCGGGGGATGCAGGATGGGACATCGAGGCTCGGTGAAGATACGAGGGGTg > 2:268245/3‑150 (MQ=255) cagcagCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCCCCAGGCGGTACTTGCGGGCGTCGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATGCAGGATGGGACGTCGAGGCTCGTCGAAGGTTCGAGGGTTGa < 1:786839/150‑1 (MQ=255) agcCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATGCAGGATGGGACGTCGAGGCTCGTCGAAGGTTCGAGGGTTGAATGt > 2:104465/1‑150 (MQ=255) gcCAGAACAGGGACCAGCCGATATAGCCGGTGGTCTTCACCACCAGGCGGTCCGTGCGGGGGATGACGTCTTTAACGCGAGGGTACCTCGAAATGCGGGGGGTGCAGGATGGGACGTCGAGGCCCGTCGAAGGGTCGAGGGTTGACTGtt > 2:89998/1‑150 (MQ=255) ccAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATGCAGGa < 2:608180/108‑1 (MQ=255) ccAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATGCAGGa > 1:608180/1‑108 (MQ=255) gAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATGCAGGATGGGACGTCGAGGCTCGTCGAAGGTTCGAGGGTTGAATGTTGCaa < 1:104465/150‑1 (MQ=255) aaCAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATGCAGGATGGGACGTCGAGGCTCGTCGAAGGTTCGAGGGTTGAATGTTGCAAg < 1:321892/150‑1 (MQ=255) gTCTCCGCCACCCGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTTGAAACGCGGGGGATGCAGGATGGGACGTCGAGGCTCGTCGAAGGTNCGAGGGTTGAATGTTGCAAGGGATTTTGTAGCTGAAACGAAGGGaaa < 1:268245/150‑1 (MQ=255) ttCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATGCAGGATGGGACGTCGAGGCTCGTCGAAGGTTCGAGGGtt < 1:141280/108‑1 (MQ=255) ttCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATGCAGGATGGGACGTCGAGGCTCGTCGAAGGTTCGAGGGtt > 2:141280/1‑108 (MQ=255) caccaGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATGCAGGATGGGACGTCGAGGCTCGTCGAAGGTTCGAGGGTTGAATGTTGCAAGGGATTTTGTAGCTGAAACGAAGGGAAAAATCCGAc > 2:716624/1‑150 (MQ=255) ccaGGCGGTACTTGCTGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATGCAGGATGGGACGTCGAGGCTCGTCGAAGGTTCGAGGGTTGAATGTTGCAAg > 2:467852/1‑113 (MQ=255) ccaGGCGGTACTTGCTGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATGCAGGATGGGACGTCGAGGCTCGTCGAAGGTTCGAGGGTTGAATGTTGCAAg > 2:467841/1‑113 (MQ=255) ccaGGCGGTACTTGCTGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATGCAGGATGGGACGTCGAGGCTCGTCGAAGGTTCGAGGGTTGAATGTTGCAAg < 1:467852/113‑1 (MQ=255) ccaGGCGGTACTTGCTGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATGCAGGATGGGACGTCGAGGCTCGTCGAAGGTTCGAGGGTTGAATGTTGCAAg < 1:467841/113‑1 (MQ=255) gTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATGCAGGATGGGACGTCGAGGCTCGTCGAAGGTTCGAGGGTTGAATGtt > 2:530203/1‑101 (MQ=255) gTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATGCAGGATGGGACGTCGAGGCTCGTCGAAGGTTCGAGGGTTGAATGtt < 1:530203/101‑1 (MQ=255) tACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATGCAGGATGGGACGTCGAGGCTCGTCGTACGTTCGAGGGtt > 1:270426/1‑93 (MQ=255) tACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATGCAGGATGGGACGTCGAGGCTCGTCGAAGGTTCGAGGGtt < 2:270426/93‑1 (MQ=255) ttGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATGCAGGATGGGACGTCGAGGCTCGTCGAAGGTTCGAGGGTTGAATGTTGCAAGGGATTTTGTAGCTGAAACGAAGGGAAAAATCCGACCGGCCTATCAGgt < 1:89998/150‑1 (MQ=255) | CCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATGCAGGATGGGACGTCGAGGCTCGTCGAAGGTTCGAGGGTTGAATGTTGCAAGGGATTTTGTAGCTGAAACGAAGGGAAAAATCCGACCGGCCTATCAGGT > NC_002947/966553‑966846 |
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: ‑ |