| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 4,082,784 | C→A | 22.2% | L229M (CTG→ATG) ‡ | PP_3592 → | RpiR family transcriptional regulator |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 4,082,784 | 0 | C | A | 22.2% | 34.4 / 4.4 | 18 | L229M (CTG→ATG) ‡ | PP_3592 | RpiR family transcriptional regulator |
| Reads supporting (aligned to +/- strand): ref base C (6/8); new base A (2/2); total (8/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01 | |||||||||||
GGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAGGCTGGCATTGCGGTCACCCTGGTGACCGACACCTTCTGTGACTGGGCGGATGCCAATGCCGACGAGGTGTTCCGCATCCCCACCGAGTTCAACCTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTC > NC_002947/4082639‑4082930 | ggTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGGGCTGGTGGTGTTCGAGGCGCCCCGTTATTCCCGGCATGCCCTGATGCTGTGCCACAAGGCGCGCTAGGCTGGCATTGCGGTCAGCCTGGt > 2:118461/1‑150 (MQ=255) aCTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAGGCTGGCATTGCGGTCACCCTGGTGACCGACACCTTCTGTGACt < 1:9761/150‑1 (MQ=255) ttCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAGGCTGGCATTGCGGTCACCCTGGTGACCGACACCTTCTGTGACTgg > 1:40364/1‑150 (MQ=255) ctgctgGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAGGCTGGCATTGCGGTCACCCTGGTGACCGACACCTTCTGTGACTGGGCGGATGCCAAt > 1:292924/1‑150 (MQ=255) tGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAGGCTGGCATTGCGGTCACCCTGGTGACCGACACCTTCTGTGACTGGGCGGATGCCAATGCCGACGAGgt < 2:311111/150‑1 (MQ=255) tggtgTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAGGCTGGCATTGCGGTCACCCTGGTGACCGACACCTTCTGTGACTGGGCGGATGCCAATGCCGACGAGGTGTTCCGCATCCCCACCGAGt < 1:256369/150‑1 (MQ=255) ttCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAGGCTGGCATTGCGGTCACCCTGGTGACCGACACCTTCTGTGACTGGGCGGATGCCAATGCCGACGAGGTGTTCCGCATCCCCACCGAGTTCAAc < 2:292924/150‑1 (MQ=255) ccGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAGGCTGGCATTGCGGTCACCCTGGTGACCGACACCTTCTGTGACTGGGCGGATGCCAATGCcgac < 1:598727/100‑1 (MQ=255) ccGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAGGCTGGCATTGCGGTCACCCTGGTGACCGACACCTTCTGTGACTGGGCGGATGCCAATGCcgac > 2:598727/1‑100 (MQ=255) cccTGATGCTGTGCCAGAAGGCGCGCGAGGCTGGCATTGCGGTCACCCTGGTGACCGACACCTTCTGTGACTGGGCGGATGCCAATGCCGACGAGGTGTTCCGCATCCCCACCGAGTTCAACCTGTTCTGGGAATCCACCTCGACGATGc > 1:477169/1‑150 (MQ=255) cgAGGCTGGCATTGCGATCCCCCTGGTGACCGACACCTTCTGTGACTGGGCGGATGCCAATGCCGACGAGGTGTTCCGCATCCCCACCGAGTTCAACCTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCaa < 1:31569/150‑1 (MQ=255) aTTGCGGTCACCCTGGTGACCGACACCTTCTGTGACTGGGCGGATGCCAATGCCGACGAGGTGTTCCGCATCCCCACCGAGTTCAACCTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGc > 2:527643/1‑150 (MQ=255) gtcGGTCACCAGGGTGACCGACACCTTCTGTGACTGGGCGGATGCCAATGCCGACGAGGTGTTCCGCATCCCCACCGAGTTCAAcc > 1:544054/3‑86 (MQ=255) gtcGGTCACCAGGGTGACCGACACCTTCTGTGACTGGGCGGATGCCAATGCCGACGAGGTGTTCCGCATCCCCACCGAGTTCAAcc < 2:544054/84‑1 (MQ=255) ggTCACCAGGGTGACCGACACCTTCTGTGACTGGGCGGATGCCAATGCCGACGAGGTGTTCCGCATCCCCACCGAGTTCAACCTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATgg > 1:584998/1‑131 (MQ=255) ggTCACCAGGGTGACCGACACCTTCTGTGACTGGGCGGATGCCAATGCCGACGAGGTGTTCCGCATCCCCACCGAGTTCAACCTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATgg < 2:584998/131‑1 (MQ=255) gTCACCCTGGTGACCGACACCTTCTGTGACTGGGCGGATGCCAATGCCGACGAGGTGTTCCGCATCCCCACCGAGTTCAACCTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCaagaag < 1:80044/150‑1 (MQ=255) aCCCTGGTGACCGACACCTTCTGTGACTGGGCGGATGCCAATGCCGACGAGGTGTTCCGCATCCCCACCGAGTTCAACCTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTc < 2:477169/150‑1 (MQ=255) | GGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAGGCTGGCATTGCGGTCACCCTGGTGACCGACACCTTCTGTGACTGGGCGGATGCCAATGCCGACGAGGTGTTCCGCATCCCCACCGAGTTCAACCTGTTCTGGGAATCCACCTCGACGATGCTGTCGTGGGTGCATCTGATGGTCAACGAGGTCTGCAAGAAGCTC > NC_002947/4082639‑4082930 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |