Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 1,874,658 | G→A | 36.2% | R81R (CGC→CGT) | ygcF ← | conserved hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 1,874,658 | 0 | G | A | 36.2% | 35.2 / 29.7 | 55 | R81R (CGC→CGT) | ygcF | conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base G (13/22); new base A (14/6); total (27/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.61e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GCAGCAAATCATGAATGCAAGGCTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTATCACTCTCTTTGGT > minE/1874594‑1874720 | gcagcaAATCATGAATGCAAGGCTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCg > 1:602503/1‑67 (MQ=255) gcagcaAATCATGAATGCAAGGCTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGa > 1:1643269/1‑68 (MQ=255) gcagcaAATCATGAATGCAAGGCTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGa > 1:396138/1‑68 (MQ=255) gcagcaAATCATGAATGCAAGGCTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGa > 1:561879/1‑68 (MQ=255) caAATCATGAATGCAAGGCTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAAt < 1:2138262/69‑1 (MQ=255) aaTCATGAATGCAAGGCTCACCACCCGTAATCACCACATGACGCGCGGTGTATCCCTGACGGCCAATGa < 1:895232/69‑1 (MQ=255) aaTCATGAATGCAAGGCTCACCACCCGTAATCACCACATGACGCGCGGTGTATCCCTGACGGCCAATGa < 1:2028702/69‑1 (MQ=255) aaTCATGAATGCAAGGCTCACCACCCGTAATCACCACATGACGCGCGGTGTATCCCTGACGGCCAATGa < 1:1606259/69‑1 (MQ=255) aaTCATGAATGCAAGGCTCACCACCCGTAATCACCACATGACGCGCGGTGTATCCCTGACGGCCAATGa < 1:1272269/69‑1 (MQ=255) aaTCATGAATGCAAGGCTCACCACCCGTAATCACCACATGACGCGCGGTGTATCCCTGACGGCCAATGa < 1:1236407/69‑1 (MQ=255) aTCATGAATGCAAGGCTCACCACCCGTAATCACCACATGACGCGCGGTGTATCCCTGACGGCCAATGa < 1:1931490/68‑1 (MQ=255) aTGAATGCAAGGCTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGaaa > 1:611790/1‑65 (MQ=255) aTGAATGCAAGGCTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGAca > 1:2993280/1‑67 (MQ=255) aTGAATGCAAGGCTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGAca > 1:1615310/1‑67 (MQ=255) aTGAATGCAAGGCTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGAca > 1:1216268/1‑67 (MQ=255) tGCAAGGCTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCa < 1:605000/69‑1 (MQ=255) tGCAAGGCTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCa < 1:2529125/69‑1 (MQ=255) tGCAAGGCTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCa < 1:2152122/69‑1 (MQ=255) tGCAAGGCTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCa < 1:1772166/69‑1 (MQ=255) tGCAAGGCTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCa < 1:1459656/69‑1 (MQ=255) tGCAAGGCTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCa < 1:1704249/69‑1 (MQ=255) gCAAGGCTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCa < 1:2570068/68‑1 (MQ=255) aGGCTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCAAATc < 1:877618/69‑1 (MQ=255) cTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCAAATCTTc < 1:1844686/69‑1 (MQ=255) cTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCAAATCTTc < 1:363965/69‑1 (MQ=255) cTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCAAATCTTc < 1:338510/69‑1 (MQ=255) cTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCAAATCTTc < 1:3194171/69‑1 (MQ=255) cTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCAAATCTTc < 1:2946219/69‑1 (MQ=255) cTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCAAATCTTc < 1:2623337/69‑1 (MQ=255) cTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCAAATCTTc < 1:1728858/69‑1 (MQ=255) accacATGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCAAATCTTCACTGCTCGCAGccccc > 1:1380270/1‑69 (MQ=255) accacATGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCAAATCTTCACTGCTCGCAGccccc > 1:1514705/1‑69 (MQ=255) ccacATGACGCGCGGTGTATCCCTGACGGCCAATGACAGCCAGCAAATCtt > 1:1575527/1‑51 (MQ=255) ccacATGACGCGCGGTGTATCCCTGACGGCCAATGACAGCCAGCAAATCTTCACTGc > 1:255814/1‑57 (MQ=255) ccacATGACGCGCGGTGTATCCCTGACGGCCAATGACAGCCAGCAAATCTTCACTGc > 1:1765361/1‑57 (MQ=255) ccacATGACGCGCGGTGTATCCCTGACGGCCAATGACAGCCAGCAAATCTTCACTGc > 1:1195847/1‑57 (MQ=255) ccacATGACGCGCGGTGTATCCCTGACGGCCAATGACAGCCAGCAAATCTTCACTGc > 1:1853926/1‑57 (MQ=255) ccacATGACGCGCGGTGTATCCCTGACGGCCAATGACAGCCAGCAAATCTTCACTGc > 1:1486934/1‑57 (MQ=255) ccacATGACGCGCGGTGTATCCCTGACGGCCAATGACAGCCAGCAAATCTTCACTGc > 1:1757782/1‑57 (MQ=255) ccacATGACGCGCGGTGTATCCCTGACGGCCAATGACAGCCAGCAAATCTTCACTGc > 1:3247514/1‑57 (MQ=255) ccacATGACGCGCGGTGTATCCCTGACGGCCAATGACAGCCAGCAAATCTTCACTGc > 1:3237963/1‑57 (MQ=255) ccacATGACGCGCGGTGTATCCCTGACGGCCAATGACAGCCAGCAAATCTTCACTGc > 1:3055221/1‑57 (MQ=255) ccacATGACGCGCGGTGTATCCCTGACGGCCAATGACAGCCAGCAAATCTTCACTGc > 1:210558/1‑57 (MQ=255) ccacATGACGCGCGGTGTATCCCTGACGGCCAATGACAGCCAGCAAATCTTCACTGc > 1:1067383/1‑57 (MQ=255) ccacATGACGCGCGGTGTATCCCTGACGGCCAATGACAGCCAGCAAATCTTCACTGc > 1:2400786/1‑57 (MQ=255) ccacATGACGCGCGGTGTATCCCTGACGGCCAATGACAGCCAGCAAATCTTCACTGc > 1:1679283/1‑57 (MQ=255) tGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTa < 1:3389335/69‑1 (MQ=255) tGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTa < 1:2093395/69‑1 (MQ=255) tGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTa < 1:318513/69‑1 (MQ=255) tGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTa < 1:601445/69‑1 (MQ=255) cgcgGTGTATCCCTGGCGACCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTATCAct < 1:360758/69‑1 (MQ=255) cgcgGTGTATCCCTGGCGACCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTATCAct < 1:433231/69‑1 (MQ=255) tCCCTGGCGCCCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCGCCCACTTATCACTCTCTTTGGt > 1:3318832/1‑69 (MQ=255) tCCCTGGCGACCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTATCACTCTCTTTGGt > 1:3325508/1‑69 (MQ=255) tCCCTGGCGACCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTATCACTCTCTTTGGt > 1:1152158/1‑69 (MQ=255) | GCAGCAAATCATGAATGCAAGGCTCACCACCCGTAATCACCACATGCCGCGCGGTGTATCCCTGGCGACCAATGACAGCCAGCAAATCTTCACTGCTCGCAGCCCCCCACTTATCACTCTCTTTGGT > minE/1874594‑1874720 |
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: ‑ |