| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_000913 | 3,620,646 | G→A | 17.4% | 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 | 17.4% | 52.3 / 4.2 | 23 | 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 (9/10); new base A (2/2); total (11/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
TAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTGACTTACCCTGCGCAACGGAAGATGCCACCGGCAGCCGGAAAACCA > NC_000913/3620612‑3620782 | gcccgtcaggcccggtggCTGACGGGCTGGAAATACGCCGGGAATATGATGAATGGGGCCGTCTGATTCAGGAAACTGcccc > 1:680945‑M2/19‑82 (MQ=255)gcccgtcaggcccggtggCTGACGGGCTGGAAATACGCCGGGAATATGATGAATGGGGCCGTCTGATTCAGGAAACTGcccc < 2:680945‑M2/64‑1 (MQ=255) cccgtcaggcccggtggCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGCCGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTg > 1:307187‑M2/18‑125 (MQ=255) cccgtcaggcccggtggCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTg < 2:307187‑M2/108‑1 (MQ=255) cgtcaggcccggtggCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTgg > 4:1007620‑M2/16‑54 (MQ=255) cgtcaggcccggtggCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTgg < 3:1007620‑M2/39‑1 (MQ=255) cgtcaggcccggtggCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTa > 4:808779‑M2/16‑106 (MQ=255) cgtcaggcccggtggCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTa < 3:808779‑M2/91‑1 (MQ=255) cgtcaggcccggtggCTGACGGGCTGGAAATACGCCGGGAATATGATGAATGGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTg > 4:276045‑M2/16‑123 (MQ=255) cgtcaggcccgggggCTGACGGGCGGGAAATACGCCGGGAATATGATGAATGGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTg < 3:276045‑M2/108‑1 (MQ=255) TGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCGGGAAACTGCCCCTGATGGCGa > 4:173027‑M2/1‑72 (MQ=255) TGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCGGGAAACTGCCCCTGACGGCGa < 3:173027‑M2/72‑1 (MQ=255) GACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGtt < 3:806904‑M2/88‑1 (MQ=255) GACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGtt > 4:806904‑M2/1‑88 (MQ=255) GACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATaa > 3:262918‑M2/1‑95 (MQ=255) GACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATaa < 4:262918‑M2/95‑1 (MQ=255) GACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTGACTTACCCTGCGCAACGGAAGATGCCACCGGCAGCCGGaaaa > 2:243401‑M2/1‑148 (MQ=255) ACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTGACTTACCCTGCGCAACGGAAGATGCCACCGGCAGCCGGAAAAc > 4:270654‑M2/1‑148 (MQ=255) CGGGCTGGAATTGCGCCGGGAATATGATGAATTGGTCCGTCTGATGCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTGACTTACCCTGCGCAACGGAAGATGCCACCGGCAGCCGGAAAAcc < 4:236934‑M2/148‑1 (MQ=255) CGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTGACTTACCCTGCGCAACGGAAGATGCCACCGGCAGCCGGAAAAcc > 3:552730‑M2/1‑148 (MQ=255) GGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTGACTTACCCTGCGTAACGGAAGATGCCACCGGCAGCCGGAAAACCa < 2:227526‑M2/148‑1 (MQ=255) GGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTGACTTACCCTGCGCAACGGAAGATGCCACCGGCAGCCGGAAAACCa < 3:803460‑M2/148‑1 (MQ=255) GGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTGACTTACCCTGCGCAACGGAAGATGCCACCGGCAGCCGGAAAACCa < 1:243401‑M2/148‑1 (MQ=255) | TAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAATTGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTATGATAATCCACACAGTGACTTACCCTGCGCAACGGAAGATGCCACCGGCAGCCGGAAAACCA > NC_000913/3620612‑3620782 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |