Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 668,267 | A→C | 12.1% | F87C (TTC→TGC) ‡ | PP_0571 ← | membrane protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 668,267 | 0 | A | C | 12.1% | 77.9 / 5.0 | 33 | F87C (TTC→TGC) ‡ | PP_0571 | membrane protein |
Reads supporting (aligned to +/- strand): ref base A (14/15); new base C (2/2); total (16/17) | |||||||||||
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.28e-01 |
TGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGCCAGCAGGTTGAACAACCCGCTCAATGC > NC_002947/668125‑668415 | tGAAAGCCTGGTGCTGCTGTGCAGCCCTGTCCCGACGCAAGGCCGCTCCCACAAAGACCTCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGcgcg < 1:389004/150‑1 (MQ=255) aaGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGAt < 1:335751/150‑1 (MQ=255) ggggCTGCTGTGCAACCCTGTCCCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGa < 1:159764/150‑1 (MQ=255) gggCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAg < 1:167309/150‑1 (MQ=255) tgctgTGCAGCCCTGTCGCCACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTAca < 1:747222/150‑1 (MQ=255) tgCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACagcagc < 1:365339/150‑1 (MQ=255) tgCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACagcagc < 1:46744/150‑1 (MQ=255) ccTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAgg < 1:194021/150‑1 (MQ=255) gcgACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGa < 2:715883/128‑1 (MQ=255) gcgACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGa > 1:715883/1‑128 (MQ=255) aaGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAggg > 1:610584/1‑150 (MQ=255) tCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATTCATGcgcgcct < 1:221835/107‑4 (MQ=255) tCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATTCATGcgcgcct > 2:221835/1‑104 (MQ=255) cacaAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGCCCGATGATCAGGGCCAGGATGCTc > 2:198838/1‑150 (MQ=255) gATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCaa > 2:599377/1‑150 (MQ=255) aTCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGc < 2:739026/145‑1 (MQ=255) aTCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGc > 1:739026/1‑145 (MQ=255) aCCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTg > 1:426242/1‑150 (MQ=255) ccAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTgcg > 2:168911/1‑150 (MQ=255) tcatcaCCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACa < 2:426242/150‑1 (MQ=255) tcaCCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGc > 1:443062/1‑150 (MQ=255) tcaCCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGc > 2:107459/1‑150 (MQ=255) aCCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCa > 1:34029/1‑150 (MQ=255) ccACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATc > 2:689814/1‑150 (MQ=255) ccACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATc < 2:631001/150‑1 (MQ=255) caGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGAt < 2:610584/150‑1 (MQ=255) ccccAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGaa < 1:107459/150‑1 (MQ=255) cGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGCcagca > 1:79669/1‑150 (MQ=255) tccgcgcATTCATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGCcagca < 2:302545/135‑1 (MQ=255) tccgcgcATTCATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGCcagca > 1:302545/3‑137 (MQ=255) gcgcATGAATGCGCGGATTGCCGAGTACATCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGCCAGCAGGTTGAACAACCCGCt > 1:800652/1‑150 (MQ=255) gcATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGCCAGCAGGTTGAACAACCCGCTCa > 2:374357/1‑150 (MQ=255) gAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGCCAGCAGGTTGAACAACCCGCTCAATGc < 1:374357/150‑1 (MQ=255) | TGAAAGCCTGGGGCTGCTGTGCAGCCCTGTCGCGACGCAAGGCCGCTCCCACAAAGACCGCGCAGATCAGTGGATCACCCAGCTCATCACCCACAGCCCCAGCACCAGCCAGATGATTCCGAGGATGATCGAGGCGCGCATGAATGCGCGGATTGCCGAGTACAGCAGCATCAGGCCGATGATCAGGGCCAGGATGCTCAGCAGCGAAGTGTCCATGCCCAAGGTGCGCGACAGGCCATCGATGAAGTTGCCACCGGCATTGGCCAGCAGGTTGAACAACCCGCTCAATGC > NC_002947/668125‑668415 |
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: ‑ |