Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 3,982,923 | A→G | 9.1% | intergenic (+112/‑35) | proM → / → aslB | tRNA‑Pro/putative AslA‑specific sulfatase‑maturating enzyme |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,982,923 | 0 | A | G | 9.1% | 130.5 / 4.2 | 44 | intergenic (+112/‑35) | proM/aslB | tRNA‑Pro/putative AslA‑specific sulfatase‑maturating enzyme |
Reads supporting (aligned to +/- strand): ref base A (23/17); new base G (2/2); total (25/19) | |||||||||||
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.00e-01 |
GGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGC > NC_000913/3982784‑3983071 | ggAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGaa > 3:295802/1‑149 (MQ=255) gCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGa < 1:48630/132‑1 (MQ=255) gCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGa > 2:48630/1‑132 (MQ=255) gACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGa < 3:430456/149‑1 (MQ=255) ggTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCa > 1:413543/1‑125 (MQ=255) ggTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCa < 2:413543/125‑1 (MQ=255) tttGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCATTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCtttt < 4:295802/149‑1 (MQ=255) ttttCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAAc > 3:418924/1‑82 (MQ=255) ttttCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAAc < 4:418924/82‑1 (MQ=255) cTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACATCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGAt < 2:408960/149‑1 (MQ=255) tgtgCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATg > 3:78265/1‑149 (MQ=255) aCCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCg > 4:198052/1‑149 (MQ=255) tccctTCGCAATAAACGCCCGTAATAACTCATTGCCCCATGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTc > 4:140004/1‑149 (MQ=255) tccctTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTc < 4:9165/149‑1 (MQ=255) tccctTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTc > 1:198754/1‑149 (MQ=255) aTAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAAt < 3:22755/149‑1 (MQ=255) cGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACtgt > 2:13371/1‑149 (MQ=255) taataaCTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGa > 1:344214/1‑149 (MQ=255) aCTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGactact > 2:425531/1‑149 (MQ=255) tCATTGCCCCACGGTATGATTTCGCCCTTAGCGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGt > 4:147374/1‑149 (MQ=255) tCATTGCCCCACGGTATGATTTCGCCCTTAGCGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGt > 2:401834/1‑149 (MQ=255) tCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGAtt < 1:131507/126‑1 (MQ=255) tCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGAtt > 2:131507/1‑126 (MQ=255) tCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGt > 3:3363/1‑149 (MQ=255) cATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGtt > 2:205202/1‑149 (MQ=255) aTTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCAtgt < 2:305095/98‑1 (MQ=255) aTTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCAtgt > 1:305095/1‑98 (MQ=255) ttGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGc > 3:267404/1‑88 (MQ=255) ttGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGc < 4:267404/88‑1 (MQ=255) ttGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGtttt < 2:344214/149‑1 (MQ=255) ccccACGGTATGATTTCGCCCTTAGCGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTAt < 1:401834/149‑1 (MQ=255) ccccACGGTATGATTTCGCCCTTAGCGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTAt < 3:147374/149‑1 (MQ=255) cccACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTAtc > 1:51721/1‑149 (MQ=255) cACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTAtctc < 4:3363/149‑1 (MQ=255) tATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCTACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGCTTGTAATCTGAACTGTGACTACTGTGTTTATCTCGaaaa > 3:390015/1‑149 (MQ=255) tATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGaaaa > 3:221935/1‑149 (MQ=255) tATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGaaaa > 1:344750/1‑149 (MQ=255) gATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACa > 2:256562/1‑149 (MQ=255) aTTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACaa > 4:8409/1‑149 (MQ=255) aTTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACaa < 4:78265/149‑1 (MQ=255) aTTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACaa > 2:472366/1‑149 (MQ=255) tttCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAAt < 3:198052/149‑1 (MQ=255) cTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTa < 1:425531/149‑1 (MQ=255) ttAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTAc < 1:13371/149‑1 (MQ=255) aCGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACcgc > 2:369210/1‑149 (MQ=255) | GGAGGTTCGAATCCTCTCTCGCCGACCAATTTTGAACCCCGCTTCGGCGGGGTTTTTTGTTTTCTGTGCATTTCGTCACCCTCCCTTCGCAATAAACGCCCGTAATAACTCATTGCCCCACGGTATGATTTCGCCCTTAACGTATTGAAGGATGACTTCAGGCAAGGAGCGACCATGCTGCAACAGGTTCCAACGCGTGCTTTTCATGTGATGGCGAAACCGAGCGGTTCCGATTGTAATCTGAACTGTGACTACTGTTTTTATCTCGAAAAACAATCCCTTTACCGC > NC_000913/3982784‑3983071 |
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: ‑ |