Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 1,607,039 | A→G | 16.7% | K14R (AAG→AGG) | lsrG → | (4S)‑4‑hydroxy‑5‑phosphonooxypentane‑2,3‑dione isomerase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,607,039 | 0 | A | G | 16.7% | 57.1 / 4.6 | 29 | K14R (AAG→AGG) | lsrG | (4S)‑4‑hydroxy‑5‑phosphonooxypentane‑2,3‑dione isomerase |
Reads supporting (aligned to +/- strand): ref base A (14/10); new base G (2/3); total (16/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.32e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.75e-01 |
GAAAGCCGTACAGGCGGTGGTTCACCATAACGAAACGGCTGATCGGGCATATGAACTCTATCTGAGTGAAAAACAGTAACTGCGGATCTAAGGAGAAGAATTATGCACGTCACACTGGTTGAAATTAACGTTCATGAAGACAAGGTTGACGAGTTTATCGAAGTTTTTCGCCAGAACCACCTGGGCTCTGTACAGGAAGAAGGCAATTTGCGCTTCGATGTCTTACAGGACCCGGAAGTGAATTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAA > NC_000913/1606897‑1607175 | gaAAGCCGTACAGGCGGTGGTTCACCATAACGAAACGGCTGATCGGGCATATGAACTCTATCTGAGTGAAAAACAGTAACTGCGGATCTAAGGAGAAGAATTATGCACGTCACACTGGTTGAAATTAACGTTCATGAAGACaa > 1:493912/1‑143 (MQ=255) tACAGGCGGTGGTTCACCATAACGAAACGGCTGATCAGGCATATGAACTCTATCTGAGTGAAAAACAGTAACTGCGGATCTAAGGAGAAGAATTATGCACGTCACACTGGTTGAAATTAACGTTCATGAAGACCAGGTTGACg > 1:167237/1‑143 (MQ=255) gtTCACCATAACGAAACGGCTGATCGGGCATATGAACTCTATCTGAGTGAAAAACAGTAACTGCGGATCTAAGGAGAAGAATTATGCACGTCACACTGGTTGAAATTAACGTTCATGAAGACAAGGTTGACGAGTTTATCGaa > 1:138766/1‑143 (MQ=255) tCACCATAACGAAACGGCTGATCGGGCATATGAACTCTATCTGAGTGAAAAACAGTAACTGCGGGTCTAAGGCGAAGAATTATGCACGTCACACTGGTTGAAATTAACGTTCATGAAGACAAGGTTGACGAGTTTATCGAAGt > 2:406180/1‑143 (MQ=255) tCACCATAACGAAACGGCTGATCCGGCATATGAACTCTATCTGAGTGAAAAACAGTAACTGCGGATCTAAGGAGAAGAATTATGCACGTCACACTGGTTGAAATTAACGTTCATGAAGACAAGGATGACGAGTGTATAGAAGt > 2:188204/1‑143 (MQ=255) aCGAAACTGCTGATCGGGCATATGAACTCTATCTGAGTGAAAAACAGTAACTGCGGATCTAAGGAGAAGAATTATGCACGTCACACTGGTTGAAATTAACGTTCATGAAGACAAGGTTGACGAGTTTATCGAAGTTTTTCGcc < 2:348385/143‑1 (MQ=255) aaaCGGTTGATCGGGCATATGAACTCTAACTGAGTGAAAAACAGTAACCGCGGATCTAAGGAGAAGAATTATGCACGTCACACTGGTTGAAATTAACGTTCATGAAGACAAGGTTGACGAGTTTATCGAAGTTTTTCGCCAGa < 1:331585/143‑1 (MQ=255) tCGGGCATATGAACTCTATCTGAGTGAAAAACAGTAACTGCGGATCTAAGGAGAAGAATTATGCACGTCACACTGGTTGAAATTAACGTTCATGAAGACAAGGTTGACGAGTTTATCGAAGTTTTTCGCCAGAACCACCTggg < 1:331190/143‑1 (MQ=255) tCGGGCATATGAACTCTATCTGAGTGAAAAACAGTAACTGCGGATCTAAGGAGAAGAATTATGCACGTCACACTGGTTGAAATTAACGTTCATGAAGACAAGGTTGACGAGTTTATCGAAGTTTTTCGCCAGAACCACCTggg < 1:327773/143‑1 (MQ=255) cGGGCATATGAACTCTGGCTGTGTGAACAACAGTGACCGCGGATCTAAGGAGAAGAATTCTGCACGTCACTCGGGTTGGAATTCGCGTTCATGAGGACAAGGTTGACGAGTTTATCGAAGTTTTTCGCCAGAACCACCTGGGc < 2:60479/143‑1 (MQ=255) atatGAACTCTATCTGAGTGAAAAACAGTAACTGCGGATCTAAGGAGAAGAATTATGCACGTCACACTGGTTGAAATTAACGTTCATGAAGACAAGGTTGACGACTTTATCGAAGTTTTTCGCCAGAACCACCTGGGCTCTGt > 2:50309/1‑143 (MQ=255) gAACTCTATCTGAGTGAAAAACAGTAACTGCGGATCTAAGGAGAAGAATTATGCACGTCACACTGGTTGAAATTAACGTTCATGACGACAAGGTTGACGAGTTTATCGAAGTTTTTCGCCAGAACCACCTGGGCTCTGTACAg < 1:326624/143‑1 (MQ=255) tCTGAGTGAAAAACAGTAACTGCGGATCTAAGGAGAAGAATTATGCACGTCACACTGGTTGAAATTAACGTTCATGAAGACAAGGTTGACGAGTTTATCGAAGTTTTTCGCCAGAACCACCTGGGCTCTGTACAGGAAGAAgg < 1:445629/143‑1 (MQ=255) aaaaaCAGTAACTGCGGATCTAAGGAGAAGAATTATGCACGTCACACTGGTTGAAATTAACGTTCATGAAGACAAGGTTGACGAGTTTATCGAAGTTTTTCGCCAGAACCACCTGGGCTCTGTACAGGAAGAAGGCAATTTgc > 1:293270/1‑143 (MQ=255) aGTGACTGCGGATCTAAGGAGAAGAATTATGCACGTCACACTGGTTGAAATTAACGTTCATGAAGACAAGGTTGACGAGTTTATCGAAGTTTTTCGCCAGAACCACCTGGGCTCTGTACAGGAAGAAGGCAATTTGCGCTTCg < 1:122700/143‑1 (MQ=255) aaCTGCGGATCTAAGGAGAAGAATTATGCACGTCACACTGGTTGAAATTAACGTTCATGAAGACAAGGTTGCCGAGTTTATCGAAGTTTTTCGCCAGAACCACCTGGGCTCTGTACAGGAGGAAGGCAATGTGCGCTTCGATg > 1:18225/1‑143 (MQ=255) aTCTAAGGAGAAGAATTATGCACGTAACACTGGTTGAAATTAACGTTCATGAAGACAAGGTTGACGAGTTTATCGAAGTTTTTCGCCAGAACCACCTGGGCTCTGTACAGGAAGAAGGCAATTTGCGCTTCGATGTCTTACAg > 1:362843/1‑143 (MQ=255) aaTTATGCACGTCACACTGGTTGAAATTAACGTTCATGAAGACAAGGTTGACGAGTTTATCGAAGTTTTTCGCCAGAACCACCTGGGCTCTGTACAGGAAGAAGGCAATTTGCGCTTCGATGTCTTACAGGACCCGGAAGTGa > 1:128238/1‑143 (MQ=255) tGCACGTCACACTGGTTGAAATTAACGTTCATGAAGACAAGGTTGACGAGTTTATCGAAGTTTTTCGCCAGAACCACCTGGGCTCTGTACAGGACGAAGGCAATTTGCGCTTCGACGTCTTACCGGCCCCGGAAGTGAATTcc > 2:215865/1‑142 (MQ=255) cGTCACACTGGTTGAAATTAACGTTCATGAAGACAAGGTTGACGAGTTTATCGAAGTTTTTCGCCAGAACCACCTGGGCTCTGTACAGGAAGAAGGCAATTTGCGCTTCGATGTCTTACAGGACCCGGAAGTGAATTCGCGCt > 1:318719/1‑143 (MQ=255) cGTCACACTGGTTGAAATTAACGTTCATGAAGACAAGGTTGACGAGTTTATCGAAGTTTTTCGCCAGAACCACCTGGGCTCTGTACAGGAAGAAGGCAATTTGCGCTTCGATGTCTTACAGGACCCGGAAGTGAATTCGCGCt < 2:293270/143‑1 (MQ=255) cacaCTGGTTGAAATGAACGGTCATGAAGACAAGGTTGACGAGTTTATCGAAGTTTTTCGCCAGAACCACCTGGGCTCTGTACAGGAAGAAGGCAATTTGCGCTTCGATGTCTTACAGGACCCGGAAGTGAATTCGCGCtttt < 1:133719/143‑1 (MQ=255) acTGGTTGAAATTAACGTTCATGAAGACAGGGTTGACGAGTTTATCGAAGTTTTTCGCCAGAACCACCTGGGCTCTGTACAGGAAGAAGGCaa < 2:190307/93‑1 (MQ=255) acTGGTTGAAATTAACGTTCATGAAGACAGGGTTGACGAGTTTATCGAAGTTTTTCGCCAGAACCACCTGGGCTCTGTACAGGAAGAAGGCaa < 2:189541/93‑1 (MQ=255) acTGGTTGAAATTAACGTTCATGAAGACAGGGTTGACGAGTTTATCGAAGTTTTTCGCCAGAACCACCTGGGCTCTGTACAGGAAGAAGGCaa > 1:189541/1‑93 (MQ=255) acTGGTTGAAATTAACGTTCATGAAGACAGGGTTGACGAGTTTATCGAAGTTTTTCGCCAGAACCACCTGGGCTCTGTACAGGAAGAAGGCaa > 1:190307/1‑93 (MQ=255) tcgTTGCAATTAACTTTCATGAAGGCCGTTTTGACGAGTTTATCGAAGTCTTTCGCCAGCACCCCCTGGGCTCTGTACAGGAAGAAGGCAATTTGCGCTTCGATGTCTTACAGGACCCGGAAGTGAATTCGCGCTTTTATATc < 2:335301/141‑1 (MQ=255) ttGAAATTAACGTTCATGAAGACAAGGTTGACGAGTTTATCGAAGTTTTTCGCCAGAACCACCTGGGCTCTGTACAGGAAGAAGGCAATTTGCGCTTCGATGTCTTACAGGACCCGGAAGTGAATTCGCGCTTTTATATCTAc > 1:314601/1‑143 (MQ=255) aaGACAAGGTTGACGAGTTTATCGAAGTTTTTCGCCAGAACCACCTGGGCTCTGTACAGGAAGAAGGCAATTTGCGCTTCGATGTCTTACAGGACCCGGAAGTGAATTCGCGCTTTTATATCTACGAAGCCTATAAAGATGaa > 1:20203/1‑143 (MQ=255) aaGACAAGGTTGACGAGTTTATCGAAGTTTTTCGCCAGAACCACCTGGGCTCTGTACAGGAAGAAGGCAATTTGCGCTTCGATGTCTTACAGGACCCGGAAGTGAATTCGCGATTTTATATCTAAGAAGCCTATAAAGATGaa > 1:12612/1‑143 (MQ=255) | GAAAGCCGTACAGGCGGTGGTTCACCATAACGAAACGGCTGATCGGGCATATGAACTCTATCTGAGTGAAAAACAGTAACTGCGGATCTAAGGAGAAGAATTATGCACGTCACACTGGTTGAAATTAACGTTCATGAAGACAAGGTTGACGAGTTTATCGAAGTTTTTCGCCAGAACCACCTGGGCTCTGTACAGGAAGAAGGCAATTTGCGCTTCGATGTCTTACAGGACCCGGAAGTGAATTCGCGCTTTTATATCTACGAAGCCTATAAAGATGAA > NC_000913/1606897‑1607175 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |