Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP009273 | 2,885,943 | A→G | 9.7% | intergenic (+5/‑73) | queD → / → ygcN | 6‑pyruvoyl tetrahydrobiopterin synthase (PTPS)/putative oxidoreductase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP009273 | 2,885,943 | 0 | A | G | 9.7% | 128.9 / 2.8 | 46 | intergenic (+5/‑73) | queD/ygcN | 6‑pyruvoyl tetrahydrobiopterin synthase (PTPS)/putative oxidoreductase |
Reads supporting (aligned to +/- strand): ref base A (20/21); new base G (2/3); total (22/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.13e-01 |
ACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGTGCCGGTATTGCAGGCACCGCTTGCGCGTTACGCTGCGCGCGAGCGGGTTTATCCGTTTTGTTACTGGAACGCGCTGAAATCCCCGGCAGCAAAAATCTTTCCGGCGGGCGGTTATATACCCATGCACTCGCGGAACTCCTCCCGCAATTTCATCTGACCGCGC > CP009273/2885686‑2886211 | aCTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAAt < 1:210793/275‑1 (MQ=255) aCTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCAACGCAGGTTGTAGCGATCGCGGCGAATGATAAGAGTGTGTCGGCGGtta > 2:373758/1‑271 (MQ=255) ggAAGTCGGTCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACg < 1:373761/275‑1 (MQ=255) aTACGGGCTGGATTATAGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGGCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACGCGTTGTGCCGCTGTTAAGGGCGGTGATGGTAAAAGAAACATGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGgtg < 1:160376/274‑1 (MQ=255) aaaGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAACGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGt > 1:733911/1‑225 (MQ=255) cggcgTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATTGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACg < 2:733905/225‑1 (MQ=255) ccTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGTTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGTg > 1:413860/1‑274 (MQ=255) tATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGTGCCg > 2:739899/1‑264 (MQ=255) tctcAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGGTAAGAGTGTGTCGGCGGTCAATTTc > 1:750281/1‑173 (MQ=255) tcAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCAGCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGTGCCGGTAt < 1:739905/264‑1 (MQ=255) aaTGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCtt < 2:750275/173‑1 (MQ=255) tGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGt > 2:679999/1‑161 (MQ=255) gAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGgtgttgtgtt > 2:3583/1‑183 (MQ=255) aaaCCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACg < 1:680005/161‑1 (MQ=255) aCCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTgg < 1:3583/183‑1 (MQ=255) aaTGGATTTGGGATCAGGTTATAGGCGTTTTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGANGACTGCGACATTATTATTATTGGTGCCGGTATTGCAGGCACCGCTTGCGCGTTACGCTtcgcg > 1:571450/1‑247 (MQ=255) ttAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGGGTGTGTCGGCGGTCAATTTccc < 1:177160/107‑1 (MQ=255) caaaCCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCGGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTaa < 1:509754/113‑1 (MQ=255) ttGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGt > 1:170693/1‑89 (MQ=255) ttGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAgacga > 2:187982/1‑159 (MQ=255) ttGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGTGc > 1:539261/1‑183 (MQ=255) gCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAAt < 2:170692/89‑1 (MQ=255) gCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGc < 1:187985/159‑1 (MQ=255) gCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGTGCCGGt < 2:539258/183‑1 (MQ=255) gCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGTGCCGGTATTGCAGGCACCGCTTGCGCGTTACGCTGCGCGCGAGCGGGTTTATCCGTTTTGTTACTGGa < 1:59366/245‑1 (MQ=255) aGTGCGGTTATGGTAAAATATACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAATTAACTCTATGTTAGGGTGTTGTGTTCTGGATATNTGGGGCATGACACGGAAGACGACTGCGACATTATTATTATTGGTGCCGGTATTGCAGGCACCGCTTGCGCGTTACGCTGCGCGCGAGCGGGTTTATCCGTTTTGTTACTGGAACGCGCTGAAATCCCCGGCAGCAAAAATCTTTCCggcggg < 1:688528/275‑1 (MQ=255) tGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGt > 1:69450/1‑83 (MQ=255) ggTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACg < 2:69450/83‑1 (MQ=255) aCCTGCACCGCAGGTTGTAACTATCGCGGCGAATGATAAGGGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGTGCCGGTATTGCAGGCACCGCTTGCGCGTTACGCTGCGCGCGAGCGGGTTTATCCGTTTTGTTACTGGAACGCGCTGAAATCCCCGGCAGCAAAAATCTTTCCGGCGGGCGGTTATATACCCATGAACTc > 1:680603/1‑275 (MQ=255) ccTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGGGTGTGTCGGCGGTCAATTTCCCTTaa > 1:57404/1‑66 (MQ=255) ccTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGTGACGGTAt > 1:741336/1‑153 (MQ=255) cACCGCAGGTTGTATCTATCGCGGCGAATGATAAGGGTGTGTCGGCGGTCAATTTCCCTTAAGTaa < 2:57404/66‑1 (MQ=255) cACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGTGCCGGTATTGCa < 2:741330/153‑1 (MQ=255) cACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGTGCCGGTATTGCAGGCACCGCTTGCGCGTTACGCTGCGCGCGAGCGGGTTTATCCGTTTTGTTACTGGGACGCGCTGAAATCCCCGGCAGCAAAAATCTTTCCggcgggc > 1:6901/1‑250 (MQ=255) atctatcGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGgtgt > 1:139062/1‑69 (MQ=255) tatcGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCtt > 1:368086/1‑44 (MQ=255) atcGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGgtgttgtg < 2:139061/69‑1 (MQ=255) tcGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACt > 1:108348/1‑103 (MQ=255) tcGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGTGCCGGTATTGCAGGCACCGCTTGCGCGTTACGCTGCGCGCGAGCGGGTTTATCCGTTTTGTTACTGGAACGCGCTGAAATCCCCGGCAGCAAAAATCTTTCCGGCGGGCGGTTATATACCAATGCACTCGCGGAACTCCTCCCGCAAttta > 1:682984/1‑273 (MQ=255) cgcgGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGTGCCGGTATTGCAGGCACCGCTTGCGCGTTACGCTgcgcgc > 1:265402/1‑162 (MQ=255) gcggcgAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGt < 2:368083/44‑1 (MQ=255) ggcgAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGa < 2:108348/103‑1 (MQ=255) gcgAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTAGGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGTGCCGGTATTGCAGGCACCGCTTGCGCGTTACGCTGCGCGCGAGc < 2:265399/162‑1 (MQ=255) aaTGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGTGCCGGTATTGCAGGCa > 2:614718/1‑131 (MQ=255) gATATGGGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGTGCCGGTATTGCAGGCACCGCTTGCGCGTTACGCTGCGCGCGAGCGGGTTTATCCGTTTTGTTACTGGAACGCGCTGAAATCCCCGGCAGCAAAAATCTTTCCGGCGGGCGGTTATATACCCATGCACTCGCGGAACTCCTCCCGCAATTTCATCTGACCgcgc < 1:53696/275‑1 (MQ=255) aTAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTGTTATTGGTGCCGGTATTGCAGGCACCGc < 1:614724/131‑1 (MQ=255) | ACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTCCAGGTCTGGAAAACCCAACCAGCGAGGTTTTAGCAAAATGGATTTGGGATCAGGTTAAACCCGTTGTGCCGCTGTTAAGTGCGGTGATGGTAAAAGAAACCTGCACCGCAGGTTGTATCTATCGCGGCGAATGATAAGAGTGTGTCGGCGGTCAATTTCCCTTAAGTAACGCTATGTTAGGGTGTTGTGTTCTGGATATCTGGGGCATGACATGGAAGACGACTGCGACATTATTATTATTGGTGCCGGTATTGCAGGCACCGCTTGCGCGTTACGCTGCGCGCGAGCGGGTTTATCCGTTTTGTTACTGGAACGCGCTGAAATCCCCGGCAGCAAAAATCTTTCCGGCGGGCGGTTATATACCCATGCACTCGCGGAACTCCTCCCGCAATTTCATCTGACCGCGC > CP009273/2885686‑2886211 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |