Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 54,435 | T→C | K90E (AAA→GAA) | surA ← | peptidyl‑prolyl cis‑trans isomerase SurA |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 54,435 | 0 | T | C | 97.1% | 122.0 / ‑6.5 | 35 | K90E (AAA→GAA) | surA | peptidyl‑prolyl cis‑trans isomerase SurA |
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (18/16); minor base G (0/1); total (18/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.86e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GCATCTGATCCAGCGTCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGT > NC_000913/54360‑54510 | gcATCTGATCCAGCGTCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTCCACTc < 1:12585/81‑1 (MQ=255) aTCTGATCCAGCGTCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTCCACTccc < 2:12586/81‑1 (MQ=255) gATCCAGCGTCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTCCACTCCCAttt > 2:12638/1‑81 (MQ=255) ccAGCGTCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTCCACTCCCATTTTCt > 2:12625/1‑81 (MQ=255) cAGCGTCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTCCACTCCCATTTTCt < 1:12587/80‑1 (MQ=255) aGCGTCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTCCACTCCCATTTTCTGc < 1:12588/81‑1 (MQ=255) aGCGTCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTCCACTCCCATTTTCTGc > 1:12636/1‑81 (MQ=255) tCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTCCACTCCCATTTTCTGCCCCa < 1:12589/81‑1 (MQ=255) tCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTCCACTCCCATTTTCTGCCCCa > 2:12643/1‑81 (MQ=255) tCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTCCACTCCCATTTTCTGCCCCa > 1:12632/1‑81 (MQ=255) cATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTCCACTCCCATTTTCTGCCCCAt > 2:12647/1‑81 (MQ=255) gttgttCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTCCACTCCCATTTTCTGCCCCATCTg > 1:12637/1‑81 (MQ=255) ttgttCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTCCACTCCCATTTTCTGCCCCATCTGc < 2:12590/81‑1 (MQ=255) tttCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTCCACTCCCATTTTCTGCCCCATCTGCAGgatgat > 1:12652/1‑81 (MQ=255) tttCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTCCACTCCCATTTTCTGCCCCATCTGCAGgatgat > 2:12631/1‑81 (MQ=255) cgATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTCCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGAt > 1:12622/1‑81 (MQ=255) cgATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTCCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGAt < 2:12591/81‑1 (MQ=255) gATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTCCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATc < 1:12592/81‑1 (MQ=255) gATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTCCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATc < 1:12593/81‑1 (MQ=255) aTGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTCCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATcc < 1:12594/81‑1 (MQ=255) gCAATCGCCTGATCCAGCTGCTCATCGGAGATTTCCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAt < 2:12595/80‑1 (MQ=255) gCAATCGCCTGATCCAGCTGCTCATCGGAGATTTCCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGAt < 2:12596/80‑1 (MQ=255) tGATCCAGCTGCTCATCGGAGATTTCCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTc > 1:12660/1‑80 (MQ=255) tGATCCAGCTGCTCATCGGAGATTTCCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTc > 1:12639/1‑80 (MQ=255) gATCCAGCTGCTCATCGGAGATTTCCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCa > 1:12645/1‑81 (MQ=255) gATCCAGCTGCTCATCGGAGATTTCCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCa > 1:12640/1‑81 (MQ=255) aTCCAGCTGCTCATCGGAGATTTCCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCAt > 1:12656/1‑81 (MQ=255) aGCTGCTCATCGGAGATTTCCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGAtt > 2:12651/1‑81 (MQ=255) ctCATCGGAGATTTGCGCTCCCATTTTCTGCGCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTtggt < 2:12597/80‑1 (MQ=255) agcGGAGATTTCCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGgcg < 1:12599/79‑1 (MQ=255) aTCGGAGATTTCCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGgcg > 2:12653/1‑81 (MQ=255) aTCGGAGATTTCCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGgc < 1:12598/80‑1 (MQ=255) cGGAGATTTCCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGgcg > 1:12644/1‑79 (MQ=255) ggAGATTTCCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCa < 2:12600/80‑1 (MQ=255) gaTTTCCACTCCCATTTTCTGCCCCATCTGCGGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAgcgt < 2:12601/81‑1 (MQ=255) | GCATCTGATCCAGCGTCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGT > NC_000913/54360‑54510 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |