Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 1,255,649 | G→A | 17.6% | A113A (GCG→GCA) | yecE → | conserved hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 1,255,649 | 0 | G | A | 17.6% | 116.2 / 16.7 | 51 | A113A (GCG→GCA) | yecE | conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base G (15/27); new base A (3/6); total (18/33) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01 |
TTGGCTCCGCGCATTGGGCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTT > minE/1255581‑1255719 | ttGGCTCCGCGCATTGGTCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCt < 1:1357052/71‑1 (MQ=255) ggCTCCGCGCATTGGGCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCttt > 1:1694014/1‑71 (MQ=255) ggCTCCGCGCATTGGGCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCttt > 1:2882525/1‑71 (MQ=255) ggCTCCGCGCATTGGGCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCtt < 1:2893194/70‑1 (MQ=255) ggCTCCGCGCATTGGGCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCtt < 1:3923784/70‑1 (MQ=255) ggCTCCGCGCATTGGGCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCtt < 1:2577941/70‑1 (MQ=255) cTCCGCGCATTGGGCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTg > 1:2270399/1‑70 (MQ=255) cATTGGGCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTtc < 1:1581803/71‑1 (MQ=255) ttGGGCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTtctc < 1:274454/71‑1 (MQ=255) cAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATtc < 1:2282760/71‑1 (MQ=255) aaTACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATtct > 1:2282387/1‑71 (MQ=255) aTACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCAttt < 1:1014249/60‑1 (MQ=255) aTACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATtctc > 1:3716368/1‑71 (MQ=255) aTACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATtctc > 1:3745330/1‑71 (MQ=255) tACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTctct < 1:789457/71‑1 (MQ=255) tACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTctct < 1:2739988/71‑1 (MQ=255) cTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTc < 1:1274923/71‑1 (MQ=255) cTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCACTTTGGCATTTTCTCg > 1:3328318/1‑62 (MQ=255) tGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTtctc < 1:1511582/60‑1 (MQ=255) gCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGa > 1:3441020/1‑60 (MQ=255) gCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATtct > 1:59007/1‑64 (MQ=255) tGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTtctc < 1:1032105/56‑1 (MQ=255) tGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCACTTTGGCATTTTCTCGATTCTCTTCCTg < 1:2447934/69‑1 (MQ=255) gCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATtctc > 1:2165122/1‑62 (MQ=255) gCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCACtt < 1:617042/43‑1 (MQ=255) cAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGa < 1:806358/71‑1 (MQ=255) aaCTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCACTTTGGCATTTTCTCGATTCTCTTCCTGGCGaa > 1:482184/1‑71 (MQ=255) cctgccACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTtct > 1:3874622/1‑47 (MQ=255) cctgccACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTa < 1:734921/70‑1 (MQ=255) tgccACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAAttt > 1:1180285/1‑67 (MQ=255) gccACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTa > 1:852514/1‑71 (MQ=255) ccACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCAtttac < 1:246899/42‑3 (MQ=255) ccACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTTGATTCTCTTCCTGGCGAATTTAATTa < 1:3889856/70‑1 (MQ=255) ccACATTCGGCCCACGGGAGCTGCCTGCACTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAAtt < 1:3209061/69‑1 (MQ=255) acaTTCGGCCCACGGGAGCTGCCTGCACTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATgg > 1:1936153/1‑71 (MQ=255) aTTCGGCCCACGGGAGCTGCCTGCACTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCg < 1:1036221/71‑1 (MQ=255) aTTCGGCCCACGGGAGCTGCCTGCACTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCg < 1:619601/71‑1 (MQ=255) ttCGGCCCACGGGAGCTGCCTGCACTTTGGCATTTTCTCGATTCTCTTCCTg < 1:1631846/52‑1 (MQ=255) ttCGG‑CCACGGGAGCAGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTa > 1:3249420/1‑60 (MQ=255) ggCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTgg < 1:3198055/70‑1 (MQ=255) ggCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTgg < 1:368080/70‑1 (MQ=255) ggCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTgg < 1:1984661/70‑1 (MQ=255) ccACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTctct < 1:1502424/41‑1 (MQ=255) cACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAAttt < 1:3150780/53‑1 (MQ=255) cACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAAttt < 1:3621757/53‑1 (MQ=255) aCGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTcc > 1:763412/1‑71 (MQ=255) cGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAAtt < 1:1903074/55‑1 (MQ=255) cGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCg > 1:3720488/1‑71 (MQ=255) gggAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAAt < 1:1637750/53‑1 (MQ=255) gggAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGc < 1:1480323/71‑1 (MQ=255) gCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATc < 1:3272407/71‑1 (MQ=255) gcgcTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGtt > 1:1371676/1‑71 (MQ=255) cgcTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTaa > 1:570281/1‑40 (MQ=255) cgcTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGttt < 1:409982/71‑1 (MQ=255) cgcTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGttt < 1:396966/71‑1 (MQ=255) gcTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGtttt > 1:3711623/1‑71 (MQ=255) gcTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGtttt < 1:769697/71‑1 (MQ=255) gcTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGtttt < 1:1457525/71‑1 (MQ=255) gcTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGtttt > 1:3224484/1‑71 (MQ=255) | TTGGCTCCGCGCATTGGGCAATACTGGCTGCAACTGCCTGCCACATTCGGCCCACGGGAGCTGCCTGCGCTTTGGCATTTTCTCGATTCTCTTCCTGGCGAATTTAATTATGGCGTGGAAGTCCGCCATCCACAGTTTT > minE/1255581‑1255719 |
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: ‑ |