Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 1,903,222 | A→G | 23.6% | Q1244R (CAG→CGG) | PP_1703 → | assimilatory nitrate reductase/sulfite reductase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 1,903,222 | 0 | A | G | 23.6% | 19.8 / 5.0 | 17 | Q1244R (CAG→CGG) | PP_1703 | assimilatory nitrate reductase/sulfite reductase |
Reads supporting (aligned to +/- strand): ref base A (6/7); new base G (2/2); total (8/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.70e-01 |
TGGGCCCTGGTACTGGCATTGCACCGTTCCGTGCGTTTCTCGAGGAGCGAGAGGCACGGGGGGCAAAAGGCGACAACTGGCTGTTCTTCGGCGAGCAGTACGCAGCGACCGATTTCTATTACCAGGAGCAGTTGCAGGCCTGGCAAGCGGCGGGGCACCTGCGGCTGGACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGT > NC_002947/1903093‑1903371 | tgagCCCTGGTACAGGCATTGCACCGTTCCGTGCGTTTCTCGAGGAGCCAGAAGCACGGGGGGCAAAAGCCGACAAGTGGCTGTTCTTCGGCGAGCAGTACGCAGCGACCGATTTCTATTATCAGGAGCAGTTGCAGGCCTGGCAAccgg > 2:17504/4‑150 (MQ=255) tGGTACTGGCATTGCACCGTTCCGTGCGTTTCTCGAGGAGCGAGAGGCACGGGGGGCAAAAGGCGACAACTGGCTGTTCTTCGGCGAGCAGTACGCAGCGACCGATTTCTATTACCAGGAGCAGTTGCAGGCCTGGCAAGCGGCGGGGCa > 1:89100/1‑150 (MQ=255) gCATTGCACCGTTCCGTGCGTTTCTCGAGGAGCGAGAGGCACGGGGGGCAAAAGGCGACAACTGGCTGTTCTTCGGCGAGCAGTACGCAGCGACCGATTTCTATTACCAGGAGCGGTTGCAGGCCTGGCAAGCGGCGGGGCACCTGCGGc > 1:127393/1‑150 (MQ=255) gCATTGCACCGTTCCGTGCGTTTCTCGAGGAGCGAGAGGCACGGGGGGCAAAAGGCGACAACTGGCTGTTCTTCGGCGAGCAGTACGCAGCGACCGATTTCTATTACCAGGAGCGGTTGCAGGCCTGGCAAGCGGCGGGGCACCTGCGGc > 1:127364/1‑150 (MQ=255) aCCGTTCCGTGCGTTTCTCGAGGAGCGAGAGGCACGGGGGGCAAAAGGCGACAACTGGCTGTTCTTCGGCGAGCAGTACGCAGCGACCGATTTCTATTACCAGGAGCAGTTGCAGGCCTGGCAAGCGGCGGGGCACCTGCGGCTGGacac < 1:360064/150‑1 (MQ=255) tttCTCGAGGAGCGAGAGGCACGGGGGGCAAAAGGCGACAACTGGCTGTTCTTCGGCGAGCAGTACGCAGCGACCGATTTCTATTACCAGGAGCGGTTGCAGGCCTGGCAAGCGGCGGGGCACCTGCGGCTGGACACGGCGTTCTCGCgg < 2:127393/150‑1 (MQ=255) tttCTCGAGGAGCGAGAGGCACGGGGGGCAAAAGGCGACAACTGGCTGTTCTTCGGCGAGCAGTACGCAGCGACCGATTTCTATTACCAGGAGCGGTTGCAGGCCTGGCAAGCGGCGGGGCACCTGCGGCTGGACACGGCGTTCTCGCgg < 2:127364/150‑1 (MQ=255) acGAGGAGCGAGAGGCACGGGGGGCAAAGCGGGACAACTGGCTGGTCAACGGCGAGCAGTACGCAGGGACCGATGTCTATTACGAGGAGCAGTTGCAGGCCTGGGAAGCGGCGGGACAACTGAGACTGGAAACGGCGTTCTCGAGGGAca > 2:184409/2‑149 (MQ=255) aCGGGGGGCAAAAGGCGACAACTGGCTGTTCTTCGGCGAGCAGTACGCAGCGACCGATTTCTATTACCAGGAGCAGTTGCAGGCCTGGCAAGCGGCGGGGCACCTGCGGCTGGACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTa < 2:89100/150‑1 (MQ=255) aCGGGGGGCAAAAGGCGACAACTGGCTGTTCTTCGGCGAGCAGTACGCAGCGAACGATTTCTATTACCAGGAGCAGTTg < 2:169756/79‑1 (MQ=255) aCGGGGGGCAAAAGGCGACAACTGGCTGTTCTTCGGCGAGCAGTACGCAGCGAACGATTTCTATTACCAGGAGCAGTTg > 1:169756/1‑79 (MQ=255) gACAACTGGCTGTTCTTCGGCGAGCAGTACGCAGCGACCGATTTCTATTACCAGGAGCAGTTGCAGGCCTGGCAAGCGGCGGGGCACCTGCGGCTGGACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATg < 2:373233/150‑1 (MQ=255) gCTGTTCTTCGGCGAGCAGTACGCAGCGACCGATTTCTATTACCAGGAGCAGTTGCAGGCCTGGCAAGCGGCGGGGCACCTGCGGCTGGACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCa < 2:270507/150‑1 (MQ=255) gCAGCGACCGATTTCTATTACCAGGAGCAGTTGCAGGCCTGGCAAGCGGCGGGGCACCTGCGGCTGGACACGGCGTTCTCGCGGGACCAGGTCGagaaga < 1:355971/100‑1 (MQ=255) gCAGCGACCGATTTCTATTACCAGGAGCAGTTGCAGGCCTGGCAAGCGGCGGGGCACCTGCGGCTGGACACGGCGTTCTCGCGGGACCAGGTCGagaaga > 2:355971/1‑100 (MQ=255) ctaGGAGCAGTTGCAGGCCTGGCAAGCGGCGGGGCACCTGCGGCTGGACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTAttt > 2:368608/3‑150 (MQ=255) gAGCAGTTGCAGGCCTGGCAAGCGGCGGGGCACCTGCGGCTGGACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTAt < 1:368608/150‑1 (MQ=255) aGTTGCAGGCCTGGCAAGCGGCGGGGCACCTGCGGCTGGACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTAtgtgt < 1:275338/150‑1 (MQ=255) | TGGGCCCTGGTACTGGCATTGCACCGTTCCGTGCGTTTCTCGAGGAGCGAGAGGCACGGGGGGCAAAAGGCGACAACTGGCTGTTCTTCGGCGAGCAGTACGCAGCGACCGATTTCTATTACCAGGAGCAGTTGCAGGCCTGGCAAGCGGCGGGGCACCTGCGGCTGGACACGGCGTTCTCGCGGGACCAGGTCGAGAAGATCTATGTGCAGCAGCGCATGCTCGAGCAGGGTGCGCTGTTGTGGCAATGGCTGGAGGCGGGGGCGTATTTCTATGTGT > NC_002947/1903093‑1903371 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |