Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 1,532,490 | G→A | M225I (ATG→ATA) | ydcC → | H repeat‑associated putative transposase YdcC |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,532,490 | 0 | G | A | 96.0% | 58.1 / ‑3.3 | 25 | M225I (ATG→ATA) | ydcC | H repeat‑associated putative transposase YdcC |
Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (15/9); total (16/9) | |||||||||||
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.87e-01 |
GATGAGAAATCTAATGAGATTACAGCTATCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGCGCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAG > NC_000913/1532251‑1532733 | gATGAGAAATCTAATGAGATTACAGCTATCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGCGCATGACAGTTACGCAATAAGTGAAAagag < 1:335625/251‑1 (MQ=18) aaaTCTAATGAGATTACAGCTATCCCAGAACTTCGTAACAAGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATAGTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTGTTCGCGGTAAAAGGAAACCAGGGGCGGCTGAATAAAGCCTTTGAGGAAAAATTTCCGCGGAAAGAATTAAATAATCCAGCGCATGATAGGTACGCAATAAGGGAAAAGAGTCACgg < 2:387912/251‑1 (MQ=11) tCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGCGCATGACAGTTACGCAATAAGTGAAAAGAGTCACGGCAGAGACGAAATCCGTCTTCat > 2:52691/1‑251 (MQ=18) ggATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGCGCATGACAGTTACGCAATAAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCtgatg > 1:88405/1‑251 (MQ=18) gATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGCGCATGACAGTTACGCAATAAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCtgatg > 1:276154/1‑250 (MQ=18) atatTAAAGCAAAAATCATCAAAACTGATGCGACGGGGTGACAGAAAGATATTGCAAAGAAGATACAAAAACCGGGAGGTGATTATTTATTCGCGGAATAAGGACACCAGGGGCGGCTAAATAAAGCCTTTGAGGCAAAATTTACGCTGAAAGAATTAAATAATCCCGCGCATGACAGCTACGCAATAAGTGAAAAGATTCACGGCAGAGAAGTACTCCGTCTTCATATTGTTTGCGATGTCCCtgatga > 2:20578/1‑250 (MQ=11) acaACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGCGCATGACAGTTACGCAATAAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTg > 1:111013/1‑250 (MQ=18) gATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGCGCATGACAGTTACGCAATAAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGa > 1:169079/1‑250 (MQ=18) gCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGCGCATGACAGTTACGCAATAAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTgaagaa < 1:52691/250‑1 (MQ=18) aGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGCGCATGACAGTTACGCAATAAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAAtt > 1:9060/1‑250 (MQ=18) aaaCAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGCGCATGACAGTTACGCAATAAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTAGTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCat < 2:203656/251‑1 (MQ=18) atttattCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGCGCATGACAGTTACGCAATAAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACaaa > 1:252904/1‑249 (MQ=18) aaaaGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGCGCATGACAGTTACGCAATAAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTGTGCGATGTCCCTGATGAACTTAGTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCTCGaaa < 1:20803/251‑1 (MQ=255) aGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGCGCATGACAGTTACGCAATAAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCTCGAAATGa < 1:296738/251‑1 (MQ=37) aaataaaGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGCGCATGACAGTTACGCAATAAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCAc > 2:408938/1‑147 (MQ=17) aaataaaGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGCGCATGACAGTTACGCAATAAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCAc < 1:408938/147‑1 (MQ=17) aaaGAATTAAATAATCCAGAGCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGc > 2:58592/1‑251 (MQ=37) ataatCCAGCGCATGACAGTTACGCAATAAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATAAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATc > 2:324182/1‑251 (MQ=255) taatCCAGCGCATGACAGTTACGCAATAAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATcc > 1:195528/1‑251 (MQ=255) ccAGCGCATGACAGTTACGCAATAAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTg > 2:301695/1‑67 (MQ=255) ccAGCGCATGACAGTTACGCAATAAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTg < 1:301695/67‑1 (MQ=255) aGCGCATGACAGTTACGCAATAAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAAcc > 1:129189/1‑251 (MQ=255) gcATGCCAGTTACGCAATAAGTGAAAAGAGTCACGGCAGAGAAGAAATACGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGACTGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCTCGAACTGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATcctataccaat > 2:323462/1‑242 (MQ=14) ttACGCAATAAGTGAAAAGAGTCACGGCAGAGAAGAACTCCGTCTTCATATTGTTTGCGATGTCCCTGACGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTgg > 2:76810/1‑251 (MQ=255) aCGCAATAAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGgag < 1:65144/251‑1 (MQ=255) | GATGAGAAATCTAATGAGATTACAGCTATCCCAGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGTTGCCAGAAAGATATTGCAGAGAAGATACAAAAACAGGGAGGTGATTATTTATTCGCGGTAAAAGGAAACCAGGGGCGGCTAAATAAAGCCTTTGAGGAAAAATTTCCGCTGAAAGAATTAAATAATCCAGCGCATGACAGTTACGCAATGAGTGAAAAGAGTCACGGCAGAGAAGAAATCCGTCTTCATATTGTTTGCGATGTCCCTGATGAACTTATTGATTTCACGTTTGAATGGAAAGGGCTGAAGAAATTATGCGTGGCAGTCTCCTTTCGGTCCATAATAGCAGAACAAAAGAAAGAGCTCGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCATGTGGAG > NC_000913/1532251‑1532733 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |