Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | USA300TCH1516_ALE | 41,171 | T→C | 6.7% | intergenic (‑38/‑62) | pbp ← / → mecR1 | Beta‑lactam‑inducible penicillin‑binding protein/Methicillin resistance mecR1 protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | USA300TCH1516_ALE | 41,171 | 0 | T | C | 6.7% | 478.1 / 16.6 | 148 | intergenic (‑38/‑62) | pbp/mecR1 | Beta‑lactam‑inducible penicillin‑binding protein/Methicillin resistance mecR1 protein |
Reads supporting (aligned to +/- strand): ref base T (69/69); new base C (7/3); total (76/72) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.28e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01 |
Only 100 of 156 total aligned reads displayed.
CCCGACAACTACAACTATTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGTTAAGTA > USA300TCH1516_ALE/41083‑41260 | cccGACAACTACAACTATTAAAATAAGTGTATCAATTTTTATCTCTTTCATCATTATTCTCCTTATATAAGACTACATTTGTAATATACTa < 1:3317151/91‑1 (MQ=255) cccGACAACTACAACTATTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTa < 1:979200/91‑1 (MQ=255) ccGACAACTACAACTATTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAATAtactac < 1:2346249/91‑1 (MQ=255) cGACAACTACAACTATTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACa > 1:2878784/1‑91 (MQ=255) cGACAACTACAACTATTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACa > 1:1356896/1‑91 (MQ=255) cGACAACTACAACTATTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACa > 1:1356811/1‑91 (MQ=255) gACAACTACAACTATTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAATATACTACaa < 1:3906517/91‑1 (MQ=255) acaactacaactATTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACaaa > 1:3891640/1‑91 (MQ=255) caactacaactaTTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAAt > 1:2433281/1‑91 (MQ=255) caactacaactaTTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAAt > 1:634998/1‑91 (MQ=255) aactacaactaTTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATg < 1:3719715/91‑1 (MQ=255) aactacaactaTTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAATATACTACAAATg < 1:1624225/91‑1 (MQ=255) actacaactaTTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATgt < 1:3125769/91‑1 (MQ=255) actacaactaTTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATgt < 1:3129705/91‑1 (MQ=255) actacaactaTTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGCATATTACAAATgt < 1:3124710/91‑1 (MQ=255) ctaTTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTAttt < 1:810549/91‑1 (MQ=255) ttAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATg > 1:12880/1‑91 (MQ=255) ttAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATg > 1:2628741/1‑91 (MQ=255) aaaTAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCa > 1:93012/1‑91 (MQ=255) aaTAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCaa < 1:2671176/91‑1 (MQ=255) aaTAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCaa > 1:1053564/1‑91 (MQ=255) aaGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAaat > 1:649159/1‑91 (MQ=255) gTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAaataa > 1:1842909/1‑91 (MQ=255) gTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAaataa > 1:2236317/1‑91 (MQ=255) gTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAaataa > 1:1842807/1‑91 (MQ=255) ggAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATg > 1:1211714/1‑91 (MQ=255) ggAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATg > 1:2285655/1‑91 (MQ=255) ggAACAATNTTTNTCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATg < 1:422027/91‑1 (MQ=255) aaCAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCATAATAATGtt < 1:2091198/91‑1 (MQ=255) aaCAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGtt < 1:1088068/91‑1 (MQ=255) aaCAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGtt < 1:1157088/91‑1 (MQ=255) aaCAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGtt < 1:869004/91‑1 (MQ=255) aaCAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGtt < 1:1541637/91‑1 (MQ=255) aaCAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGtt < 1:3543157/91‑1 (MQ=255) cAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTtat < 1:1514244/91‑1 (MQ=255) cAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTtat < 1:304695/91‑1 (MQ=255) cAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTtat < 1:1784841/91‑1 (MQ=255) aaTTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTtata < 1:3139337/91‑1 (MQ=255) aaTTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTtata > 1:3800018/1‑91 (MQ=255) tttATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAAttt > 1:1026460/1‑91 (MQ=255) ttATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAAtttt > 1:1878572/1‑91 (MQ=255) tATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAAttttt < 1:1094713/91‑1 (MQ=255) tATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAAttttt > 1:2860091/1‑91 (MQ=255) tctcTTTCATCAATATACTCCTTAAATAAGACTACATTTGTAGTATATTACAAATGTAGGATCCATGTCAAAATAATGTTATAATTTTtgt < 1:2976150/91‑1 (MQ=255) tCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTtgt < 1:2806776/91‑1 (MQ=255) tCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTtgt > 1:2750313/1‑91 (MQ=255) tCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTtgt < 1:986049/91‑1 (MQ=255) tCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTtgt > 1:2709413/1‑91 (MQ=255) tttCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGata > 1:2646537/1‑91 (MQ=255) tttCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGata < 1:2021960/91‑1 (MQ=255) ttCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGatat < 1:117452/91‑1 (MQ=255) ttCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGatat < 1:770793/91‑1 (MQ=255) ttCATCAATATACTCCTCATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGatat < 1:3556564/91‑1 (MQ=255) tcaATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAg < 1:2694215/91‑1 (MQ=255) tcaATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAg < 1:595764/91‑1 (MQ=255) aaTATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGgt > 1:1414036/1‑91 (MQ=255) aaTATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGgt > 1:85107/1‑91 (MQ=255) aaTATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGgt < 1:2634061/91‑1 (MQ=255) atatACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGgtg > 1:3781354/1‑91 (MQ=255) atatACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGgtg > 1:2034120/1‑91 (MQ=255) atatACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGgtg > 1:3303407/1‑91 (MQ=255) atatACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGgtg < 1:3238268/91‑1 (MQ=255) tCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAg > 1:2311903/1‑91 (MQ=255) tCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAg < 1:1498097/91‑1 (MQ=255) tatataAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGgtgt > 1:907967/1‑91 (MQ=255) tatataAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGgtgt > 1:2158678/1‑91 (MQ=255) atataAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGtt > 1:172600/1‑91 (MQ=255) atataAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGtt < 1:1591117/91‑1 (MQ=255) tataAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTa < 1:562934/91‑1 (MQ=255) tataAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTa > 1:2406334/1‑91 (MQ=255) tataAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTa > 1:2461675/1‑91 (MQ=255) tataAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTa < 1:1449401/91‑1 (MQ=255) tataAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTa > 1:2405431/1‑91 (MQ=255) ataAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTat < 1:2633670/91‑1 (MQ=255) taAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTatc > 1:2088175/1‑91 (MQ=255) taAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTatc < 1:4082468/91‑1 (MQ=255) gACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATGATTTTTGTGATATGGAGGTGTAGAAGGTGTTatcatc < 1:3506623/91‑1 (MQ=255) gACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTatcatc < 1:2054786/91‑1 (MQ=255) gACTACATTTGTAGTATATTACAAATGTAGTATGTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTatcatc < 1:3681032/91‑1 (MQ=255) tACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCttt < 1:3291291/91‑1 (MQ=255) tACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCttt < 1:3291175/91‑1 (MQ=255) aCATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCtttt > 1:974680/1‑91 (MQ=255) aCATTTGTAATATACTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCtttt > 1:2776684/1‑91 (MQ=255) aCATTTGTAATATACTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCtttt > 1:406530/1‑91 (MQ=255) aTTTGTAGTATATTACAGATGTAGTATTTATGGCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCtttttt < 1:3546833/91‑1 (MQ=255) aTTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCtttttt > 1:3567834/1‑91 (MQ=255) aTTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCtttttt < 1:554697/91‑1 (MQ=255) tttGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTa > 1:2563806/1‑91 (MQ=255) ttGTAATATACTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTaa > 1:4101680/1‑91 (MQ=255) tGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAAt > 1:1233240/1‑91 (MQ=255) gtagtaTATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATg < 1:4008565/91‑1 (MQ=255) tagtaTATTACGAATGTGGGATTTATGTCGAAATAATGTTAGAATTTTTGTGGTATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGt < 1:3333042/91‑1 (MQ=255) taatataCTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGt > 1:2034848/4‑91 (MQ=255) agtaTATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGtt < 1:345111/91‑1 (MQ=255) aatataCTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGtt > 1:2430913/3‑91 (MQ=255) gtaTATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGTTa > 1:2174465/1‑91 (MQ=255) atatTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGTTAAg < 1:3957117/91‑1 (MQ=255) atatTACAAATGGAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGTTAAg > 1:3200796/1‑91 (MQ=255) tatTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGTTAAGt < 1:237308/91‑1 (MQ=255) atTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGTTAAGta < 1:3114916/91‑1 (MQ=255) | CCCGACAACTACAACTATTAAAATAAGTGGAACAATTTTTATCTTTTTCATCAATATACTCCTTATATAAGACTACATTTGTAGTATATTACAAATGTAGTATTTATGTCAAAATAATGTTATAATTTTTGTGATATGGAGGTGTAGAAGGTGTTATCATCTTTTTTAATGTTAAGTA > USA300TCH1516_ALE/41083‑41260 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |