| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,216,702 | T→C | 47.5% | Q677Q (CAA→CAG) | prc ← | carboxy‑terminal protease for penicillin‑binding protein 3 |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,216,702 | 0 | T | C | 47.5% | ‑4.7 / 55.2 | 40 | Q677Q (CAA→CAG) | prc | carboxy‑terminal protease for penicillin‑binding protein 3 |
| Reads supporting (aligned to +/- strand): ref base T (7/14); new base C (6/13); total (13/27) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.53e-02 | |||||||||||
CGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGAT > minE/1216641‑1216769 | cGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGCTGTTCcgcg > 1:940825/1‑71 (MQ=255) gCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGCTGTTCCGCgg < 1:3590293/71‑1 (MQ=255) ttAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGt < 1:1536877/69‑1 (MQ=255) ttAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGt < 1:772533/69‑1 (MQ=255) ttAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTc < 1:2238637/70‑1 (MQ=255) ttAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGCTGTTCCGCGGGTc < 1:1796495/70‑1 (MQ=255) tAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTg > 1:1481072/1‑71 (MQ=255) aaaaaTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGc < 1:776386/71‑1 (MQ=255) aaaTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGCTGTTCCGCGGGTCTGGCtt > 1:651032/1‑71 (MQ=255) aTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCtttt < 1:43951/71‑1 (MQ=255) aTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCtttt < 1:2059974/71‑1 (MQ=255) aTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGCTGTTCCGCGGGTCTGGCtttt < 1:561759/71‑1 (MQ=255) aTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGCTGTTCCGCGGGTCCGGCtttt < 1:708869/71‑1 (MQ=255) tCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGCTGTTCCGCGGGTCTGGCttttt < 1:84240/71‑1 (MQ=255) tCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGCTGTTCCGCGGGTCTGGCtttt > 1:2922181/1‑70 (MQ=255) aGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGCTGTTCCGCGGGTCTGGCTTTTTCa < 1:558096/71‑1 (MQ=255) acaATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCt < 1:2962992/71‑1 (MQ=255) acaATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCt > 1:1227180/1‑71 (MQ=255) aaTTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTc > 1:2442573/1‑71 (MQ=255) aTTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTc < 1:3251056/53‑1 (MQ=255) tttCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGc > 1:1789675/1‑71 (MQ=255) ttCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGCTGTTCCGCGGGTCTGGCTTTTTCa < 1:1622715/62‑1 (MQ=255) aattGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCcg > 1:3532698/3‑44 (MQ=255) ttGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTTTGGCTTTTTCAAGCTTCGCCa < 1:1992505/69‑1 (MQ=255) ttccTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATc < 1:1477415/68‑1 (MQ=255) ccTGATTGATATTACTTGACGGGAGCGGGCTGTTCCGCGGGTc < 1:191076/43‑1 (MQ=255) ccTGATTGATATTACTTGACGGGAGCGGGCTGTTCCGCGGGTc < 1:2703118/43‑1 (MQ=255) attgatATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTg < 1:1005589/69‑1 (MQ=255) ttgatATTACTTGACGGGAGCGGGCTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGa > 1:1038577/1‑71 (MQ=255) atTACTTGACGGGAGCGGGCTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATAtt < 1:236918/70‑1 (MQ=255) gACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGtct > 1:1087673/1‑71 (MQ=255) gACGGGAGCGGGCTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGtc < 1:2477758/70‑1 (MQ=255) aCGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATc < 1:3252647/50‑1 (MQ=255) aCGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGtctc < 1:601337/71‑1 (MQ=255) cGGGAGCGGGCTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCa < 1:2880456/71‑1 (MQ=255) gggAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGa > 1:2639630/1‑50 (MQ=255) gAGCGGGCTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATc > 1:2484004/1‑70 (MQ=255) gAGCGGGCTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATc > 1:679612/1‑70 (MQ=255) aGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATcc < 1:2417351/70‑1 (MQ=255) gggCTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGAt < 1:2483511/71‑1 (MQ=255) | CGCTGTTAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGAT > minE/1216641‑1216769 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |