Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 86,245 | A→C | 20.9% | M158R (ATG→AGG) | PP_0075 ← | choline sulfate transporter |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 86,245 | 0 | A | C | 20.9% | 33.3 / 9.9 | 24 | M158R (ATG→AGG) | PP_0075 | choline sulfate transporter |
Reads supporting (aligned to +/- strand): ref base A (11/8); new base C (3/2); total (14/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.46e-01 |
TGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATGCCGCCGAGGAAGGTC > NC_002947/86109‑86394 | tGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCATGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGtt > 2:571791/1‑150 (MQ=255) tGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGtt < 2:286517/150‑1 (MQ=255) agagTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCt > 1:5336/1‑150 (MQ=255) tttGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGcc > 2:24312/1‑150 (MQ=255) gAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGAccggc < 2:264417/150‑4 (MQ=255) cACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGGGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGa < 1:438852/150‑1 (MQ=255) cGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGa > 2:169403/1‑150 (MQ=255) cACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATggc < 1:169403/150‑1 (MQ=255) ccAGCAGCAAGGAGGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGc > 2:753554/1‑150 (MQ=255) gcagcaAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGGCAGGTCCTGCCCGGAGATCAGCGCGGTGCCCGGTCTGCGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGa < 1:60750/150‑1 (MQ=255) caAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCATGGCCCTGCCca < 1:324693/89‑1 (MQ=255) caAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCATGGCCCTGCCca > 2:324693/1‑89 (MQ=255) gACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGAc > 1:284553/1‑150 (MQ=255) tGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGAc > 1:254928/1‑150 (MQ=255) cTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCaa > 2:233098/1‑150 (MQ=255) cGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGga > 2:20584/1‑115 (MQ=255) cGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGga < 1:20584/115‑1 (MQ=255) ggTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGGGCGCCGAAGCGCAATAGccc > 1:26715/1‑150 (MQ=255) aGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGagcagc < 2:5336/150‑1 (MQ=255) gggctGGCCTGGCTGCGCAGGGACATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGagcagcagc < 2:180695/146‑1 (MQ=255) tgctggGCATGGCCCTGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAgcgc < 2:113542/121‑1 (MQ=255) tgctggGCATGGCCCTGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAgcgc > 1:113542/2‑122 (MQ=255) gctggGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATgcc > 2:163231/1‑150 (MQ=255) ctggGCATGGCCCTGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATgccg > 2:646916/1‑150 (MQ=255) aTGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATGCCGCCGAGGAAGGTc > 1:289886/1‑150 (MQ=255) | TGGGCGCCAGAGTTTGAAGGCCGCACCGATCACGACGGTGGCCAGGCCCAGCAGCAAGGAAGGTAGGTCGACCTCGCCTGCGTGGCTGGCCAGGTCCTGCACGGTCATCAGCGCGGTGGCCTGGCTGGGCAGGGCCATGCCCAGCAGGTTGGGTAGCTGGCCCAAGGCAATGACGATGGCGGCACCCAGCGTGAAGCCGAGCACCACGGAATGGGAGACGAAGTTGACCAGCGCGCCGAAGCGCAATAGCCCGAGCAGCAGCTGGAAAATGCCGCCGAGGAAGGTC > NC_002947/86109‑86394 |
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: ‑ |