Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 254,223 | A→T | 10.0% | intergenic (+21/+36) | prfH → / ← pepD | pseudogene, RF‑1 domain family; probable peptide chain release factor/aminoacyl‑histidine dipeptidase (peptidase D) |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 254,223 | 0 | A | T | 10.0% | 126.2 / 6.5 | 40 | intergenic (+21/+36) | prfH/pepD | pseudogene, RF‑1 domain family; probable peptide chain release factor/aminoacyl‑histidine dipeptidase (peptidase D) |
Reads supporting (aligned to +/- strand): ref base A (15/21); new base T (2/2); total (17/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.38e-01 |
GAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGG > NC_000913/254076‑254356 | gAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGaa < 1:865319/148‑1 (MQ=255) aaGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAg < 2:220354/148‑1 (MQ=255) cTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCg < 4:182687/148‑1 (MQ=255) cTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCg < 3:849162/148‑1 (MQ=255) gcaacagcaacaGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTaa < 4:642316/148‑1 (MQ=255) caacagcaacagGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTaaa < 3:342790/148‑1 (MQ=255) aacagcaacagGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAAc < 1:923149/148‑1 (MQ=255) cagcaacagGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGa > 3:461519/1‑140 (MQ=255) cagcaacagGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGa < 4:461519/140‑1 (MQ=255) cagcaacagGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGa < 4:577874/140‑1 (MQ=255) cagcaacagGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGa > 3:577874/1‑140 (MQ=255) cagcaacagGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTaa > 4:108475/1‑144 (MQ=255) cagcaacagGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTaa < 3:108475/144‑1 (MQ=255) acagGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATc < 3:331422/147‑1 (MQ=255) acagGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATc > 4:331422/1‑147 (MQ=255) gAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgc < 4:780692/147‑1 (MQ=255) gAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgc > 3:780692/1‑147 (MQ=255) aaaTCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGa < 2:89604/113‑1 (MQ=255) aaaTCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGa > 1:89604/1‑113 (MQ=255) cGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaaatcaa > 1:388128/1‑135 (MQ=255) cGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaaatcaa < 2:388128/135‑1 (MQ=255) aTCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTaat < 3:136624/96‑2 (MQ=255) aTCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTaat > 4:136624/1‑95 (MQ=255) ccAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaaatcaaa < 1:1022495/114‑1 (MQ=255) ccAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaaatcaaa > 2:1022495/1‑114 (MQ=255) ttGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTaaa < 3:1034647/87‑1 (MQ=255) ttGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTaaa > 4:1034647/1‑87 (MQ=255) aCGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaa < 2:959639/99‑1 (MQ=255) aCGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGcaa > 1:959639/1‑99 (MQ=255) aaCCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGccttc > 1:897143/1‑80 (MQ=37) aaCCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGATGGCGTTTAATCGCCTTc < 2:897143/83‑1 (MQ=38) cAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTa > 3:183110/1‑85 (MQ=255) cAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTa < 4:183110/85‑1 (MQ=255) tGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCg < 2:580112/113‑1 (MQ=255) tGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCg > 1:580112/1‑113 (MQ=255) tGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGa > 4:969253/1‑114 (MQ=255) tGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGa < 3:969253/114‑1 (MQ=255) gCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGaa < 4:470931/113‑1 (MQ=255) gCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGaa > 3:470931/1‑113 (MQ=255) cGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATg > 4:542539/1‑91 (MQ=255) cGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATg < 3:542539/91‑1 (MQ=255) aaGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTa > 3:931845/1‑40 (MQ=255) aaGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTa < 4:931845/40‑1 (MQ=255) aaGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCgg > 4:628701/1‑135 (MQ=255) aaGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCgg < 3:628701/135‑1 (MQ=255) | GAAGCTGGAGCAACAGCAACAGGAAAATAGCGCGGCGCTGAAATCGCAGCGGCGAATGTTCCATCACCAGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTACTTCGCCGGAATTTCTTTCAGCAGTTCAGTCAGCAGTGTCCAGTAATGACCTACGCTTTCGATGTGAACTTGCTCATCCGGAGAGTGTGGACCGG > NC_000913/254076‑254356 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |