| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 5,563 | G→A | 28.9% | intergenic (+33/+120) | 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,563 | 0 | G | A | 28.9% | 53.4 / 23.7 | 38 | intergenic (+33/+120) | yaaX/yaaA | hypothetical protein/conserved hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base G (18/9); new base A (2/9); total (20/18) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.13e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCT > minE/5501‑5632 | ggTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAg < 1:2602466/63‑1 (MQ=255) tCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGc < 1:3380927/71‑1 (MQ=255) aTGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGa > 1:3615515/1‑71 (MQ=255) gTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGc < 1:3169001/71‑1 (MQ=255) ccAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTg > 1:3252330/1‑71 (MQ=255) ggCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGgcg > 1:1573287/1‑71 (MQ=255) ggCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGgc > 1:2778058/1‑70 (MQ=255) gCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTg < 1:3062188/67‑1 (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) 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) | GGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCT > minE/5501‑5632 |
| 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: ‑ |