| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 5,537,075 | G→C | 24.5% | intergenic (+52/+51) | nadE → / ← PP_4870 | NH(3)‑dependent NAD(+) synthetase/azurin |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 5,537,075 | 0 | G | C | 24.5% | 65.5 / 14.4 | 32 | intergenic (+52/+51) | nadE/PP_4870 | NH(3)‑dependent NAD(+) synthetase/azurin |
| Reads supporting (aligned to +/- strand): ref base G (18/6); new base C (2/5); total (21/11) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.65e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.24e-01 | |||||||||||
GCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTCG > NC_002947/5536930‑5537222 | gCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAACACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACg > 1:137475/1‑149 (MQ=255) cATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAACACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACg < 2:137475/148‑1 (MQ=255) gCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGC‑gg < 2:392863/149‑1 (MQ=25) aacGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGC‑gg < 2:420479/147‑1 (MQ=25) cGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGC‑GGGt < 2:288619/149‑1 (MQ=255) gTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGC‑GGGTGa < 5:234790/149‑1 (MQ=255) aGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGGACCCACCg < 8:265565/126‑1 (MQ=255) aGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGGACCCACCg > 7:265565/1‑126 (MQ=255) aaGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCTGCGAATACGAGGGTGAACCCACCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGc > 5:329927/1‑149 (MQ=255) gAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGaagaag > 3:174401/1‑149 (MQ=255) gAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGaagaag > 7:293510/1‑149 (MQ=255) gggTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGCCGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGt > 3:430088/1‑149 (MQ=255) ttCACCCTCGTATTCGCGGGTGAACCCACCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATAGAGATGTGGACCGGGAAGGTTCTGAAGATGCTGTAg > 8:208436/1‑148 (MQ=255) ttCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGt > 6:398464/1‑149 (MQ=255) tCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGC‑GTGTTAACCCGCTCCTGCAGGGTTCGCCGTCATCCTTTGATTACTTGACGATTACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAATAAGCt < 5:264230/143‑1 (MQ=18) tCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAg < 3:100732/129‑1 (MQ=18) tCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAg > 4:100732/1‑129 (MQ=18) tCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCt > 6:270122/1‑143 (MQ=21) tCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCt < 5:270122/143‑1 (MQ=21) tCACCCGCGAATGCGTCGGTGGGTTCACAATCGTATTCGC‑GGGTGAACCAGCTCCTGCAGGGTTAGCCGTCAGCCTTTGATTTCTTGACGATGACCTTGCCTTTCATCATCGAGATGTGTCACGGGAAGGTGCAGAAGAAGATgatatc > 6:264230/1‑144 (MQ=18) ccGCGAATACGAGGGTGAACCCTCCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCatcatc < 6:44511/107‑1 (MQ=255) ccGCGAATACGAGGGTGAACCCACCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCatcatc > 5:44511/1‑107 (MQ=255) ccGCGAATACGAGGGTGAACCCACCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcacca > 7:201674/1‑149 (MQ=255) ccGCGAATACGAGGGTGAACCCACCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcacca > 3:300219/1‑149 (MQ=255) ccGCGAATACGAGGGTGAACCCACCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcacca > 1:230628/1‑149 (MQ=255) ccGCGAATACGAGGGTGAACCCACCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCATTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcacca > 5:34267/1‑149 (MQ=255) ccGCGAATACGAGGGTGAACCCACCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGAAGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcacca > 2:378736/1‑149 (MQ=255) tACGAGGGTGAACCCACCGACGCATTCGCAGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCt > 1:34540/1‑149 (MQ=255) tcggTGAACCCACCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAg > 4:427812/3‑149 (MQ=255) tGAACCCACCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTg > 5:359063/1‑149 (MQ=255) gAACCCACCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTg > 4:6415/1‑148 (MQ=255) aCCCACCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGa < 6:34267/148‑1 (MQ=255) aCCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAgtgt > 2:67754/1‑149 (MQ=255) ccGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCCTGGAAGTGTc < 6:329927/149‑1 (MQ=255) gACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTCg > 2:301854/1‑148 (MQ=255) | GCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTCG > NC_002947/5536930‑5537222 |
| 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: ‑ |