Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP009273 | 2,145,170 | G→T | 26.4% | intergenic (+43/+22) | yegJ → / ← yegK | uncharacterized protein/ser/thr phosphatase‑related protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP009273 | 2,145,170 | 0 | G | T | 26.4% | 58.0 / 6.6 | 51 | intergenic (+43/+22) | yegJ/yegK | uncharacterized protein/ser/thr phosphatase‑related protein |
Reads supporting (aligned to +/- strand): ref base G (17/19); new base T (9/5); total (27/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.52e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01 |
GGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑TTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGC > CP009273/2145033‑2145317 | gggATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑TTGAATTGTTGTCTTCAACCGcc < 1:471068/150‑1 (MQ=255) ggATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑TTGAATTGTTGTCTTCAACCGCCa < 2:736608/150‑1 (MQ=255) cgtttATACATTTGATCTTTTATTGGCCCCCAACCAAAAAGAAAAAGCCGACGAGCAAAGAGAAAAATTATTCCCCTAATCAGCCGTTATTTTTTTTTTTTAAACGGGCGT‑TTAAATTGTTGTCTTCAACCGCCATGAGAAACTCActct < 2:606328/147‑1 (MQ=255) gTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑TTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGc > 1:681267/1‑150 (MQ=255) ttATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑TTGAATTTTTTTTTTCAACCGCCCTGAGAAAATCACTCTGcc > 2:971839/1‑150 (MQ=255) tatcCCATTGATCTTTTATTGGCCCCCTACCAAAAAAAAGAAGCCGACGAGCTAAGAGAAAAATTATTTCGCTAACCAGCCGTTATTTTTATTGTTTAACCGGGCGG‑TTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGcc < 1:971839/149‑1 (MQ=255) ttGATCTTTTATTGGTCACCTAACAAAAAAAAGAAGCCGACGAGAAAAGAAAAAAATTAGTTCGATAATCAGCCTTTATATTTATTTTTAAAACGGGCGG‑TTAATTGGTTGTCTTAAACCGCCATGAGAAACTCACTCTGCCCACAgcg < 2:48610/149‑1 (MQ=255) gATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑TTTAATTGTTTTCTTTAACCGCCATGAGAAACTCACTCTGCCCCCAgcgc > 2:101589/1‑148 (MQ=255) gATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑TTGAATTTTTGTTTTCAACCGCCATGAGAAACTCACTCTGCCCACAacgc > 2:831952/1‑148 (MQ=255) tcttttATTGCTCCCCTATAAAAAAGAAGAAGCCGACGAGCTAAGAAAAAAATTTTTCCGCTAATCAGCCGTTATATTTTTTGTTAAACCTGGCGT‑TTAAATTTTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAg < 2:1075648/146‑1 (MQ=255) tttATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑TTGAATTGTTTTCTTCAACCGCCCTGAGAAACTCACTCTGCCCACAGCGCCAGGGcc > 2:409637/1‑150 (MQ=255) tttttGGCCACCTACCCAAAAAAAGAACCCGAGAAGATAAGAGAAAAATTATTTCCTAAATCACCCTTTATATTTATTTTTAAAACGGGCGT‑TTGAATTTTTTTCTTAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTg < 1:233038/144‑1 (MQ=255) ttATTGGCCCCCTATCAAAAAGAAGAACCCGAGGAGCTAAGAAAAAAATTATTTCCCAAATCGGCCTTTATTTTTATTGTTTAACCTGGCGG‑TTAAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAg < 1:24345/145‑1 (MQ=255) ttGGTCACCTATACAAAAGAAGAAGCCGACGAGCTAAGAGAAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTAAAACTGGGGG‑TTGAATTGTTGTCTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAAc < 2:873441/149‑1 (MQ=255) tGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑TTGGATTGTTTTCTTCAACCGCCATGAGAAAATCACTCTGCCCACAGCGCCAGCGGCAACGt > 2:425191/1‑150 (MQ=255) cTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑TTGAATTGTTGTCTTCAACCGCCCTGAGAAACTCACTCTGCCCCCAGCGCCAGTGCCCACGGTTTttca > 2:640454/1‑150 (MQ=255) tATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑TTGAATTGTTGTCTTTAACCGCCATGAGAAAATCACTCTGCCCACAGCGCCAGGGCCAACGGTTTatcat > 1:519617/1‑150 (MQ=255) tCCAAAAGAAGAAGCCGACGAGATAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTAAAACGGGCGG‑TTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTatcatca < 1:179021/150‑1 (MQ=255) caaaaaaaAGAAGCCGACGAGCTAAGAGCAAATTTAGTTCGCTAATCAGCCTTTATATTTTTTTTTTAAACGGGCGG‑TTGATTTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTatcatcat < 1:28315/148‑1 (MQ=255) cAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑TTGAATTGTTGTCTTCAACCGCCATGAGAAAATCACTCTGCCCACAGCGCCAGTGCCAACGTTTTatcatcatc > 1:72564/1‑150 (MQ=255) aaaaaGAAGAAGCCGACGAGAAAAGAGAAAAATTAGTTCGTTAATCAGCCTTTATATTTTTTTTTTAAACGGGCGG‑TTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTatcatcat < 1:101370/148‑1 (MQ=255) aaaaaaGAACCCGACGAGCAAAGAACAAATTTATTCCCCTATTCACCTTTTTTTTTTTTTTTTTAACCGGGCGT‑TTAATTTTTTTTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCg < 2:584942/150‑1 (MQ=255) cgacgaGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑TTGAATTGTTTTCTTCAACCGCCATGAGAAAATCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGttt > 2:276557/1‑150 (MQ=255) cgacgaGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑TTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAACGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGtt > 2:304610/1‑149 (MQ=255) gacgaGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑TTGAATTGTTTTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTa > 1:846463/1‑150 (MQ=255) ggggCTAAGAGAAAATTTAGTTCGCTAATCAGCCTTTATTTTTATTGTTTAAACGGGCGG‑TTGAATTGTTGTCTCAACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTAcc < 2:1067666/147‑1 (MQ=255) cgaGGTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTTTTTAAACTGGCGG‑TTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTAcc < 2:277560/150‑1 (MQ=255) gcaaaGACCAAATTTATTCCGCTAATCAG‑CGTTTTTTTTTTTTTTTAAACCGGGGGTTTGAATTGTTTTCTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTAccgcc < 1:811712/147‑1 (MQ=255) aaGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑TTGGATTGTTGTCTTCAACCGCCATGAGAAAATCACTCTGCCCACAGCGCCAGGGCCAACGGTTTATCATCATCCGGACGCTCGCTTACCGCCggg > 1:306772/1‑150 (MQ=255) agCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑TTTAATTGTTTTCTTCAACCGCCATTAGAAACTCACTCTGCCCACAGCGCCAGTGCCCACGTTTTATCATCATCCGTACGCTCGTTTACCCCCGGGCTg > 1:713606/1‑150 (MQ=255) cAAAATTGTTTCGCTAATCAGCCTTTATATTTTTTGTTTAACCTGGGGG‑TTGATTTTTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGa < 1:353457/150‑1 (MQ=255) cAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑TTGAATTGTTTTCTTCAACCGCCCTGAGAAACTCACTCTTCCCACAGCGCCAGTGCCAACGTTTTTTCATCATCCCTACGCTCGTTTACCCCCGGGCTggg > 2:480661/1‑149 (MQ=255) aaaaaTTATTTCGTTAACCAGCTTTTATATTTATTTTTAAAACTGGCGT‑TTGATTTTTTGTCTTAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGa < 1:101589/149‑1 (MQ=255) aTTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑TTGAATTGTTGTCTTCAAACCCCAAGAGAAAATCACTCTGCCCACAGCGCCAATGCCAACGTTTTATTATTATTCGTACGCCCGTTTTCCGCCGGGCTGGGCaa > 1:1102004/1‑149 (MQ=255) tAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑GTGAAATGTTGTCTTCAACCGCCATGAGAAAATCACTCTGCCCACAGCGCCAGTGCCAAAGGTTTATTATCATCCGTAAGCCCGGTTAACGCCCGGCTGGGCaaaa > 2:1076385/1‑149 (MQ=255) ttCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑TTGAATTGTTGTCTTCAAACGCCATGAGAAACTCACTCTGCCCACAACGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAAATg > 1:387521/1‑150 (MQ=255) tCAGCCGTTATTTTTATTTTTTAAACGGGCGG‑TTGATTTGTTTTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACa < 2:681267/150‑1 (MQ=255) cGTTATATTTATTGTTTAAACTGGCGG‑ATGAATTTTTGTCTTCAACAGCCCTGAGAAACACACTCTCCCCCCTGCGCCAGGGCCAACGATTTTTCATAAATCGGACGCCCCTTTACCGCCGGGCGGGGGAAAAAACGCttttaaaaat > 1:821851/1‑141 (MQ=255) ttATATTTATTGTTTAAACTGGCGG‑TTGAATTGCTCTCTTCCACCGTCATTAGAAACTCCCTCTGCCCACACCGCCACAGTCCACGTTTTATCATCATCCGTTCCCTCGTTTACCCCCCGGCTGGACAAAAACCGCTTTAACAATTGAgg > 1:718967/1‑150 (MQ=255) atttattGTTTAAACTGGCGG‑TTGAATTTTTGTCTTCAACCCCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCCACGTTTTATCATCATCCGTACGCTCGTTTTCCCCCGGGCTGGACAAAAACTGCTTTTACAATTCCGGgaa > 2:441306/1‑149 (MQ=255) ttttttttTAAAACGGGGGT‑TTAATTTGTTGTTTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGtaat < 1:34444/149‑1 (MQ=255) tttattGTTTAAACTGGCGG‑TTGAATTTTTTTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTTTCATCATCCGTACGCTCGTTTACCGCCCGGCTGGGCAAAAACTGCTTTAACAATTCAGGgaat > 2:48136/1‑149 (MQ=255) tAAACTGGCGG‑TTGAATTGTTGTCTTTAACCGCCATGAGAAACTCACTCCGCCCCCAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAgtt > 2:681448/1‑150 (MQ=255) tAAACTGGCGG‑TTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATTGAgtt > 1:758875/1‑150 (MQ=255) cgggggg‑TTGATTTGTTGTCTCAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTc < 1:640454/145‑1 (MQ=255) cTGGCGG‑TTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTc > 2:641670/1‑150 (MQ=255) ggCGG‑TTAAATTGTTGTGTTAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCt < 2:158804/150‑1 (MQ=255) tGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTtg > 2:710685/1‑150 (MQ=255) aaTTTTTTTCTTAAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTtgc < 1:871558/149‑1 (MQ=255) aaTTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTtgct > 1:627623/1‑150 (MQ=255) aaTTGTTGTCTCCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTtgct < 1:216619/150‑1 (MQ=255) tttgtttTCTCCACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTtgctg < 2:306772/149‑1 (MQ=255) ttgtCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTgccgc > 1:1041453/1‑150 (MQ=255) | GGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGG‑TTGAATTGTTGTCTTCAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATCATCATCCGTACGCTCGTTTACCGCCGGGCTGGACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAGTTGTTCCTGCGTTGCTGCCGC > CP009273/2145033‑2145317 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |