Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP000730 | 1,117,630 | C→T | 8.0% | intergenic (+3416/‑104) | USA300HOU_1046 → / → USA300HOU_1047 | hypothetical protein/GTP‑binding protein TypA |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP000730 | 1,117,630 | 0 | C | T | 8.0% | 166.5 / 3.7 | 62 | intergenic (+3416/‑104) | USA300HOU_1046/USA300HOU_1047 | hypothetical protein/GTP‑binding protein TypA |
Reads supporting (aligned to +/- strand): ref base C (27/30); new base T (2/3); total (29/33) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01 |
CCAAAGGTTGCAACGAGTAATATAAATGTTAAAATCATCATTGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTATTTATATGACTAATAAAAGAGAAGATGTCCGCAATATAGCAA > CP000730/1117490‑1117770 | ccAAAGGTTGCAACGAGTAATATAAATGTTAAAATCATCATTGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGA‑TTTTTTTTGTTTc > 1:168775/1‑141 (MQ=255) aaaGGTTGCAACGAGTAATATAAATGTTAAAATCATCATTGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGA‑TTTTTTTTGTTTCAt > 1:397813/1‑141 (MQ=255) aaaGGTGGCAACGAGTAATATAAATGTTAAAATCATCATTGGTACATTTGATGCACCAAAAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGA‑TTTTTTTTGTTTCAt < 2:938322/141‑1 (MQ=255) aaGGTTGCAACGAGTAATATAAATGTTAAAATCATCATTGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGATTTTTTTTTGTTTCAt > 1:506953/1‑141 (MQ=255) cAACGAGTAATATAAATGTTAAAATCATCATTGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATACTTTAAGA‑TTTTTTTTGTTTCATCCTATccc > 2:50732/1‑141 (MQ=255) gAGTAATATAAATGTTAAAATCATCATTGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGA‑TTTTTTTTGTTTCATCCTATCCCCAcc > 1:786897/1‑141 (MQ=255) aGTAATATAAATGTTAAAATCATCATTGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGA‑TTTTTATTGTTTCATCCTATCCCCACCt < 1:694197/141‑1 (MQ=255) aGTAATATAAATGTTAAAATCATCATTGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCAAATATAGAAAATACCCAAAAGATATTTTTAGA‑TTTTTTTTGTTTCATCCTATCCCCACCt < 1:103816/141‑1 (MQ=255) tAATATAAATGTTAAAATCATCATTGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGA‑TTTTTTTTGTTTCATCCTATCCCCACCttt < 2:75254/141‑1 (MQ=255) aaTATAAATGTTAAAATCATCATTGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGA‑TTTTTTTTGTTTCATCCTATCCCCACCtttt > 2:502824/1‑141 (MQ=255) ataAACGTTAAAATCATCATTAGCACACTTCATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCt < 1:611056/141‑1 (MQ=255) aggTTAAAATAATCATTGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATc < 2:367725/139‑1 (MQ=255) aTGTTAAAATCATCATTGGTAAATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATc < 2:343501/141‑1 (MQ=255) ttAAAATCATCATTGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCttt < 1:542924/141‑1 (MQ=255) aaaTCATCATTTGTACATTTGATGCACCAACTGCAAAACTAAAAAATACTAAGAACACAACAGCTAATATAGAAAATACCCAAACGATATTTTTAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTtct < 2:892912/141‑1 (MQ=255) catcatTGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCCAAAGATAATTTTAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAAt < 1:743578/141‑1 (MQ=255) catcatTGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGA‑TTTTTTTTGTTTTATCCTATCCCCACCTTTTGCTTATCTTTTCTCAAt < 1:387729/141‑1 (MQ=255) catcatTGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGA‑TTTTTTTTGTTTTATCCTATCCCCACCTTTTGCTTATCTTTTCTCAAt < 1:454704/141‑1 (MQ=255) catcatTGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGA‑TTTTTTTTGTTTTATCCTATCCCCACCTTTTGCTTATCTTTTCTCAAt < 1:710259/141‑1 (MQ=255) tGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTTTCTTTTCTCAATTATATg > 2:929182/1‑141 (MQ=255) tGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTatata > 1:185358/1‑140 (MQ=255) tGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATg > 1:623245/1‑141 (MQ=255) tGATGCACCAACAGCAAAACTAAATAATAATAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATAtaataa < 1:791581/141‑1 (MQ=255) accaccAGCAAAACTAAACCATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAATATATTTTTAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTGAATTATATGATATAATAAAAAAGt < 1:601390/141‑1 (MQ=255) ccAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATTAAAAAGtt > 2:925831/1‑141 (MQ=255) aCAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAATTACCAAAAAGAAATTTTTAGA‑ATTTTTTTGTTTCTTTCTATACCCACCGTTTGCTTATCTTTACACAATAATATGATATAATACAAAAGTTgaa > 2:53802/1‑139 (MQ=255) cAGGAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTaa < 2:376775/141‑1 (MQ=255) aaCTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTaaaa < 2:929085/141‑1 (MQ=255) aaCTAAATAAGACTAAAAAAACAACAGCTAATATAGAAAATACCCAAAAGAAATTTTTAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTaaaa < 1:194937/141‑1 (MQ=255) aaaTACTACTAACAAAACAACAGCGAATATAGAAAAAACCCAAAAGATATTTTGAGA‑TTTTTTTTGTTTCATCCTAGCCCCACCTTTTGCTTATCTTTTCTCAATTATAGGATATAATAAAAAAGTTGTAATTAAAAGTgg < 2:530431/141‑1 (MQ=255) aataCTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGATTTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGAt < 2:506953/141‑1 (MQ=255) aaaaaaaGAACAGCGAATATAGAAAATATCCAAAAGATATTTTTAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACtt < 1:874547/138‑1 (MQ=255) aaGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACtt > 1:618466/1‑141 (MQ=255) aGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTa < 1:240486/141‑1 (MQ=255) aacaacAGCTAATATAGAAAATACCCAAAAGATATTTTTAGA‑TTATTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAGAGTGGGATTTTACTTaagaa > 1:533906/1‑141 (MQ=255) acaacaGCTAATATAGAAAATACCCAAAAGATATTTTTAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTaagaaa < 1:539457/141‑1 (MQ=255) aacaGCTAATATAGAAAATACCCAAAAGATATTTTTAGA‑TTTTTTTTGTTTCATACTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTACTTAAAAGTGGGATTTGACTTtagaaaga > 2:267138/1‑141 (MQ=255) aaGCTAATATAGAAAATACCCAAAAGATATTTTTAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAaagaag < 1:917082/140‑1 (MQ=255) aaTATAGAAAATACCCAAAAGATATTTTTAGA‑TTTTTTTTTTTTTTTTTTATCCCCACCTTTTGCTTATCTTTTCTCAATTTTATGATATAATAAAAAAGTTGTAATTAACAGTTGGATTCTACTTAAGAAAGAAAGAAAc > 2:774399/1‑141 (MQ=255) tatataAAATACCCAAAAGATATTGTTAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTTACTTATATGGTTTGATAAAAGAGGTGTAATTTAAAGTGGGATTTTACTTAAGGATGGAGGAAACTa < 2:199125/141‑1 (MQ=255) tataGAAAATACCCAAAAGATATTTTTAGA‑‑TTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTAt < 1:592481/141‑1 (MQ=255) taGAAAATACCCAAAAGATATTTTTAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTAtt < 1:606180/141‑1 (MQ=255) aaaaTACCCAAAAGATATTTTTAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGACATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTATTtat < 2:403067/141‑1 (MQ=255) aaaTACCCAAAAGATATTTTTAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTATTtata > 1:874051/1‑141 (MQ=255) tACCCAAAAGATATTTTTAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTATTTATATGa < 2:168775/141‑1 (MQ=255) aaCCCAAAAGATATTTATAGA‑TTTTTCTTGTTTCATCCTATCCCCACCTTTTGCTTAACATATCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTATTTATATGa < 1:925831/140‑1 (MQ=255) aaaaGATATTTTTAGA‑TTTTTTTTTTTTCTTCTTATCCCCACATTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGGCGGATTTTACTTAAGAAAGAAGAAAACTATTTATATGACtaat > 2:813010/1‑141 (MQ=255) aaaaGATATTTTTAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTATTTATATGaataat > 1:553047/1‑141 (MQ=255) aaaaGATATTTTTAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTATTTATATGACtaat < 2:191342/141‑1 (MQ=255) aaaaGATATTTTTAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTATTTATATGACtaat < 2:786897/141‑1 (MQ=255) aaaaGATATTTTTAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTATTAATATGACtaat > 1:564748/1‑141 (MQ=255) gATATTTTTAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTATTTATATGACTAATaaaa > 1:335741/1‑141 (MQ=255) atTTTTAGA‑TTTTTTTTTTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTATAAGTGGGATTTTACTTAATAAAGAAGGATATTATTTATATGTCttaaataagct > 2:174423/1‑131 (MQ=255) atTTTTAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTATTTATATGCCTAATaaaaaaa > 2:495514/1‑138 (MQ=255) atTTTTAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAACTTGTAATTAAAAGTGGGATTTTACTTAACAAAGAAGGAAACTATTTATATAACTAATAAAagaa > 2:125890/1‑140 (MQ=255) tttttAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTATTTATATGACTAATAAAAgaga > 1:864091/1‑141 (MQ=255) tttttAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTATTTATATGACTAATAAAAgaga > 1:161931/1‑141 (MQ=255) tttttAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTATTTATATGACTAATAAAAgaga > 1:847188/1‑141 (MQ=255) tttttAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAAATATTTATATGACTAGTAAAAgaga > 1:161795/1‑141 (MQ=255) tttttAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGAGGGATTTTACTTAAGAAAGAAGGAAACTATTTATATGACTAATAAAAgaga > 1:491048/1‑141 (MQ=255) tttAGA‑TTTTTTTTTTTTTTTTCTATCACCAACTTTTGCTTGTCTTTTCTCAATTATATGATATAATAAAAAAGTTGCAATTCAAAGTCGGATTTTACTTAAGAAAGATTGAAACTTTTTATATGACTAATAAAAGagaaa > 2:493494/1‑140 (MQ=255) gTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATACAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTATTTATATAGCTACTAAAAAAGAAGATATCCGCacaata > 2:122034/1‑136 (MQ=255) tttCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTATTTATATGACTAATAAAAGAGAAGATGTTCGCAATATAg < 2:638276/141‑1 (MQ=255) cATCCTATCCCCACCTTTTGCTTATCTTTTCTAAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTATTTATATGACTAATAAAAGAGAAGATGTCCGCAATATAGCaa < 2:87142/141‑1 (MQ=255) | CCAAAGGTTGCAACGAGTAATATAAATGTTAAAATCATCATTGGTACATTTGATGCACCAACAGCAAAACTAAATAATACTAAGAAAACAACAGCTAATATAGAAAATACCCAAAAGATATTTTTAGA‑TTTTTTTTGTTTCATCCTATCCCCACCTTTTGCTTATCTTTTCTCAATTATATGATATAATAAAAAAGTTGTAATTAAAAGTGGGATTTTACTTAAGAAAGAAGGAAACTATTTATATGACTAATAAAAGAGAAGATGTCCGCAATATAGCAA > CP000730/1117490‑1117770 |
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: ‑ |