Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 1,491,389 | A→G | K163R (AAG→AGG) | cybB → | cytochrome b561 |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,491,389 | 0 | A | G | 93.5% | 88.7 / ‑2.2 | 31 | K163R (AAG→AGG) | cybB | cytochrome b561 |
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (12/17); total (13/18) | |||||||||||
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.74e-01 |
CGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCTATT > NC_000913/1491254‑1491514 | cGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGa < 1:642752/139‑1 (MQ=255) aCGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTaaa < 2:1023721/139‑1 (MQ=255) cTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGTACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGa < 2:323871/96‑1 (MQ=255) cTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGTACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGa > 1:323871/1‑96 (MQ=255) cTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGt < 2:7820/130‑1 (MQ=255) cTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGt > 1:7820/1‑130 (MQ=255) cTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGAGATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAgg > 2:849320/1‑139 (MQ=255) tAAAGTCGTTGCATGAAACGCTGGCGAATCTGGGATCTTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGCCGACACACTTATACGCATGATGCCGCGTAAACGTCCCTGAAGGat > 2:720898/1‑139 (MQ=255) tAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACACTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGat > 2:857427/1‑139 (MQ=255) aGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATAtt > 1:550277/1‑139 (MQ=255) tGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCATACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAg > 2:40741/1‑139 (MQ=255) tGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAg > 2:58257/1‑139 (MQ=255) tGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAg > 2:838449/1‑139 (MQ=255) tGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAg > 2:100750/1‑139 (MQ=255) aTGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGaaaa < 1:849320/139‑1 (MQ=255) cTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCg < 1:100750/139‑1 (MQ=255) tatTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGt < 2:101414/139‑1 (MQ=255) ttGTCATCGGGTTGCACGCTGCGGCTGCACTGGCATACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTa < 1:40741/139‑1 (MQ=255) ttGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTa < 2:1050956/139‑1 (MQ=255) tCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAacaca > 2:1080161/1‑59 (MQ=255) tCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAacaca < 1:1080161/59‑1 (MQ=255) tCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAg < 1:58257/139‑1 (MQ=255) tGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATAc < 2:574656/134‑1 (MQ=255) tGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATAc > 1:574656/1‑134 (MQ=255) cGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTa > 2:208231/1‑118 (MQ=255) cGGCTGCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTa < 1:208231/118‑1 (MQ=255) gCACTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCg < 1:857427/139‑1 (MQ=255) aCTGGCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGa < 1:113760/139‑1 (MQ=255) ggCACACCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAgc < 2:776273/139‑1 (MQ=255) acacCACTATTTCTGGAGGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGcc < 2:550277/139‑1 (MQ=255) cacGATTTCTGGAGGGACAACACACCTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGACTACGGGTCTAAATCGGAAGCGCCtatt < 1:836728/139‑1 (MQ=255) | CGATGCCTTACGCTTCAGAGGCCAATTTCGAACGGGTAGATAGCTTAAAGTCGTGGCATGAAACGCTGGCGAATCTGGGATATTTTGTCATCGGGTTGCACGCTGCGGCTGCACTGGCACACCACTATTTCTGGAAGGACAACACACTTCTACGCATGATGCCGCGTAAACGTTCCTGAAGGATATTTAAAGAAAACGCCTGTACTAAAACCGACCCGTGGTACAGGCGAAGAATACGGGTCTACATCGGAAGCGCCTATT > NC_000913/1491254‑1491514 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |