| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 158,183 | A→G | 28.2% | N41S (AAC→AGC) | dgt → | deoxyguanosine triphosphate triphosphohydrolase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 158,183 | 0 | A | G | 28.2% | 56.6 / 26.5 | 39 | N41S (AAC→AGC) | dgt | deoxyguanosine triphosphate triphosphohydrolase |
| Reads supporting (aligned to +/- strand): ref base A (11/17); new base G (7/4); total (18/21) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.85e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCG > minE/158118‑158252 | gCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAGCtctc < 1:17707/71‑1 (MQ=255) cAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTcc > 1:1193855/1‑71 (MQ=255) cGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAGCTCTCCGGCaa > 1:1321608/1‑71 (MQ=255) aaaaCCGAACATGAGATCCTTCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtcg > 1:1050887/1‑71 (MQ=255) aaaaCCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtcg > 1:1912732/1‑71 (MQ=255) aaaaCCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtc < 1:1113937/70‑1 (MQ=255) aaaaCCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtc < 1:1528069/70‑1 (MQ=255) aaaCGGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtcg < 1:621771/70‑1 (MQ=255) aaaCCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtcgt < 1:599572/71‑1 (MQ=255) ccGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGc < 1:1485988/61‑1 (MQ=255) ccGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTcgtcgt > 1:1272675/1‑71 (MQ=255) cGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAGCTCTCCGGCAATTCgtcgtc < 1:75058/71‑1 (MQ=255) cATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGc > 1:386287/1‑70 (MQ=255) tCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAAt > 1:939073/1‑52 (MQ=255) ccTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCCGCTCTCCGGCAATTCGTCGTCTGCAACAAAAGa > 1:1680405/1‑71 (MQ=255) ccTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAGCTCTCCGGCAATTCGTCGTCTGCAACaaaa > 1:1415500/1‑69 (MQ=255) ccTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAGCTCTCCGGCAATTCGTCGTCTGCAACaaaa > 1:56876/1‑69 (MQ=255) ccTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAGCTCTCCGGCAATTCGTCGTCTGCAACaaaa > 1:1922429/1‑69 (MQ=255) ccTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGcaaca < 1:454686/66‑1 (MQ=255) ccTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGcaaca < 1:1944197/66‑1 (MQ=255) ccTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGcaaca < 1:419556/66‑1 (MQ=255) ccTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACaaaa > 1:360874/1‑69 (MQ=255) ccTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGa > 1:1782474/1‑71 (MQ=255) cGGATCTTCGAGAGCGATCGCGGGCGTATCATCAGCTCTCCGGCAATTCGTCGTCTGCAACAAAAGAcc > 1:1258005/1‑69 (MQ=255) ttCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtcgt < 1:2007086/45‑1 (MQ=255) ttCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtcgt < 1:1816160/45‑1 (MQ=255) cGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGttttt < 1:401322/70‑1 (MQ=255) cGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGtttt < 1:428768/69‑1 (MQ=255) cGAGAGCGATAGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGt < 1:2190895/66‑1 (MQ=255) agATCGCGGGCGTATCATCAACTCTCCAGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTg < 1:1387091/69‑1 (MQ=255) gcgGGCGTATCATCAGCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAgcgc < 1:30240/71‑1 (MQ=255) gcgGGCGTATCATCAGCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAgcgc < 1:379985/71‑1 (MQ=255) gcgGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCGAGGTTTTTCCACTGGAgcgc < 1:1433608/71‑1 (MQ=255) gcgGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAgcgc < 1:1962600/71‑1 (MQ=255) tcatcaGCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATgccgcc > 1:473124/1‑71 (MQ=255) catcaACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGccgccg > 1:1442829/1‑71 (MQ=255) atcaACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGt > 1:1158261/1‑71 (MQ=255) atcaACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGt > 1:1551578/1‑71 (MQ=255) tcaACTCTCCGGCAATTCGGCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTg < 1:500195/71‑1 (MQ=255) aaCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTgcg < 1:928591/71‑1 (MQ=255) | GCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCG > minE/158118‑158252 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |