| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 158,183 | A→G | 8.9% | 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 | 8.9% | 210.9 / 10.0 | 78 | N41S (AAC→AGC) | dgt | deoxyguanosine triphosphate triphosphohydrolase |
| Reads supporting (aligned to +/- strand): ref base A (34/37); new base G (2/5); total (36/42) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.42e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGC > minE/158114‑158253 | cACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCaa > 1:2269278/1‑70 (MQ=255)cACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAAc > 1:3311004/1‑71 (MQ=255)cACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAAc > 1:2702138/1‑71 (MQ=255)cACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAAc > 1:581969/1‑71 (MQ=255)cACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAAc > 1:2161235/1‑71 (MQ=255)cACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAAc > 1:3811614/1‑71 (MQ=255)cACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAAc > 1:3257464/1‑71 (MQ=255) gCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACtctc > 1:468239/1‑71 (MQ=255) cAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTcc > 1:1490393/1‑71 (MQ=255) ggCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGc < 1:2620094/71‑1 (MQ=255) ggCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGc < 1:1829469/71‑1 (MQ=255) ggCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGc < 1:3762994/71‑1 (MQ=255) aaaaCCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtcg > 1:3946766/1‑71 (MQ=255) aaaaCCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtcg > 1:2015898/1‑71 (MQ=255) aaCCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGc < 1:4041535/63‑1 (MQ=255) aaCCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtcgtc < 1:2719718/71‑1 (MQ=255) aaCCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtcgtc < 1:2685141/71‑1 (MQ=255) aaCCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtcgtc < 1:3533103/71‑1 (MQ=255) aaCCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtcgtc < 1:1809292/71‑1 (MQ=255) aaCCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtcgtc < 1:3600398/71‑1 (MQ=255) aaCCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtcgtc < 1:1728757/71‑1 (MQ=255) aCCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAGCTCTCCGGCAATtcgtcg > 1:3099450/1‑71 (MQ=255) aCCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAGCTCTCCGGCAATtcgtcg > 1:1042433/1‑71 (MQ=255) aaCATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTg > 1:3893644/1‑71 (MQ=255) tGAGATCCTGCGGATCTTCGCGAGCGATCGCGGGCGTATCATCAACTCTCCGGc < 1:2049889/54‑1 (MQ=255) tGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGc < 1:1237388/54‑1 (MQ=255) aTCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCa < 1:193947/51‑1 (MQ=255) aTCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCa < 1:3085379/51‑1 (MQ=255) aTCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCa < 1:930974/51‑1 (MQ=255) aTCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTcgtcgt < 1:3476611/60‑1 (MQ=255) aTCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTcgtcgt < 1:406188/60‑1 (MQ=255) aTCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTcgtcgt < 1:747974/60‑1 (MQ=255) aTCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTcgtcgt < 1:1046902/60‑1 (MQ=255) ccTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGc < 1:935678/48‑1 (MQ=255) ccTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAAt > 1:1727326/1‑51 (MQ=255) ccTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAAt > 1:608909/1‑51 (MQ=255) ccTGCGGATCTTCGAGAGAGATCGCGGGCGTATCATCAACTCTCCGGCAAt > 1:1408253/1‑51 (MQ=255) cTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCaa > 1:1061837/1‑38 (MQ=255) cTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCaa > 1:2005104/1‑38 (MQ=255) cTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCaa > 1:3539702/1‑38 (MQ=255) cGGATCTTCGAGAGCGATCGCGGGCGTATCATCaa > 1:3062284/1‑35 (MQ=255) cGGATCTTCGAGAGCGATCGCGGGCGTATCATCaa > 1:369617/1‑35 (MQ=255) cGGATCTTCGAGAGCGATCGCGGGCGTATCATCaa > 1:314346/1‑35 (MQ=255) cGGATCTTCGAGAGCGATCGCGGGCGTATCATCaa > 1:3083353/1‑35 (MQ=255) cGGATCTTCGAGAGCGATCGCGGGCGTATCATCaa > 1:2910672/1‑35 (MQ=255) cGGATCTTCGAGAGCGATCGCGGGCGTATCATCaa > 1:765517/1‑35 (MQ=255) cGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGc < 1:3958429/44‑1 (MQ=255) cGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGc < 1:3392006/44‑1 (MQ=255) gATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCaa > 1:4120019/1‑44 (MQ=255) gATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCaa > 1:1124383/1‑44 (MQ=255) gATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCaa > 1:775966/1‑44 (MQ=255) cGAGAGCGATCGCGGGCGTATCATCAGCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGtttt < 1:4046756/69‑1 (MQ=255) cGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACaaaa > 1:2625315/1‑57 (MQ=255) cGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACaaaa > 1:140325/1‑57 (MQ=255) cGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACaaaa > 1:3954607/1‑57 (MQ=255) cGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACaaaa > 1:416013/1‑57 (MQ=255) cGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACaaaa > 1:69997/1‑57 (MQ=255) cGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACaaaa > 1:2107359/1‑57 (MQ=255) cGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACaaaa > 1:2137144/1‑57 (MQ=255) cGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTc > 1:4001880/1‑71 (MQ=255) gagaGCGATCGCGGGCGTATCATCAGCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTcc < 1:436470/71‑1 (MQ=255) gagaGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTcc < 1:3120657/71‑1 (MQ=255) agagCGATCGCGGGCGTATCATCAGCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTcc < 1:3054493/70‑1 (MQ=255) gagCGATCGCGGGCGTATCATCAGCTCTCCGGCaa < 1:113001/35‑1 (MQ=255) cGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAAGAAAAGACCCAGGTTTTTCCACTgg < 1:546949/71‑1 (MQ=255) cgcgGGCGTATCATCAACTCTCCGGCAATtcgtcg > 1:1603737/1‑35 (MQ=255) cgcgGGCGTATCATCAACTCTCCGGCAATtcgtcg > 1:1788553/1‑35 (MQ=255) cgcgGGCGTATCATCAACTCTCCGGCAATtcgtcg > 1:2435522/1‑35 (MQ=255) cgcgGGCGTATCATCAACTCTCCGGCAATtcgtcg > 1:43485/1‑35 (MQ=255) gcgGGCGTATCATCAGCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAgcgc < 1:1030411/71‑1 (MQ=255) gcgGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAgcgc < 1:3164541/71‑1 (MQ=255) gcgGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAgcgc < 1:3287801/71‑1 (MQ=255) gcgGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAgcgc < 1:2958982/71‑1 (MQ=255) gcgGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAgcgc < 1:2827919/71‑1 (MQ=255) gcgGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAgcgc < 1:1384188/71‑1 (MQ=255) tATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGtttt < 1:3700442/51‑1 (MQ=255) tcatcaACTCTCCTGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAgcgc < 1:1664011/62‑1 (MQ=255) tcatcaACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAgcgc < 1:774682/62‑1 (MQ=255) tcatcaACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATgccgcc < 1:2623891/71‑1 (MQ=255) tcatcaACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATgccgcc > 1:3528180/1‑71 (MQ=255) tcatcaACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATgccgcc < 1:2797784/71‑1 (MQ=255) tcatcaACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATgccgcc > 1:2645643/1‑71 (MQ=255) tcatcaACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATgccgcc < 1:3804302/71‑1 (MQ=255) tcatcaACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATgccgcc > 1:1956349/1‑71 (MQ=255) tcatcaACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATgccgcc < 1:3357043/71‑1 (MQ=255) atcaACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGa < 1:2553225/56‑1 (MQ=255) caACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAgc > 1:693641/1‑56 (MQ=255) caACTCTCCGGCAATTCGTCCTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAgc > 1:1400503/1‑56 (MQ=255) aaCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTgcg > 1:2229491/1‑71 (MQ=255) aaCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTgc > 1:2613748/1‑70 (MQ=255) aCTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTgcgc < 1:2350061/71‑1 (MQ=255) | CACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCGTGCGC > minE/158114‑158253 |
| 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: ‑ |