Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,013,258 | A→T | 31.1% | intergenic (+83/+48) | prfC → / ← potF‑I | peptide chain release factor 3/putrescine‑binding protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,013,258 | 0 | A | T | 31.1% | 36.6 / 16.9 | 24 | intergenic (+83/+48) | prfC/potF‑I | peptide chain release factor 3/putrescine‑binding protein |
Reads supporting (aligned to +/- strand): ref base A (6/10); new base T (3/4); total (9/15) | |||||||||||
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.54e-01 |
AACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCC‑ATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGCTTGCCCATCACCTCGTCG > NC_002947/1013110‑1013405 | aaCCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATAc < 2:350012/150‑1 (MQ=255) ggCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCGGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCtgtg < 2:273086/150‑1 (MQ=21) gcgcTGATGGAAGGGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTgg < 2:577752/150‑1 (MQ=21) atcgAAGTGCGGTGCCCGGACATCAACTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCgg < 2:278155/147‑1 (MQ=38) tCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCtgtg > 2:350486/1‑121 (MQ=17) tCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCtgtg < 1:350486/121‑1 (MQ=17) cAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAg < 1:742555/150‑1 (MQ=255) ttCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAggcgg > 1:586663/1‑150 (MQ=255) gcgcTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATa > 1:214998/1‑106 (MQ=255) gcgcTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATa < 2:214998/106‑1 (MQ=255) tGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCC‑ATTTCAATTTCCGGCTTTGATACTGGt > 2:62433/1‑150 (MQ=255) cGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCC‑ATTTCAATTTCCGGCTTTGATACTGGTCCAGAt < 2:504308/150‑1 (MQ=255) cTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCC‑ATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTc < 1:62433/150‑1 (MQ=255) ggCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCC‑ATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTc < 2:756493/150‑1 (MQ=255) gCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCC‑ATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCg > 1:447298/1‑150 (MQ=255) ccACAGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCC‑ATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTa > 2:227015/1‑119 (MQ=17) ccACAGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCC‑ATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTa < 1:227015/119‑1 (MQ=17) ccACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCCTAGCGGccctt < 2:628685/80‑3 (MQ=25) cacaGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCC‑ATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTc > 2:705455/1‑150 (MQ=255) ggAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCC‑ATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAgc < 2:13538/150‑1 (MQ=255) tCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCC‑ATTTCAATTTCCGGCTTTGATACTGGTCCAGa < 2:111175/102‑1 (MQ=255) tCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCC‑ATTTCAATTTCCGGCTTTGATACTGGTCCAGa > 1:111175/1‑102 (MQ=255) aGACCTGGAGGTCATGTAAAGGCTGGGGGAGCTGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCTC‑ATTTCCATTTCCCGCTTTGATACTTGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGc < 1:472095/150‑1 (MQ=255) ccaGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCTGCAAAGCAGCCCCGATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGAAAGATCTTCGCCGGCATCGCCTCCAGCGCATACAGCTTGcc > 1:225057/3‑150 (MQ=18) tCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCC‑ATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGCTTGCCCATCAc < 2:586663/150‑1 (MQ=255) aTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCC‑ATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGCTTGCCCATCACCtcgtcg > 2:206687/1‑150 (MQ=255) | AACCTGGCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCC‑ATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGCTTGCCCATCACCTCGTCG > NC_002947/1013110‑1013405 |
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: ‑ |