| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,330,586 | A→G | 18.2% | Y334H (TAC→CAC) | cpsG ← | phosphomannomutase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,330,586 | 0 | A | G | 18.2% | 62.4 / 11.6 | 33 | Y334H (TAC→CAC) | cpsG | phosphomannomutase |
| Reads supporting (aligned to +/- strand): ref base A (14/13); new base G (4/2); total (18/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.65e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.36e-01 | |||||||||||
GCACACCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGT > minE/1330518‑1330654 | gCACACCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAAt < 1:1080165/71‑1 (MQ=255) acacCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGTAAGTAAtgg > 1:1676453/1‑71 (MQ=255) acacCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTa > 1:1949175/1‑67 (MQ=255) acacCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAAtgg > 1:90483/1‑71 (MQ=255) acacCAGTTCGGCGACCAGAAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAAtgg > 1:1057068/1‑71 (MQ=255) cacCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAAtggt < 1:1706253/71‑1 (MQ=255) gCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCAt > 1:1862939/1‑71 (MQ=255) aCCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTc < 1:1300078/71‑1 (MQ=255) cagcagCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGc > 1:1606354/1‑71 (MQ=255) cagcagCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGc > 1:1239670/1‑71 (MQ=255) cagcCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGccac < 1:1482400/71‑1 (MQ=255) cATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGt > 1:1026122/1‑71 (MQ=255) aTGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAAtgg > 1:2194528/1‑39 (MQ=255) aTGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAAtgg > 1:1153025/1‑39 (MQ=255) tGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCt > 1:2123886/1‑71 (MQ=255) ccGCTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTc < 1:825262/45‑1 (MQ=255) gCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCAt < 1:723361/50‑1 (MQ=255) gCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCa < 1:1288507/71‑1 (MQ=255) aGTAAGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATa > 1:1662217/1‑71 (MQ=255) tAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATa < 1:334117/47‑1 (MQ=255) aGCGAAATCACGGAAGTGATGGTGGTCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGtt > 1:677160/1‑71 (MQ=255) aGCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGtt > 1:724798/1‑71 (MQ=255) aGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGa > 1:111136/1‑47 (MQ=255) aGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGtt > 1:374738/1‑71 (MQ=255) aGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGt > 1:1584070/1‑70 (MQ=255) aGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGt > 1:1452252/1‑70 (MQ=255) gCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATAGATGGc < 1:201654/50‑1 (MQ=255) gCGAAATCACGGAAGTGATGGTGGGCGCTCATTTCGCCACCATAGATGGc < 1:1134821/50‑1 (MQ=255) gAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCtt < 1:234199/71‑1 (MQ=255) gAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCtt < 1:918263/71‑1 (MQ=255) aCGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTtaata < 1:1557702/70‑1 (MQ=255) gTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCttt < 1:588290/60‑1 (MQ=255) gTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCttt < 1:199592/60‑1 (MQ=255) gTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCgt > 1:2004594/1‑71 (MQ=255) aaTGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATa > 1:526729/1‑50 (MQ=255) aaTGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATa > 1:1074365/1‑50 (MQ=255) aaTGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATa > 1:969095/1‑50 (MQ=255) | GCACACCAGTTCGGCGACCAGCAGCCACGGGATCATGCCGCTGTCGCAGTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGT > minE/1330518‑1330654 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |