Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP000730 | 39,378 | C→A | 100% | D586Y (GAT→TAT) | mecA ← | penicillin binding protein 2 prime |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP000730 | 39,378 | 0 | C | A | 84.6% | 25.0 / 0.0 | 13 | D586Y (GAT→TAT) | mecA | penicillin binding protein 2 prime |
Reads supporting (aligned to +/- strand): ref base C (0/0); major base A (5/6); minor base T (1/1); total (6/7) | |||||||||||
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.39e-01 | |||||||||||
Rejected as polymorphism: E-value score below prediction cutoff. | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
GTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATATTGCC > CP000730/39256‑39510 | gTTTGGATTGTCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCACTTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTtttatt > 2:10421/1‑141 (MQ=255) atGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCAt > 1:156010/1‑141 (MQ=255) ttCTACTTGTTTCATTTTGAGTTCTGCAAAGCCGGGTTTGCCACTTACGTGTGCATAAGAGCTATAAAAAGATGCTTTATGTGTTTTATTTGCGACGTGTTGCATACATTCAGTTACTAGAGTGATGTTGTCTTTGGAgat > 2:152024/1‑141 (MQ=255) tGTTTCTTTTCGAGTGCTGCCGTGCCGGATGTGCCACTTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATAttt < 2:57052/141‑1 (MQ=255) tCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTttct < 2:44552/141‑1 (MQ=255) tCTGCAGTACCGTATTTGCCAATTAAGTTTGCATAAGAACTATAAATTTATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTg < 1:244297/141‑1 (MQ=255) cAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCAAAAGTTCGtctt > 2:108915/1‑141 (MQ=255) gTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATTTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTc < 1:115421/141‑1 (MQ=255) cGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCgtgt > 1:49547/1‑141 (MQ=255) tCTATAAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATAtt < 2:160328/141‑1 (MQ=255) tataAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATATTGc < 2:5211/141‑1 (MQ=255) ataAATATATTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATATTGcc > 2:106392/1‑141 (MQ=255) aaaaaTATATTCCTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATAATTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAAGTGCGTTAATATTGcc < 1:95523/139‑1 (MQ=255) | GTTTGGATTATCTTTATCATATGATATAAACCACCCAATTTGTCTGCCAGTTTCTCCTTGTTTCATTTTGAGTTCTGCAGTACCGGATTTGCCAATTAAGTTTGCATAAGATCTATAAATATCTTCTTTATGTGTTTTATTTACGACTTGTTGCATACCATCAGTTAATAGATTGATATTTTCTTTGGAAATAATATTTTTCTTCCAAACTTTGTTTTTCGTGTCTTTTAATAAGTGAGGTGCGTTAATATTGCC > CP000730/39256‑39510 |
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: ‑ |