Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | USA300TCH1516_ALE | 1,081,439 | G→T | 6.8% | R47I (AGA→ATA) | purN → | Phosphoribosylglycinamide formyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | USA300TCH1516_ALE | 1,081,439 | 0 | G | T | 6.8% | 163.2 / 4.1 | 60 | R47I (AGA→ATA) | purN | Phosphoribosylglycinamide formyltransferase |
Reads supporting (aligned to +/- strand): ref base G (20/35); new base T (2/2); total (22/38) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.24e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-02 |
ATGGTTAAAATTGCGATTTTTGCATCAGGTTCAGGAAGTAACTTTGAAAATATAGTTGAGCATGTTGAATCAGGAAAACTTGAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGATATTCCTGTTTATATTAATGAACCAAAACAATTTGATTCAAAAGCAGCGTATGAACAACATTTAGTAACATTATTAAATGAAGATAAGGTAGAGTGGATTATTCTAGCTGGC > USA300TCH1516_ALE/1081300‑1081566 | aTGGTTAAAATTGCGAGTTTTGCATCAGGTTCAGGAAGTAACTTTGAAAATATAGTTGAGCATGTTGAATCAGGAAAACTTGAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTatagatag < 1:489743/140‑1 (MQ=255) tAAAATTGCGATTTTTGCATCAGGTTCAGGAAGTAACTTTGAAAATATAGTTGAGCATGTTGAATCAGGAAAACTTGAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑aa > 1:405121/1‑140 (MQ=255) tAAAATTGCGATTTTTGCATCAGGTTCAGGAAGTAACTTTGAAAATATAGTTGAGCATGTTGAATCAGGAAAACTTGAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGAGAGAGC‑aa > 1:343505/1‑140 (MQ=255) tttttGCATCAGGTTCAGGAAGTAACTTTGAAAATATAGTTGAGCATGTTGAATCAGGAAAACTTGAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGat < 2:201812/140‑1 (MQ=255) tttttGCATCAGGTACAGGAAGTAACTTTGAAAATATAGTTGAGCATTTAGAATCAGGAAAACTTGAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGat < 1:432468/140‑1 (MQ=255) ttttgGCATCAGGTTCAGGAAGTAACTTTGAAAATATAGTTGAGCATGTTGAATCAGGAAAACTTGAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGat < 2:80427/140‑1 (MQ=255) ttttGCATCAGGTTCAGGAAGTAACTTTGAAAATATAGTTGAGCATGTTGAATCAGGAAAACTTGAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGata < 1:410122/140‑1 (MQ=255) tttGCATCAGGTTCATGAAGTAACTTTGAAAATATAGTTGAGCATGTTGAATCAGGAAAACTTGAAAATATTGAAGTTATGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGatat < 2:324653/140‑1 (MQ=255) tttGCAGCAGGGTCAGGAAGTAACTTTGAAAATATAGTTGAGCATGTTGAATCAGGAAAACTTGAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGatat < 2:654395/140‑1 (MQ=255) tGCATCAGGTTCAGGAAGTAACTTTGAAAATATAGTTGAGCATGTTGAATCAGGAAAACTTGAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGATATTc < 1:435358/140‑1 (MQ=255) gCATCAGGTTCAGGAAGTAACTTTGAAAATATAGTTCAGCATGTTGAATCAGGAAAACTTGAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGATATTcc < 1:263915/140‑1 (MQ=255) cATCAGGTTCAGGAAGTAACTTTGAAAATATAGTTGAGCATGTTGAATCAGGAAAACTTGAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGATATTCCt < 2:593670/140‑1 (MQ=255) gTTCAGGATGTAACTTGTAAAATATAGTTGAGCATGTGGTATCAGGAAAACTTGAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGATATTCCTGTTtat < 2:576598/140‑1 (MQ=255) cgaggaaCTTTGAAAATATAGTTGAGCATGTTGAATCAGGAAAACTTGAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGATATTCCTGTTTATATTAAt < 2:526106/135‑1 (MQ=255) accTTTGAAAATATAGTTGAGGATGTTGAATCGGTAAAGACGGCGTCTGGTGAGTGTACGGCGCT‑ATATATGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGATATTCCTGTTTATATTAATGAAcc < 2:339154/138‑1 (MQ=255) aacgtggAAAATATAGTTGAGCATGTGGAATCAGGAAAACTTGAAAATATTGAAGTTAAGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGATATTCCTGTTTATATTAATGAAcc < 1:436123/134‑1 (MQ=255) aCTTTGAAAATATAGTTGAGCATGTTGAATCAGGAAAACTGGAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGATATTCCTGTTTATATTAATGAACCa < 1:196533/140‑1 (MQ=255) ccttGAAAATATAGTTGAGCATGTTGCATCAGGAAAACTTGAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGATATTCCTGTTTATATTAATGAACCaa < 1:330581/138‑1 (MQ=255) cTTTGAAAATATAGTTGAGCATGTTGAATCAGGAAAACTTGAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGATATTCCTTTTTATATTAATGAACCaa < 1:685451/140‑1 (MQ=255) tttGAAAATATAGTTGAGCATGTTGAATCAGGAAAACTTGAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGATATTCCTGTTTATATTAATGAACCaaa > 1:128203/1‑140 (MQ=255) tGAAAATATAGTTGAGCATGTTGAATCAGGAAAACTTGAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGATATTCCTGTTTATATTAATGAACCAAAAc > 1:657604/1‑140 (MQ=255) aaTATAGTTGAGCATGTTGAATCAGGAAAACTTGAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGATATTCCTGTTTATATTAATGAACCAAAACAAtt > 1:416290/1‑140 (MQ=255) atatAGTTGAGCATGTTGAATCAGGAAAACTCGAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGGATAGATAGAGCAAAAAAAGCACGATATTCCTGTTTAAATTAATGAACCAAAACAAtt > 2:624425/1‑140 (MQ=255) atatAGTGGAGCATGTTGAATCAGGAAAACTTGAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGATATTCCTGTTTATATTAATGAACCAAAACAAttt < 1:458618/140‑1 (MQ=255) taGTTGAGCATGTTGAATCAGGAAAACTTGAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGATATTCCTGTTTATATTAATGAACCAAAACAATTTGAt < 1:526769/140‑1 (MQ=255) taGTTGAGCATGTTGAATCAGGAAAACTTGAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGATATTCCTGTTTATATTAATGAACCAAAACAATTTGAt > 2:584859/1‑140 (MQ=255) gTTTAGCATGTTGAATTAGGAAAACTTGTAAATATTAAACAT‑GGCCGCTGATTTACGGAT‑CATCAAAATGCGTTTTGTATATATAGAGC‑AAAAAAGCACGATATTCCTGTTTATATTAATGAACCAAAACAATTTGATTc < 2:159037/140‑1 (MQ=255) tttggCATGTTGAATCAGGAAAACTTGAAGATATGGAAGTTACGGCGCG‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGATACTCCTGTTTATATTAATGAACCAAAACAATTTGATTCa < 1:62328/136‑1 (MQ=255) tGAGCATGTTGAATCAGGAAAACTTGAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGATATTCCTGTTTATATTAATGAACCAAAACAATTTGATTCaa > 1:338586/1‑140 (MQ=255) aggaggatGAAGCATGAAAGCGGGAAAATATTGAAGTTACGGCGCT‑ATA‑CTGAATACAGCAGAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCTTGATACTCCTTGGTATATTAATGAACCAAAACAATTTGATTCaaaa < 2:305158/133‑1 (MQ=255) aTGTTGAATCAGGAAAACTTGAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGATATTCCTGTTTATATTAATGAACCAAAACAATTTGATTCAAAagca > 1:407122/1‑140 (MQ=255) tGAATCAGGAAAACTTGAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGATATTCCTGTTTATATTAATGAACCAAAACAATTTGATTCAAAAGCAGCGt < 1:743025/140‑1 (MQ=255) tGAATCAGGAAAACTTGAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGATATTCCTGTTTATATTAATGAACCAAAACAATTTGATTCAAAAGCAGCGt < 1:59345/140‑1 (MQ=255) gAATCAGGAAAACTTGAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATATAGC‑AAAAAAGCACGATATTCCTGTTTATATTAATGAACCAAAACAATTTGATTCAAAAGCAGCGTa > 2:206230/1‑140 (MQ=255) gAATCAGGAAAACTTGAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATATAGC‑AAAAAAGCACGATATTCCTGTTTATATTAATGAACCAAAACAATTTGATTCAAAAGCAGCGTa > 2:205437/1‑140 (MQ=255) cAGGAAAACTTGAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGATATTCCTGTTTATATTAATGAACCAAAACAATTTGATTCAAAAGCAGCGTATGaa < 1:425553/140‑1 (MQ=255) gAAAACTTGAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCCAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGATATTCCTGTTTATATTAATGCACCAAAACAATTTGCTTCAAAACCAGCGTGTGgacaa > 1:362980/1‑140 (MQ=255) gAAAACTTGAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGATATTCCTGTTTATATTAATGAACCAAAACAATTTGATTCAAAAGCAGCGTATGaacaa > 1:377268/1‑140 (MQ=255) aaaCTTGAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGATATTCCTGTTTATATTAATGAACCAAAACAATTTGATTCAAAAGCAGCGTATGAacaaca > 2:590698/1‑140 (MQ=255) aaCTTGAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGATATTCCTGTTTATATTAATGAACCAAAACAATTTGATTCAAAAGCAGCGTATGAACAACAt > 1:550304/1‑140 (MQ=255) cTTGAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAACGCCTTTTTTATAGATAGAGC‑AAAAAAGCACGATATTCCTGTTTATATTAATGAACCAAAACAATTTGATTCAAAAGCAGCGTATGAACAACAttt < 1:321636/140‑1 (MQ=255) tttacaaTCTTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATACGTTTTGTATAGATAGAGC‑AAAAAAGCATGATATTCCCTTTTATATCAATGAACCAAAACAATTTGATTCAAAAGCAGCGTATGAACAACATTTa < 1:201259/135‑1 (MQ=255) gAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATATAGC‑AAAAAAGCACGATATTCCTGTTTATATTAATGAACCAAAACAATTTGATTCAAAAGCAGCGTATGAACAACATTTAGt < 1:205437/140‑1 (MQ=255) gAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATATAGC‑AAAAAAGCACGATATTCCTGTTTATATTAATGAACCAAAACAATTTGATTCAAAAGCAGCGTATGAACAACATTTAGt < 1:206230/140‑1 (MQ=255) aaTATTGCAGTTACCGCGCT‑ATATACGTAT‑CATCAAAATGCGTTGTGTATACATACCGC‑AAAAAAGCACGATATTCCTGTTTATATTAATGAACCAAAACAATTTGATTCAAAAGCAGCGTATGAACAACATTTAGTAAc < 1:657678/140‑1 (MQ=255) aaTATTGAAGGTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGATATTCCTGTTTATATTAATGAACCAAAACAATTTGATTCAAAAGCAGCGTATGAACAACATTTAGTAAc < 1:247037/140‑1 (MQ=255) atTGAAGTTACTGCGCT‑ATATACTGAT‑CATGAATATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGATATTCCTGTTTATATTAATGAACCAAAACAATTTGATTCAAAAGCAGCGTATGAACAACATTTAGTAACatt < 2:254602/140‑1 (MQ=255) ttGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGATATTCCTGTTTATATTAATGAACCAAAACAATTTGATTCAAAAGCAGCGTATGAACAACATTTAGTAACatta > 1:368675/1‑140 (MQ=255) gAAGTGACGGCGGT‑ATATGCGGAT‑CAGGAAAATGCGTTTTGTAGAGATAGAGG‑AAAAAAGCACGATATTCCTGTTTATATTAATGAACCAAAACAATTTGATTCAAAAGCAGCGTATGAACAACATTTAGTAACattatt < 1:584703/140‑1 (MQ=255) cGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGATATTCCTGTTTATATTAATGAACCAAAACAATTTGATTCAAAAGCAGCGTATGAACAACATTTAGTAACATTATTAAATGaa > 1:642355/1‑140 (MQ=255) cGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGATATTCCTGTTTATATTAATGAACCAAAACAATTTGATTCAAAAGCAGCGTATGAACAACATTTAGTAACATTATTAAATGaa > 1:593855/1‑140 (MQ=255) ggCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGATATTCCTGTTTATATTAATGAACCAAAACAATTTGATTCAAAAGCAGCGTATGAACAACATTTAGTAACATTATTAAATGAAg < 2:498179/140‑1 (MQ=255) ggCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGATATTCCTGTTTATATTAATGAACCAAAACAATTTGATTCAAAAGCAGCGTATGAACAACATTTAGTAACATGATTAAATGAAg > 2:627513/1‑140 (MQ=255) gcgcT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGATATTCCTGTTTATATTAATGAACCAAAACAATTTGATTCAAAAGCAGCGTATGAACAACATTTAGTAACATTATTAAATGAAGa > 1:609313/1‑140 (MQ=255) atataCGTAT‑CATCATAATGCGTTTCGTATAGATAGAGA‑AAAAAAGCACGATATTCCTGTTTATATTAATGAACCAAAACAATTTCATTCAAAAGCAGCGTATGAACAACATTTAGTAACATTATTAAATGAAGATAAgg < 2:242019/140‑1 (MQ=255) ggAT‑CATCAAAATGCGTTTTGTATAGAGAGAAC‑AAAAAAGCACGATATTCCTGTTTATATTAATGAACCAAAACAATTTGATTCAAAAGCAGCGTATGAACAACATTTAGTAACATTATTAAATGAAGATAAGGTAGAGt < 2:110674/140‑1 (MQ=255) t‑catcaAATTGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGAGATTCCTGTTTATATTAATGAACCAAAACAATTTGATTCAAAAGCAGCGTATGAACAACATTTAGTAACATTATTAAATGAAGATAAGGTAGAGTGGa < 2:688756/140‑1 (MQ=255) catcgAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGATAGTCCTGTTTATATTAATGAAACAAAACAATTTGATTCAACACCAGCGTATGAACACCATTTAGTAACAGTATTAAATGGAGAGATGGTCGAGTGGaa > 2:734154/1‑139 (MQ=255) catcaAAATGCGTGTTGTATAGATAGAGC‑AAAAAAGCACGATATTCCTGTTTATATTAATGAACCAAAACAATTAGATTCAAAAGCAGCGTATGAACAACATTTAGTAACATTATTAAATGAAGATAAGGTAGAGTGGat > 2:49310/1‑140 (MQ=255) gtattGTAGAGATAGAGC‑AAAAAAGCACGATAGTCCTGTTTATATTAATGAACCAAAACAATTTGATTCAAAAGCAGCGTATGAACAACATTTAGTAACATTATTAAATGAAGATAAGGTAGAGTGGATTATTCTAGCTg < 1:359386/137‑1 (MQ=255) tttGGATAGATAGAGC‑AAAAAAGCACTATATTCCTGGTTATATTAATGAACCAAAACAATTTGATTCAAAAGCAGCGGATGAACAACATTTAGTAACATTATTAAATGAAGATAAGGTAGAGTGGATTATTCTAGCTGGc < 2:763263/140‑1 (MQ=255) | ATGGTTAAAATTGCGATTTTTGCATCAGGTTCAGGAAGTAACTTTGAAAATATAGTTGAGCATGTTGAATCAGGAAAACTTGAAAATATTGAAGTTACGGCGCT‑ATATACGGAT‑CATCAAAATGCGTTTTGTATAGATAGAGC‑AAAAAAGCACGATATTCCTGTTTATATTAATGAACCAAAACAATTTGATTCAAAAGCAGCGTATGAACAACATTTAGTAACATTATTAAATGAAGATAAGGTAGAGTGGATTATTCTAGCTGGC > USA300TCH1516_ALE/1081300‑1081566 |
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: ‑ |