Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | USA300TCH1516_ALE | 2,373,809 | A→T | 15.4% | intergenic (‑34/+79) | pbuG ← / ← topB | Guanine/hypoxanthine permease PbuG/DNA topoisomerase 3 |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | USA300TCH1516_ALE | 2,373,809 | 0 | A | T | 15.4% | 137.3 / 13.2 | 62 | intergenic (‑34/+79) | pbuG/topB | Guanine/hypoxanthine permease PbuG/DNA topoisomerase 3 |
Reads supporting (aligned to +/- strand): ref base A (28/24); new base T (7/3); total (35/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.91e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01 |
AACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCAAGTTCTTTAATGCATCTTTAAACGGATTATTGTCTAAACCTTCA > USA300TCH1516_ALE/2373674‑2373942 | aaCTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGttt < 1:160234/140‑1 (MQ=255) cggCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTCTCGCACTGCAAATAGTTTTTCACGAATGGATGCTCCTTTTTGTAAAATACCTTAATAGTTTaa < 1:100140/138‑1 (MQ=255) gCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAata < 1:589589/140‑1 (MQ=255) tAAAATATATGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATc < 2:233639/140‑1 (MQ=255) atatAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTc > 1:35124/1‑140 (MQ=255) atatAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTc > 1:367697/1‑140 (MQ=255) aGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTtaa > 1:571178/1‑140 (MQ=255) aGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGATCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTtaa > 1:666867/1‑140 (MQ=255) cATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATaaa > 1:228412/1‑140 (MQ=255) cATAGATAAGAAAGTTGTGATACAGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACATTAATAGTTTAATATCATTCTtaataca > 2:392450/1‑138 (MQ=255) ataataAAGTTGTGATACCGCCTAAGATTTCTATTTTAAAGTTTGTTTCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACa < 2:472976/140‑1 (MQ=255) tAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATATCTTAATAGTTTAATATCATTCTTAATaaaataaaa > 2:338860/1‑135 (MQ=255) aGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTACTTGTTTAATATCTTTCTTCATAATAAACAATa > 2:344101/1‑140 (MQ=255) aaaGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAACATACCTGAATTTTTTACTATAATTCTGATTAATAAACAGTaac > 2:244011/1‑138 (MQ=255) gTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCaaa > 2:434505/1‑140 (MQ=255) gTTGTAATACCGACTAAGATTTCTATTTTATAGTTTGTTCCATATATATCGAAATGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCaaa < 2:392371/140‑1 (MQ=255) ttGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCaaaa > 2:619200/1‑140 (MQ=255) gATACAGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACTAAAAGCGa > 2:580314/1‑140 (MQ=255) caccTAATATTTCGATTTTAAAGTTTGTTCCATATTTATCCAAATGGAAAATAACAATCACGAATGGCTGCTCCGCGATTTAAAATACCTTAATAGTTTAATATAATTCTTAATAAAAAACAATACCAAAAGCGAACttt < 2:315321/138‑1 (MQ=255) ccTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTTAAATACCTTAATAGTTTAATATCATTCTTAATAAAAACCAATACCAAAAGCGTCCTTTTc > 1:589258/1‑140 (MQ=255) ccTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTATAATACCTTAATTGTTTAATTTCATTCTTTATTAAAATCAATACCAAAAGCATTCTTTTc > 2:410627/1‑140 (MQ=255) cggaGAATTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGGTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCGTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCg < 2:593930/137‑1 (MQ=255) tAAGATTTCTCTTTTAAAGTTTGTGCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGt < 1:37963/140‑1 (MQ=255) aGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATTGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGttt < 1:675147/140‑1 (MQ=255) aGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGttt < 1:434505/140‑1 (MQ=255) tctcTTTTAAAGTTTGTTACATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATTGTTTAATATCATTCTTAATAAAAAACAATACAAAAAGAGAACTTTTTGTATATCTg > 2:366404/1‑140 (MQ=255) ttttAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACTAAAAGCGAACTTTTCGTTTTTTTGATTc > 2:612515/1‑140 (MQ=255) tttAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCa > 1:600316/1‑140 (MQ=255) aaaGTGTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATTGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACt < 2:589525/140‑1 (MQ=255) ttGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACTTTAATTGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGtt > 2:588408/1‑140 (MQ=255) ttGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATTGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCTAACTTTTCTTTTTTCTGATTCAACtttttt > 2:119526/1‑140 (MQ=255) ccATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGAtt > 2:2178/1‑140 (MQ=255) ccATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGAtt > 1:210779/1‑140 (MQ=255) tATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCta > 1:253496/1‑140 (MQ=255) cTGGAAATAGCTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTc < 2:94033/140‑1 (MQ=255) tGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCa > 1:193439/1‑140 (MQ=255) tGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCa > 1:499657/1‑140 (MQ=255) ggAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTTTTTTTTTTTAATAGTTTAATATCATTCTTAATTATTTACAATACCAAAAGCGTACTTTTCGGTTTTCTGATTCATCTTTGTTCGTTTTTTTCTTTAATTTTaa > 2:31509/1‑140 (MQ=255) ggAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATTCCTTAATATTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGCTTCAACTTTGTTCGATTTTATCTATAAATTCaa > 1:187458/1‑140 (MQ=255) ggAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTCTCTATAAATTCaa > 2:453441/1‑140 (MQ=255) ggAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCaa > 1:231508/1‑140 (MQ=255) gAAATATGTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTAAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCAAg < 2:52714/140‑1 (MQ=255) gAAATAGTTTTTGACGCTTGGATGCTACGTTTTTTAAAATAACTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCAAg < 2:144278/140‑1 (MQ=255) gAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCAAg < 1:385886/140‑1 (MQ=255) gAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCAAg > 1:704171/1‑140 (MQ=255) aaaaTAATTTTTCACGAAGAGATGATCATTTTTTTAAAATAACTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCAAg < 2:338711/139‑1 (MQ=255) aaaTAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAATTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATTTATAATTTCaatt > 2:329976/1‑138 (MQ=255) aTAGTTTTTCACGAATGGATGCT‑CTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAATTTCAAGTTCt > 2:324764/1‑140 (MQ=255) tAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCAAGTTCt < 1:2178/140‑1 (MQ=255) ttttCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCAAGTTCTTTaa < 1:303350/140‑1 (MQ=255) cGAATGGATGCTCCTTTTTTTAAAATACTTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGAGAACTTTTCGTTTTTCTTATTCAACTTTGTTCCATTCTATCTATAAACTCAAGTTCTTTTAtctatt > 2:457009/1‑135 (MQ=255) aaTGGATGCTCCTTTTTTTAAAATACTTTATTAGTTTAATATCATTTTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTACCTGCGTTGTTATCTATATATTCAAGTTCTTTAATGTATCtt > 2:311990/1‑140 (MQ=255) aTGGATGCTCCTTTTTTTAAAATACCTTAATTGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCAAGTTCTTTAATGCATCttt < 1:561651/140‑1 (MQ=255) gATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCAAGTTCTTTAATGCATCTTTaaa < 2:571178/140‑1 (MQ=255) gATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCAAGTTCTTTAATGCATCTTTaaa < 2:666867/140‑1 (MQ=255) gATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCAAGTTCTTTAATGCATCTTTaaa < 2:704171/140‑1 (MQ=255) cTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCAAGTTCTTTAATGCATCTTTAAACGGa < 1:619200/140‑1 (MQ=255) ttttAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCAAGTTCTTTAATGCATCTTTAAACGGATTATTGt < 2:69644/140‑1 (MQ=255) ttttttttAATAGTTTAATATCATTCTTAATACATTTCAATACCATAAGCGTACTTTTCGTTTTTCTGAATCCACTTTTTTTTTTTTTATTTATAACTTCAACTTCTTTAATTCTTCTTTAAACGGATTATTTTCTAAAc > 2:568235/7‑140 (MQ=255) aTACCTTAATTGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCAAGTTTTTTAATGTATCTTTAAATGGATTATTGTCTAAAcc > 2:141795/1‑140 (MQ=255) ccTTGAGAGTTTAATAGCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTGTCTGATTCAACTTTGTTCGATTTTATCTATCAATTCAAGTTCTTTAAGGCATCTCTAAACGGATTATTGTCTAAACCTTc < 2:13424/140‑1 (MQ=255) taataTAGTTTAATACCAGTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCAAGTTCTTTAATGCATCTTTAAACGGATTATTGTCTAAACCTTCa < 2:517448/136‑1 (MQ=255) | AACTGCTAAAATATAGGCCATAGATAAGAAAGTTGTGATACCGCCTAAGATTTCTCTTTTAAAGTTTGTTCCATATTTATCGAACTGGAAATAGTTTTTCACGAATGGATGCTCCTTTTTTTAAAATACCTTAATAGTTTAATATCATTCTTAATAAAAAACAATACCAAAAGCGAACTTTTCGTTTTTCTGATTCAACTTTGTTCGATTTTATCTATAAATTCAAGTTCTTTAATGCATCTTTAAACGGATTATTGTCTAAACCTTCA > USA300TCH1516_ALE/2373674‑2373942 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |