| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 617,556 | T→A | 100% | D198E (GAT→GAA) | ycaL → | predicted peptidase with chaperone function |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 617,556 | 0 | T | A | 86.4% | 57.5 / 3.0 | 22 | D198E (GAT→GAA) | ycaL | predicted peptidase with chaperone function |
| Reads supporting (aligned to +/- strand): ref base T (0/3); new base A (19/0); total (19/3) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.49e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01 | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
CTATCAATGCTAAATACTCCCGTGATAAAGAGTCCGAAGCAGATGATTTCTCCTTTGATCTGTTGAAGAAACGTGGC > minE/617510‑617586 | cTATCAATGCTAAATACTCCCGTGATAAAGAGTCCGAAGCAGATGATTTCTCCTTTGATCTGTTgaa < 1:993110/67‑1 (MQ=255)cTATCAATGCTAAATACTCCCGTGATAAAGAGTCCGAAGCAGATGATTTCTCCTTTGATCTGTTgaa < 1:821171/67‑1 (MQ=255)cTATCAATGCTAAATACTCCCGTGATAAAGAGTCCGAAGCAGATGATTTCTCCTTTGATCTGTTgaa < 1:1368766/67‑1 (MQ=255) aaaTACTCCCGTGATAAAGAGTCCGAAGCAGATGAATTCTCCTTTGATCTGTTGAAGAAACGTGGc > 1:1938560/1‑66 (MQ=255) aaaTACTCCCGTGATAAAGAGTCCGAAGCAGATGAATTCTCCTTTGATCTGTTGAAGAAACGTGGc > 1:823712/1‑66 (MQ=255) aaaTACTCCCGTGATAAAGAGTCCGAAGCAGATGAATTCTCCTTTGATCTGTTGAAGAAACGTGGc > 1:285726/1‑66 (MQ=255) aaaTACTCCCGTGATAAAGAGTCCGAAGCAGATGAATTCTCCTTTGATCTGTTGAAGAAACGTGGc > 1:2091175/1‑66 (MQ=255) aaaTACTCCCGTGATAAAGAGTCCGAAGCAGATGAATTCTCCTTTGATCTGTTGAAGAAACGTGGc > 1:2074938/1‑66 (MQ=255) aaaTACTCCCGTGATAAAGAGTCCGAAGCAGATGAATTCTCCTTTGATCTGTTGAAGAAACGTGGc > 1:1990742/1‑66 (MQ=255) aaaTACTCCCGTGATAAAGAGTCCGAAGCAGATGAATTCTCCTTTGATCTGTTGAAGAAACGTGGc > 1:1967839/1‑66 (MQ=255) aaaTACTCCCGTGATAAAGAGTCCGAAGCAGATGAATTCTCCTTTGATCTGTTGAAGAAACGTGGc > 1:1102114/1‑66 (MQ=255) aaaTACTCCCGTGATAAAGAGTCCGAAGCAGATGAATTCTCCTTTGATCTGTTGAAGAAACGTGGc > 1:1916056/1‑66 (MQ=255) aaaTACTCCCGTGATAAAGAGTCCGAAGCAGATGAATTCTCCTTTGATCTGTTGAAGAAACGTGGc > 1:177870/1‑66 (MQ=255) aaaTACTCCCGTGATAAAGAGTCCGAAGCAGATGAATTCTCCTTTGATCTGTTGAAGAAACGTGGc > 1:1689714/1‑66 (MQ=255) aaaTACTCCCGTGATAAAGAGTCCGAAGCAGATGAATTCTCCTTTGATCTGTTGAAGAAACGTGGc > 1:1609770/1‑66 (MQ=255) aaaTACTCCCGTGATAAAGAGTCCGAAGCAGATGAATTCTCCTTTGATCTGTTGAAGAAACGTGGc > 1:1484677/1‑66 (MQ=255) aaaTACTCCCGTGATAAAGAGTCCGAAGCAGATGAATTCTCCTTTGATCTGTTGAAGAAACGTGGc > 1:1422592/1‑66 (MQ=255) aaaTACTCCCGTGATAAAGAGTCCGAAGCAGATGAATTCTCCTTTGATCTGTTGAAGAAACGTGGc > 1:1214750/1‑66 (MQ=255) aaaTACTCCCGTGATAAAGAGTCCGAAGCAGATGAATTCTCCTTTGATCTGTTGAAGAAACGTGGc > 1:1170407/1‑66 (MQ=255) aaaTACTCCCGTGATAAAGAGTCCGAAGCAGATGAATTCTCCTTTGATCTGTTGAAGAAACGTGGc > 1:1108946/1‑66 (MQ=255) aaaTACTCCCGTGATAAAGAGTCCGAAGCAGATGAATTCTCCTTTGATCGGTTGAAGAAACGTGGc > 1:796116/1‑66 (MQ=255) aaaTACTCCCGTGATAAAGAGTCCGAAGCAGATGAATTCACCTTTGATCTGTTGAAGAAACGTGGc > 1:652938/1‑66 (MQ=255) | CTATCAATGCTAAATACTCCCGTGATAAAGAGTCCGAAGCAGATGATTTCTCCTTTGATCTGTTGAAGAAACGTGGC > minE/617510‑617586 |
| 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: ‑ |