Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 3,156,591 | T→G | 23.4% | intergenic (+73/‑32) | yqhD → / → dkgA | NADPH‑dependent aldehyde reductase YqhD/methylglyoxal reductase DkgA |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,156,591 | 0 | T | G | 23.4% | 39.5 / 11.8 | 33 | intergenic (+73/‑32) | yqhD/dkgA | NADPH‑dependent aldehyde reductase YqhD/methylglyoxal reductase DkgA |
Reads supporting (aligned to +/- strand): ref base T (13/12); new base G (4/4); total (17/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.27e-01 |
TGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAG > NC_000913/3156457‑3156731 | tGGGCGAAATTCATTACATTACGTTGGATGTCAGCCGCCGTATATCGCCACGCTCCCGCTAAGTTTTTGTCTCCTCAAACGTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACAACTCACCGGAGCAGGCTCCGGt < 1:548728/143‑1 (MQ=255) gCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAg < 2:323928/143‑1 (MQ=255) aaaaTCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCCTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACAACTCACCGGAGCAGGCTCCGGTGAggtg < 1:519918/143‑4 (MQ=255) aaTCATGACATTACGTTGTATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCAGGCTCCGGTGAggtgtt < 1:152134/143‑6 (MQ=255) cgcATGACATTCCGTTGGATGTCGGCCGCCCTATATACGAAGCCGCTCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCata < 2:423554/141‑1 (MQ=255) gtaCATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATaaa < 2:151165/141‑1 (MQ=255) aTGTCAGCCGCCGTATATACGATGCTGGCCGCTTAGCTTTTTTCTTCTTATTGTTTCGTTTTCGGGCAATACGTCCAGACTGAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATgg < 1:131399/143‑1 (MQ=255) catcccccgTATATACTATGCCGCCCGCTAAGCTTTTTACGCCTCTATCGTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATCAAGGAGGAACGTATGGCTaa < 2:538516/137‑1 (MQ=255) ggcgccgTAGATGCGACGCCGTCTGCTAGGCCTTTTACGCCTCCAACTTTCGCTTTCGGGCAGTTCGTCCAGACTTAAGTGCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATc < 2:530936/141‑1 (MQ=255) ccgccgTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTACGTTTTCGGGCATTTCGTCCAGAGTTAAGCTCACAACACCTCACCGGAGCCTGCTCCGGGGAGTTCATATAAAGGAGGAACGTATGACTAATcc > 2:452265/1‑143 (MQ=255) ccgccgTATATACGAAACCGCCCGCTAAGCTTTTTACACCTCAAACTTTCGTTTTCGGGCATTTCGTAAAGACTTAAGGTAACAACACCTCACCGACGCATGCTCCGCAGAGTTCATATAAAGGAGGAACGAACAGCTAATca > 2:59770/1‑142 (MQ=255) tatataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTAGGGCATTTCGTCCAGACTTAAGTTCACAACACATCACCGGAGCCTGCTCCGGTGAGATCATATAACGGAAGAACGTATGGCTAATACAAcgat > 2:430193/1‑140 (MQ=255) tataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGAAGCCGGCTCCGGTGAGTTCATAAAAAAGAGGAACGTATGACTAATCCAACCGtta > 1:481632/1‑143 (MQ=255) tataCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCCTCCATACTTAAGTTCACAACACCTCACCAAAAACTGATCCGGTGAATTCATATAAAGGAGGAACGTATGGCTAATCCAACCGtta > 1:481646/1‑143 (MQ=255) ataCGAAGCCGCCCGCTAAGCTTTTTTCGCCTCAAACTTTCGTTTTGGGGCATTTCGTCCAGACTTAAGTTCACAAGACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGttat < 1:382624/143‑1 (MQ=255) taCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCattatt > 2:257533/1‑143 (MQ=255) aCGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTCAGTTCACAACACCTCACCGAAGCCTGCTCCGGTGAGTTCATTAAAAGGGGGAACGTCTGACTAATCCACTCGttatta > 2:110889/1‑143 (MQ=255) cccgcTAAGCTTTTTACGCCTCACATTTTCATGTGCGGGCATTTCGTCCAGACTTAACTTCACAACACCTCACCGGAGCAGGCTCCGGTGAGGTAATATAAAGGAGGAGGGTGGGGTGAATCCAACCGATATTAAGCTAccag > 2:343325/1‑140 (MQ=255) ccTACGCTCTTGACGGCGCACACTTGTGTTTGCGGGCATTTCGTCNAGACTTATGATCACAACACATAGCCGGAGCCTGCTCTGATCAGTTCATATCAAGTAGGAACGTATGGCTAATCCAACCGTTATTAAACTACAGGATg < 2:129616/142‑1 (MQ=255) gCTTTTTTCGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCATACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAAt < 1:465484/143‑1 (MQ=255) cTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTAATATAAAGGAGGAACGTATGGCTAATCCAAACGTTATTAAGCTACAGGATGGCAATg > 1:267339/1‑143 (MQ=255) tgtgttCGTCTGCGGTCATTTCCGCCAGATCTAAGTTCCAGACACCTCATCGGCGTGGGCTTGGGTGAGTTAATAGGAGGTAGACACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTg < 2:565626/139‑1 (MQ=255) ttttCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATAGCAATGTCATGCCCCAGCTGCGCCTggg > 2:258157/1‑143 (MQ=255) tCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGAGt > 1:515537/1‑143 (MQ=255) ggCATTTCGTCCATACTTTATTTCACAACACCTCTCCGGAGCTTGCTCCTGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTg < 1:532657/143‑1 (MQ=255) gCATTTCGTCCAGACTTAAATTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGAGTGGGCGCCTgg > 2:235389/1‑143 (MQ=255) gttcataCTTAAGTTCACATCACCTCTCCGTAGCCTGCTCCGGTTAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGgcaagcaa < 1:476419/137‑1 (MQ=255) tgagatCACAGCGCCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGTTATTCCAGTAATGatg < 2:309753/138‑3 (MQ=255) gTTCACAACACCTCACCGCAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGAATGGCAATGTCATGCCCCCGCTGGGACTGGGCGTCTGGCAAGCAAGTACTGAGGacg > 2:265586/1‑141 (MQ=255) cTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGCCTGGGCGTCTGGCAAGCAAGTCATGAGGAAGTAATCaccacc > 1:579418/1‑143 (MQ=255) aCCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATACCCCAACTGGAACTGAGCGTATGGCAAGCAAGTAATGAGGAAGTAATAACCGCCAtt > 1:132656/1‑143 (MQ=255) aCCGGAGCAGGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACAGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGGACGTAGTGAGGTAGTAATCAcagatatt > 1:240374/1‑136 (MQ=255) ccTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGTAGTAATCACCGGCATCGAAAAAg > 2:575930/1‑143 (MQ=255) | TGGGCGAAAATCATGACATTACGTTGGATGTCAGCCGCCGTATATACGAAGCCGCCCGCTAAGCTTTTTACGCCTCAAACTTTCGTTTTCGGGCATTTCGTCCAGACTTAAGTTCACAACACCTCACCGGAGCCTGCTCCGGTGAGTTCATATAAAGGAGGAACGTATGGCTAATCCAACCGTTATTAAGCTACAGGATGGCAATGTCATGCCCCAGCTGGGACTGGGCGTCTGGCAAGCAAGTAATGAGGAAGTAATCACCGCCATTCAAAAAG > NC_000913/3156457‑3156731 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |