Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 2,173,821 | A→C | 11.5% | A83A (GCT→GCG) ‡ | PP_1926 ← | phosphotyrosine‑protein phosphatase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 2,173,821 | 0 | A | C | 11.5% | 79.9 / 5.1 | 35 | A83A (GCT→GCG) ‡ | PP_1926 | phosphotyrosine‑protein phosphatase |
Reads supporting (aligned to +/- strand): ref base A (16/15); new base C (2/2); total (18/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 = 8.89e-01 |
TCAGCAGACGGGCAGCGAACAGCCCGGTGCGGCCGGAGCCGCCCTTGCAGTGGATTGCCAGGGCCTTGTTGGAAGCGAGCAAAGCGTCCAACTGGTGGCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAGCGCTACGGCACCGGCCTGCTTGAGCGTATCC > NC_002947/2173685‑2173968 | tCAGCAGACGGGCAGCGAACAGCCCGGTGCGGCCGGAGCCGCCCTTGCAGTGGATTGCCAGGGCCTTGTTGGAAGCGAGCAAAGCGTCCAACTGGTGGCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGCTCCTGatcaa < 1:35666/147‑2 (MQ=255) tCAGCAGACGGGCAGCGAACAGCCCGGTGCGGCCGGAGCCGCCCTTGCAGTGGATTGCCAGGGCCTTGTTGGAAGCGAGCAAAGCGTCCAACTGGTGGCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGCTCCTGatcaa > 2:35666/1‑146 (MQ=255) tCAGCAGACGGGCAGCGAACAGCCCGGTGCGGCCGGAGCCGCCCTTGCAGTGGATTGCCAGGGCCTTGTTGGAAGCGAGCAAAGCGTCCAACTGGTGGCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCt > 1:167394/1‑150 (MQ=255) tCAGCAGACGGGCACCGAACAGCGCGGTGCGGCCGGAGCCGCCCTTGCAGTGGATTGCCAGGGCCTTGTTGGAAGCGAGCAAAGCGTCCAACTGGTGGCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCt > 2:70739/1‑150 (MQ=255) agACGGGCAGCGAACAGCCCGGTGCGGCCGGAGCCGCCCTTGCAGTGGATTGCCAGGGCCTTGTTGGAAGCGAGCAAAGCGTCCAACTGGTGGCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATc < 2:557600/150‑1 (MQ=255) aCAGCCCGGTGCGGCCGGAGCCGCCCTTGCAGTGGATTGCCAGGGCCTTGTGGGAAGCGAGCAAAGCGTCCAACTGGTGGCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAccc > 2:98445/1‑150 (MQ=255) gCGGCCGGAGCAGCCCTTGCAGTGGATTGCCAGGGCCTTGTTGGAAGCGAGCAAAGCGTCCAACTGGTGGCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAAccc < 1:164367/150‑1 (MQ=255) gccCTTGCAGTGGATTGCCAGGGCCTTGTTGGAAGCGAGCAAAGCGTCCAACTGGTGGCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGCTCCTGatcaaaa < 1:423949/109‑4 (MQ=255) gccCTTGCAGTGGATTGCCAGGGCCTTGTTGGAAGCGAGCAAAGCGTCCAACTGGTGGCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGCTCCTGatcaaaa > 2:423949/1‑106 (MQ=255) gccCTTGCAGTGGATTGCCAGGGCCTTGTTGGAAGCGAGCAAAGCGTCCAACTGGTGGCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGatcatc > 1:430516/1‑108 (MQ=255) gccCTTGCAGTGGATTGCCAGGGCCTTGTTGGAAGCGAGCAAAGCGTCCAACTGGTGGCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGatcatc < 2:430516/108‑1 (MQ=255) gccCTTGCAGTGGATTGCCAGGGCCTTGTTGGAAGCGAGCAAAGCGTCCAACTGGTGGCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCt < 2:366309/110‑1 (MQ=255) gccCTTGCAGTGGATTGCCAGGGCCTTGTTGGAAGCGAGCAAAGCGTCCAACTGGTGGCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCt > 1:366309/1‑110 (MQ=255) gCAGTGGATTGCCAGGGCCTTGTTGGAAGCGAGCAAAGCGTCCAACTGGTGGCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCCCTCCAACCCATGCTGCTTGGCCGCGAt < 2:167394/150‑1 (MQ=255) ggAAGCGAGCAAAGCGTCCAACTGGTGGCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTc > 2:473882/1‑150 (MQ=255) aaGCGAGCAAAGCGTCCAACTGGTGGCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCt < 1:473882/150‑1 (MQ=255) aaCTGGTGGCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATcc < 2:590086/60‑1 (MQ=255) aaCTGGTGGCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATcc > 1:590086/1‑60 (MQ=255) ggtggCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTcc > 1:592423/1‑143 (MQ=255) ggtggCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTcc < 2:592423/143‑1 (MQ=255) tggCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTc < 1:747349/125‑1 (MQ=255) tggCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTc > 2:747349/1‑125 (MQ=255) tggCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGACCGCGATGCCAAGCTGGGCAACGCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGc < 1:98445/150‑1 (MQ=255) gAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCaa > 1:627353/1‑76 (MQ=255) gAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCaa < 2:627353/76‑1 (MQ=255) tctctcCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACaa > 1:688382/1‑150 (MQ=255) ctctcCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAg > 1:757244/1‑150 (MQ=255) cccAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAAcc < 2:193730/70‑1 (MQ=255) cccAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAAcc > 1:193730/1‑70 (MQ=255) cTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAGCGCTACGGCACCGGc > 1:62049/1‑150 (MQ=255) cGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAGCGCTACGGCACCGGCCt < 2:757244/150‑1 (MQ=255) gTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAGCGCTACGGCACCGGCCTGCtt < 2:62049/150‑1 (MQ=255) ttGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAGCGCTACGGCACCGGCCTGCTTGAg < 1:55765/150‑1 (MQ=255) tCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAGCGCTACGGCACCGGCCTGCTTGAGCGTa > 1:193496/1‑150 (MQ=255) ggAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAGCGCTACGGCACCGGCCTGCTTGAGCGTATcc > 2:604991/1‑150 (MQ=255) | TCAGCAGACGGGCAGCGAACAGCCCGGTGCGGCCGGAGCCGCCCTTGCAGTGGATTGCCAGGGCCTTGTTGGAAGCGAGCAAAGCGTCCAACTGGTGGCGAATCTCTCCCAGCCGAGCCTCGAAGTTTTGATCAGGAGCCTGATCATCCTCGATCGGCAGCTGGAACCACTCCAACCCATGCTGCTTGGCCGCGATGCCAAGCTGGGCAACTCCGTTTTCCTGAAGCTCTGCGTCCGACATCAGCGTAACAAGCGCTACGGCACCGGCCTGCTTGAGCGTATCC > NC_002947/2173685‑2173968 |
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: ‑ |