| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | CP000730 | 1,076,774 | T→C | 19.2% | W51R (TGG→CGG) | purL → | phosphoribosylformylglycinamidine synthase subunit II |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | CP000730 | 1,076,774 | 0 | T | C | 19.2% | 59.5 / 8.5 | 32 | W51R (TGG→CGG) | purL | phosphoribosylformylglycinamidine synthase subunit II |
| Reads supporting (aligned to +/- strand): ref base T (12/13); new base C (3/3); total (16/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 = 9.85e-01 | |||||||||||
ATCGAACCAAGCGTTGAAGAAATTAAACTTGAAAAAGTATATCAAGATATGGGATTAAGTGATCAAGAATATGAAAAAGTTTGCGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAATTTCCTACGTCAGGTGACCATGTGCTTATGGGGCCTGGTGAAGGTGCAGGGGTAGTCGATATAGGTGATAATCAAGCCGTAGTATTTAAA > CP000730/1076636‑1076911 | aTCGAACCAAACGATGAAGAAACTAAACTTGAAAAAGTAGATCCAGATATGGGATTAAGTGAACAAGAATATGAAAAAGATTGCGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTgg < 2:285593/141‑1 (MQ=255) gaagaaATTAAACTTGAAAAAGTATATCAAGATATGGGATTAAGTGATCAAGAATATGAAAAAGTTTGCGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTATTCTGTTATGTGGAGTGAACATTGctct > 1:249074/1‑141 (MQ=255) aaCTTGAAAAAGTATATCAAGATATGGGATTAAGTGATCAAGAATATGAAAAAGTTTGCGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTACTCTTATAAACAtt > 1:33111/1‑141 (MQ=255) aCTTGAAAAAGTATATCAAGATATGGGATTAAGTGATCAAGAATATGAAAAAGTTTGCGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCGCTTATAAACATTc > 1:284293/1‑141 (MQ=255) aaaaGTATATCAAGATATGGGATTAAGTGATCAAGAATATGAAAAAGTTTGCGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAAcc < 2:326103/141‑1 (MQ=255) aaaaGTATATCAAGATATGGGATTAAGTGATCAAGAATATGAAAAAGTTTGCGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGCGGAGTGAACATTGCTCTTATAAACATTCTAAAcc > 2:79487/1‑141 (MQ=255) aaaaGTATATCAAGATATGGGATTAAGTGATCAAGAATATGAAAAAGTTTGCGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGCGGAGTGAACATTGCTCTTATAAACATTCTAAAcc > 2:320932/1‑141 (MQ=255) aaaaGTATATCAAGATATGGGATTAAGTGATCAAGAATATGAAAAAGTTTGCGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGCGGAGTGAACATTGCTCTTATAAACATTCTAAAcc > 2:28420/1‑141 (MQ=255) tataTCAAGATATGGGATTAAGTGATCAAGAATATGAAAAAGTTTGCGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGtttt > 1:652/1‑141 (MQ=255) cAAGATATGGGATTAAGTGATCAAGAATATGAAAAAGTTTGCGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAg < 2:3971/141‑1 (MQ=255) aaaGATATGGGATTAAGTGATCAAGAATATGAAAAAGTTTGCGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAg < 2:3975/140‑1 (MQ=255) tatGGGATTAAGTGATCAAGAATATGAAAAAGTTTGCGATATTTTAGTCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAAtt < 2:354207/141‑1 (MQ=255) atGGGATTAAGTGATCAAGAATATGAAAAAGTTTGCGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAAttt < 2:130164/141‑1 (MQ=255) atGGGATTAAGTGATCAAGAATATGAAAAAGTTTGCGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAAttt < 1:215333/141‑1 (MQ=255) tGGGACTAAGTGATCAAGAATATGAAAAAGTTTGCGATATTTTAGGCAGACAACCTAACTTTACACAATCAGGTATCTTTTCTGATATCGGCAGTGAACAAAGCTCTTATAAACATTCTAAACCGCTTTTCACGCAATTTc > 2:247063/1‑141 (MQ=255) aGTGAACAAGAAGATGAAAAAGTTTGCGATAGTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAATTTCCTACGTCa < 1:146348/141‑1 (MQ=255) cAAGAATATGAAAAAGTTTGCGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAATTTCCTACGTCAGGTGAc > 1:213828/1‑141 (MQ=255) atGAAAAAGTTTGCGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAATTTCCTACGTCAGGTGACCATGTGc > 2:159550/1‑141 (MQ=255) atGAAAAAGTTTGCGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAATTTCCTACGTCAGGTGACCATGTGc < 1:120089/141‑1 (MQ=255) aGTTTGCGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAATTTCCTACGTCAGGTGACCATGTGCTTATggg > 2:90999/1‑141 (MQ=255) cGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTTAATTTATTTCCTACGTCAGGTGGCCATGTGCTTATGGGGCCTgg > 2:89716/1‑141 (MQ=255) tatTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAATTTCCTACGTCAGGTGACCATGTGCTTATGGGGCCTGGTGa < 2:139700/141‑1 (MQ=255) aGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAATTTCCTACGTCAGGTGACCATGTGCTTCTGGGGCCTGGTGAAGGTGCAggg < 2:256221/141‑1 (MQ=255) aaCCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGCGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAATTTCCTACGTCAGGTGACCATGTGCTTATGGGGCCTGGTGAAGGTGCAGGGGTAg < 1:28420/141‑1 (MQ=255) aaCCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGCGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAATTTCCTACGTCAGGTGACCATGTGCTTATGGGGCCTGGTGAAGGTGCAGGGGTAg < 1:79487/141‑1 (MQ=255) aaCCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGCGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAATTCCCTACGTCAGGTGACCATGTGCTTATGGGGCCTGGTGAAGGTGCAGGGGTAg < 1:320932/141‑1 (MQ=255) cTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAATTTCCTACGTCAGGTGACCATGTGCTTATGGGGCCTGGTGAAGGTGCAGGGGTAGTCg > 2:115992/1‑141 (MQ=255) cAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAATTTCCTACGTCAAGTGACCATGTGCTTATGGGGCCTGGTGAAGGTGCAGGGGTATTCGATATAGGTGatcatc > 2:152505/1‑141 (MQ=255) aaGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAATTTCCTACGTCAGGTGACCATGTGCTTATGGGGCCTGGTGAAGGTGCAGGGGTAGTCGATATAGGTGATAATc < 1:308289/140‑1 (MQ=255) cTTTTCTGTTATGTGGAGTGAACGTTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAATTTCCTACGTCAGGTGACCATGTGCTTATGGGGCCTGGTGAAGGTGCAGGGGTAGTCGATATAGGTGATAATCAAGCCgg > 2:256472/1‑140 (MQ=255) ttttCTTTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAATTTCCTACGTCAGGTGACCATGTGCTTATGGGGCCTGGTGAAGGTGCAGGGGTAGTCGATATAGGTGATAATCAAGCCgta < 2:325688/141‑1 (MQ=255) aTGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAATTTCCTACGTCAGGTGACCATGTGCTTATGGGGCCTGGTGAAGGTGCAGGGGTAGTCGATATAGGTGATAATCAAGCCGTAGTATTTaaa > 1:337946/1‑141 (MQ=255) | ATCGAACCAAGCGTTGAAGAAATTAAACTTGAAAAAGTATATCAAGATATGGGATTAAGTGATCAAGAATATGAAAAAGTTTGCGATATTTTAGGCAGACAACCTAACTTTACAGAAACAGGTATCTTTTCTGTTATGTGGAGTGAACATTGCTCTTATAAACATTCTAAACCGTTTTTAAAGCAATTTCCTACGTCAGGTGACCATGTGCTTATGGGGCCTGGTGAAGGTGCAGGGGTAGTCGATATAGGTGATAATCAAGCCGTAGTATTTAAA > CP000730/1076636‑1076911 |
| 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: ‑ |