Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP000730 | 2,455,698 | 2 bp→AT | 100% | coding (1216‑1217/1260 nt) | USA300HOU_2315 → | hypothetical membrane protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP000730 | 2,455,698 | 0 | C | A | 93.3% | 82.1 / 1.5 | 30 | coding (1216/1260 nt) | USA300HOU_2315 | hypothetical membrane protein |
Reads supporting (aligned to +/- strand): ref base C (0/0); major base A (13/14); minor base . (1/1); total (15/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.61e-01 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch. | |||||||||||
* | CP000730 | 2,455,699 | 0 | A | T | 97.8% | 77.0 / ‑6.4 | 30 | H406L (CAT→CTT) | USA300HOU_2315 | hypothetical membrane protein |
Reads supporting (aligned to +/- strand): ref base A (2/0); major base T (13/13); minor base G (0/2); total (15/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.84e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Frequency below/above cutoff threshold. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
GTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTT > CP000730/2455569‑2455825 || gTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCCaaaaaa > 2:265718/1‑141 (MQ=255) cAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTgaag < 1:120061/141‑1 (MQ=255) aGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCGTCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTgaagaa < 1:134458/141‑1 (MQ=255) caccacAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACa > 2:138643/1‑141 (MQ=255) tGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGAt < 2:49282/141‑1 (MQ=255) ggACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAAAAATTAAGAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGAt < 1:27568/140‑1 (MQ=255) aCGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAAATATCGAAAAAACTGAAGAACAATCTTCAGATaa > 2:195076/1‑141 (MQ=255) ctaTTAAACAACAGATTCAAATGGAGACGAAAGAGAGCAATCATTAAAAGAACCTGAGACTTACAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGGAGAACAATCTTCAGATaaa < 2:119511/139‑1 (MQ=255) cAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAAGTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAg < 1:247716/141‑1 (MQ=255) aaaTGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAgc < 1:119010/141‑1 (MQ=255) aGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAACAAATA‑ATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTAGc > 2:320753/1‑141 (MQ=255) gagCAATCATTAATAGAACCTGCGACTTATAAAGAAGACAGACGTGCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAgcatcatc < 2:270533/141‑8 (MQ=255) cAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTaa > 1:191168/1‑141 (MQ=255) tAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAgatg > 2:128580/1‑141 (MQ=255) aaaaGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAgatga < 1:150015/141‑1 (MQ=255) aGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGc > 1:354283/1‑141 (MQ=255) aGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATATTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTAACCCTTAAAAAGATGATGc > 2:368/1‑141 (MQ=255) gAACCTGGGACGTATAGAGAAGACAGACGTTCATCAGTTGTAAGTGATGCA‑AAAAAAGTATCGAAAAGACTGAAGAACAATGTTCAGATAAAAACAAATAATAGCAGAGCGCTTGCTAGTACCTCTTGAAAAGATGATGCt < 2:352733/141‑1 (MQ=255) aCCTGCGACTTATAAAGAAGACAGACGTCCATTAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAg < 1:8217/141‑1 (MQ=255) ccTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAAAAAAA‑TTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGc > 1:238741/1‑141 (MQ=255) ccTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTCCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGc > 1:125752/1‑141 (MQ=255) cTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGACGATGCTAGCAAGCGCt < 1:212172/141‑1 (MQ=255) ataAAGAAGACAGACGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCtttt > 2:67898/1‑141 (MQ=255) acagacGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCACTTAAAAAGATGATGCTAGCAAGCGCTTTTctatactac > 1:8574/1‑140 (MQ=255) agacGTTCATCAGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACtatat > 2:118728/1‑141 (MQ=255) acGTTCATCAGTTGTAATTGATGCA‑AACAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTatatat > 1:263157/1‑141 (MQ=255) tcaGTTGTAATTGATGCAAAAAAA‑TTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTAttttt < 1:138643/141‑1 (MQ=255) caGTTGTAATTGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTc > 1:171336/1‑141 (MQ=255) tGTAATTGATG‑AAAAAAAATTATCGAAAAAACTGAATAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTTGCAAGCGCTTTTCTATACTATATATTATTTTTCTaaa < 2:84742/141‑1 (MQ=255) ttGATGCA‑AAAAAATTATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAAtttt < 2:356822/141‑1 (MQ=255) || GTTAAACAAGCACCACAAAGTGACGATTTAACAAACGATTCAAATGAAGATGAAATAGAGCAATCATTAAAAGAACCTGCGACTTATAAAGAAGACAGACGTTCATCAGTTGTAATTGATGCAGAAAAACATATCGAAAAAACTGAAGAACAATCTTCAGATAAAAACAAATAATAACAAAGCGCTTGCTAGTACCTCTTAAAAAGATGATGCTAGCAAGCGCTTTTCTATACTATATATTATTTTTCTAAAATTTT > CP000730/2455569‑2455825 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |