Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 1,633,561 | T→C | 60.5% | G521G (GGT→GGC) | narQ → | sensory histidine kinase in two‑component regulatory system with NarP (NarL) |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 1,633,561 | 0 | T | C | 60.5% | 21.8 / 29.1 | 48 | G521G (GGT→GGC) | narQ | sensory histidine kinase in two‑component regulatory system with NarP (NarL) |
Reads supporting (aligned to +/- strand): ref base T (5/14); new base C (23/6); total (28/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.82e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
ACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTGACTTTTT > minE/1633503‑1633622 | acacGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGt > 1:645201/1‑68 (MQ=255) acacGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGt > 1:3014736/1‑68 (MQ=255) acacGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGt > 1:2768475/1‑68 (MQ=255) acacGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGt > 1:1406361/1‑68 (MQ=255) ttATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAatat < 1:1063339/69‑1 (MQ=255) ttATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAatat < 1:3069985/69‑1 (MQ=255) ttATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAatat < 1:252171/69‑1 (MQ=255) ttATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAatat < 1:1837728/69‑1 (MQ=255) ttATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAatat < 1:1787130/69‑1 (MQ=255) ttATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAatat < 1:1306956/69‑1 (MQ=255) ttATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAatat < 1:1088279/69‑1 (MQ=255) tataTCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTg < 1:914314/63‑1 (MQ=255) tataTCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTg < 1:239041/63‑1 (MQ=255) tCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATgc > 1:1239349/1‑69 (MQ=255) cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata > 1:335600/1‑53 (MQ=255) cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata > 1:159683/1‑53 (MQ=255) cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata > 1:1173032/1‑53 (MQ=255) cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata > 1:3016622/1‑53 (MQ=255) cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata > 1:3034686/1‑53 (MQ=255) cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata > 1:3095366/1‑53 (MQ=255) cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata > 1:3244194/1‑53 (MQ=255) cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata > 1:3006792/1‑53 (MQ=255) cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata > 1:404004/1‑53 (MQ=255) cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata > 1:492753/1‑53 (MQ=255) cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata > 1:767058/1‑53 (MQ=255) cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata > 1:90068/1‑53 (MQ=255) cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata > 1:946550/1‑53 (MQ=255) cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata > 1:994286/1‑53 (MQ=255) cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata > 1:2165866/1‑53 (MQ=255) cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata > 1:1621864/1‑53 (MQ=255) cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata > 1:1808971/1‑53 (MQ=255) cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata > 1:1857189/1‑53 (MQ=255) cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata > 1:2985394/1‑53 (MQ=255) cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata > 1:2253647/1‑53 (MQ=255) cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata > 1:2347877/1‑53 (MQ=255) cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGTCTGAata > 1:246209/1‑53 (MQ=255) cGGGATTGGCATCGGTGAACCGAAAGAACCCGAAGGCCATTATGGGCTGAat > 1:1904761/1‑52 (MQ=255) gATTGGTATAGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTg < 1:3358390/46‑1 (MQ=255) aTCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTa < 1:2452269/69‑1 (MQ=255) aTCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTa < 1:2880799/69‑1 (MQ=255) tCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTa < 1:1402288/68‑1 (MQ=255) aCCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTa < 1:2842367/52‑1 (MQ=255) ccGAATGCCATTATGGTCTGAATATCATGCGCGAACGCGCGGAAAGATTGGGTGGGACGCTGACttttt < 1:767594/69‑1 (MQ=255) ccGAAGGCCATTATGGTCTGAATATCATGCGCGAACGCGCGGAAAGATTGGGTGGGACTCTGACttttt < 1:1006933/69‑1 (MQ=255) ccGAAGGCCATTATGGTCTGAATATCATGCGCGAACGCGCGGAAAGATTGGGTGGGACGCTGACttttt < 1:3187553/69‑1 (MQ=255) ccGAAGGCCATTATGGTCTGAATATCATGCGCGAACGCGCGGAAAGATTGGGTGGGACGCTGACttttt < 1:210106/69‑1 (MQ=255) ccGAAGGCCATTATGGTCTGAATATCATGCGCGAACGCGCGGAAAGATTGGGTGGGACGCTGACttttt < 1:702739/69‑1 (MQ=255) ccGAAGGCCATTATGGTCTGAATATCATGCGCGAACGCGCGGAAAGATTGGGTGGGACGCTGACttttt < 1:1175865/69‑1 (MQ=255) | ACACGGTTTATATCCGTGATAACGGGATTGGTATCGGTGAACCGAAAGAACCCGAAGGTCATTATGGTCTGAATATCATGCGCGAACGCGCGGAACGGCTAGGTGGGACGCTGACTTTTT > minE/1633503‑1633622 |
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: ‑ |