Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 3,342,964 | G→A | 9.0% | Q15* (CAG→TAG) | yefM ← | YoeB‑YefM toxin‑antitoxin system antitoxin/DNA‑binding transcriptional repressor |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 3,342,964 | 0 | G | A | 9.0% | 120.1 / 2.9 | 45 | Q15* (CAG→TAG) | yefM | YoeB‑YefM toxin‑antitoxin system antitoxin/DNA‑binding transcriptional repressor |
Reads supporting (aligned to +/- strand): ref base G (19/21); new base A (2/2); total (21/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01 |
AGCGCGAAGCGTTTTTGGACGATCCCAACAGGTGAATGGTCTCGTTCATCCCGTTGTAGTCCTCCAGAGACATCATCACTACGGGTTCGCCACGCTGTCGCGTGATTACGGCGGGCTCATGGTCCCGGCAAACATCATCCATTGTCTGCTTCAGTTCGGCGCGAGCCTGGCTAAAAGTCAGTACGTGCATTTGTCCATCCTCCAGGCTCTGCCCGGGTAAGGGTTGCTCTTTGGACCGCCACAGGTGACGCCCAAAGTTCTAAGAAAGTGAAGCCTGTCGACCCGCACGGCG > NC_002947/3342817‑3343108 | aGCGCGAAGCGTTTTTGGACGATCCCAACAGGTGAATGGTCTCGTTCATCCCGTTGTAGTCCTCCAGAGACATCATCACTACGGGTTCGCCACGCTGTCGCGTGATTACGGCGGGCTCATGGTCCCGGCAAACATCATCCATTGTCTGc < 2:85176/149‑1 (MQ=255) gcgcGAAGCGTTTTTGGACGATCCCAACAGGTGAATGGTCTCGTTCATCCCGTTGTAGTCCTCCAGAGACATCATCACTACGGGTTCGCCACGCTGTCGCGTGATTACGGCGGGCTCATGGTCCCGGCAAACATCATCCATTGTCTGCt > 5:53987/1‑149 (MQ=255) gACGATCCCAACAGGTGAATGGTCTCGTTCATCCCGTTGTAGTCCTCCAGAGACATCATCACTACGGGTTCGCCACGCTGTCGCGTGATTACGGCGGGCTCATGGTCCCGGCAAACATCATCCATTGTCTGCTTCAGTTCGGCGCGAGc > 4:52424/1‑149 (MQ=255) gATCCCAACAGGTGAATGGTCTCGTTCATCCCGTTGTAGTCCTCCAGAGACATCATCACTACGGGTTCGCCACGCTGTCGCGTGATTACGGCGGGCTCATGGTCCCGGCAAACATCATCCATTGTCTACTTCAGTTCGGCGCGAGCCTg > 4:87384/1‑149 (MQ=255) gATCCCAACAGGTGAATGGTCTCGTTCATCCCGTTGTAGTCCTCCAGAGACATCATCACTACGGGTTCGCCACGCTGTCGCGTGATTACGGCGGGCTCATGGTCCCGGCAAACATCATCCATTGTCTACTTCAGTTCGGCGCGAGCCTg > 2:202163/1‑149 (MQ=255) aTCCCAACAGGTGAATGGTCTCGTTCATCCCGTTGTAGTCCTCCAGAGACATCATCACTACGGGTTCGCCACGCTGTCGCGTGATTACGGCGGGCTCATGGTCCCGGCAAACATCATCCATTGTCTGCTTCAGTTCGGCGCGAGCCTgg < 1:107392/149‑1 (MQ=255) aCAGGTGAATGGTCTCGTTCATCCCGTTGTAGTCCTCCAGAGACATCATCACTACGGGTTCGCCACGCTGTCGCGTGATTACGGCGGGCTCATGGTCCCGGCAAACATCATCCATTGTCTGCCTCAGTTCGGCGCGAGCCTGGCTaaaa < 3:101581/149‑1 (MQ=255) aCAGGTGAATGGTCTCGTTCATCCCGTTGTAGTCCTCCAGAGACATCATCACTACGGGTTCGCCACGCTGTCGCGTGATTACGGCGGGCTCATGGTCCCGGCAAACATCATCCATTGTCTACTTCAGTTCGGCGCGAGCCTGGCTaaaa < 3:87384/149‑1 (MQ=255) aCAGGTGAATGGTCTCGTTCATCCCGTTGTAGTCCTCCAGAGACATCATCACTACGGGTTCGCCACGCTGTCGCGTGATTACGGCGGGCTCATGGTCCCGGCAAACATCATCCATTGTCTACTTCAGTTCGGCGCGAGCCTGGCTaaaa < 1:202163/149‑1 (MQ=255) gTGAATGGTCTCGTTCATCCCGTTGTAGTGCTCCAGAGACATCATCACTACGGGTTCGCCACGCTGTCGCGTGATTACGGCGGGCTCATGGTCCCGGCAAACATCATCCATTGTCTGCTTCAGTTCGGCGCGAGCCTGGCTAAAAGTCa < 5:72624/149‑1 (MQ=255) gAATGGTCTCGTTCATCCCGTTGTAGTCCTCCAGAGACATCATCACTACGGGTTCGCCACGCTGTCGCGTGATTACGGCGGGCTCATGGTCCCGGCAAACATCATCCATTGTCTGCTTCAGTTCGGCGCGAGCCTGGCTAAAAga > 6:8168/1‑144 (MQ=255) gAATGGTCTCGTTCATCCCGTTGTAGTCCTCCAGAGACATCATCACTACGGGTTCGCCACGCTGTCGCGTGATTACGGCGGGCTCATGGTCCCGGCAAACATCATCCATTGTCTGCTTCAGTTCGGCGCGAGCCTGGCTAAAAga < 5:8168/145‑2 (MQ=255) gTCTCGTTCATCCCGTTGTAGTCCTCCAGAGACATCATCACTGCGGGTTCGCCACGCTGTAGCGTGATTACGGCGGGCTCATGGTCCCGGCAAACATCATCCATTGTCTGCTTCAGTTAGGCGCGAGCCTGGCTAAAAGTCAGTACGTg > 2:218584/1‑149 (MQ=255) tcGTTCATCCCGTTGTAGTCCTCCAGAGACATCATCACTACGGGTTCGCCACGCTGTCGCGTGATTACGGCGGGCTCATGGTCCCGGCAAACATCATCCATTGTCTGCTTCAGTTCGGCGCGAGCCTGGCTAAAAGTCAGt < 6:177142/141‑1 (MQ=255) tcGTTCATCCCGTTGTAGTCCTCCAGAGACATCATCACTACGGGTTCGCCACGCTGTCGCGTGATTACGGCGGGCTCATGGTCCCGGCAAACATCATCCATTGTCTGCTTCAGTTCGGCGCGAGCCTGGCTAAAAGTCAGt > 5:177142/1‑141 (MQ=255) gTTCATCCCGTTGTAGTCCTCCAGAGACATCATCACTACGGGTTCGCCACGCTGTCGCGTGATTACGGCGGGCTCATGGTCCCGGCAAACATCATCCATTGTCTGCTTCAGTTCGGCGCGAGCCTGGCTAAAAGTCAGTACGTGCAt < 2:267548/147‑1 (MQ=255) gatgtCCCGTTGTAGTCCTCCAGAGACATGATCACTACGGTTTCGCCTCGCTGTCGCGTGATTACGGCGGGCTCATGGTCCCGGCAAAGATCATCCATTGTATGCTTCAGTTCGGCGCTAGCCTGGCTAAAAGTCAGTACGTGCATTTg < 5:63745/145‑1 (MQ=255) tcctccAGAGACATCATCACTACGGGTTCGCCACGCTGTCGCGTGATTACGGCGGGCTCATGGTCCCGGCAAACATCATCCATTGTCTGCTTCAGTTCGGCGCGAGCCTGGCTAAAAGTCAGTACGTGCATTTGTCCATCCTCCAGGct > 4:19625/1‑149 (MQ=255) ctccAGAGACATCATCACTACGGGTTCGCCACGCTGTCGCGTGATTACGGCGGGCTCATGGTCCCGGCAAACATCATCCATTGTCTGCTTCAGTTCGGCGCGAGCCTGGCTAAAAGTCAGTACGTGCATTTGTCCATCCTCCAGGctct < 3:52424/149‑1 (MQ=255) ctccAGAGACATCATCACTACGGGTTCGCCACGCTGTCGCGTGATTACGGCGGGCTCATGGTCCCGGCAAACATCATCCATTGTCTGCTTCAGTTCGGCGCGAGCCTGGCTAAAAGTCAGTACGTGCATTTGTCCATCCGCCAGGctct > 2:257032/1‑149 (MQ=255) tccAGAGACATCATCACTACGGGTTCGCCACGCTGTCGCGTGATTACGGCGGGCTCATGGTCCCGGCAAACATCATCCATTGTCTGCTTCAGTTCGGCGCGAGCCTGGCTAAAAGTCAGTACGTGCATTTGTCCATCCTCCAGGCTCTg > 5:207455/1‑149 (MQ=255) tccAGAGACAACAACACTACGCGTTCGCCACCCTGACGCTTGATTACGGCGGGCTCATGGTCCCGGCAAAAATCATCCATTGTCTGCTTCAGTTCGGCGCGACCCTGGCTAAAAGTCAGTACGTGCATTTGGCCATCCTCCAGCctcag > 6:266415/1‑147 (MQ=255) ccAGAGACATCATCACTACGGGTTCGCCACGCTGTCGCGTGATTACGGCGGGCTCATGGTCCCGGCAAACATCATCCATTGTCTGCt > 3:176607/1‑87 (MQ=255) ccAGAGACATCATCACTACGGGTTCGCCACGCTGTCGCGTGATTACGGCGGGCTCATGGTCCCGGCAAACATCATCCATTGTCTGCt < 4:176607/87‑1 (MQ=255) gaCATCATCACTACGGGTTCGCCACGCTGTCGCGTGATTACGGCGGGCTCATGGTCCCGGCAAACATCATCCATTGTCTGCTTCAGTTCGGCGCGAGCCTGGCTAAAAGTCAGTACGTGCATTTGTCCa < 8:130900/129‑1 (MQ=255) gaCATCATCACTACGGGTTCGCCACGCTGTCGCGTGATTACGGCGGGCTCATGGTCCCGGCAAACATCATCCATTGTCTGCTTCAGTTCGGCGCGAGCCTGGCTAAAAGTCAGTACGTGCATTTGTCCa > 7:130900/1‑129 (MQ=255) atcaCTACGGGTTCGCCACGCTGTCGCGTGATTACGGCGGGCTCATGGTCCCGGCAAACATCATCCATTGTCTGCTTCAGTTCGGCGCGAGCCTGGCTAAAAGTCAGTACGTGCATTTGTCCATCCTCCAGGCTCTGCCCGGGTAAgg < 3:19625/148‑1 (MQ=255) cACGCTGTCGCGTGATTACGGCGGGCTCATGGTCCCGGCAAACATCATCCATTGTCTGCTTCAGTTCGGCGCGAGCCTGGCTAAAAGTCAGTACGTGCATTTGTCCATCCTCCAGGCTCTGCCCGGGTaa > 5:221575/1‑130 (MQ=255) cACGCTGTCGCGTGATTACGGCGGGCTCATGGTCCCGGCAAACATCATCCATTGTCTGCTTCAGTTCGGCGCGAGCCTGGCTAAAAGTCAGTACGTGCATTTGTCCATCCTCCAGGCTCTGCCCGGGTaa < 6:221575/130‑1 (MQ=255) tGTCGCGTGATTACGGCGGGCTCATGGTCCCGGCAAACATCATCCATTGTCTGCTTCAGTTCGGCGCGAGCCTGGCTAAAAGTCAGTACGTGCATTTGTCCATCCTCCAGGCTCTGCCCGGGTAAGGGTTGCTCTTTGGACCGCCACAg > 3:190618/1‑149 (MQ=255) gcgTGATTACGGCGGGCTCATGGTCCCGGCAAACATCATCCATTGTCTGCTTCAGTTCGGCGCGAGCCTGGCTAAAAGTCAGTACGTGCATTTGTCCATCCTCCAGGCTCTGCCCGGGTAAGGGTTGCTCTTTGGACCGCCACAGGTGa < 6:53987/149‑1 (MQ=255) aTTACGGCGGGCTCATGGTCCCGGCAAACATCATCCATTGTCTGCTTCAGTTCGGCGCGAGCCTGGCTAAAAGTCAGTACGTGCATTTGTCCATCCTACAGGCTCTGCCCGGGTAAGGGTTGCTCTTTGGACCGCAACAGGTGACGcc > 7:252207/1‑148 (MQ=255) ggcggGCTCATGGTCCCGGCAAACATCATCCCTTGTCTGCTTCAGTTCGGCGCGAGCCTGGCTAAAAGTCAGTACGTGCATTTGTCCATCCTCCAGGCTCTGCCCGGGTAAGGGTTGCTCTTTGGACCGCCACAGGTGACGCCCAAAGt < 6:182917/149‑1 (MQ=255) gcggGCTCATGGTCCCGGCAAACTTCATCCATTGTCTTCTTCAGTTCGGCGCGAGCCTGTCTAAAAGTCAGTACGTGCATTTGTCCATCCTCCAGGCTCTGCCCGGGTAAGGGTTGCTCTTTGGACCGCCACAGGTGACGCCCAAAGtt < 5:266415/149‑1 (MQ=255) gCTCATGGTCCCGGCAAACATCATCCATTGTCTGCTTCAGTTCGGCGCGAGCCTGGCTAAAAGTCAGTACGTGCATTTGTCCATCCGCCAGGCTCTGCCCGGGTAAGGGTTGCTCTTTGGACCGCCACAGGTGACGCCCAAAGTTCTaa < 1:257032/149‑1 (MQ=255) tCATGGTCCCGGCAAACATCATCCATTGTCTGCTTCAGTTCGGCGCGAGCCTGGCTAAAAGTCAGTACGTGCATTTGTCCATCCTCCAGGCTCTGCCCGGGTAAGGGTTGCTCTTTGGACCGCCACAGGTGACGCCCAAAGTTCTAAg > 8:174989/1‑148 (MQ=255) aTGGTCCCGGCAAACATCATCCATTGTCTGCTTCAGTTCGGCGCGAGCCTGGCTAAAAGTCAGTACGTGCATTTGTCCATCCTCCAGGCTCTGCCCGGGTAAGGGTTGCTCTTTGGACCGCCACAGGTGACGCCCAAAGTTCTAAGa < 6:207455/147‑1 (MQ=255) tCCCGGCAAACATCATCCATTGTCTGCTTCAGTTCGGCGCGAGCCTGGCTAAAAGTCAGTACGTGCATTTGTCCATCCTCCAGGCTCTGCCCGGGTAAGGGTTGCTCTTTGGACCGCCACAGGTGACGCCCAAAGTTCTAAGAAAGTGa < 7:174989/149‑1 (MQ=255) cGGCAAACATCATCCATTGTCTGCTTCAGTTCGGCGCGAGCCTGGCTAAAAGTCAGTACGTGCATTTGTCCATCCTCCAGGCTCTGCCCGGGTAAGGGTTGCTCTTTGGACCGCCACAGGTGACGCCCAAAGTTCTAAGAAAGTGaa < 3:195981/147‑1 (MQ=255) cGGCAAACATCATCCATTGTCTGCTTCAGTTCGGCGCGAGCCTGGCTAAAAGTCAGTACGTGCATTTGTCCATCCTCCAGGCTCTGCCCGGGTAAGGGTTGCTCTTTGGACCGCCACAGGTGACGCCCAAAGTTCTAAGAAAGTGaa > 4:195981/1‑147 (MQ=255) cGGCAAACATCATCCATTGTCTGCTTCAGTTCGGCGCGAGCCTGGCTAAAAGTCAGTACGTGCATTTGTCCATCCTCCAGGCTCTGCCCGGGTAAGGGTTGCTCTTTGGACCGCCACAGGTGACGCCCAAAGTTCTAAGAAAGTGAAGc > 4:48341/1‑149 (MQ=255) cAAACATCAACCAGTGTGCGCTTCACTTCGGCGCGACCCGGGCTAAAACTCACTACGTCCATTTGTCCATCCTCCAGGCTCTGCCCGGGTAAGGGTTGCTCTTTGGACCGCCACAGGTCACGCCCAAAGTTCTAAGAAAGTGAAGCCTg < 4:9650/149‑1 (MQ=255) taatcCATTGTCTGCTTCAGTTCGTCGCGAGCCTGGCTAAAAGTCAGAACGTGCATTTGTCCATCCTACAGGCTCTGCCCGGGTAAGGGTTGCTCTTTGGACCGCCACAGGTGACGCCCAAAGTTCTAAGAAAGTGAAGCCTGTCGAcc < 8:252207/147‑1 (MQ=255) tcCATTGTCTGCTTCAGTTCGGCGCGAGCCTGGCTAAAAGTCAGTACGTGCATTTGTCCATCCTCCAGGCTCTGCCCGGGTAAGGGTTGCTCTTTGGACCGCCACAGGTGACGCCCAAAGTTCTAAGAAAGTGAAGCCTGTCGACCCg > 1:199489/1‑148 (MQ=255) gTCTGCTTCAGTTCGGCGCGAGCCTGGCTAAAAGTCAGTACGTGCATTTGTCCATCCTCCAGGCTCTGCCCGGGTAAGGGTTGCTCTTTGGACCGCCACAGGTGACGCCCAAAGTTCTAAGAAAGTGAAGCCTGTCGACCCGCAcggcg > 7:68555/1‑149 (MQ=255) | AGCGCGAAGCGTTTTTGGACGATCCCAACAGGTGAATGGTCTCGTTCATCCCGTTGTAGTCCTCCAGAGACATCATCACTACGGGTTCGCCACGCTGTCGCGTGATTACGGCGGGCTCATGGTCCCGGCAAACATCATCCATTGTCTGCTTCAGTTCGGCGCGAGCCTGGCTAAAAGTCAGTACGTGCATTTGTCCATCCTCCAGGCTCTGCCCGGGTAAGGGTTGCTCTTTGGACCGCCACAGGTGACGCCCAAAGTTCTAAGAAAGTGAAGCCTGTCGACCCGCACGGCG > NC_002947/3342817‑3343108 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |