Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 3,620,646 | G→A | 28.6% | L485L (TTG→TTA) ‡ | rhsB → | Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,620,646 | 0 | G | A | 28.6% | 27.4 / 12.6 | 21 | L485L (TTG→TTA) ‡ | rhsB | Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor |
Reads supporting (aligned to +/- strand): ref base G (11/4); new base A (4/2); total (15/6) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01 |
AACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACAC > NC_000913/3620559‑3620733 | taaaacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTgc < 4:485061‑M2/66‑1 (MQ=255) taaaacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTgc > 3:485061‑M2/24‑89 (MQ=255) aaaacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGCGCTGGAATTGCCCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCAc > 4:659786‑M2/23‑149 (MQ=255) aacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAAATACGCCGGGAATatga > 1:463617‑M2/21‑99 (MQ=255) aacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAAATACGCCGGGAATatga < 2:463617‑M2/79‑1 (MQ=255) aacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAAATACGCCGGGAATATGATGAATGGGGCCGTCTGATTCAGGAAACTGcccc > 2:556310‑M2/21‑132 (MQ=255) aacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAAATACGCCGGGAATATGATGAATGGGGCCGTCTGATTCAGGAAACTGcccc < 1:556310‑M2/112‑1 (MQ=255) aacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAAATACGCCGGGAATATGATGAATGGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCAccc > 4:369770‑M2/21‑149 (MQ=255) aacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAAATACGCCGGGAATATGATGAATGGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCAccc > 2:199640‑M2/21‑149 (MQ=255) aacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGccc < 4:617683‑M2/111‑1 (MQ=255) aacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGccc > 3:617683‑M2/21‑131 (MQ=255) aacgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCAccc > 3:873854‑M2/21‑149 (MQ=255) acgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCg > 4:906702‑M2/20‑149 (MQ=255) acgccgatgccctgccatcCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCg > 2:841402‑M2/20‑149 (MQ=255) tttaatTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATcc > 3:231805‑M2/7‑149 (MQ=255) CGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCca > 4:112863‑M2/1‑149 (MQ=255) CGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCca > 1:148553‑M2/1‑149 (MQ=255) CGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTCGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCca > 2:35178‑M2/1‑149 (MQ=255) GGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCcac > 2:751485‑M2/1‑149 (MQ=255) GCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTAt < 4:164573‑M2/138‑1 (MQ=255) GCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTAt > 3:164573‑M2/1‑138 (MQ=255) GCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCcaca < 3:814223‑M2/149‑1 (MQ=255) CGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCcacac < 3:562155‑M2/149‑1 (MQ=255) | AACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAACCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACAC > NC_000913/3620559‑3620733 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |