| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 5,581 | G→A | 12.2% | intergenic (+51/+102) | yaaX → / ← yaaA | hypothetical protein/conserved hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 5,581 | 0 | G | A | 12.2% | 127.5 / 8.6 | 49 | intergenic (+51/+102) | yaaX/yaaA | hypothetical protein/conserved hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base G (22/21); new base A (2/4); total (24/25) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.67e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01 | |||||||||||
GGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGC > minE/5513‑5651 | ggCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGgcg > 1:1573287/1‑71 (MQ=255)ggCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGgc > 1:2778058/1‑70 (MQ=255) aCATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCt > 1:2055083/1‑71 (MQ=255) cACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCtt > 1:2002606/1‑68 (MQ=255) aaaTGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTa > 1:158555/1‑71 (MQ=255) aaTGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTa < 1:1861934/70‑1 (MQ=255) aaTGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAc < 1:1659446/71‑1 (MQ=255) aCAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGtt > 1:3717708/1‑70 (MQ=255) aaaTGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTaa < 1:979416/70‑1 (MQ=255) aTGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTc > 1:1329133/1‑71 (MQ=255) tAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAatat < 1:739098/71‑1 (MQ=255) cgcatccGGTATGTAACGCCTGATGCGACGCTGACGCGTCTTATCTGGCCTACacgctgcgattttgtagg < 1:2440899‑M1/68‑19 (MQ=255) cgcatccGGTATGTAACGCCTGATGCGACGCTGACGCGTCTTATCTGGCCTACacgctgcgattttgtagg < 1:2281548‑M1/68‑19 (MQ=255) cgcatccGGCATTCAACGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACgcgctgcgattttgtagg < 1:3492425‑M1/68‑19 (MQ=255) aCAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAatata > 1:3280752/1‑70 (MQ=255) gcatcCGGTATGTAACGCCTGATGCGACGCTGACGCGTCTTATCTGGCCTACacgctgcgattttgtaggc < 1:759492‑M1/69‑20 (MQ=255) cAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCa > 1:1480437/1‑64 (MQ=255) catCCGGCATTCAACGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACacgctgcgattttgtaggcc < 1:719446‑M1/70‑21 (MQ=255) aTCCGGCATTCAACGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACgcgctgcgatt < 1:315886‑M1/61‑12 (MQ=255) aTCCGGCATTCAACGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACacgctgcgattt < 1:392510‑M1/62‑13 (MQ=255) tCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAc > 1:3737090/1‑49 (MQ=255) tCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTTATTCTGCAATATATTGaa > 1:1590653/1‑71 (MQ=255) tCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAatata < 1:3868539/66‑1 (MQ=255) ccGGCATTCAACGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACgcgctgcgattttgtaggccgga > 1:2010629‑M1/1‑48 (MQ=255) ggCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCt < 1:2009738/71‑1 (MQ=255) ggCATTCAACGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCt < 1:167075/44‑1 (MQ=255) gCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTg > 1:1148074/1‑71 (MQ=255) gCATTCAACGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACacgctgcgattttgtaggccggat > 1:579630‑M1/1‑45 (MQ=255) aTTCAACGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACgcgctgcgattttgtaggccggataagc < 1:746257‑M1/71‑29 (MQ=255) tCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATg > 1:2321021/1‑71 (MQ=255) caacgcctgATGCGACGCTGGCGCGTCTTATCAGGCCTACgcgctgcgattttgtaggccggataagc > 1:317321‑M1/4‑40 (MQ=255) cAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTTATTCTGCAATATATTGAATCTGCATGc > 1:772863/1‑71 (MQ=255) cAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGc > 1:1961496/1‑71 (MQ=255) aGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCt > 1:3159980/1‑71 (MQ=255) gcgcCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGc < 1:2370522/65‑1 (MQ=255) gcCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCtttt > 1:1413923/1‑71 (MQ=255) ccTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCtttt < 1:1291823/70‑1 (MQ=255) aTGCGACGCTGGCGCGTCTTATCAGGCCTACacgctgcgattttgtaggccggataagcaaagcgcatc < 1:11235‑M1/69‑39 (MQ=255) aTGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTc < 1:1888829/39‑1 (MQ=255) aTGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAgg < 1:2561009/71‑1 (MQ=255) aTGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAg < 1:376284/70‑1 (MQ=255) aTGCGACGCTGACGCGTCTTATCTGGCCTACacgctgcgattttgtaggccggataagcaaagcgcatccg > 1:1996752‑M1/1‑31 (MQ=255) tGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAgg > 1:3251649/1‑70 (MQ=255) gACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAgg < 1:2569631/67‑1 (MQ=255) gACGCTGACGCGTCTTATCTGGCCTACacgctgcgattttgtaggccggataagc < 1:1042922‑M1/55‑29 (MQ=255) gACGCTGACGCGTCTTATCTGGCCTACacgctgcgatt > 1:1632913‑M1/1‑27 (MQ=255) cGCTGGCGCGTCTTATCAGGCCTACacgctgcgattttgtaggccggataagcaaagcgcatccggcacg > 1:1249894‑M1/1‑25 (MQ=255) cGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGa > 1:3668316/1‑70 (MQ=255) cGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGAt < 1:2723348/71‑1 (MQ=255) gCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGAt > 1:3482613/1‑70 (MQ=255) gCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATa < 1:1389208/71‑1 (MQ=255) ggCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGa < 1:841959/66‑1 (MQ=255) gcgcGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGc > 1:3199098/1‑71 (MQ=255) | GGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGC > minE/5513‑5651 |
| 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: ‑ |