| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 1,834,094 | A→G | 13.9% | D65G (GAT→GGT) | ygbM → | conserved hypothetical protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 1,834,094 | 0 | A | G | 13.9% | 81.9 / 8.2 | 36 | D65G (GAT→GGT) | ygbM | conserved hypothetical protein |
| Reads supporting (aligned to +/- strand): ref base A (14/17); new base G (3/2); total (17/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.50e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAG > minE/1834027‑1834147 | cTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGatat > 1:1423492/1‑71 (MQ=255) cAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATAtt < 1:196913/69‑1 (MQ=255) cAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTaa < 1:1362657/71‑1 (MQ=255) aaaTCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATAtt < 1:1314852/68‑1 (MQ=255) aaTCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATAtt < 1:1284517/67‑1 (MQ=255) tCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCg > 1:1313821/1‑71 (MQ=255) tCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCg > 1:1587310/1‑71 (MQ=255) ccAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCgg > 1:199377/1‑71 (MQ=255) ccAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCgg > 1:2152040/1‑71 (MQ=255) ccAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCgg > 1:1823232/1‑71 (MQ=255) cAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCgg < 1:1159377/70‑1 (MQ=255) cAAAAGCAACTGGAGCAAAACCATCTGACACTGGAGCTGTTCAATACCGCTCCCGGAGATATTAATGCggg < 1:1420321/71‑1 (MQ=255) aaaGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTa < 1:1779775/62‑1 (MQ=255) aaaGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCgggg > 1:1584621/1‑70 (MQ=255) gCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAAt < 1:519617/71‑1 (MQ=255) aCTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATAtt < 1:1441226/55‑1 (MQ=255) aCTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATAtt < 1:22223/55‑1 (MQ=255) tGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGGTATTAATGCGTGGGAATGGGGa > 1:1822591/1‑71 (MQ=255) tGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGa > 1:1617240/1‑71 (MQ=255) tGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGa > 1:1015414/1‑71 (MQ=255) gAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGGTATTAATGCGGGGGAATGGGGa < 1:2210405/69‑1 (MQ=255) cAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATc < 1:244371/71‑1 (MQ=255) aaCCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATg < 1:859102/59‑1 (MQ=255) aaCCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTa < 1:747944/66‑1 (MQ=255) aaCCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGc > 1:1977077/1‑71 (MQ=255) aCCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGcc > 1:1250511/1‑71 (MQ=255) tGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCgg > 1:1113127/1‑45 (MQ=255) tGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCgg > 1:1228485/1‑45 (MQ=255) gACACTGGCGCTGTTCAATACCGCTCCCGGAGGTAtt < 1:421386/37‑1 (MQ=255) gcTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGTATGGGGATTATCCGCCCTTCCCGGACGTGAGc < 1:1809023/71‑1 (MQ=255) ttCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGAc < 1:1829137/61‑1 (MQ=255) ttCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACg < 1:1618823/70‑1 (MQ=255) cAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAg > 1:1786723/1‑71 (MQ=255) cAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAg > 1:1729028/1‑71 (MQ=255) aaTACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAg < 1:917342/70‑1 (MQ=255) ccGCTCCCGGAGGTATTAATGCGGGGGAATGGGGATTATc > 1:1316465/1‑40 (MQ=255) ccGCTCCCGGAGGTATTAATGCGGGGGAATGGGGATTATc > 1:914483/1‑40 (MQ=255) | CTGCAAATCCAAAAGCAACTGGAGCAAAACCATCTGACACTGGCGCTGTTCAATACCGCTCCCGGAGATATTAATGCGGGGGAATGGGGATTATCCGCCCTTCCCGGACGTGAGCACGAAG > minE/1834027‑1834147 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |