Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 1,603,289 | T→C | 67.4% | V105V (GTA→GTG) | yfeX ← | conserved hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 1,603,289 | 0 | T | C | 67.4% | 47.0 / 27.5 | 58 | V105V (GTA→GTG) | yfeX | conserved hypothetical protein |
Reads supporting (aligned to +/- strand): ref base T (5/14); new base C (28/11); total (33/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.67e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.52e-01 |
ACCAAAGGCTTCCATCGCCGCCTGGGCGACAGAGAAGTTTACGTCGTGACGCAGAGAAAGAATGTGGATCAACACATCGAACTGGGTCGTCGGCGCAAGGCCTT > minE/1603250‑1603353 | aCCAAAGGCTTCCATCGCCGCCTGGGCGACAGAGAAGTTTACGTCGTGACGCAGagaaagaa < 1:3390710/62‑1 (MQ=255) aCCAAAGGCTTCCATCGCCGCCTGGGCGACAGAGAAGTTTACGTCGTGACGCAGagaaagaa < 1:2532979/62‑1 (MQ=255) aCCAAAGGCTTCCATCGCCGCCTGGGCGACAGAGAAGTTTACGTCGTGACGCAGagaaagaa < 1:590872/62‑1 (MQ=255) aCCAAAGGCTTCCATCGCCGCCTGGGCGACAGAGAAGTTTACGTCGTGACGCAGAGAAAGAATGTgg > 1:223445/1‑67 (MQ=255) aCCAAAGGCTTCCATCGCCGCCTGGGCGACAGAGAAGTTTACGTCGTGACGCAGAGAAAGAATGTgg > 1:1775274/1‑67 (MQ=255) aCCAAAGGCTTCCATCGCCGCCTGGGCGACAGAGAAGTTTACGTCGTGACGCAGAGAAAGAATGTGGa > 1:63077/1‑68 (MQ=255) aCCAAAGGCTTCCATCGCCGCCTGGGCGACAGAGAAGTTTACGTCGTGACGCAGAGAAAGAATGTGGa > 1:3151881/1‑68 (MQ=255) aCCAAAGGCTTCCATCGCCGCCTGGGCGACAGAGAAGTTTACATCGTGACGCAGAGAAAGAATGTGGa > 1:600134/1‑68 (MQ=255) cAAAGGCTTCCATCGCCGCCTGGGCGACAGAGAAGTTTACGTCGTGACGCAGAGAAAGAATGTGGATc < 1:15764/68‑1 (MQ=255) aaGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGAGAAAGACTGTGGa < 1:2816/64‑1 (MQ=255) aaGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGAGAAAGAATGTGGa < 1:932031/64‑1 (MQ=255) aaGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGAGAAAGAATGTGGa < 1:864872/64‑1 (MQ=255) aaGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGAGAAAGAATGTGGa < 1:2999388/64‑1 (MQ=255) aaGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGAGAAAGAATGTGGa < 1:3059031/64‑1 (MQ=255) aaGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGAGAAAGAATGTGGa < 1:3331374/64‑1 (MQ=255) aaGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGAGAAAGAATGTGGa < 1:484490/64‑1 (MQ=255) aaGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGAGAAAGAATGTGGa < 1:155709/64‑1 (MQ=255) aaGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGAGAAAGAATGTGGa < 1:1271022/64‑1 (MQ=255) aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa > 1:1035698/1‑57 (MQ=255) aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa > 1:3291628/1‑57 (MQ=255) aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa > 1:3287824/1‑57 (MQ=255) aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa > 1:3269829/1‑57 (MQ=255) aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa > 1:417791/1‑57 (MQ=255) aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa > 1:304938/1‑57 (MQ=255) aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa > 1:878064/1‑57 (MQ=255) aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa > 1:2980380/1‑57 (MQ=255) aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa > 1:2722180/1‑57 (MQ=255) aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa > 1:618212/1‑57 (MQ=255) aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa > 1:2268555/1‑57 (MQ=255) aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa > 1:177270/1‑57 (MQ=255) aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa > 1:1447878/1‑57 (MQ=255) aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa > 1:2158382/1‑57 (MQ=255) aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa > 1:1444097/1‑57 (MQ=255) aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa > 1:2212989/1‑57 (MQ=255) aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa > 1:1142863/1‑57 (MQ=255) aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa > 1:2224051/1‑57 (MQ=255) aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa > 1:2278924/1‑57 (MQ=255) aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa > 1:2313071/1‑57 (MQ=255) aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa > 1:111971/1‑57 (MQ=255) aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa > 1:2329453/1‑57 (MQ=255) aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa > 1:1114636/1‑57 (MQ=255) aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa > 1:1007172/1‑57 (MQ=255) aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa > 1:249508/1‑57 (MQ=255) aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaaga > 1:562614/1‑56 (MQ=255) aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGAGAAAGAATGTGGa < 1:2618757/63‑1 (MQ=255) aGGCTTCCATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGAAGCagag > 1:1714011/1‑51 (MQ=255) aGGCTTACATCGCCGCCTGAGCGACAGAGAAGTTCACGTCGTGACGCAGagaaagaa > 1:425651/1‑57 (MQ=255) aTCGCCGCCTGGGCGACAGAGAAGTTTACGTCGTGACGCAGAGAAAGAATGTGGATCAACACATCGa < 1:1788709/67‑1 (MQ=255) aTCGCCGCCTGGGCGACAGAGAAGTTTACGTCGTGACGCAGAGAAAGAATGTGGATCAACACATCGa < 1:2039998/67‑1 (MQ=255) aTCGCCGCCTGGGCGACAGAGAAGTTTACGTCGTGACGCAGAGAAAGAATGTGGATCAACACATCGa < 1:2098250/67‑1 (MQ=255) aTCGCCGCCTGGGCGACAGAGAAGTTTACGTCGTGACGCAGAGAAAGAATGTGGATCAACACATCGa < 1:215028/67‑1 (MQ=255) aTCGCCGCCTGGGCGACAGAGAAGTTTACGTCGTGACGCAGAGAAAGAATGTGGATCAACACATCGa < 1:3223093/67‑1 (MQ=255) aTCGCCGCCTGGGCGACAGAGAAGTTTACGTCGTGACGCAGAGAAAGAATGTGGATCAACACATCGa < 1:1236697/67‑1 (MQ=255) aTCGCCGCCTGGGCGACAGAGAAGTTTACGTCGTGACGCAGAGAAAGAATGTGGATCAACACATCGa < 1:3068108/67‑1 (MQ=255) aTCGCCGCCTGGGCGACAGAGAAGTTTACGTCGTGACGCAGAGAAAGAATGTGGATCAACACATCGa < 1:2675735/67‑1 (MQ=255) ccgccTGAGCGACAGAGAAGTTCACGTCGTGACGCAGAGAAAGAATGTGGa < 1:2904597/51‑1 (MQ=255) agaAGTTTACGTCGTGACGCAGAGAAAGAATGTGGATCAACACATCGAACTGGGTCGTCGGCGCAAGGc < 1:674330/69‑1 (MQ=255) aGTTTACGTCGTGACGCAGAGAAAGAATGTGGATCAACACATCGAACTGGGTCGTCGGCGCAAGGCCtt < 1:2427565/69‑1 (MQ=255) | ACCAAAGGCTTCCATCGCCGCCTGGGCGACAGAGAAGTTTACGTCGTGACGCAGAGAAAGAATGTGGATCAACACATCGAACTGGGTCGTCGGCGCAAGGCCTT > minE/1603250‑1603353 |
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: ‑ |